Commit 7d8e1113 authored by Chad Hanna's avatar Chad Hanna Committed by Chad Hanna

gstlal_inspiral_create_prior_diststats: make the signal model depend on bandwidth


(cherry picked from commit 79dbac43)
parent c4d47e58
Pipeline #78824 passed with stages
in 35 minutes and 25 seconds
......@@ -119,8 +119,7 @@ def parse_command_line():
bandwidths += [templates.bandwidth(row.mass1, row.mass2, row.spin1z, row.spin2z, f_min = 10.0, f_max = row.f_final, delta_f = 0.25, psd = psd[ifo])]
horizon_factors.update(bank.horizon_factors)
if options.df == "bandwidth":
# NOTE the 4000 is tuned by looking at real data distributions
options.df = int(4000. / min(bandwidths))
options.df = int(min(bandwidths))
return options, process_params, filenames, template_ids, horizon_factors
......@@ -160,10 +159,10 @@ process = ligolw_process.register_to_xmldoc(xmldoc, u"gstlal_inspiral_create_pri
rankingstat = far.RankingStat(template_ids = template_ids, instruments = options.instrument, delta_t = options.coincidence_threshold, min_instruments = options.min_instruments, population_model_file = options.mass_model_file, dtdphi_file = options.dtdphi_file, horizon_factors = horizon_factors)
if options.background_prior > 0:
rankingstat.denominator.add_noise_model(number_of_events = options.background_prior, df = int(options.df))
rankingstat.denominator.add_noise_model(number_of_events = options.background_prior)
# Add the numerator
rankingstat.numerator.add_signal_model()
rankingstat.numerator.add_signal_model(df = int(options.df))
#
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment