From 7d8e1113e6457da351ce4a5ad8ed92bc8221e625 Mon Sep 17 00:00:00 2001 From: Chad Hanna <chad.hanna@comp-hd-002.gwave.ics.psu.edu> Date: Sun, 8 Sep 2019 20:33:13 +0000 Subject: [PATCH] gstlal_inspiral_create_prior_diststats: make the signal model depend on bandwidth (cherry picked from commit 79dbac43d2e796ce0048a37d892362a21bb75b5a) --- gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats b/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats index b0e3350443..5630377c7c 100755 --- a/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats +++ b/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats @@ -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)) # -- GitLab