add pairty violation / birefringence pycbc waveform plugin

birefringence-waveform/waveform.py

+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