All results included here have been obtained using approximant IMRPhenomPv2.
PSD Comparison
- lalinference version: 6.50
- bilby version: 0.3.5
Comparison between PSD generated by lalinference and bilby. Note the roll-off at high frequency in Bilby's PSD.
Comparison between Lalinference and bilby
Short summary of settings:
- approximant = IMRPhenomPv2
- sample_rate = 4096
- reference_frequency = 50
- npoints = 512
- no calibration!
Based on the default
prior_files/GW150914.prior,GW150914_new.priorhas a modified range for the component masses. This was changed to match the default priors used in lalinference for the masses. The full prior setting used are listed below for completeness.
mass_1 = Uniform(name='mass_1', minimum=20, maximum=50, unit='$M_{\\odot}$')
mass_2 = Uniform(name='mass_2', minimum=20, maximum=50, unit='$M_{\\odot}$')
a_1 = Uniform(name='a_1', minimum=0, maximum=0.9)
a_2 = Uniform(name='a_2', minimum=0, maximum=0.9)
tilt_1 = Sine(name='tilt_1')
tilt_2 = Sine(name='tilt_2')
phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi)
phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi)
luminosity_distance = bilby.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=0, maximum=2e3, unit='Mpc')
dec = Cosine(name='dec', minimum=-np.pi/2, maximum=np.pi/2)
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi)
iota = Sine(name='iota', minimum=0, maximum=np.pi)
psi = Uniform(name='psi', minimum=0, maximum=np.pi)
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi)
geocent_time = Uniform(1126259462.322, 1126259462.522, name='geocent_time', unit='$s$')-> For future need to use PowerLaw prior for distance:
luminosity_distance = bilby.gw.prior.PowerLaw(name='luminosity_distance', minimum=0, maximum=2e3, unit='Mpc')
Version comparison (v 0.3.1/0.3.4/0.3.5)
