... | ... | @@ -317,9 +317,7 @@ Nathan and David sign off on !254, as of 10th March 2023, with the reference out |
|
|
|
|
|
## checks against an independent script
|
|
|
|
|
|
We compared against ... from Nathan's [independent script](https://ldas-jobs.ligo.caltech.edu/~nathan.johnson-mcdaniel/o3-postprocessing/apply_wf_fits.py) (which is wrapped with some [Condor parallelization](https://ldas-jobs.ligo.caltech.edu/~nathan.johnson-mcdaniel/o3-postprocessing/make_wf_fit_application_pipeline.py)) on ...
|
|
|
|
|
|
The comparison was carried out using this script:
|
|
|
We compared against the results (CIT:/home/nathan.johnson-mcdaniel/o3-postprocessing/SEOBNRv4PHM_run_Cartesian_spins/SEOBNRv4PHM_final_state_samples.dat and CIT:/home/nathan.johnson-mcdaniel/GW190521p1/IMRPhenomPv3HM_final_state_check_final_spin_with_chi_p.dat) from Nathan's [independent script](https://ldas-jobs.ligo.caltech.edu/~nathan.johnson-mcdaniel/o3-postprocessing/apply_wf_fits.py) (which is wrapped with some [Condor parallelization](https://ldas-jobs.ligo.caltech.edu/~nathan.johnson-mcdaniel/o3-postprocessing/make_wf_fit_application_pipeline.py)). The comparison was carried out using this script:
|
|
|
|
|
|
```python
|
|
|
from pesummary.gw.file.read import read
|
... | ... | @@ -338,7 +336,7 @@ print("max difference for 'final_spin': {}".format(np.max(np.abs(difference)))) |
|
|
|
|
|
print("-------- SEOBNRv4PHM --------")
|
|
|
f = read("SEOBNRv4PHM_pesummary.dat")
|
|
|
nathan = np.genfromtxt("SEOBNRv4PHM_final_state_samples_bugfix.dat", names=True)
|
|
|
nathan = np.genfromtxt("SEOBNRv4PHM_final_state_samples.dat", names=True)
|
|
|
fm_n = nathan["Mf_SEOBNRv4PHM"]
|
|
|
fs_n = nathan["af_SEOBNRv4PHM"]
|
|
|
|
... | ... | |