From 89985e322e44bec7b53e0e62e6ea2e349335321e Mon Sep 17 00:00:00 2001
From: Stephanie Brown <stephanie.brown@condor1.atlas.local>
Date: Tue, 10 Aug 2021 17:31:51 +0000
Subject: [PATCH] add new file
---
waveform/parity_XPHM.py | 23 +++++++++++++++++++++++
1 file changed, 23 insertions(+)
create mode 100644 waveform/parity_XPHM.py
diff --git a/waveform/parity_XPHM.py b/waveform/parity_XPHM.py
new file mode 100644
index 0000000..7894074
--- /dev/null
+++ b/waveform/parity_XPHM.py
@@ -0,0 +1,23 @@
+def gen(**kwds):
+ from pycbc.waveform import get_fd_waveform
+ import lal
+ #print(kwds)
+
+ if 'approximant' in kwds:
+ kwds.pop("approximant")
+ hp, hc = get_fd_waveform(approximant="IMRPhenomXPHM", **kwds)
+
+ temp = kwds['parity_amu'] / 1e9 / lal.QE_SI * lal.H_SI * lal.PI * lal.PI / lal.H0_SI
+
+ hp_parity = hp + hc * temp * hp.sample_frequencies **2
+ hc_parity = hc - hp * temp * hp.sample_frequencies **2
+
+ return hp_parity, hc_parity
+
+
+#pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)
+
+def add_me(**kwds):
+ kwds['cpu_fd']['IMRPhenomXPHMParity'] = gen
+ kwds['filter_time_lengths']['IMRPhenomXPHMParity'] = kwds['filter_time_lengths']['IMRPhenomXPHM']
+
--
GitLab