Functions
	lookfornuhist (flukafile)

Variables
	parser = argparse.ArgumentParser(description="Extract interacting neutrinos")

	type

	dest

	default

	help

	options = parser.parse_args()

	geniefile = r.TFile.Open(options.geniepath, "UPDATE")

	flukafile = r.TFile.Open(options.flukapath, "READ")

	hfluka_nuE = lookfornuhist(flukafile)

	nbins = hfluka_nuE.GetNbinsX()

	minE = hfluka_nuE.GetXaxis().GetBinLowEdge(1)

	maxE = hfluka_nuE.GetXaxis().GetBinUpEdge(nbins)

	df = r.RDataFrame("gst", geniefile)

	hgenie_nuE

	c1 = r.TCanvas()

	c2 = r.TCanvas()

	hnuratio

	c3 = r.TCanvas()

	maxratio = hnuratio.GetMaximum()

	flukatree = flukafile.Get("t")

	output_nutree = flukatree.CloneTree(0)

	nneutrinos = flukatree.GetEntries()

	Ekin = nuevent.Ekin

	randomratio = r.gRandom.Uniform(0, maxratio)

	matchingratio = hnuratio.GetBinContent(hnuratio.FindBin(Ekin))

Detailed Description

from the neutrinos produced by the FLUKA simulation, extract a subsample of interacting according to the GENIE cross section

Function Documentation

◆ lookfornuhist()

extract_interacting_neutrinos.lookfornuhist ( flukafile )

the actual name of the histograms depend on the neutrino we are reading,
i loop over all possible candidates

Definition at line 30 of file extract_interacting_neutrinos.py.

def lookfornuhist(flukafile):
    """the actual name of the histograms depend on the neutrino we are reading,
    i loop over all possible candidates"""
    nup_histnames = ["1012", "1014", "1016", "2012", "2014", "2016"]
    myhist = r.TH1F()
    for name in nup_histnames:
        temphist = flukafile.Get(name)
        if temphist:
            myhist = temphist
    return myhist
 
 

Variable Documentation

◆ c1

extract_interacting_neutrinos.c1 = r.TCanvas()

Definition at line 62 of file extract_interacting_neutrinos.py.

◆ c2

extract_interacting_neutrinos.c2 = r.TCanvas()

Definition at line 65 of file extract_interacting_neutrinos.py.

◆ c3

extract_interacting_neutrinos.c3 = r.TCanvas()

Definition at line 84 of file extract_interacting_neutrinos.py.

◆ default

extract_interacting_neutrinos.default

Definition at line 11 of file extract_interacting_neutrinos.py.

◆ dest

extract_interacting_neutrinos.dest

Definition at line 10 of file extract_interacting_neutrinos.py.

◆ df

extract_interacting_neutrinos.df = r.RDataFrame("gst", geniefile)

Definition at line 48 of file extract_interacting_neutrinos.py.

◆ Ekin

extract_interacting_neutrinos.Ekin = nuevent.Ekin

Definition at line 100 of file extract_interacting_neutrinos.py.

◆ flukafile

extract_interacting_neutrinos.flukafile = r.TFile.Open(options.flukapath, "READ")

Definition at line 26 of file extract_interacting_neutrinos.py.

◆ flukatree

extract_interacting_neutrinos.flukatree = flukafile.Get("t")

Definition at line 91 of file extract_interacting_neutrinos.py.

◆ geniefile

extract_interacting_neutrinos.geniefile = r.TFile.Open(options.geniepath, "UPDATE")

Definition at line 25 of file extract_interacting_neutrinos.py.

◆ help

extract_interacting_neutrinos.help

Definition at line 12 of file extract_interacting_neutrinos.py.

◆ hfluka_nuE

extract_interacting_neutrinos.hfluka_nuE = lookfornuhist(flukafile)

Definition at line 42 of file extract_interacting_neutrinos.py.

◆ hgenie_nuE

extract_interacting_neutrinos.hgenie_nuE

Initial value:

=  df.Histo1D(
    (
        "hgenie_nuE",
        "Energy of muon neutrinos from genie spectrum;E[GeV]",
        nbins,
        minE,
        maxE,
    ),
    "pzv",
)

Definition at line 49 of file extract_interacting_neutrinos.py.

◆ hnuratio

extract_interacting_neutrinos.hnuratio

Initial value:

=  r.TH1D(
    "hnuratio",
    "ratio between Genie and Fluka energies, proportional to the cross section",
    nbins,
    minE,
    maxE,
)

Definition at line 74 of file extract_interacting_neutrinos.py.

◆ matchingratio

extract_interacting_neutrinos.matchingratio = hnuratio.GetBinContent(hnuratio.FindBin(Ekin))

Definition at line 104 of file extract_interacting_neutrinos.py.

◆ maxE

extract_interacting_neutrinos.maxE = hfluka_nuE.GetXaxis().GetBinUpEdge(nbins)

Definition at line 45 of file extract_interacting_neutrinos.py.

◆ maxratio

extract_interacting_neutrinos.maxratio = hnuratio.GetMaximum()

Definition at line 88 of file extract_interacting_neutrinos.py.

◆ minE

extract_interacting_neutrinos.minE = hfluka_nuE.GetXaxis().GetBinLowEdge(1)

Definition at line 44 of file extract_interacting_neutrinos.py.

◆ nbins

extract_interacting_neutrinos.nbins = hfluka_nuE.GetNbinsX()

Definition at line 43 of file extract_interacting_neutrinos.py.

◆ nneutrinos

extract_interacting_neutrinos.nneutrinos = flukatree.GetEntries()

Definition at line 96 of file extract_interacting_neutrinos.py.

◆ options

extract_interacting_neutrinos.options = parser.parse_args()

Definition at line 22 of file extract_interacting_neutrinos.py.

◆ output_nutree

extract_interacting_neutrinos.output_nutree = flukatree.CloneTree(0)

Definition at line 92 of file extract_interacting_neutrinos.py.

◆ parser

extract_interacting_neutrinos.parser = argparse.ArgumentParser(description="Extract interacting neutrinos")

Definition at line 6 of file extract_interacting_neutrinos.py.

◆ randomratio

extract_interacting_neutrinos.randomratio = r.gRandom.Uniform(0, maxratio)

Definition at line 102 of file extract_interacting_neutrinos.py.

◆ type

extract_interacting_neutrinos.type

Definition at line 9 of file extract_interacting_neutrinos.py.

Functions

Variables