Skip to content
Snippets Groups Projects
Commit 3baf8672 authored by Patrick Godwin's avatar Patrick Godwin
Browse files

inspiral_lr.py: enable iDQ only for singles

parent 93b19a9d
No related branches found
No related merge requests found
Pipeline #155435 failed
......@@ -394,14 +394,16 @@ class LnSignalDensity(LnLRDensity):
# instruments)
interp = self.interps["snr_chi"]
# Evaluate the IDQ glitch probability # FIXME put in denominator
for ifo, seg in segments.items():
# only proceed if we have a trigger from this ifo
if ifo in snrs:
# FIXME don't just use the last segment, somehow include whole template duration?
t = float(seg[1])
# NOTE choose max over +-1 seconds because the sampling is only at 1 Hz.
lnP -= max(self.idq_glitch_lnl[ifo]([t-1., t, t+1.]))
# Evaluate the IDQ glitch probability (singles only)
# FIXME put in denominator
if len(snrs) == 1:
for ifo, seg in segments.items():
# only proceed if we have a trigger from this ifo
if ifo in snrs:
# FIXME don't just use the last segment, somehow include whole template duration?
t = float(seg[1])
# NOTE choose max over +-1 seconds because the sampling is only at 1 Hz.
lnP -= max(self.idq_glitch_lnl[ifo]([t-1., t, t+1.]))
return lnP + sum(interp(snrs[instrument], chi2_over_snr2) for instrument, chi2_over_snr2 in chi2s_over_snr2s.items())
......
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