... | @@ -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
|
... | | ... | |