def medianLabel(config, emin, emax):
if config[:2] == 'IC':
median, sigL, sigR = readDist_IC(config, emin, emax)
if config[:2] == 'IT':
median, sigL, sigR = readDist_IT(config, emin, emax)
median = 10**(9+median)
median = round(median, -int(math.floor(math.log10(median))) + 1)
if median >= 1e15:
return '%.1fPeV' % (median / 1e15)
if median >= 1e12:
if median >= 1e13:
return '%iTeV' % (median / 1e12)
return '%.1fTeV' % (median / 1e12)
configs, files = configs[csort], np.asarray(files)[csort]
x = range(len(configs))
smooth = 5
if args.files == "energy":
files = getEnergyMaps("IC")
ebins = [float(i) for i in getEbins()]
nbins = len(ebins) - 1
x = [readDist_IC("IC", ebins[i], ebins[i + 1])[0] for i in range(nbins)]
xL = [readDist_IC("IC", ebins[i], ebins[i + 1])[1] for i in range(nbins)]
xR = [readDist_IC("IC", ebins[i], ebins[i + 1])[2] for i in range(nbins)]
# x = getMids(ebins, infvalue=100)
smooth = 20
# Custom treatment of IT map
files += [["{}/IT_24H_sid_STA8.fits".format(my.ani_maps)]]
# x = np.append(x, np.log10(2*10**6))
median, sigL, sigR = readDist_IT("IT", 8, 100)
x += [median]
xL += [sigL]
xR += [sigR]
if args.files == "comp":
files = glob.glob("%s/IT_24H_sid_STA8_?.fits" % my.ani_maps)
comps = [f[-6] for f in files]
# comp2e = {'p':1.,'h':2.,'o':8.,'f':26.}
comp2e = {"p": 1.0, "h": 4.0, "o": 16.0, "f": 56.0}
x = [1.0 / comp2e[comp] for comp in comps]
smooth = 20
if args.files == "solar":
files = glob.glob("%s/IC_24H_solar.fits" % my.ani_maps)
smooth = 5
x = ["solar"]