def gtlike_analysis(roi, name, plotdir):
gtlike = Gtlike(roi)
like = gtlike
print_summary = lambda: print_summary(like, maxdist=10)
print_summary()
paranoid_gtlike_fit(like)
print_summary()
results = sourcedict(like, name)
results['powerlaw_upper_limit'] = powerlaw_upper_limit(
like,
name,
powerlaw_index=2.1,
emin=emin,
emax=emax,
cl=.95,
delta_log_like_limits=10)
sed = GtlikeSED(like, name, always_upper_limit=True)
sed.plot('%s/sed_gtlike_%s.png' % (plotdir, name))
sed.save('%s/sed_gtlike_%s.yaml' % (plotdir, name))
return results
def pointlike_analysis(roi, name, plotdir):
print_summary = lambda: roi.print_summary(galactic=True, maxdist=10)
print_summary()
fit_prefactor(roi, name)
roi.fit()
print_summary()
results = sourcedict(roi, name)
results['powerlaw_upper_limit'] = powerlaw_upper_limit(roi, name, powerlaw_index=2.1,
emin=emin, emax=emax, cl=.95)
roi.plot_sed(which=name, filename='%s/sed_pointlike_%s.png' % (plotdir,name), use_ergs=True)
return results
def gtlike_analysis(roi, name, plotdir):
gtlike = Gtlike(roi)
like = gtlike
print_summary = lambda: print_summary(like, maxdist=10)
print_summary()
paranoid_gtlike_fit(like)
print_summary()
results=sourcedict(like, name)
results['powerlaw_upper_limit'] = powerlaw_upper_limit(like, name, powerlaw_index=2.1,
emin=emin, emax=emax, cl=.95, delta_log_like_limits=10)
sed = GtlikeSED(like, name, always_upper_limit=True)
sed.plot('%s/sed_gtlike_%s.png' % (plotdir,name))
sed.save('%s/sed_gtlike_%s.yaml' % (plotdir,name))
return results
def pointlike_analysis(roi, name, plotdir):
print_summary = lambda: roi.print_summary(galactic=True, maxdist=10)
print_summary()
fit_prefactor(roi, name)
roi.fit()
print_summary()
results = sourcedict(roi, name)
results['powerlaw_upper_limit'] = powerlaw_upper_limit(roi,
name,
powerlaw_index=2.1,
emin=emin,
emax=emax,
cl=.95)
roi.plot_sed(which=name,
filename='%s/sed_pointlike_%s.png' % (plotdir, name),
use_ergs=True)
return results
