Python get_final_freq Examples

Example #1

File: waveform.py Project: bhooshan-gadre/my_pycbc_old

def seobnrrom_final_frequency(**kwds):
    from pycbc.pnutils import get_final_freq
    return get_final_freq("SEOBNRv2", kwds['mass1'], kwds['mass2'],
                          kwds['spin1z'], kwds['spin2z'])

Example #2

def seobnrv2_final_frequency(**kwds):
    return pnutils.get_final_freq("SEOBNRv2", kwds['mass1'], kwds['mass2'],
                                  kwds['spin1z'], kwds['spin2z'])

Example #3

File: waveform.py Project: bhooshan-gadre/pycbc

def seobnrv2_final_frequency(**kwds):
    return pnutils.get_final_freq("SEOBNRv2", kwds['mass1'], kwds['mass2'],
                                  kwds['spin1z'], kwds['spin2z'])

Example #4

File: waveform.py Project: sfairhur/pycbc

def seobnrrom_final_frequency(**kwds):
    from pycbc.pnutils import get_final_freq

    return get_final_freq("SEOBNRv2", kwds["mass1"], kwds["mass2"], kwds["spin1z"], kwds["spin2z"])

Example #5

File: gatech_cat_spectra.py Project: astroclark/numrel_bursts

#

print "Extracting and generating waveform"

mass = simulations.simulations[sim_number]['Mmin30Hz']

# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Common params

mass1, mass2 = pnutils.mtotal_eta_to_mass1_mass2(mass,
        simulations.simulations[sim_number]['eta'])

# Estimate ffinal 
chi = pnutils.phenomb_chi(mass1, mass2,
        simulations.simulations[sim_number]['spin1z'],simulations.simulations[sim_number]['spin2z'])
ffinal = pnutils.get_final_freq(approx, mass1, mass2, 
        simulations.simulations[sim_number]['spin1z'],simulations.simulations[sim_number]['spin2z'])
#upp_bound = ffinal
#upp_bound = 1.5*ffinal
upp_bound = 0.5/delta_t

# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# NUMERICAL RELATIVITY

# --- Generate the polarisations
hplus_NR, hcross_NR = nrbu.get_wf_pols(
       simulations.simulations[sim_number]['wavefile'], mass, inclination=inc,
       delta_t=delta_t, f_lower=30.0001, distance=distance)


# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# APPROXIMANT 

Example #6

File: gatech_cat_matches.py Project: astroclark/numrel_bursts

matches = []

f, ax = pl.subplots(nrows=len(masses), ncols=2, figsize=(15, 15))

for m, mass in enumerate(masses):
    "Extracting and generating mass %d of %d (%.2f)" % (m, len(masses), mass)

    # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    # Common params
    #    mass = masses[-1]

    mass1, mass2 = pnutils.mtotal_eta_to_mass1_mass2(mass, simulations.simulations[0]["eta"])

    # Estimate ffinal
    ffinal = pnutils.get_final_freq(
        "SEOBNRv2", mass1, mass2, simulations.simulations[0]["spin1z"], simulations.simulations[0]["spin2z"]
    )

    # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    # NUMERICAL RELATIVITY

    if errors_file is None:

        # --- Generate the polarisations from hdf5
        hplus_NR, hcross_NR = nrbu.get_wf_pols(
            simulations.simulations[0]["wavefile"],
            mass,
            inclination=inc,
            delta_t=delta_t,
            f_lower=30.0001 * min(masses) / mass,
            distance=distance,

Example #7

File: waveform.py Project: gayathrigcc/pycbc

def seobnrrom_final_frequency(**kwds):
    from pycbc.pnutils import get_final_freq
    return get_final_freq("SEOBNRv2", kwds['mass1'], kwds['mass2'],
                   kwds['spin1z'], kwds['spin2z'])