From 8f8d7a5c28c74e70e439518087a279d31cf65fa7 Mon Sep 17 00:00:00 2001 From: Kipp Cannon <kcannon@cita.utoronto.ca> Date: Mon, 4 Dec 2017 02:09:22 +0900 Subject: [PATCH] inspiral_lr: edit a comment --- gstlal-inspiral/python/stats/inspiral_lr.py | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) diff --git a/gstlal-inspiral/python/stats/inspiral_lr.py b/gstlal-inspiral/python/stats/inspiral_lr.py index cc0f58259e..33870a6ec8 100644 --- a/gstlal-inspiral/python/stats/inspiral_lr.py +++ b/gstlal-inspiral/python/stats/inspiral_lr.py @@ -563,15 +563,15 @@ class LnNoiseDensity(LnLRDensity): # sanity check rates assert all(triggers_per_second_per_template[instrument] for instrument in snrs), "impossible candidate in %s at %s when rates were %s triggers/s/template" % (", ".join(sorted(snrs)), ", ".join("%s s in %s" % (str(seg[1]), instrument) for instrument, seg in sorted(segments.items())), str(triggers_per_second_per_template)) - # P(t | noise) = (candidates per unit time @ t) / livetime. - # by not normalizing by the livetime the return value can - # only ever be proportional to the probability density, but - # we avoid the problem of the ranking statistic definition - # changing on-the-fly while running online, allowing - # candidates collected later to have their ranking - # statistics compared meaningfully to the values assigned - # to candidates collected earlier, when the total livetime - # was smaller. + # P(t | noise) = (candidates per unit time @ t) / total + # candidates. by not normalizing by the total candidates + # the return value can only ever be proportional to the + # probability density, but we avoid the problem of the + # ranking statistic definition changing on-the-fly while + # running online, allowing candidates collected later to + # have their ranking statistics compared meaningfully to + # the values assigned to candidates collected earlier, when + # the total number of candidates was smaller. lnP = math.log(sum(self.coinc_rates.strict_coinc_rates(**triggers_per_second_per_template).values()) * self.num_templates) # P(instruments | t, noise) -- GitLab