for final_map_key in currentmaps.keys():
outdir = "Plots/Reference%s/NH/%s/"%(reference,init_map_key.split("_maps")[0])
currentmap = currentmaps[final_map_key]
bestmap = NH_vals[str(int(reference))][init_map_key][final_map_key]
currentmapobj = {}
currentmapobj['ebins'] = ebins
currentmapobj['czbins'] = czbins
currentmapobj['map'] = np.array(currentmap)
bestmapobj = {}
bestmapobj['ebins'] = ebins
bestmapobj['czbins'] = czbins
bestmapobj['map'] = np.array(bestmap)
diffmapobj = delta_map(currentmapobj, bestmapobj)
ratiomapobj = ratio_map(diffmapobj, bestmapobj)
ratiomax[init_map_key][final_map_key].append(np.amax(ratiomapobj['map']))
if init_map_key == 'nue_maps':
inittitle = r'$\nu_e$'
if init_map_key == 'nue_bar_maps':
inittitle = r'$\bar{\nu}_e$'
if init_map_key == 'numu_maps':
inittitle = r'$\nu_{\mu}$'
if init_map_key == 'numu_bar_maps':
inittitle = r'$\bar{\nu}_{\mu}$'
if final_map_key == 'nue':
finaltitle = r'$\nu_e$'
if final_map_key == 'nue_bar':
plt.suptitle(plottitle, fontsize=18)
plt.savefig('%s_honda_flux_maps.png'%(prim))
plt.close()
IPmapobj = {}
IPmapobj['ebins'] = IPfh[prim]['ebins']
IPmapobj['czbins'] = IPfh[prim]['czbins']
IPmapobj['map'] = np.array(IPfh[prim]['map'])
SImapobj = {}
SImapobj['ebins'] = SIfh[prim]['ebins']
SImapobj['czbins'] = SIfh[prim]['czbins']
SImapobj['map'] = np.array(SIfh[prim]['map'])
diffmapobj = delta_map(IPmapobj, SImapobj)
ratiomapobj = ratio_map(diffmapobj, SImapobj)
IPlogmapobj = {}
IPlogmapobj['ebins'] = IPfh[prim]['ebins']
IPlogmapobj['czbins'] = IPfh[prim]['czbins']
IPlogmapobj['map'] = np.log10(IPfh[prim]['map'])
SIlogmapobj = {}
SIlogmapobj['ebins'] = SIfh[prim]['ebins']
SIlogmapobj['czbins'] = SIfh[prim]['czbins']
SIlogmapobj['map'] = np.log10(SIfh[prim]['map'])
difflogmapobj = delta_map(IPlogmapobj, SIlogmapobj)
ratiologmapobj = ratio_map(difflogmapobj, SIlogmapobj)
for binning in cake_dict['binning']['dimensions']:
if binning['name'] == 'energy':
cake_map['ebins'] = binning['bin_edges']
if binning['name'] == 'coszen':
cake_map['czbins'] = binning['bin_edges']
if binning['name'] == 'true_energy':
cake_map['ebins'] = binning['bin_edges']
if binning['name'] == 'true_coszen':
cake_map['czbins'] = binning['bin_edges']
if binning['name'] == 'reco_energy':
cake_map['ebins'] = binning['bin_edges']
if binning['name'] == 'reco_coszen':
cake_map['czbins'] = binning['bin_edges']
RatioMapObj = ratio_map(cake_map, pisa_map)
DiffMapObj = delta_map(pisa_map, cake_map)
DiffRatioMapObj = ratio_map(DiffMapObj, pisa_map)
plt.figure(figsize = (20,5))
plt.subplot(1,5,1)
show_map(pisa_map)
plt.xlabel(r'$\cos\theta_Z$')
plt.ylabel(r'Energy [GeV]')
plt.title('%s %s PISA V2'%(titles[cake_dict['name']],args.stage))
plt.subplot(1,5,2)
show_map(cake_map)
plt.xlabel(r'$\cos\theta_Z$')
plt.ylabel(r'Energy [GeV]')
plt.title('%s %s PISA V3'%(titles[cake_dict['name']],args.stage))
for binning in cake_dict2["binning"]["dimensions"]:
if binning["name"] == "energy":
cake_map2["ebins"] = binning["bin_edges"]
if binning["name"] == "coszen":
cake_map2["czbins"] = binning["bin_edges"]
if binning["name"] == "true_energy":
cake_map2["ebins"] = binning["bin_edges"]
if binning["name"] == "true_coszen":
cake_map2["czbins"] = binning["bin_edges"]
if binning["name"] == "reco_energy":
cake_map2["ebins"] = binning["bin_edges"]
if binning["name"] == "reco_coszen":
cake_map2["czbins"] = binning["bin_edges"]
RatioMapObj = ratio_map(cake_map1, cake_map2)
DiffMapObj = delta_map(cake_map1, cake_map2)
DiffRatioMapObj = ratio_map(DiffMapObj, cake_map1)
plt.figure(figsize=(20, 5))
plt.subplot(1, 5, 1)
show_map(cake_map1)
plt.xlabel(r"$\cos\theta_Z$")
plt.ylabel(r"Energy [GeV]")
plt.title("%s %s PISA V3(1)" % (titles[cake_dict1["name"]], args.description1))
plt.subplot(1, 5, 2)
show_map(cake_map2)
plt.xlabel(r"$\cos\theta_Z$")
plt.ylabel(r"Energy [GeV]")
plt.title("%s %s PISA V3(2)" % (titles[cake_dict1["name"]], args.description2))