analytic waveforms

Direct check of waveform calling from LALSuite. Call the same waveform parameters and approximants with Bilby and PyCBC to ensure calling of LALSuite is correct.

Comparison of IMRPhenomPv2 waveform: input parameters:

injection_parameters = dict(
    mass_1=36., mass_2=32., a_1=0.2, a_2=0.4, tilt_1=0.4, tilt_2=0.7,
    phi_12=1.7, phi_jl=0.3, luminosity_distance=4000., theta_jn=0.4, psi=2.659,
    phase= 1.3 + np.pi/2., geocent_time=1126259642.413, ra=1.375, dec=0.2108)

With spin parameters transformed to lal/pycbc ones using: bilby.gw.conversion.bilby_to_lalsimulation_spins conversion function.

Script used to generate the plots/calculate the overlap: test.py

Plots:

Overlaps: not normalised with power spectral density

cross polarizion: 0.9999950422312026
plus polarizon: 0.9999958113444587

Notes

• There is a discrepancy between the solar mass maybe other constants in bilby and lal/pycbc
• The waveforms outputted between bilby and pycbc have a phase difference of 2pi. The phase here refers to (np.angle(hf)). No idea why.

Comparison of IMRPhenomD_NRTidal waveform, PyCBC does not have IMRPhenom_Pv2_NRTidal: input parameters:

injection_parameters = dict(
    mass_1=1.5, mass_2=1.3, a_1=0.2, a_2=0.4, tilt_1=0.0, tilt_2=0.0,
    phi_12=1.7,phi_jl=0.3, luminosity_distance=40., theta_jn=0.4, psi=2.659,
    phase=1.3 + np.pi/2., geocent_time=1126259642.413, ra=1.375, dec=-1.2108, lambda_1=400, lambda_2=450)

Similar process, plots generated using script test_tidal.py

Overlaps: not normalised with power spectral density for IMRPhenomD_NRTidal

cross polarizion: 1.0000000000000038
plus polarizon: 0.9999999999999961