def gtlike_analysis(roi, hypothesis, upper_limit=False, cutoff=False):
print 'Performing Gtlike crosscheck for %s' % hypothesis
gtlike=Gtlike(roi)
like=gtlike.like
like.fit(covar=True)
r=sourcedict(like, name)
if upper_limit:
r['upper_limit'] = powerlaw_upper_limit(like, name, emin=emin, emax=emax, cl=.95)
if cutoff:
r['test_cutoff']=test_cutoff(like,name)
for kind, kwargs in [['4bpd',dict(bin_edges=np.logspace(2,5,13))],
['1bpd',dict(bin_edges=np.logspace(2,5,4))]]:
print 'Making %s SED' % kind
sed = SED(like, name, **kwargs)
sed.plot('sed_gtlike_%s_%s_%s.png' % (kind,hypothesis,name))
sed.verbosity=True
sed.save('sed_gtlike_%s_%s_%s.dat' % (kind,hypothesis,name))
return r
def gtlike_analysis(roi, hypothesis, upper_limit=False, cutoff=False):
print "Performing Gtlike crosscheck for %s" % hypothesis
gtlike = Gtlike(roi)
like = gtlike.like
like.fit(covar=True)
r = sourcedict(like, name)
if upper_limit:
r["upper_limit"] = powerlaw_upper_limit(like, name, emin=emin, emax=emax, cl=0.95)
if cutoff:
r["test_cutoff"] = test_cutoff(like, name)
for kind, kwargs in [
["4bpd", dict(bin_edges=np.logspace(2, 5, 13))],
["1bpd", dict(bin_edges=np.logspace(2, 5, 4))],
]:
print "Making %s SED" % kind
sed = SED(like, name, **kwargs)
sed.plot("sed_gtlike_%s_%s_%s.png" % (kind, hypothesis, name))
sed.verbosity = True
sed.save("sed_gtlike_%s_%s_%s.dat" % (kind, hypothesis, name))
return r
