Changes

Page history
Update SummaryReview script authored by Nicola De Lillo's avatar Nicola De Lillo
Show whitespace changes
Inline Side-by-side
SummaryReview-script.md
View page @ 8d36fab6
......@@ -72,125 +72,43 @@ for index in params:
which outputs
```
luminosity_distance [194.16 282.01]
phi_12 [0.33 5.94]
#### Obtained using PESummary
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]
chi_p [0.01 0.11]
a_1 [0. 0.11]
luminosity_distance [194.16 282.01]
theta_jn [0.63 2.39]
L1_matched_filter_snr_abs [21.99 22.32]
H1_matched_filter_snr_abs [10.51 10.76]
V1_matched_filter_snr_abs [3.69 4.49]
```
These numbers have been compared to the reviewed results for GW190814 `combinedPHM` at the bottom of this page [1]
* upper value `mass ratio`
* lower value `total_mass_source`
* lower value `mass_1_source`
* lower limit `theta_jn `
* both limits for `L1_matched_filter_snr_abs ` (0.05 difference for the lower limit)
* Same for `H1_matched_filter_snr_abs `
* Same for `V1_matched_filter_snr_abs` , with a 0.13 difference.
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
parameter_dict = {
"chirp_mass_source": "Chirp mass $\mathcal{M}/M_{\odot}$",
"total_mass_source": "Total mass $M/M_{\odot}$",
"mass_ratio": "Mass ratio $q$",
"mass_1_source": "Primary mass $m_{1}/M_{\odot}$",
"mass_2_source": "Secondary mass ${m_{2}/M_{\odot}}$",
"chi_eff": "Effective inspiral spin parameter $\chi_{eff}$",
"chi_p": "Effective precession parameter $\chi_{p}$",
"a_1": "Dimensionless primary spin magnitude $a_{1}$",
"luminosity_distance": "Luminosity Distance $D_{L}/\text{Mpc}$",
"redshift": "Source redshift $z$",
"theta_jn": "Inclination angle $\\theta_{\mathrm{JN}}$",
"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}}$",
}
waveform_model = "combinedPHM"
print("---------- {} ----------\n".format(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(
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.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`. Plus I have added two lines that print the intervals in the form [min, max]
The output of the script is:
These numbers have been compared to the reviewed results for GW190814 `combinedPHM` at the bottom of this page [1], which correspond to numbers cited in the discovery paper [2], which we report her:
```
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
#### GW190814 reviewed numbers
chirp_mass_source = 6.09^0.06_0.06 | [6.03 6.15]
total_mass_source = 25.52^1.21_1.22 | [24.31 26.73]
mass_ratio = 0.11^0.01_0.01 | [0.1 0.13]
mass_1_source = 22.91^1.32_1.32 | [21.58 24.23]
mass_2_source = 2.62^0.11_0.12 | [2.5 2.73]
chi_eff = -0.02^0.06_0.07 | [-0.09 0.04]
chi_p = 0.05^0.06_0.04 | [0.01 0.11]
a_1 = 0.03^0.08_0.03 | [0.0 0.11]
luminosity_distance = 237.97^44.04_43.81 | [194.16 282.01]
redshift = 0.05^0.01_0.01 | [0.04 0.06]
theta_jn = 0.9^1.49_0.28 | [0.63 2.39]
L1_matched_filter_snr_abs = 22.2^0.13_0.21 | [21.99 22.32]
H1_matched_filter_snr_abs = 10.66^0.1_0.15 | [10.51 10.76]
V1_matched_filter_snr_abs = 4.27^0.22_0.58 | [3.69 4.49]
network_matched_filter_snr = 24.98^0.14_0.22 | [24.76 25.12]
```
[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
[2] [https://dcc.ligo.org/LIGO-P190814-v2](https://dcc.ligo.org/LIGO-P190814-v2)
\ No newline at end of file