Skip to content
Snippets Groups Projects
Commit fc0839f9 authored by Ryan Lang's avatar Ryan Lang
Browse files

gstlal_inspiral_injection_snr: Added functionality for NR and IMBH

parent ef1b6b6a
No related branches found
No related tags found
No related merge requests found
......@@ -93,6 +93,16 @@ def calc_expected_snr(inj):
# FIXME have better scheme for calculating the needed sample_rate
sample_rate = 16384.0
approximant = lalsimulation.GetApproximantFromString(str(inj.waveform))
if approximant == lalsimulation.NR_hdf5:
LALparams = lal.CreateDict()
lalsimulation.SimInspiralWaveformParamsInsertNumRelData(LALparams, str(inj.numrel_data))
f_min = inj.f_lower
else:
LALparams = None
f_min = options.flow
h_plus, h_cross = lalsimulation.SimInspiralTD(
m1 = inj.mass1*lal.MSUN_SI,
m2 = inj.mass2*lal.MSUN_SI,
......@@ -109,10 +119,10 @@ def calc_expected_snr(inj):
eccentricity = 0.0,
meanPerAno = 0.0,
deltaT = 1.0 / sample_rate,
f_min = options.flow,
f_min = f_min,
f_ref = 0.0,
LALparams = None,
approximant = lalsimulation.GetApproximantFromString(str(inj.waveform))
LALparams = LALparams,
approximant = approximant
)
h_plus.epoch += injtime
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment