From 1874312e5c297a8abe84fc53014fc335c5cf3e03 Mon Sep 17 00:00:00 2001
From: Yifan Wang <yifan.wang@aei.mpg.de>
Date: Mon, 31 May 2021 13:49:36 +0200
Subject: [PATCH] add pairty violation / birefringence pycbc waveform plugin
---
birefringence-waveform/waveform.py | 19 +++++++++++++++++++
1 file changed, 19 insertions(+)
create mode 100644 birefringence-waveform/waveform.py
diff --git a/birefringence-waveform/waveform.py b/birefringence-waveform/waveform.py
new file mode 100644
index 0000000..4cd1a86
--- /dev/null
+++ b/birefringence-waveform/waveform.py
@@ -0,0 +1,19 @@
+from pycbc.waveform import get_fd_waveform
+import numpy as np
+import pycbc.conversions
+import lal
+
+def IMRPhenomParity(parity_lambdatilt, **kwds):
+
+ if 'approximant' in kwds:
+ kwds.pop("approximant")
+ hp, hc = get_fd_waveform(approximant="IMRPhenomXPHM", **kwds)
+
+ tempVal = lal.H_SI * lal.PI * lal.PI * r / parity_lambdatilt;
+
+ deltaPhi1 = tempVal * f * f;
+
+ hp_parity = hp + hc * deltaPhi1;
+ hc_parity = hc - hp * deltaPhi1;
+
+ return hpd,hcd
--
GitLab