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"]