Skip to content
Snippets Groups Projects

GW170817 in open_data_examples

Merged Virginia d'Emilio requested to merge virginia.demilio/bilby:GW170817_example into master
All threads resolved!
Files
2
+ 21
0
# Example of priors for GW170817
# mass_1 = Constraint(name='mass_1', minimum=1.2, maximum=1.99)
# mass_2 = Constraint(name='mass_2', minimum=1., maximum=1.4)
chirp_mass = Uniform(name='chirp_mass', minimum=0.87, maximum=1.74, unit='$M_{\odot}$', boundary=None)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1.0, boundary=None)
a_1 = Uniform(name='a_1', minimum=0, maximum=0.05, boundary='reflective')
a_2 = Uniform(name='a_2', minimum=0, maximum=0.05, boundary='reflective')
tilt_1 = Sine(name='tilt_1', boundary='reflective')
tilt_2 = Sine(name='tilt_2', boundary='reflective')
phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=10, maximum=100, unit='Mpc', boundary=None)
dec = Cosine(name='dec', boundary='reflective')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
theta_jn = Sine(name='theta_jn', boundary='reflective')
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
# lambda_1 = Constraint(name='lambda_1', minimum=0, maximum=5000)
# lambda_2 = Constraint(name='lambda_2', minimum=0, maximum=5000)
lambda_tilde = Uniform(name='lambda_tilde', minimum=0, maximum=1000, boundary=None)
delta_lambda_tilde = Uniform(name='delta_lambda_tilde', minimum=-5000, maximum=5000, boundary=None)
Loading