... | ... | @@ -736,11 +736,36 @@ assert "final_mass" in data.dtype.names |
|
|
To confirm that the remnant quantities calculated for non-spinning samples are named `{}` rather than `{}_non_evolved`, we run the following command line:
|
|
|
|
|
|
```bash
|
|
|
rsync -avz --rsh="gsissh -p 2222" ldas-pcdev1.ligo.caltech.edu:/home/anuradha.gupta/public_html/spin_evolution_test/pe_runs_160719/non_spinning_ns_run/1135136363.571641445-0/H1L1/posterior_samples.dat ./non_spinning.dat
|
|
|
|
|
|
rsync -avz --rsh="gsissh -p 2222" ldas-pcdev1.ligo.caltech.edu:/home/anuradha.gupta/spin_evolution_test/pe_runs_160719/non_spinning_ns_run/run/posterior_samples/posterior_H1L1_1135136363.571641445-0.hdf5 ./non_spinning.hdf5
|
|
|
|
|
|
summarypages --webdir /home/charlie.hoy/public_html/LVC/projects/pesummary_review/v0.5.6/non_spinning_remnant --samples /home/charlie.hoy/projects/pesummary_review/remnant/non_spinning/non_spinning.hdf5 --gw --labels non_spinning --multi_process 10
|
|
|
```
|
|
|
|
|
|
To check the samples, we compare to the `cbcBayesPostproc` code with the following script:
|
|
|
|
|
|
```python
|
|
|
import numpy as np
|
|
|
|
|
|
posplots = np.genfromtxt("/home/charlie.hoy/projects/pesummary_review/remnant/non_spinning/non_spinning.dat", names=True)
|
|
|
pesummary = np.genfromtxt("/home/charlie.hoy/public_html/LVC/projects/pesummary_review/v0.5.6/non_spinning_remnant/samples/non_spinning_pesummary.dat", names=True)
|
|
|
|
|
|
remnant_posplots = ["mf_evol", "af_evol", "mf_source_evol", "e_rad_evol", "l_peak_evol"]
|
|
|
remnant_pesummary = ["final_mass", "final_spin", "final_mass_source", "radiated_energy", "peak_luminosity"]
|
|
|
for _pos_param, _pes_param in zip(remnant_posplots, remnant_pesummary):
|
|
|
print("{} - max difference: {}".format(
|
|
|
_pos_param, np.max(np.abs(posplots[_pos_param] - pesummary[_pes_param]))
|
|
|
))
|
|
|
```
|
|
|
```bash
|
|
|
mf_evol - max difference: 8.526512829121202e-14
|
|
|
af_evol - max difference: 2.0794477251229182e-13
|
|
|
mf_source_evol - max difference: 0.012373066036786895
|
|
|
e_rad_evol - max difference: 0.0006290199782483796
|
|
|
l_peak_evol - max difference: 0.00025579293297539607
|
|
|
```
|
|
|
|
|
|
## Review comments
|
|
|
|
|
|
## sign-off |
|
|
\ No newline at end of file |