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))