| ... | ... | @@ -48,3 +48,69 @@ The single summary pages can be found at: |
|
|
|
In `pesummary` the credible interval for 1D histograms is calculated in a way that includes the 90% of the probability: [median - 5% percentile, median + 95%], as you can see from [plot.py#L106](https://git.ligo.org/lscsoft/pesummary/-/blob/master/pesummary/core/plots/plot.py#L106) and [plot.py#L107](https://git.ligo.org/lscsoft/pesummary/-/blob/master/pesummary/core/plots/plot.py#L107)
|
|
|
|
in the source code.
|
|
|
|
|
|
|
|
To review the credible interval produced, we use
|
|
|
|
|
|
|
|
|
|
|
|
We have run this script that prints the credible interval
|
|
|
|
|
|
|
|
|
|
|
|
```python
|
|
|
|
import numpy as np
|
|
|
|
import pesummary
|
|
|
|
|
|
|
|
from pesummary.gw.file.read import read
|
|
|
|
|
|
|
|
# grabbing the data
|
|
|
|
f = read('/home/nicola.delillo/pesummary-review/summaryreview/combined_PHM.dat')
|
|
|
|
|
|
|
|
# grabbing parameters name
|
|
|
|
params = f.samples_dict.keys()
|
|
|
|
|
|
|
|
#calculating and printing credible intervals for each parameter
|
|
|
|
for index in params:
|
|
|
|
c_intervals = f.samples_dict[index].confidence_interval([5, 95])
|
|
|
|
print(index + ' ' + str(np.round(c_intervals,2)))
|
|
|
|
```
|
|
|
|
|
|
|
|
which outputs
|
|
|
|
|
|
|
|
```
|
|
|
|
luminosity_distance [194.16 282.01]
|
|
|
|
phi_12 [0.33 5.94]
|
|
|
|
chirp_mass_source [6.03 6.15]
|
|
|
|
cos_tilt_1 [-0.97 0.85]
|
|
|
|
phi_jl [0.36 5.83]
|
|
|
|
total_mass_source [24.31 26.73]
|
|
|
|
mass_ratio [0.1 0.13]
|
|
|
|
chirp_mass [6.39 6.43]
|
|
|
|
comoving_distance [186.24 265.82]
|
|
|
|
V1_matched_filter_snr_abs [3.69 4.49]
|
|
|
|
H1_matched_filter_snr_angle [-0.11 0.1 ]
|
|
|
|
mass_2 [2.64 2.86]
|
|
|
|
tilt_2 [0.62 2.64]
|
|
|
|
tilt_1 [0.55 2.89]
|
|
|
|
chi_p [0.01 0.11]
|
|
|
|
mass_1 [22.68 25.56]
|
|
|
|
L1_matched_filter_snr_abs [21.99 22.32]
|
|
|
|
geocent_time [1.24985226e+09 1.24985226e+09]
|
|
|
|
mass_1_source [21.58 24.23]
|
|
|
|
symmetric_mass_ratio [0.08 0.1 ]
|
|
|
|
mass_2_source [2.5 2.73]
|
|
|
|
cos_tilt_2 [-0.88 0.81]
|
|
|
|
phi_2 [0.36 5.95]
|
|
|
|
psi [0.06 3.05]
|
|
|
|
V1_matched_filter_snr_angle [-0.31 0.39]
|
|
|
|
theta_jn [0.63 2.39]
|
|
|
|
ra [0.2 0.4]
|
|
|
|
dec [-0.57 -0.4 ]
|
|
|
|
cos_theta_jn [-0.73 0.81]
|
|
|
|
a_2 [0.06 0.94]
|
|
|
|
total_mass [25.54 28.2 ]
|
|
|
|
phase [0.4 5.88]
|
|
|
|
H1_matched_filter_snr_abs [10.51 10.76]
|
|
|
|
chi_eff [-0.09 0.04]
|
|
|
|
log_likelihood [-45774.32 302.36]
|
|
|
|
phi_1 [0.34 5.93]
|
|
|
|
L1_matched_filter_snr_angle [-0.07 0.06]
|
|
|
|
a_1 [0. 0.11]
|
|
|
|
``` |
