예제 #1
0
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,