diff --git a/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats b/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats index b0e3350443caa30e0cef44a0a4aa39632d9a2008..5630377c7cbda7d5709105837d28f8ad70c6b773 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)) #