|
|
|
# GW150914 IMRPhenomPv2 (lalinference) |
|
|
\ No newline at end of file |
|
|
|
# Setup etc...
|
|
|
|
|
|
|
|
| bilby | parallel bilby |
|
|
|
|
| ------ | ------ |
|
|
|
|
| [ini file](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_review/calibration_review/bilby/run_GW150914_with_calibration.ini) | cell |
|
|
|
|
| cell | cell |
|
|
|
|
|
|
|
|
# Data files
|
|
|
|
|
|
|
|
| H1 | L1 |
|
|
|
|
| ------ | ------ |
|
|
|
|
| [cal envelope](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_review/calibration_review/bilby/GWTC1_GW150914_H_CalEnv.txt) | [cal envelope](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_review/calibration_review/bilby/GWTC1_GW150914_L_CalEnv.txt) |
|
|
|
|
| [psd](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_review/calibration_review/bilby/h1_psd.dat) | [psd](https://ldas-jobs.ligo.caltech.edu/~gregory.ashton/bilby_review/calibration_review/bilby/l1_psd.dat) |
|
|
|
|
|
|
|
|
# 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')
|
|
|
|
```
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
# GW150914 IMRPhenomPv2 (parallel bilby)
|
|
|
|
|
|
|
|
Run settings are same as in regular bilby run, but the workflow is executed slightly differently and certain args are passed to the run scripts in a slightly different way.
|
|
|
|
|
|
|
|
## Workflow:
|
|
|
|
1. `python get_gwf.py`
|
|
|
|
2. `parallel_bilby_generation run_GW150914_with_cal.ini `
|
|
|
|
3. `sbatch slurm_submit.sh`
|
|
|
|
|
|
|
|
The ini file does not include the sampler args. These are passed to parallel_bilby_analysis within the slurm submit file:
|
|
|
|
|
|
|
|
`mpirun parallel_bilby_analysis --nlive 1000 --dynesty-walks 500` |
|
|
\ No newline at end of file |
