diff --git a/gstlal-inspiral/bin/gstlal_inspiral_mass_model b/gstlal-inspiral/bin/gstlal_inspiral_mass_model index 3ddd3f0c9d0c23d624a42af78ab7cf0292570679..fc73d3561a7236d2f9b32939bd55cf86d42b8a18 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")