Skip to content

Instantly share code, notes, and snippets.

@avivajpeyi
Created November 3, 2022 00:18
Show Gist options
  • Save avivajpeyi/22bc160785976a906613d3caaa4622c4 to your computer and use it in GitHub Desktop.
Save avivajpeyi/22bc160785976a906613d3caaa4622c4 to your computer and use it in GitHub Desktop.
Example GW150914 ini
ntasks-per-node = 48
time = 10:30:00
nodes = 4
mem-per-cpu="2GB"
extra-lines = "module --force purge && module load git/2.18.0 git-lfs/2.4.0 gcc/9.2.0 openmpi/4.0.2 numpy/1.19.2-python-3.8.5 mpi4py/3.0.3-python-3.8.5 && module unload zlib && source /fred/oz980/avajpeyi/envs/sstar_venv/bin/activate"
# The following are the same as bilby_pipe
label = GW150914
outdir = outdir_GW150914
trigger-time = 1126259462.391
psd-dict = {H1=GWTC1_H1_PSD.dat,L1=GWTC1_L1_PSD.dat}
channel-dict = {H1:DCS-CALIB_STRAIN_C02, L1:DCS-CALIB_STRAIN_C02}
data-dict = {H1=H1_data.gwf,L1=L1_data.gwf}
frequency-domain-source-model = lal_binary_black_hole
prior-dict = {
chirp_mass = bilby.gw.prior.UniformInComponentsChirpMass(name='chirp_mass', minimum=9, maximum=69, unit='$M_{\odot}$'),
mass_ratio = bilby.gw.prior.UniformInComponentsMassRatio(name='mass_ratio', minimum=0.125, maximum=1),
mass_1 = Constraint(name='mass_1', 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_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, 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$'),
dec = Cosine(name='dec'),
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic'),
theta_jn = Sine(name='theta_jn'),
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic'),
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic'),
}
gaussian-noise = False
detectors = [H1, L1]
duration = 8
coherence-test=False
minimum-frequency = 20
sampling-frequency = 2048
distance-marginalization=True
phase-marginalization = False
time-marginalization = False
jitter-time = False
reference-frequency = 20.0
waveform-approximant = IMRPhenomXPHM
nlive = 512
nact = 10
walks=200
maxmcmc=15000
check-point-deltaT=3600
nestcheck=False
n-parallel=1
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment