Skip to content
Snippets Groups Projects
Commit 9b5a8472 authored by chad.hanna's avatar chad.hanna
Browse files

gstlal_inspiral_mass_model: tweak max ligo mass

parent ae29b4e5
No related branches found
No related tags found
No related merge requests found
...@@ -72,8 +72,8 @@ prob = {} ...@@ -72,8 +72,8 @@ prob = {}
mchirps = {} mchirps = {}
ligo_min = 3 ligo_min = 3
ligo_max = 30 ligo_max = 40
ligo_peak = 30 ligo_peak = 40
ligo_alpha = -1 ligo_alpha = -1
ligonorm = schechter_norm(ligo_min, ligo_max, ligo_peak, ligo_alpha) ligonorm = schechter_norm(ligo_min, ligo_max, ligo_peak, ligo_alpha)
...@@ -150,14 +150,14 @@ coefficients = numpy.zeros((1, 1, max(ids)+1), dtype=float) ...@@ -150,14 +150,14 @@ coefficients = numpy.zeros((1, 1, max(ids)+1), dtype=float)
for tid in ids: for tid in ids:
coefficients[0,0,tid] = numpy.log(prob[tid]) - numpy.log(norm) coefficients[0,0,tid] = numpy.log(prob[tid]) - numpy.log(norm)
import matplotlib #import matplotlib
matplotlib.use('agg') #matplotlib.use('agg')
from matplotlib import pyplot #from matplotlib import pyplot
#pyplot.semilogx(chirpmasses, coefficients[0,0,:] + numpy.log((chirpmasses)**(15./6.)), "*") ##pyplot.semilogx(chirpmasses, coefficients[0,0,:] + numpy.log((chirpmasses)**(15./6.)), "*")
pyplot.semilogx(chirpmasses, coefficients[0,0,:], "*") #pyplot.semilogx(chirpmasses, coefficients[0,0,:], "*")
pyplot.grid() #pyplot.grid()
pyplot.ylim([-20, 0]) #pyplot.ylim([-20, 0])
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