SummaryReview-script.md

@@ -126,8 +126,7 @@ These numbers have been compared to the reviewed results for GW190814 `combinedP
 
 My guess is that this is due to different versions of `np.percentile` (pesummary uses ). To cross-check this, I have personally run the script at the page linked [1]:
 
-```
-
+```python
 import numpy as np
 from pesummary.gw.file.read import read
 
@@ -143,32 +142,55 @@ parameter_dict = {
         "luminosity_distance": "Luminosity Distance $D_{L}/\text{Mpc}$",
         "redshift": "Source redshift $z$",
         "theta_jn": "Inclination angle $\\theta_{\mathrm{JN}}$",
-        "L1_matched_filter_snr": "Signal to Noise ratio in \\llo $\\rho_{\mathrm{L}}$",
-        "H1_matched_filter_snr": "Signal to Noise ratio in \\lho $\\rho_{\mathrm{H}}$",
-        "V1_matched_filter_snr": "Signal to Noise ratio in \\virgo $\\rho_{\mathrm{V}}$",
+        "L1_matched_filter_snr_abs": "Signal to Noise ratio in \\llo $\\rho_{\mathrm{L}}$",
+        "H1_matched_filter_snr_abs": "Signal to Noise ratio in \\lho $\\rho_{\mathrm{H}}$",
+        "V1_matched_filter_snr_abs": "Signal to Noise ratio in \\virgo $\\rho_{\mathrm{V}}$",
         "network_matched_filter_snr": "Network Signal to Noise ratio $\\rho_{\mathrm{HLV}}$",
     }
 
-interested = ["PhenomPv3HM", "SEOBNRv4PHM", "combinedPHM"]
-for waveform_model in interested:
+
+waveform_model =  "combinedPHM"
+
 print("---------- {} ----------\n".format(waveform_model))
-    f = read(result_file_map[waveform_model])
+f = read('/home/nicola.delillo/pesummary-review/summaryreview/combined_PHM.dat')
 f.generate_all_posterior_samples()
 params = parameter_dict.keys()
 
 for param in params:
     samples = f.samples_dict[param]
     print(
-            "{} = {}^{}_{}".format(
+        "{} = {}^{}_{} [{} {}]".format(
             param, np.round(np.median(samples), 2),
             np.round(np.percentile(samples, 95) - np.median(samples), 2),
-                np.round(np.median(samples) - np.percentile(samples, 5), 2)
+            np.round(np.median(samples) - np.percentile(samples, 5), 2),
+            np.round(np.percentile(samples, 5), 2),
+            np.round(np.percentile(samples, 95), 2)
+            
         )
     )
-
 ```
 
-I have adapted the script to run only on `combinedPHM` data. I had also to change `matched_filter_snr` tag in `matched_filter_snr_abs`.
+I have adapted the script to run only on `combinedPHM` data. I had also to change `matched_filter_snr` tag in `matched_filter_snr_abs`. Plus I have added  two lines that print the intervals in the form [min, max] 
+
 
+The output of the script is:
+
+```
+chirp_mass_source = 6.09^0.06_0.06
+total_mass_source = 25.52^1.21_1.22
+mass_ratio = 0.11^0.01_0.01
+mass_1_source = 22.91^1.32_1.32
+mass_2_source = 2.62^0.11_0.12
+chi_eff = -0.02^0.06_0.07
+chi_p = 0.05^0.06_0.04
+a_1 = 0.03^0.08_0.03
+luminosity_distance = 237.97^44.04_43.81
+redshift = 0.05^0.01_0.01
+theta_jn = 0.9^1.49_0.28
+L1_matched_filter_snr_abs = 22.2^0.13_0.21
+H1_matched_filter_snr_abs = 10.66^0.1_0.15
+V1_matched_filter_snr_abs = 4.27^0.22_0.58
+network_matched_filter_snr = 24.98^0.14_0.22
+```
 
 [1]: [https://git.ligo.org/publications/gw190814/gw190814-discovery/-/wikis/PE-result-table](https://git.ligo.org/publications/gw190814/gw190814-discovery/-/wikis/PE-result-table) . There are sum difference in the following parameters
\ No newline at end of file