Skip to content
Snippets Groups Projects
Select Git revision
  • 5eb8713e6c788894db565a6b45b2c8d32b1fc05e
  • master default protected
  • develop-GA
  • timeFstatmap
  • add-higher-spindown-components
  • develop-DK
  • adds-header-to-grid-search
  • v1.2
  • v1.1.2
  • v1.1.0
  • v1.0.1
11 results

generate_data.py

Blame
  • Forked from Gregory Ashton / PyFstat
    Source project has a limited visibility.
    waveform.py 497 B
    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