| ... | ... | @@ -94,18 +94,47 @@ $ tar -xvf GW190521.tar |
|
|
|
```python
|
|
|
|
from pesummary.gw.fetch import fetch_open_data
|
|
|
|
from pesummary.io import read
|
|
|
|
import numpy as np
|
|
|
|
|
|
|
|
for event in ["GW190412", "GW190521"]:
|
|
|
|
f = fetch_open_data(event, unpack=True, read_file=True, path=f"{event}.h5")
|
|
|
|
g = read(f"{event}/{event}.h5")
|
|
|
|
assert sorted(f.labels) == sorted(g.labels)
|
|
|
|
np.testing.assert_almost_equal(np.array(f.samples), np.array(g.samples))
|
|
|
|
assert f.parameters == g.parameters
|
|
|
|
samples_f = f.samples_dict
|
|
|
|
samples_g = g.samples_dict
|
|
|
|
for num, label in enumerate(f.labels):
|
|
|
|
for param in f.parameters[num]:
|
|
|
|
np.testing.assert_almost_equal(
|
|
|
|
samples_f[label][param], samples_g[label][param]
|
|
|
|
)
|
|
|
|
```
|
|
|
|
|
|
|
|
with output:
|
|
|
|
|
|
|
|
```bash
|
|
|
|
|
|
|
|
```
|
|
|
|
Downloading https://dcc.ligo.org/public/0169/P2000223/005/GW190412.tar
|
|
|
|
|===============================================================| 1.6G/1.6G (100.00%) 3m35s
|
|
|
|
2020-11-20 11:07:53 PESummary WARNING : Could not find reference_frequency in input file. Using 20Hz as default
|
|
|
|
2020-11-20 11:08:03 PESummary INFO : Averaging the final spin from the following fits: function bbh_final_spin_precessing_projected_UIB2016 at 0x2b1562563268, function bbh_final_spin_precessing_projected_Healyetal at 0x2b15625631e0, function bbh_final_spin_precessing_HBR2016 at 0x2b15625632f0
|
|
|
|
2020-11-20 11:08:03 PESummary INFO : Averaging the peak luminosity from the following fits: function bbh_peak_luminosity_non_precessing_UIB2016 at 0x2b1562563400, function bbh_peak_luminosity_non_precessing_Healyetal at 0x2b1562563488
|
|
|
|
2020-11-20 11:08:10 PESummary INFO : Averaging the final mass from the following fits: function bbh_final_mass_non_precessing_UIB2016 at 0x2b1562560d08, function bbh_final_mass_non_precessing_Healyetal at 0x2b1562560b70
|
|
|
|
2020-11-20 11:08:10 PESummary WARNING : Could not find reference_frequency in input file. Using 20Hz as default
|
|
|
|
2020-11-20 11:08:20 PESummary INFO : Averaging the final spin from the following fits: function bbh_final_spin_precessing_projected_UIB2016 at 0x2b1562563268, function bbh_final_spin_precessing_projected_Healyetal at 0x2b15625631e0, function bbh_final_spin_precessing_HBR2016 at 0x2b15625632f0
|
|
|
|
2020-11-20 11:08:20 PESummary INFO : Averaging the peak luminosity from the following fits: function bbh_peak_luminosity_non_precessing_UIB2016 at 0x2b1562563400, function bbh_peak_luminosity_non_precessing_Healyetal at 0x2b1562563488
|
|
|
|
2020-11-20 11:08:28 PESummary INFO : Averaging the final mass from the following fits: function bbh_final_mass_non_precessing_UIB2016 at 0x2b1562560d08, function bbh_final_mass_non_precessing_Healyetal at 0x2b1562560b70
|
|
|
|
2020-11-20 11:09:04 PESummary WARNING : Could not find reference_frequency in input file. Using 20Hz as default
|
|
|
|
2020-11-20 11:09:14 PESummary INFO : Averaging the final spin from the following fits: function bbh_final_spin_precessing_projected_UIB2016 at 0x2b1562563268, function bbh_final_spin_precessing_projected_Healyetal at 0x2b15625631e0, function bbh_final_spin_precessing_HBR2016 at 0x2b15625632f0
|
|
|
|
2020-11-20 11:09:14 PESummary INFO : Averaging the peak luminosity from the following fits: function bbh_peak_luminosity_non_precessing_UIB2016 at 0x2b1562563400, function bbh_peak_luminosity_non_precessing_Healyetal at 0x2b1562563488
|
|
|
|
2020-11-20 11:09:21 PESummary INFO : Averaging the final mass from the following fits: function bbh_final_mass_non_precessing_UIB2016 at 0x2b1562560d08, function bbh_final_mass_non_precessing_Healyetal at 0x2b1562560b70
|
|
|
|
2020-11-20 11:09:22 PESummary WARNING : Could not find reference_frequency in input file. Using 20Hz as default
|
|
|
|
2020-11-20 11:09:32 PESummary INFO : Averaging the final spin from the following fits: function bbh_final_spin_precessing_projected_UIB2016 at 0x2b1562563268, function bbh_final_spin_precessing_projected_Healyetal at 0x2b15625631e0, function bbh_final_spin_precessing_HBR2016 at 0x2b15625632f0
|
|
|
|
2020-11-20 11:09:32 PESummary INFO : Averaging the peak luminosity from the following fits: function bbh_peak_luminosity_non_precessing_UIB2016 at 0x2b1562563400, function bbh_peak_luminosity_non_precessing_Healyetal at 0x2b1562563488
|
|
|
|
2020-11-20 11:09:39 PESummary INFO : Averaging the final mass from the following fits: function bbh_final_mass_non_precessing_UIB2016 at 0x2b1562560d08, function bbh_final_mass_non_precessing_Healyetal at 0x2b1562560b70
|
|
|
|
Downloading https://dcc.ligo.org/public/0169/P2000223/005/GW190521.tar
|
|
|
|
|===============================================================| 252M/252M (100.00%) 2m51s
|
|
|
|
````
|
|
|
|
|
|
|
|
Meaning now assertion errors are raised.
|
|
|
|
|
|
|
|
## Reviewer comments |
|
|
\ No newline at end of file |
