Commit 3487481c authored by GstLAL CBC Offline's avatar GstLAL CBC Offline Committed by Chad Hanna make the snr power law so that it is same as exp(-snr**2/2) at SNR 8

(cherry picked from commit 2ec02a5d)
parent c39711d0
Pipeline #80454 passed with stages
in 32 minutes and 32 seconds
......@@ -857,7 +857,9 @@ class LnNoiseDensity(LnLRDensity):
prcoss = numpy.ones(len(rcoss))
# This adds a faint power law that falls off faster than GWs
psnr = snr**-12
# but not exponential (to avoid numerical errors). It has
# approximately the same value at SNR 8 as does exp(-8**2 / 2.)
psnr = snr**-15
psnr = numpy.outer(psnr, numpy.ones(len(rcoss)))
arr[snrindices, rcossindices] = psnr
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment