From 9b5c8a8498f408fd754c375c9698f6a986687054 Mon Sep 17 00:00:00 2001 From: Chad Hanna <chad.hanna@ligo.org> Date: Thu, 6 Jun 2019 13:28:55 -0400 Subject: [PATCH] inspiral_lr.py: add missing snr power law term --- gstlal-inspiral/python/stats/inspiral_lr.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/gstlal-inspiral/python/stats/inspiral_lr.py b/gstlal-inspiral/python/stats/inspiral_lr.py index b6b2739f79..90ac676333 100644 --- a/gstlal-inspiral/python/stats/inspiral_lr.py +++ b/gstlal-inspiral/python/stats/inspiral_lr.py @@ -813,8 +813,9 @@ class LnNoiseDensity(LnLRDensity): prcoss = numpy.ones(len(rcoss)) psnr = 1e-8 * snr**-6 #(1. + 10**6) / (1. + snr**6) + psnr = numpy.outer(psnr, numpy.ones(len(rcoss))) psnrdcoss = numpy.outer(numpy.exp(-(snr - 2**.5)**2/ 2.) * dsnr, numpy.exp(-(rcoss - .05)**2 / .00015*2) * drcoss) - arr[snrindices, rcossindices] = psnrdcoss + arr[snrindices, rcossindices] = psnrdcoss + psnr # normalize to the requested count. give 99% of the # requested events to this portion of the model -- GitLab