add working waveform for parity violation

.gitignore

+__pycache__
+.ipynb_checkpoints

birefringence-waveform/__init__.py

+import .birefringence_waveform

birefringence-waveform/birefringence_waveform.py

@@ -18,4 +18,4 @@ def IMRPhenomParity(parity_Amu, **kwds):
     return hp_parity, hc_parity
 
 
-pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)
\ No newline at end of file
+pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)

birefringence-waveform/waveform.py

-from pycbc.waveform import get_fd_waveform
-import numpy as np
-import pycbc.conversions
-import lal
-
-def IMRPhenomParity(parity_Amu, **kwds):
-
-    if 'approximant' in kwds:
-        kwds.pop("approximant")
-    hp, hc = get_fd_waveform(approximant="IMRPhenomXPHM", **kwds)
-
-    temp =  parity_Amu * 1e9 * lal.QE_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