diff --git a/gstlal-inspiral/python/stats/inspiral_lr.py b/gstlal-inspiral/python/stats/inspiral_lr.py
index b6b2739f791a32d4da4bd19b7739855ac2024713..90ac676333d133553ad867ce1073b46106e55a98 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