O3-review/GW150914_IMRPhenomPv2_lalinference.md

@@ -79,3 +79,75 @@ dataL1 = TimeSeries.get('L1:DCS-CALIB_STRAIN_C02', start_time, end_time)
 TimeSeries.write(dataH1, target='dataH1.gwf',format='gwf')
 TimeSeries.write(dataL1, target='dataL1.gwf',format='gwf')
 ```
+
+# `run_GW150914_with_cal.ini`
+
+```
+label = GW150914_pbilby_master_data_dump_1500live
+outdir = /fred/oz117/rsmith/projects/pbilby_review/150914/outdir
+
+trigger_time = 1126259462.3910
+duration = 4
+
+waveform_approximant = IMRPhenomPv2
+
+prior-file = /fred/oz117/rsmith/projects/pbilby_review/150914/GW150914.prior
+
+psd_dict = {H1=h1_psd.dat, L1=l1_psd.dat}
+channel_dict = {H1:DCS-CALIB_STRAIN_C02, L1:DCS-CALIB_STRAIN_C02}
+data_dict = {H1=dataH1.gwf, L1=dataL1.gwf}
+
+minimum-frequency=20.0
+maximum-frequency=1024
+
+calibration-model=CubicSpline
+spline-calibration-envelope-dict = {H1:GWTC1_GW150914_H_CalEnv.txt, L1:GWTC1_GW150914_L_CalEnv.txt}
+spline_calibration-nodes = 10
+distance-marginalization=True
+phase-marginalization=True
+time-marginalization=True
+```
+
+# `GW150914.prior`
+
+```
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=9.0, maximum=69.9, unit='$M_{\odot}$')
+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')
+```
+
+# `slurm_submit.sh`
+
+```
+#!/bin/bash
+#
+#SBATCH --job-name=150914
+#
+#SBATCH --nodes=40
+#SBATCH --ntasks-per-node=16
+
+#SBATCH --time=10:00:00
+#SBATCH --mem-per-cpu=300
+#SBATCH --output=logs/%A_%a.out
+
+source /home/rsmith/.bashrc
+source activate parallel_bilby
+export MKL_NUM_THREADS="1"
+export MKL_DYNAMIC="FALSE"
+export OMP_NUM_THREADS=1
+export MPI_PER_NODE=16
+mpirun parallel_bilby_analysis /fred/oz117/rsmith/projects/pbilby_review/150914/outdir/GW150914_pbilby_master_data_dump_1500live_data_dump.pickle --nlive 1500
+```
\ No newline at end of file