add pairty violation waveform

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