Skip to content
Snippets Groups Projects
Commit 58ff77ff authored by Wanting Niu's avatar Wanting Niu :speech_balloon:
Browse files

directly read spin1 and spin2 from sgnl_inspiral_table

parent 966ad8ca
No related branches found
No related tags found
1 merge request!369Document the template_duration in svd_bank match to the actual waveform
Pipeline #482686 failed
......@@ -276,9 +276,9 @@ for outputrows in banks_subbins:
row.template_duration = chirptime.imr_time(options.f_low,
lal.MSUN_SI * row.mass1,
lal.MSUN_SI * row.mass2,
numpy.dot(row.spin1x**2 + row.spin1y**2 + row.spin1z**2, row.spin1x**2 + row.spin1y**2 + row.spin1z**2)**.5,
numpy.dot(row.spin2x**2 + row.spin2y**2 + row.spin2z**2, row.spin2x**2 + row.spin2y**2 + row.spin2z**2)**.5,
f_max = min(row.f_final, 2 * chirptime.ringf(lal.MSUN_SI * row.mass1 + lal.MSUN_SI * row.mass2, chirptime.overestimate_j_from_chi(max(numpy.dot(row.spin1x**2 + row.spin1y**2 + row.spin1z**2, row.spin1x**2 + row.spin1y**2 + row.spin1z**2)**.5, numpy.dot(row.spin2x**2 + row.spin2y**2 + row.spin2z**2, row.spin2x**2 + row.spin2y**2 + row.spin2z**2)**.5))) if approximant in templates.gstlal_IMR_approximants else spawaveform.ffinal(row.mass1, row.mass2, 'bkl_isco')))
numpy.dot(row.spin1, row.spin1)**.5,
numpy.dot(row.spin2, row.spin2)**.5,
f_max = min(row.f_final, 2 * chirptime.ringf(lal.MSUN_SI * row.mass1 + lal.MSUN_SI * row.mass2, chirptime.overestimate_j_from_chi(max(numpy.dot(row.spin1, row.spin1)**.5, numpy.dot(row.spin2, row.spin2)**.5))) if approximant in templates.gstlal_IMR_approximants else spawaveform.ffinal(row.mass1, row.mass2, 'bkl_isco')))
# template_duration may be overwritten by svd_bank if the flow is increased in order to satisfy the target maximum duration
outputrows[ind] = (approximant, rows)
num_banks += len(outputrows)
......
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