Skip to content
Snippets Groups Projects
Commit 9b5c8a84 authored by Chad Hanna's avatar Chad Hanna
Browse files

inspiral_lr.py: add missing snr power law term

parent 255b559d
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
gstlal-inspiral/python/stats/inspiral_lr.py
+ 2
− 1
View file @ 9b5c8a84
...@@ -813,8 +813,9 @@ class LnNoiseDensity(LnLRDensity): ...@@ -813,8 +813,9 @@ class LnNoiseDensity(LnLRDensity):
prcoss = numpy.ones(len(rcoss)) prcoss = numpy.ones(len(rcoss))
psnr = 1e-8 * snr**-6 #(1. + 10**6) / (1. + snr**6) 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) 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 # normalize to the requested count. give 99% of the
# requested events to this portion of the model # requested events to this portion of the model
......
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