GW170809

# An ini file for running H1L1V1 searches for GW170809
gracedb = G296853
label = GW170809
detectors = [H1, L1, V1]
coherence-test = True
duration = 16
psd-duration = 1024
outdir = results_GW170809
create-summary = True
sampler = dynesty
sampler-kwargs = {'nlive': 2048}
channel_names = [H1:DCH-CLEAN_STRAIN_C02, L1:DCH-CLEAN_STRAIN_C02, V1:Hrec_hoft_V1O2Repro2A_16384Hz]
prior_file = GW170809.prior
reference-frequency = 20
query-types = [H1_CLEANED_HOFT_C02, L1_CLEANED_HOFT_C02, V1O2Repro2A]
accounting = ligo.dev.o3.cbc.pe.lalinference
create-summary = True
create-output = True

# These are the priors for analysing GW170809.
mass_1 = Uniform(name='mass_1', minimum=14.3, maximum=88.6, unit='$M_{\\odot}$')
mass_2 = Uniform(name='mass_2', minimum=5.51, maximum=49.9, unit='$M_{\\odot}$')
chirp_mass =  Constraint(name='chirp_mass', minimum=12.3, maximum=43.5)
mass_ratio =  Constraint(name='mass_ratio', minimum=0.125, maximum=1)
total_mass = Constraint(name='total_mass', minimum=28.3, maximum=100)
a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
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 = PowerLaw(alpha=2, name='luminosity_distance', minimum=168, maximum=4000, unit='Mpc', latex_label='$d_L$')
dec =  Cosine(name='dec')
ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi)
theta_jn =  Sine(name='theta_jn')
psi =  Uniform(name='psi', minimum=0, maximum=np.pi)
phase =  Uniform(name='phase', minimum=0, maximum=2 * np.pi)