From 2c0edbc3fc61d179b1e26f6f6484cd1224889459 Mon Sep 17 00:00:00 2001
From: Yifan Wang <yifan.wang@aei.mpg.de>
Date: Wed, 28 Jul 2021 17:54:24 +0200
Subject: [PATCH] add IMRPhenomD_NRTitdal_Parity
---
IMRPDT_waveform/.DS_Store | Bin 0 -> 6148 bytes
IMRPDT_waveform/IMRPhenomDNRTparity.py | 22 ++++++++++++++++++++++
2 files changed, 22 insertions(+)
create mode 100644 IMRPDT_waveform/.DS_Store
create mode 100644 IMRPDT_waveform/IMRPhenomDNRTparity.py
diff --git a/IMRPDT_waveform/.DS_Store b/IMRPDT_waveform/.DS_Store
new file mode 100644
index 0000000000000000000000000000000000000000..5008ddfcf53c02e82d7eee2e57c38e5672ef89f6
GIT binary patch
literal 6148
zcmZQzU|@7AO)+F(5MW?n;9!8z45|!R0Z1N%F(jFgL>QrFAPJ2!M?+vV1V%$(Gz3ON
zU^D~<VF)ln+{D2Rp-0Kl5Eu=C(GY-#0H}OW0QD6Z7#JL&bOVG2Nii@oFo3%Nj0_Ac
zFio(203!nfNGnJUNGpg2X=PvpvA|}4wK6b5wK9UcAq)(R;4TS>25V<v1ltVagS9g-
zf^BACV1#IAV1(Mt2<@RTf_gL{^C8+97{Ru~TsKOOhQMeDz(Rl-!Vmz}|E>%SxcdJP
zRior+2#kinunYl47MEZbCs3t{!+W4QHvuXKVuPw;Mo^s$(F3lEVT}ML$bg~*R5_@+
b2Uo?6kTwK}57Iu`5P${HC_Nei0}uiLNUI8I
literal 0
HcmV?d00001
diff --git a/IMRPDT_waveform/IMRPhenomDNRTparity.py b/IMRPDT_waveform/IMRPhenomDNRTparity.py
new file mode 100644
index 0000000..5878ee4
--- /dev/null
+++ b/IMRPDT_waveform/IMRPhenomDNRTparity.py
@@ -0,0 +1,22 @@
+def gen(**kwds):
+ from pycbc.waveform import get_fd_waveform
+ import lal
+
+ if 'approximant' in kwds:
+ kwds.pop("approximant")
+ hp, hc = get_fd_waveform(approximant="IMRPhenomD_NRTidal", **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']['IMRPhenomD_NRTidal_Parity'] = gen
+ kwds['filter_time_lengths']['IMRPhenomD_NRTidal_Parity'] = kwds['filter_time_lengths']['IMRPhenomD']
+
--
GitLab