extract_interacting_neutrinos.py File Reference

Go to the source code of this file.

Namespaces
namespace	extract_interacting_neutrinos

Functions
	extract_interacting_neutrinos.lookfornuhist (flukafile)

Variables
	extract_interacting_neutrinos.parser = argparse.ArgumentParser(description="Extract interacting neutrinos")
	extract_interacting_neutrinos.type
	extract_interacting_neutrinos.dest
	extract_interacting_neutrinos.default
	extract_interacting_neutrinos.help
	extract_interacting_neutrinos.options = parser.parse_args()
	extract_interacting_neutrinos.geniefile = r.TFile.Open(options.geniepath, "UPDATE")
	extract_interacting_neutrinos.flukafile = r.TFile.Open(options.flukapath, "READ")
	extract_interacting_neutrinos.hfluka_nuE = lookfornuhist(flukafile)
	extract_interacting_neutrinos.nbins = hfluka_nuE.GetNbinsX()
	extract_interacting_neutrinos.minE = hfluka_nuE.GetXaxis().GetBinLowEdge(1)
	extract_interacting_neutrinos.maxE = hfluka_nuE.GetXaxis().GetBinUpEdge(nbins)
	extract_interacting_neutrinos.df = r.RDataFrame("gst", geniefile)
	extract_interacting_neutrinos.hgenie_nuE
	extract_interacting_neutrinos.c1 = r.TCanvas()
	extract_interacting_neutrinos.c2 = r.TCanvas()
	extract_interacting_neutrinos.hnuratio
	extract_interacting_neutrinos.c3 = r.TCanvas()
	extract_interacting_neutrinos.maxratio = hnuratio.GetMaximum()
	extract_interacting_neutrinos.flukatree = flukafile.Get("t")
	extract_interacting_neutrinos.output_nutree = flukatree.CloneTree(0)
	extract_interacting_neutrinos.nneutrinos = flukatree.GetEntries()
	extract_interacting_neutrinos.Ekin = nuevent.Ekin
	extract_interacting_neutrinos.randomratio = r.gRandom.Uniform(0, maxratio)
	extract_interacting_neutrinos.matchingratio = hnuratio.GetBinContent(hnuratio.FindBin(Ekin))