Skip to content

Changes

Page history
Update Remnant fits authored by David Keitel's avatar David Keitel
Show whitespace changes
Inline Side-by-side
Remnant-fits.md
View page @ 5e3e28c9
...@@ -6,16 +6,21 @@ To calculate the non-evolved remnant fits, we call the functions in `lalinferenc ...@@ -6,16 +6,21 @@ To calculate the non-evolved remnant fits, we call the functions in `lalinferenc
The conversions are handled in by the [`pesummary.gw.file.conversions` module](https://git.ligo.org/lscsoft/pesummary/-/blob/master/pesummary/gw/file/conversions.py) and specifically the [`pesummary.gw.file.conversion._Conversion` class](https://git.ligo.org/lscsoft/pesummary/-/blob/master/pesummary/gw/file/conversions.py#L534). The conversions are handled in by the [`pesummary.gw.file.conversions` module](https://git.ligo.org/lscsoft/pesummary/-/blob/master/pesummary/gw/file/conversions.py) and specifically the [`pesummary.gw.file.conversion._Conversion` class](https://git.ligo.org/lscsoft/pesummary/-/blob/master/pesummary/gw/file/conversions.py#L534).
To check the samples produced by PESummary are correct, we compared them to cbcBayesPostProc. We copied the following files from HAWK to CIT: ## Related MRs
* !241
* !242
## example outputs
* https://ldas-jobs.ligo.caltech.edu/~charlie.hoy/projects/pesummary_review/remnant_fits/non_evolved_average_precessing/comparison_to_master/html/PhenomPv3HM_PhenomPv3HM.html
* To check the samples produced by PESummary are correct, we compared them to cbcBayesPostProc. We copied the following files from HAWK to CIT:
```bash ```bash
$ rsync -avz --rsh="gsissh" ligo.gravity.cf.ac.uk:/home/sebastian.khan/public_html/LVC/projects/o3/s190521g/run2-8s/1242442967.4472656-333631/H1L1V1/posterior_samples.dat . $ rsync -avz --rsh="gsissh" ligo.gravity.cf.ac.uk:/home/sebastian.khan/public_html/LVC/projects/o3/s190521g/run2-8s/1242442967.4472656-333631/H1L1V1/posterior_samples.dat .
$ rsync -avz --rsh="gsissh" ligo.gravity.cf.ac.uk:/home/sebastian.khan/projects/o3/s190521g/paper-runs/run2-8s/run/posterior_samples/posterior_H1L1V1_1242442967.4472656-333631.hdf5 . $ rsync -avz --rsh="gsissh" ligo.gravity.cf.ac.uk:/home/sebastian.khan/projects/o3/s190521g/paper-runs/run2-8s/run/posterior_samples/posterior_H1L1V1_1242442967.4472656-333631.hdf5 .
``` ```
We then ran the following script: We then ran the following script:
```python ```python
from pesummary.gw.file.read import read from pesummary.gw.file.read import read
from pesummary.gw.file.standard_names import lalinference_map from pesummary.gw.file.standard_names import lalinference_map
...@@ -63,9 +68,7 @@ if len(not_included): ...@@ -63,9 +68,7 @@ if len(not_included):
print("Parameters not included in this analysis: {}".format(", ".join(not_included))) print("Parameters not included in this analysis: {}".format(", ".join(not_included)))
print("--------------------") print("--------------------")
``` ```
With output: With output:
```bash ```bash
$ python compare_output.py $ python compare_output.py
Maximum difference for l_peak_nonevol: 1.0045297926808416e-12 Maximum difference for l_peak_nonevol: 1.0045297926808416e-12
...@@ -73,12 +76,6 @@ Maximum difference for af_nonevol: 1.3800072196090696e-13 ...@@ -73,12 +76,6 @@ Maximum difference for af_nonevol: 1.3800072196090696e-13
Maximum difference for mf_nonevol: 7.958078640513122e-13 Maximum difference for mf_nonevol: 7.958078640513122e-13
Maximum difference for e_rad_nonevol: 4.1033842990145786e-13 Maximum difference for e_rad_nonevol: 4.1033842990145786e-13
``` ```
## Related MRs
* !241
* !242
## Review comments ## Review comments
### David Keitel 2020.03.04 ### David Keitel 2020.03.04
......