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

update minimal 150914 example

parent 4e56722c
No related branches found
No related tags found
1 merge request!59Updating priors
......@@ -7,7 +7,7 @@ stimation on GW150914 using open data.
import tupak
t0 = tupak.gw.utils.get_event_time("GW150914")
prior = dict(geocent_time=tupak.core.prior.Uniform(t0 - 0.1, t0 + 0.1, name='geocent_time'))
prior = tupak.gw.prior.BBHPriorSet(filename='GW150914.prior')
interferometers = tupak.gw.detector.get_event_data("GW150914")
likelihood = tupak.gw.likelihood.get_binary_black_hole_likelihood(interferometers)
result = tupak.run_sampler(likelihood, prior, label='GW150914')
......
# These are the default priors for analysing GW150914.
mass_1 = Uniform(name='mass_1', minimum=30, maximum=50)
mass_2 = Uniform(name='mass_2', minimum=20, maximum=40)
a_1 = Uniform(name='a_1', minimum=0, maximum=0.8)
a_2 = Uniform(name='a_2', minimum=0, maximum=0.8)
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 = tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=1e3)
dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi)
iota = Sine(name='iota')
psi = Uniform(name='psi', minimum=0, maximum=2 * np.pi)
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi)
geocent_time = Uniform(1126259462.322, 1126259462.522, name='geocent_time')
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