flavor=nu_type,sample=sample_name) plot_bin_slices(true_val, cnn_val, energy_truth=cnn_Z, weights=true_weights, old_reco = retro_val,old_reco_truth=retro_true_val, reco_energy_truth=retro_Z, old_reco_weights=retro_weights,use_fraction = True, bins=syst_bin, min_val=-550, max_val=-150,\ save=save, savefolder=save_folder_name,\ variable=plot_name, units="(m)", reco_name=reco_name, xvariable="Reconstructed Z", flavor=nu_type,sample=sample_name) if variable == "energy": if plot_others: plot_2D_prediction_fraction(true_val, cnn_val,weights=true_weights,\ save=save, savefolder=save_folder_name,bins=bins,\ xminval=minval, xmaxval=maxval, yminval=-2,ymaxval=2.,\ variable=plot_name, units=plot_units, reco_name="CNN") plot_2D_prediction_fraction(retro_true_val, retro_val, weights=retro_weights, save=save, savefolder=save_folder_name,bins=bins,\ xminval=minval, xmaxval=maxval, yminval=-2,ymaxval=2.,\ variable=plot_name, units=plot_units, reco_name=reco_name) plot_2D_prediction_fraction(true_val, cnn_val,weights=true_weights,\ save=save, savefolder=save_folder_name,bins=bins,\ variable=plot_name, units=plot_units, reco_name="CNN") plot_2D_prediction_fraction(retro_true_val, retro_val, weights=retro_weights, save=save, savefolder=save_folder_name,bins=bins,\ variable=plot_name, units=plot_units, reco_name=reco_name) if plot_main: if nu_type == "NuMu" or nu_type == "numu": dist_title = r'$\nu_\mu$ '
plot_single_resolution(Y_test_use[:,true_index]*maxabs_factor,\ Y_test_predicted[:,NN_index]*maxabs_factor,\ use_old_reco = True, old_reco = reco_test_use[:,true_index],\ use_fraction=True,bins=bins,maxaxis=2.,minaxis=-2.,\ save=save,savefolder=save_folder_name,\ variable=plot_name,units=plot_units) plot_bin_slices(Y_test_use[:,true_index]*maxabs_factor, Y_test_predicted[:,NN_index]*maxabs_factor,\ old_reco = reco_test_use[:,true_index],\ use_fraction = use_frac,\ bins=10,min_val=minval,max_val=maxval,\ save=True,savefolder=save_folder_name,\ variable=plot_name,units=plot_units, epochs = epoch,reco_name=reco_name) if first_var == "energy" and num == 0: plot_2D_prediction_fraction(Y_test_use[:,true_index]*maxabs_factor,\ Y_test_predicted[:,NN_index]*maxabs_factor,\ save,save_folder_name,bins=bins,axis_square=True,\ minval=-2.,maxval=2,\ variable=plot_name,units=plot_units) plot_bin_slices(Y_test_use[:,true_index]*maxabs_factor, Y_test_predicted[:,NN_index]*maxabs_factor,\ use_fraction = use_frac,\ bins=10,min_val=minval,max_val=maxval,\ save=True,savefolder=save_folder_name,\ variable=plot_name,units=plot_units) #plot_length_energy(Y_test_use, Y_test_predicted[:,NN_index]*maxabs_factor,\ # use_fraction=True,ebins=20,tbins=20,\ # emin=minvals[first_var_index], emax=maxvals[first_var_index],\ # tmin=0.,tmax=450.,tfactor=maxabs_factors[2],\ # savefolder=save_folder_name) if num > 0 or first_var == "zenith": plot_bin_slices(Y_test_use[:,true_index], Y_test_predicted[:,NN_index], \ energy_truth=Y_test_use[:,0]*max_energy, \
variable=plot_name,units=plot_units, reco_name=reco_name) plot_single_resolution(Y_test_use[:,true_index]*maxabs_factor,\ reco_test_use[:,NN_index],\ minaxis=-maxval,maxaxis=maxval,bins=bins, save=save,savefolder=save_folder_name,\ variable=plot_name,units=plot_units) plot_bin_slices(Y_test_use[:,true_index]*maxabs_factor, reco_test_use[:,NN_index],\ use_fraction = False,\ bins=10,min_val=minval,max_val=maxval,\ save=True,savefolder=save_folder_name,\ variable=plot_name,units=plot_units) if num ==0: plot_2D_prediction_fraction(Y_test_use[:,true_index]*maxabs_factor,\ reco_test_use[:,NN_index],\ save,save_folder_name,bins=bins,\ minval=0,maxval=2,\ variable=plot_name,units=plot_units) plot_single_resolution(Y_test_use[:,true_index]*maxabs_factor,\ reco_test_use[:,NN_index]*maxabs_factor,\ use_fraction=True,\ minaxis=-2.,maxaxis=2., save=save,savefolder=save_folder_name,\ variable=plot_name,units=plot_units) plot_bin_slices(Y_test_use[:,true_index]*maxabs_factor, reco_test_use[:,NN_index],\ use_fraction = use_frac,\ bins=10,min_val=minval,max_val=maxval,\ save=True,savefolder=save_folder_name,\ variable=plot_name,units=plot_units) #plot_length_energy(Y_test_use, reco_test_use[:,NN_index]*maxabs_factor,\ # use_fraction=True,ebins=20,tbins=20,\