Update SummaryReview script authored by Nicola De Lillo's avatar Nicola De Lillo
...@@ -12,7 +12,7 @@ Can be seen here. ...@@ -12,7 +12,7 @@ Can be seen here.
#### Step 1 (Run) #### Step 1 (Run)
We run `pesummary` and `cbcBayesPostProc` on the posterior file `.hdf5` indicated at [PESummary-Review#Notes](https://git.ligo.org/lscsoft/pesummary/-/wikis/PESummary-Review#notes) We run `pesummary` and `cbcBayesPostProc` on the posterior file `.hdf5` indicated at [PESummary-Review#Notes](https://git.ligo.org/lscsoft/pesummary/-/wikis/PESummary-Review#notes)
``` ```bash
summaryreview -w /home/nicola.delillo/public_html/pesummary-review/summaryreview -s /home/nicola.delillo/pesummary-review/samples/PROD0_posterior_samples.hdf5 summaryreview -w /home/nicola.delillo/public_html/pesummary-review/summaryreview -s /home/nicola.delillo/pesummary-review/samples/PROD0_posterior_samples.hdf5
``` ```
Where `-w` indicates the output directory. Where `-w` indicates the output directory.
...@@ -23,11 +23,19 @@ The `.hdf5` file was produced with a different version of lalsimulation. This re ...@@ -23,11 +23,19 @@ The `.hdf5` file was produced with a different version of lalsimulation. This re
The webpage produced for the comparison can be seen at: [https://ligo.gravity.cf.ac.uk/~nicola.delillo/pesummary-review/summaryreview/home.html](https://ligo.gravity.cf.ac.uk/~nicola.delillo/pesummary-review/summaryreview/home.html) The webpage produced for the comparison can be seen at: [https://ligo.gravity.cf.ac.uk/~nicola.delillo/pesummary-review/summaryreview/home.html](https://ligo.gravity.cf.ac.uk/~nicola.delillo/pesummary-review/summaryreview/home.html)
Our current version of lalsimulation gives: Our current version of lalsimulation gives:
``` ```bash
$ lalsim.GetApproximantFromString('IMRPhenomPv2_NRtidal') $ lalsim.GetApproximantFromString('IMRPhenomPv2_NRtidal')
79 79
``` ```
We edit the file running
```bash
f = h5py.File('PROD0_posterior_samples.hdf5', 'a')
f['lalinference/lalinference_nest/posterior_samples'].attrs['LAL_APPROXIMANT'] = 79
f.close()
```
## To do list ## To do list
* [ ] check 90% credible intervals * [ ] check 90% credible intervals
* [ ] bayes factors * [ ] bayes factors
... ...
......