i = args.i istr = "%05d" % i np.random.seed(i) difftype = args.difftype time = args.time position = args.position emin = args.emin emax = args.emax diffdir = "/afs/slac/g/glast/groups/diffuse/rings/2year/" iso = get_background(join(diffdir, "isotrop_2year_P76_source_v0.txt")) ring = get_background(join(diffdir, "ring_2year_P76_v0.fits")) sreekumar = get_sreekumar(diff_factor=1) if difftype == "galactic": for p in iso.smodel.param_names: iso.smodel.freeze(p) ring.smodel = Constant() diffuse_sources = [iso, ring] elif difftype == "isotropic": diffuse_sources = [iso] elif difftype == "sreekumar": diffuse_sources = [sreekumar] tempdir = mkdtemp(prefix="/scratch/") if position == "highlat": while True:
np.random.shuffle(extensions) # formula to interpolate from the flux at lowest to highest extension flux_mc = lambda extension: np.exp(np.log(min_flux_mc) + (np.log(max_flux_mc) - np.log(min_flux_mc))*\ (extension-min(extensions))/(max(extensions)-min(extensions))) index_mc = args.index emin=1e2 emax=1e5 irf="P7SOURCE_V6" skydir_mc = SkyDir() bg = get_sreekumar() ft2 = dict2fgl['ft2'] ltcube = dict2fgl['ltcube'] results = [] for extension_mc in extensions: print 'Looping over extension_mc=%g' % extension_mc model_mc = PowerLaw(index=index_mc) model_mc.set_flux(flux_mc(extension_mc), emin, emax) r = dict( type = args.type,
if source == 'W44': skydir = SkyDir(284.005,1.345) gtifile = '/u/gl/funk/data/GLAST/ExtendedSources/NewAnalysis/gtlike/W44/sel_W44_60_100000.fits' elif source == 'IC443': skydir = SkyDir(94.310,22.580) gtifile = '/u/gl/funk/data/GLAST/ExtendedSources/NewAnalysis/gtlike/IC443/sel_IC443_60_100000.fits' if args.diffuse == 'galactic': ds = get_default_diffuse( diffdir="/nfs/slac/g/ki/ki03/lande/fermi/diffuse/", gfile="gal_2yearp7v6_v0.fits", ifile="iso_p7v6source.txt") elif args.diffuse == 'sreekumar': ds = [ get_sreekumar() ] elif args.diffuse == 'nobackground': ds = [] elif args.diffuse == 'extrapolated': ds = get_background( "/afs/slac/g/glast/groups/scienceTools/bugs/Likelihood/energy_dispersion/Simulations/galdif_extended.fits", '/nfs/slac/g/ki/ki03/lande/fermi/diffuse/iso_p7v6source_extrapolated.txt') else: raise Exception("...") ps,ds = catalog.merge_lists(skydir, radius=15, user_diffuse_list=ds) d = {s.name:s for s in ds} if source == 'W44':