GW170818

# An ini file for running H1, L1, V1, and H1L1V1 searches for GW170818
trigger-time = 1187058327.1
label = cpnest_GW170818_lalinf_psd
maximum-frequency=1024
minimum-frequency=16
sampling-frequency=4096
detectors = [H1, L1, V1]
coherence-test = False
duration = 4
outdir = .
sampler = cpnest
request-cpus = 12
sampler-kwargs = {nlive: 2000, nthreads=12, 'maxmcmc': 1000, n_periodic_checkpoint=1000}
reference-frequency = 20
channel-dict = {H1:DCH-CLEAN_STRAIN_C02, L1:DCH-CLEAN_STRAIN_C02, V1:Hrec_hoft_V1O2Repro2A_16384Hz}
prior_file = GW170818.prior
accounting = ligo.dev.o3.cbc.pe.lalinference
time-marginalization=True
distance-marginalization=True
phase-marginalization=True
create-plots=True
local-generation = True
psd-dict = {H1:GW170818/psd/BayesWave_median_PSD_H1.dat, L1:GW170818/psd/BayesWave_median_PSD_L1.dat, V1:GW170818/psd/BayesWave_median_PSD_V1.dat}

# These are the priors for analysing GW170818.
mass_1 = Uniform(name='mass_1', minimum=14.5, maximum=159, unit='$M_{\\odot}$')
mass_2 = Uniform(name='mass_2', minimum=5.65, maximum=51.5, unit='$M_{\\odot}$')
chirp_mass =  Constraint(name='chirp_mass', minimum=12.3, maximum=45)
mass_ratio =  Constraint(name='mass_ratio', minimum=0.125, maximum=1)
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=100, maximum=2000, 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)