From 255b559d7221299433a64e219e900ffe615321b1 Mon Sep 17 00:00:00 2001
From: Chad Hanna <chad.hanna@ligo.org>
Date: Thu, 6 Jun 2019 13:27:57 -0400
Subject: [PATCH] gstlal_inspiral_create_prior_diststats: allow user to specify
df
---
gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats b/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats
index 2363e5bf52..e6925009d3 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats
+++ b/gstlal-inspiral/bin/gstlal_inspiral_create_prior_diststats
@@ -75,6 +75,7 @@ def parse_command_line():
parser.add_option("--write-likelihood", metavar = "filename", help = "Write merged raw likelihood data to this file.")
parser.add_option("--instrument", action = "append", help = "Append to a list of instruments to create dist stats for. List must be whatever instruments you intend to analyze.")
parser.add_option("-p", "--background-prior", metavar = "N", default = 1, type = "float", help = "Include an exponential background prior with the maximum bin count = N, default 1")
+ parser.add_option("--df", metavar = "N", default = 40, type = "int", help = "set the degrees of freedom for the background chisq prior: default 40")
parser.add_option("--svd-file", metavar = "filename", help = "The SVD file to read the template ids from")
parser.add_option("--mass-model-file", metavar = "filename", help = "The mass model file to read from (hdf5 format)")
options, filenames = parser.parse_args()
@@ -137,7 +138,7 @@ 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, horizon_factors = horizon_factors)
if options.background_prior > 0:
- rankingstat.denominator.add_noise_model(number_of_events = options.background_prior)
+ rankingstat.denominator.add_noise_model(number_of_events = options.background_prior, df = options.df)
if options.synthesize_numerator:
rankingstat.numerator.add_signal_model()
--
GitLab