Skip to content
Snippets Groups Projects
Commit 79dbac43 authored by Chad Hanna's avatar Chad Hanna
Browse files

gstlal_inspiral_create_prior_diststats: make the signal model depend on bandwidth

parent 233217f7
No related branches found
No related tags found
No related merge requests found
......@@ -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))
#
......
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