From fb38b02b079cea92d24c0df357e8cc37881837b2 Mon Sep 17 00:00:00 2001
From: Chad Hanna <chad.hanna@ligo.org>
Date: Tue, 27 Nov 2018 08:29:52 -0800
Subject: [PATCH] gstlal_inspiral_mass_model: choose a minimum number of bins
---
gstlal-inspiral/bin/gstlal_inspiral_mass_model | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/gstlal-inspiral/bin/gstlal_inspiral_mass_model b/gstlal-inspiral/bin/gstlal_inspiral_mass_model
index 7ccb8a9a5d..c77c70b53c 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_mass_model
+++ b/gstlal-inspiral/bin/gstlal_inspiral_mass_model
@@ -55,7 +55,7 @@ sngl_inspiral_table = lsctables.SnglInspiralTable.get_table(xmldoc)
mchirps_of_tmps = chirpmass(sngl_inspiral_table.get_column("mass1"), sngl_inspiral_table.get_column("mass2"))
num_templates = len(mchirps_of_tmps)
-num_bins = num_templates / 500
+num_bins = max(5, num_templates / 500)
massBA = rate.BinnedDensity(rate.NDBins((rate.LogarithmicBins(min(mchirps_of_tmps)-1e-6, max(mchirps_of_tmps)+1e-6, num_bins),)))
for m in mchirps_of_tmps:
massBA.count[(m,)] += 1
--
GitLab