Skip to content
Snippets Groups Projects

Update IMRPhenomTPHMParity

Merged Stephanie Alex Brown requested to merge alexbrown6291/3ogc-projects:master into master
2 files
+ 32
23
Compare changes
  • Side-by-side
  • Inline
Files
2
+ 0
23
def genIMRPhenomParity(**kwds):
from pycbc.waveform import get_fd_waveform
import lal
#print(kwds)
if 'approximant' in kwds:
kwds.pop("approximant")
hp, hc = get_fd_waveform(approximant="IMRPhenomXPHM", **kwds)
temp = kwds['parity_amu'] / 1e9 / lal.QE_SI * lal.H_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
#pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)
def add_me(**kwds):
kwds['cpu_fd']['IMRPhenomParity'] = genIMRPhenomParity
kwds['filter_time_lengths']['IMRPhenomParity'] = kwds['filter_time_lengths']['IMRPhenomXPHM']
Loading