From 57cfbe528c53d79bfbb9606fda539b65b0958f92 Mon Sep 17 00:00:00 2001
From: Duncan Meacher <duncan.meacher@ligo.org>
Date: Sun, 15 Jul 2018 13:44:35 -0700
Subject: [PATCH] Lowered single det SNR threshold
---
gstlal-inspiral/python/lloidparts.py | 5 ++++-
gstlal-inspiral/python/stats/inspiral_extrinsics.py | 2 +-
gstlal-inspiral/python/stats/inspiral_lr.py | 2 +-
3 files changed, 6 insertions(+), 3 deletions(-)
diff --git a/gstlal-inspiral/python/lloidparts.py b/gstlal-inspiral/python/lloidparts.py
index f16a7c00db..ba3f91d0e3 100644
--- a/gstlal-inspiral/python/lloidparts.py
+++ b/gstlal-inspiral/python/lloidparts.py
@@ -1107,7 +1107,10 @@ def mkLLOIDmulti(pipeline, detectors, banks, psd, psd_fft_length = 32, ht_gate_t
if chisq_type == 'autochisq':
# FIXME don't hardcode
# peak finding window (n) in samples is 1 second at max rate, ie max(rates)
- head = pipeparts.mkitac(pipeline, snr, 1 * max(bank.get_rates()), bank.template_bank_filename, autocorrelation_matrix = bank.autocorrelation_bank, mask_matrix = bank.autocorrelation_mask, snr_thresh = bank.snr_threshold, sigmasq = bank.sigmasq)
+ if instrument == 'H1' or instrument == 'L1':
+ head = pipeparts.mkitac(pipeline, snr, 1 * max(bank.get_rates()), bank.template_bank_filename, autocorrelation_matrix = bank.autocorrelation_bank, mask_matrix = bank.autocorrelation_mask, snr_thresh = 4, sigmasq = bank.sigmasq)
+ else:
+ head = pipeparts.mkitac(pipeline, snr, 1 * max(bank.get_rates()), bank.template_bank_filename, autocorrelation_matrix = bank.autocorrelation_bank, mask_matrix = bank.autocorrelation_mask, snr_thresh = bank.snr_threshold, sigmasq = bank.sigmasq)
if verbose:
head = pipeparts.mkprogressreport(pipeline, head, "progress_xml_%s" % suffix)
triggersrcs[instrument].add(head)
diff --git a/gstlal-inspiral/python/stats/inspiral_extrinsics.py b/gstlal-inspiral/python/stats/inspiral_extrinsics.py
index a39e460464..14d68e0bdf 100644
--- a/gstlal-inspiral/python/stats/inspiral_extrinsics.py
+++ b/gstlal-inspiral/python/stats/inspiral_extrinsics.py
@@ -901,7 +901,7 @@ class NumeratorSNRCHIPDF(rate.BinnedLnPDF):
self.norm.shape = (len(self.norm), 1)
@staticmethod
- def add_signal_model(lnpdf, n, prefactors_range, df, inv_snr_pow = 4., snr_min = 4., progressbar = None):
+ def add_signal_model(lnpdf, n, prefactors_range, df, inv_snr_pow = 4., snr_min = 3.5, progressbar = None):
if df <= 0.:
raise ValueError("require df >= 0: %s" % repr(df))
pfs = numpy.linspace(prefactors_range[0], prefactors_range[1], 100)
diff --git a/gstlal-inspiral/python/stats/inspiral_lr.py b/gstlal-inspiral/python/stats/inspiral_lr.py
index f4642f1de9..fadd518136 100644
--- a/gstlal-inspiral/python/stats/inspiral_lr.py
+++ b/gstlal-inspiral/python/stats/inspiral_lr.py
@@ -76,7 +76,7 @@ __all__ = [
class LnLRDensity(snglcoinc.LnLRDensity):
# range of SNRs covered by this object
- snr_min = 4.
+ snr_min = 3.5
# SNR, \chi^2 binning definition
snr_chi_binning = rate.NDBins((rate.ATanLogarithmicBins(2.6, 26., 300), rate.ATanLogarithmicBins(.001, 0.2, 280)))
--
GitLab