from argparse import ArgumentParser import yaml import numpy as np from lande_toolbag import tolist from snr_setup import get_snr, setup_roi from snr_helper import pointlike_analysis, gtlike_analysis, plots from likelihood_tools import force_gradient force_gradient(use_gradient=False) parser = ArgumentParser() parser.add_argument("--superfile", required=True) parser.add_argument("--name", required=True) parser.add_argument("--emin", type=float, default=1e4) parser.add_argument("--emax", type=float, default=1e5) parser.add_argument("--no-plots", action='store_true', default=False) parser.add_argument("--no-gtlike", action='store_true', default=False) args=parser.parse_args() do_plots = not args.no_plots do_gtlike = not args.no_gtlike name=args.name superfile=args.superfile
parser.add_argument("--max-free", default=10, type=float) args=parser.parse_args() initial=args.initialization do_at_tev = not args.no_at_tev do_point = not args.no_point do_extended = not args.no_extended do_gtlike = not args.no_gtlike do_plots = not args.no_plots do_upper_limits = not args.no_upper_limits do_seds = not args.no_seds force_gradient(use_gradient=args.use_gradient) name=args.name emin=args.emin emax=args.emax ft1 = yaml.load(open(args.tevsources))[name]['ft1'] ft2 = yaml.load(open(args.tevsources))[name]['ft2'] ltcube = yaml.load(open(args.tevsources))[name]['ltcube'] from modify_tev import * def get_roi(**kwargs): print 'Building the ROI' roi=setup_tev(name,args.tevsources, ft1, ft2, ltcube, free_radius=5, max_free=args.max_free, fit_emin=emin, fit_emax=emax,