SummaryReview-script.md

@@ -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]
+```