|
|
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](uploads/2649b62e8d339b235402bd486bb53100/test.py)
|
|
|
|
|
|
Plots:
|
|
|
Plot of the waveform
|
|
|
![fd_waveform](uploads/933eb5eaefc27020810a712b799dd9e4/fd_waveform.png)
|
|
|
![abs_fd_waveform](uploads/f36a64504ebbee84cdd7a0f719e50e70/abs_fd_waveform.png)
|
|
|
Plot of the phase and strain difference
|
|
|
![phasedifference](uploads/5537981b3128fff21a5f9c2726481325/phasedifference.png)
|
|
|
![straindifference](uploads/f1264213786f3a1e6cbd745abda045cd/straindifference.png)
|
|
|
|
... | ... | @@ -21,4 +28,4 @@ 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. |
|
|
* The waveforms outputted between bilby and pycbc have a phase difference of 2pi. The phase here refers to (np.angle(hf)). No idea why. |
|
|
\ No newline at end of file |