from argparse import ArgumentParser
from setup_pointlike import setup_pointlike
import yaml
parser = ArgumentParser()
parser.add_argument("-l", "--pwnlist", required=True, help="List of all yaml sources")
parser.add_argument("-n", "--name", required=True, help="Name of the pulsar")
args = parser.parse_args()
fit = lambda: roi.fit(method="minuit")
name = args.name
pwnlist = args.pwnlist
roi = setup_pointlike(name, pwnlist)
model = roi.get_model(which=name)
model.setp(1, 2)
model.freeze(1)
roi.modify(which=name, model=model)
roi.print_summary()
fit()
roi.print_summary()
ul = roi.upper_limit(which=args.name)
ts = roi.TS(which=name, quick=False, quiet=True)
from uw.like.roi_tsmap import TSCalc,TSCalcPySkyFunction
from argparse import ArgumentParser
from skymaps import SkyImage,SkyDir
import yaml
from uw.like.sed_plotter import plot_sed
parser = ArgumentParser()
parser.add_argument("-l", "--list", required=True, help="List of all yaml sources")
parser.add_argument("-n", "--name", required=True, help="Name of the pulsar")
args=parser.parse_args()
name=args.name
roi=setup_pointlike(name,args.list)
# First, calculate upper limit
fit=lambda: roi.fit(method='minuit')
source=roi.get_source(which=name)
pulsar_position=source.skydir
fit()
ts_at_pulsar = roi.TS(which=name,quick=False,quiet=True)
p,p_err=source.model.statistical(absolute=True)
flux_at_pulsar,flux_at_pulsar_err=source.model.i_flux(100,100000,error=True)
index_at_pulsar,index_at_pulsar_err=p[1],p_err[1]