Update GW150914_IMRPhenomPv2_lalinference authored by Rory Smith's avatar Rory Smith
...@@ -40,11 +40,11 @@ The ini file does not include the sampler args. These are passed to parallel_bil ...@@ -40,11 +40,11 @@ The ini file does not include the sampler args. These are passed to parallel_bil
The files `slurm_submit.sh`, `get_gwf.py` and `run_GW150914_with_cal.ini` can be found here on CIT `/home/rory.smith/projects/pbilby_review/150914` The files `slurm_submit.sh`, `get_gwf.py` and `run_GW150914_with_cal.ini` can be found here on CIT `/home/rory.smith/projects/pbilby_review/150914`
## How we get data ## Scripts
The first step in the pipeline `python get_gwf.py` uses `gwpy` to get strain data in gwf format. The `.gwf` file is then passed into the subsequent stages of the pipeline. An example of how the data retrieval works is shown below: The first step in the pipeline `python get_gwf.py` uses `gwpy` to get strain data in gwf format. The `.gwf` file is then passed into the subsequent stages of the pipeline. An example of how the data retrieval works is shown below:
# `get_gwf.py` ## `get_gwf.py`
``` ```
from gwpy.timeseries import TimeSeries from gwpy.timeseries import TimeSeries
...@@ -61,7 +61,7 @@ TimeSeries.write(dataH1, target='dataH1.gwf',format='gwf') ...@@ -61,7 +61,7 @@ TimeSeries.write(dataH1, target='dataH1.gwf',format='gwf')
TimeSeries.write(dataL1, target='dataL1.gwf',format='gwf') TimeSeries.write(dataL1, target='dataL1.gwf',format='gwf')
``` ```
# `run_GW150914_with_cal.ini` ## `run_GW150914_with_cal.ini`
``` ```
label = GW150914_pbilby_master_data_dump_1500live label = GW150914_pbilby_master_data_dump_1500live
...@@ -89,7 +89,7 @@ phase-marginalization=True ...@@ -89,7 +89,7 @@ phase-marginalization=True
time-marginalization=True time-marginalization=True
``` ```
# `GW150914.prior` ## `GW150914.prior`
``` ```
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1) mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
...@@ -110,7 +110,7 @@ psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic') ...@@ -110,7 +110,7 @@ psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic') phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
``` ```
# `slurm_submit.sh` ## `slurm_submit.sh`
``` ```
#!/bin/bash #!/bin/bash
... ...
......