From 740214b6471897e6acb322e02cbf37f1ac1630f0 Mon Sep 17 00:00:00 2001
From: Yifan Wang <yifan.wang@aei.mpg.de>
Date: Mon, 31 May 2021 16:20:39 +0200
Subject: [PATCH] add pairty violation waveform
---
birefringence-waveform/waveform.py | 12 +++++-------
1 file changed, 5 insertions(+), 7 deletions(-)
diff --git a/birefringence-waveform/waveform.py b/birefringence-waveform/waveform.py
index 4cd1a86..41ec79b 100644
--- a/birefringence-waveform/waveform.py
+++ b/birefringence-waveform/waveform.py
@@ -3,17 +3,15 @@ import numpy as np
import pycbc.conversions
import lal
-def IMRPhenomParity(parity_lambdatilt, **kwds):
+def IMRPhenomParity(parity_Amu, **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;
+ temp = parity_Amu * 1e9 * lal.QE_SI * lal.PI * lal.PI / lal.H0_SI
- deltaPhi1 = tempVal * f * f;
+ hp_parity = hp + hc * temp * hp.sample_frequencies **2
+ hc_parity = hc - hp * temp * hp.sample_frequencies **2
- hp_parity = hp + hc * deltaPhi1;
- hc_parity = hc - hp * deltaPhi1;
-
- return hpd,hcd
+ return hp_parity, hc_parity
--
GitLab