From a4f455ef3e15b85b0f276b1096c6c21e990732ab Mon Sep 17 00:00:00 2001
From: Chad Hanna <chad.hanna@ligo.org>
Date: Fri, 23 Nov 2018 17:42:55 -0800
Subject: [PATCH] gstlal_inspiral_mass_model: tweak
---
gstlal-inspiral/bin/gstlal_inspiral_mass_model | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/gstlal-inspiral/bin/gstlal_inspiral_mass_model b/gstlal-inspiral/bin/gstlal_inspiral_mass_model
index 3ddd3f0c9d..fc73d3561a 100755
--- a/gstlal-inspiral/bin/gstlal_inspiral_mass_model
+++ b/gstlal-inspiral/bin/gstlal_inspiral_mass_model
@@ -67,7 +67,7 @@ tmplt_ids = []
mchirps = []
probs = []
bbhnorm = schechter_norm(5., 45, 45., -1.6)
-othernorm = schechter_norm(min(mchirps_of_tmps), max(mchirps_of_tmps), 400., -2.3)
+othernorm = schechter_norm(1, 200, 200., -2.35)
for row in sngl_inspiral_table:
assert row.template_id not in ids
@@ -97,7 +97,7 @@ for row in sngl_inspiral_table:
# Other: made up
#
- otherprob = schechter(mchirp, 400., -2.35) / othernorm
+ otherprob = schechter(mchirp, 200., -2.35) / othernorm
#
# Combined
@@ -118,13 +118,13 @@ coefficients = numpy.zeros((1, 1, max(ids)+1), dtype=float)
for tid in ids:
coefficients[0,0,tid] = numpy.log(ids[tid])
-#print coefficients.min()
+#print coefficients.max()
#import matplotlib
#matplotlib.use('agg')
#from matplotlib import pyplot
#pyplot.loglog(mchirps, probs, "*")
#pyplot.savefig("blah.png")
-
+#
# Write it out
f = h5py.File(options.output, "w")
--
GitLab