Remnant-fits.md

@@ -314,7 +314,48 @@ Nathan and David sign off on !254, as of 10th March 2023, with the reference out
 
 ## example outputs
 
-🚧
+
+## 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:
+
+```python
+from pesummary.gw.file.read import read
+import numpy as np
+
+print("-------- IMRPhenomPv3HM --------")
+f = read("PhenomPv3HM/samples/PhemomPv3HM_pesummary.dat")
+nathan = np.genfromtxt("IMRPhenomPv3HM_final_state_check_final_spin_with_chi_p.dat", names=True)
+fm_n = nathan["Mf_IMRPhenomPv3HM"]
+fs_n = nathan["af_IMRPhenomPv3HM"]
+
+difference = f.samples_dict["final_mass"] - fm_n
+print("max difference for 'final_mass': {}".format(np.max(np.abs(difference))))
+difference = f.samples_dict["final_spin"] - fs_n
+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.dat", names=True)
+fm_n = nathan["Mf_SEOBNRv4PHM"]
+fs_n = nathan["af_SEOBNRv4PHM"]
+
+difference = f.samples_dict["final_mass"] - fm_n
+print(np.max(np.abs(difference)))
+difference = f.samples_dict["final_spin"] - fs_n
+print(np.max(np.abs(difference)))
+```
+
+With output:
+
+```bash
+-------- IMRPhenomPv3HM --------
+max difference for 'final_mass': 6.66773303237278e-11
+max difference for 'final_spin': 2.536082455151245e-12
+-------- SEOBNRv4PHM --------
+```
 
 ## Review comments