Skip to content
Snippets Groups Projects
Commit fb582da9 authored by Colm Talbot's avatar Colm Talbot
Browse files

Merge branch '615-update-example-priors' into 'master'

Resolve "Update example priors"

Closes #615

See merge request !1075
parents e516bbce d1f7246c
No related branches found
No related tags found
1 merge request!1075Resolve "Update example priors"
Pipeline #355284 passed
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=31) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=25, maximum=31)
mass_1 = Constraint(name='mass_1', minimum=10, maximum=80) mass_1 = Constraint(name='mass_1', minimum=10, maximum=80)
mass_2 = Constraint(name='mass_2', minimum=10, maximum=80) mass_2 = Constraint(name='mass_2', minimum=10, maximum=80)
a_1 = Uniform(name='a_1', minimum=0, maximum=0.99) a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
a_2 = Uniform(name='a_2', minimum=0, maximum=0.99) a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
tilt_1 = Sine(name='tilt_1', boundary='reflective') tilt_1 = Sine(name='tilt_1')
tilt_2 = Sine(name='tilt_2', boundary='reflective') tilt_2 = Sine(name='tilt_2')
phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic') 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') phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
luminosity_distance = PowerLaw(alpha=2, name='luminosity_distance', minimum=50, maximum=2000, unit='Mpc', latex_label='$d_L$') luminosity_distance = PowerLaw(alpha=2, name='luminosity_distance', minimum=50, maximum=2000, unit='Mpc', latex_label='$d_L$')
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
theta_jn = Sine(name='theta_jn', boundary='reflective') theta_jn = Sine(name='theta_jn')
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic') psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic') phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
# Example of priors for GW170817 # Example of priors for GW170817
# mass_1 = Constraint(name='mass_1', minimum=1.2, maximum=1.99) # mass_1 = Constraint(name='mass_1', minimum=1.2, maximum=1.99)
# mass_2 = Constraint(name='mass_2', minimum=1., maximum=1.4) # 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) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=1.197, maximum=1.198, unit='$M_{\odot}$', boundary=None)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1.0, boundary=None) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1.0, boundary=None)
a_1 = Uniform(name='a_1', minimum=0, maximum=0.05) a_1 = Uniform(name='a_1', minimum=0, maximum=0.05)
a_2 = Uniform(name='a_2', minimum=0, maximum=0.05) a_2 = Uniform(name='a_2', minimum=0, maximum=0.05)
tilt_1 = Sine(name='tilt_1') tilt_1 = Sine(name='tilt_1')
...@@ -15,7 +15,7 @@ ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ...@@ -15,7 +15,7 @@ ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
theta_jn = Sine(name='theta_jn') theta_jn = Sine(name='theta_jn')
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic') psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * 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_1 = Constraint(name='lambda_1', minimum=0, maximum=5000)
# lambda_2 = Constraint(name='lambda_2', 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) 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) delta_lambda_tilde = Uniform(name='delta_lambda_tilde', minimum=-5000, maximum=5000, boundary=None)
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=5, maximum=100) mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100) mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=25, maximum=100)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=5, maximum=100) mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100) mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=25, maximum=100)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5) mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5) mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=0.4, maximum=4.4)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5) mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5) mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=0.4, maximum=4.4)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=5, maximum=100) mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100) mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=25, maximum=100)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=5, maximum=100) mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100) mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=25, maximum=100)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5) mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5) mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=0.4, maximum=4.4)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
...@@ -4,8 +4,8 @@ ...@@ -4,8 +4,8 @@
# Lines beginning "#" are ignored. # Lines beginning "#" are ignored.
mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5) mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5) mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4) chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=0.4, maximum=4.4)
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3) luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec') dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic') ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment