diff --git a/gstlal-inspiral/bin/gstlal_inspiral_injection_snr b/gstlal-inspiral/bin/gstlal_inspiral_injection_snr index 5b6cefd23c4d7cbc8cc5b0e2652dae6c8a76a888..6ba29c86b85fad0fc7b83f8da2c196b8dd95b8fb 100755 --- a/gstlal-inspiral/bin/gstlal_inspiral_injection_snr +++ b/gstlal-inspiral/bin/gstlal_inspiral_injection_snr @@ -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