Skip to content
Snippets Groups Projects
Commit a4f455ef authored by Chad Hanna's avatar Chad Hanna
Browse files

gstlal_inspiral_mass_model: tweak

parent c8868acd
No related branches found
No related tags found
No related merge requests found
...@@ -67,7 +67,7 @@ tmplt_ids = [] ...@@ -67,7 +67,7 @@ tmplt_ids = []
mchirps = [] mchirps = []
probs = [] probs = []
bbhnorm = schechter_norm(5., 45, 45., -1.6) 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: for row in sngl_inspiral_table:
assert row.template_id not in ids assert row.template_id not in ids
...@@ -97,7 +97,7 @@ for row in sngl_inspiral_table: ...@@ -97,7 +97,7 @@ for row in sngl_inspiral_table:
# Other: made up # Other: made up
# #
otherprob = schechter(mchirp, 400., -2.35) / othernorm otherprob = schechter(mchirp, 200., -2.35) / othernorm
# #
# Combined # Combined
...@@ -118,13 +118,13 @@ coefficients = numpy.zeros((1, 1, max(ids)+1), dtype=float) ...@@ -118,13 +118,13 @@ coefficients = numpy.zeros((1, 1, max(ids)+1), dtype=float)
for tid in ids: for tid in ids:
coefficients[0,0,tid] = numpy.log(ids[tid]) coefficients[0,0,tid] = numpy.log(ids[tid])
#print coefficients.min() #print coefficients.max()
#import matplotlib #import matplotlib
#matplotlib.use('agg') #matplotlib.use('agg')
#from matplotlib import pyplot #from matplotlib import pyplot
#pyplot.loglog(mchirps, probs, "*") #pyplot.loglog(mchirps, probs, "*")
#pyplot.savefig("blah.png") #pyplot.savefig("blah.png")
#
# Write it out # Write it out
f = h5py.File(options.output, "w") f = h5py.File(options.output, "w")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment