From 8d438d4a7d5d9820c711c18517d0702f94ca09f9 Mon Sep 17 00:00:00 2001
From: Kipp Cannon <kipp.cannon@ligo.org>
Date: Sun, 26 Nov 2017 22:08:29 -0600
Subject: [PATCH] inspiral_lr: disable P(t) in numerator
- mucks things up somehow, don't understand it
---
gstlal-inspiral/python/stats/inspiral_lr.py | 11 +++++++----
1 file changed, 7 insertions(+), 4 deletions(-)
diff --git a/gstlal-inspiral/python/stats/inspiral_lr.py b/gstlal-inspiral/python/stats/inspiral_lr.py
index b566a8e625..cc0f58259e 100644
--- a/gstlal-inspiral/python/stats/inspiral_lr.py
+++ b/gstlal-inspiral/python/stats/inspiral_lr.py
@@ -269,10 +269,11 @@ class LnSignalDensity(LnLRDensity):
# overall scale: ln L = 0 is not a special value, as it
# would be if the numerator and denominator were both
# normalized properly.
- horizon = sorted(horizons.values())[-self.min_instruments] / 150.
- if not horizon:
- return float("-inf")
- lnP += 3. * math.log(sorted(horizons.values())[-self.min_instruments] / 150.)
+ # FIXME: disabled for now
+ #horizon = sorted(horizons.values())[-self.min_instruments] / 150.
+ #if not horizon:
+ # return float("-inf")
+ #lnP += 3. * math.log(sorted(horizons.values())[-self.min_instruments] / 150.)
# evaluate dt and dphi parameters
lnP += inspiral_extrinsics.lnP_dt_dphi_signal(snrs, phase, dt, horizons, self.delta_t)
@@ -320,6 +321,8 @@ class LnSignalDensity(LnLRDensity):
parameter sets the nominal signal rate in units of Gpc^-3
a^-1.
"""
+ # FIXME: disabled for now
+ return 1.
# FIXME: this needs to understand a mass distribution
# model and what part of the mass space this numerator PDF
# is for
--
GitLab