8.1. Gravitational waveform analysis (GW)

In this tutorial, you will learn how to extract the tidal deformability from GW analysis.

---

Import the libraries that will be employed in this tutorial.

# Import numpy
import numpy as np
# Import matplotlib
import matplotlib.pyplot as plt
# Import nucleardatapy package
import nucleardatapy as nuda

You can simply print out the properties of the nuda’s function that we will use:

# Explore the nucleardatapy module to find the correct attribute
print(dir(nuda.astro.setupGW))

['__class__', '__delattr__', '__dict__', '__dir__', '__doc__', '__eq__', '__format__', '__ge__', '__getattribute__', '__getstate__', '__gt__', '__hash__', '__init__', '__init_subclass__', '__le__', '__lt__', '__module__', '__ne__', '__new__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', 'init_self', 'print_latex', 'print_output']

---

Complete list of available sources:

sources, sources_small = nuda.astro.gw_sources( )
print('sources:',sources)

sources: ['GW170817', 'GW190425']

source = 'GW170817'
hyps = nuda.astro.gw_hyps( source = source )
print(f'source: {source}, hyps: {hyps}')

source: GW170817, hyps: [1, 2, 3, 4, 5]

for source in sources:
    hyps = nuda.astro.gw_hyps( source = source )
    print(f'source: {source}, hyps: {hyps}')
    for hyp in hyps:
        gw = nuda.astro.setupGW( source = source, hyp = hyp )
    gwav = nuda.astro.setupGWAverage( source = source )

source: GW170817, hyps: [1, 2, 3, 4, 5]
source: GW190425, hyps: [1, 2]

Figure:

nuda.fig.astro_setupGW_fig( None, sources )

Plot name: None
obss: [1, 2, 3, 4, 5]
- table: GW170817 & 1 & $1.1880^{0.004}_{-0.002}$ & $[0.7:1.0]$ & ${400.00}^{400.0}_{-400.0}$ & \cite{BPAbbott:2017} \\\\
- table: GW170817 & 2 & $1.1880^{0.004}_{-0.002}$ & $[0.4:1.0]$ & ${350.00}^{350.0}_{-350.0}$ & \cite{BPAbbott:2017} \\\\
- table: GW170817 & 3 & $1.1867^{0.000119}_{-9.5e-05}$ & $[0.639:0.982]$ & ${222.00}^{420.0}_{-138.0}$ & \cite{SDe:2018} \\\\
- table: GW170817 & 4 & $1.1860^{0.001}_{-0.001}$ & $[0.73:1.0]$ & ${300.00}^{420.0}_{-230.0}$ & \cite{BPAbbott:2019} \\\\
- table: GW170817 & 5 & $1.1860^{0.001}_{-0.001}$ & $[0.53:1.0]$ & ${315.00}^{315.0}_{-315.0}$ & \cite{BPAbbott:2019} \\\\
- table: GW170817 &  &  &  &  &  \\\\
obss: [1, 2]
- table: GW190425 & 1 & $1.4400^{0.02}_{-0.02}$ & $[0.8:1.0]$ & ${300.00}^{300.0}_{-300.0}$ & \cite{BPAbbott:2020} \\\\
- table: GW190425 & 2 & $1.4400^{0.02}_{-0.02}$ & $[0.4:1.0]$ & ${550.00}^{550.0}_{-550.0}$ & \cite{BPAbbott:2020} \\\\
- table: GW190425 &  &  &  &  &  \\\\