|
| | truth_ana.tm = ROOT.TreeManager("./structuredtree.root") |
| | central garana tree manager it links to all of the individual tree accessors given a file with garana trees produced in the default garana format More...
|
| |
| | truth_ana.bt = ROOT.Backtracker(tm) |
| |
| | truth_ana.pdgDB = ROOT.TDatabasePDG() |
| |
| dictionary | truth_ana.pdgToName |
| |
| | truth_ana.gen = tm.GetGenTree() |
| |
| | truth_ana.g4 = tm.GetG4Tree() |
| |
| | truth_ana.hg4nus = TH2F("hg4nus","g4Tree;N^{0} #nu's;GENIE scatter code",3,0,3,30,0,30) |
| |
| | truth_ana.hnureg = TH1F("hnureg",";#nu vertex region;",10,0,10) |
| |
| | truth_ana.hngen = TH1F("hngen","genTree;N^{0} MCTruth;",5,0,5) |
| |
| | truth_ana.hnprim = TH1F("hnprim","g4Tree;N^{0} primaries;",40,0,40) |
| |
| | truth_ana.hprim_pdg = TH1F("hprim_pdg","g4Tree; MCParticle PDG code;",3400,-400,3000) |
| |
| | truth_ana.hprim_mumom = TH1F("hprim_mumom","g4Tree;true momentum [GeV/c];",50,0,5) |
| |
| | truth_ana.hprim_pmom = TH1F("hprim_pmom","g4Tree;true momentum [GeV/c];",50,0,5) |
| |
| | truth_ana.hprim_mumult = TH1F("hprim_mumult","g4Tree;primary particle multiplicity;",19,0,19) |
| |
| | truth_ana.hprim_pmult = TH1F("hprim_pmult","g4Tree;primary particle multiplicity;",19,0,19) |
| |
| | truth_ana.hprim_stopmult_tpc = TH1F("hprim_stopmult_tpc","g4Tree; N^{0} primaries stopping in active volume",10,0,10) |
| |
| | truth_ana.hprim_stopmult_calo = TH1F("hprim_stopmult_calo","g4Tree; N^{0} primaries stopping in active volume",10,0,10) |
| |
| | truth_ana.hprim_chgmult_tpc = TH1F("hprim_chgmult_tpc","g4Tree;N^{0} charged primaries;",19,0,19) |
| |
| | truth_ana.hprim_chgmult_calo = TH1F("hprim_chgmult_calo","g4Tree;N^{0} charged primaries;",19,0,19) |
| |
| | truth_ana.hprim_theta_all = TH1F("hprim_theta_all","g4Tree: primaries;#theta [deg]",90,0,180) |
| |
| | truth_ana.hprim_theta_mu = TH1F("hprim_theta_mu", "g4Tree: primaries;#theta [deg]",90,0,180) |
| |
| | truth_ana.hprim_theta_p = TH1F("hprim_theta_p", "g4Tree: primaries;#theta [deg]",90,0,180) |
| |
| | truth_ana.hprim_theta_pic = TH1F("hprim_theta_pic","g4Tree: primaries;#theta [deg]",90,0,180) |
| |
| | truth_ana.hprim_t0 = TH1F("hprim_t0", "g4Tree; time [ns];",100,0,1e4) |
| |
| | truth_ana.hprim_tf = TH1F("hprim_tf", "g4Tree; time [ns];",100,0,1e4) |
| |
| | truth_ana.hall_t0 = TH1F("hall_t0", "g4Tree; time [ns];",100,0,1e4) |
| |
| | truth_ana.hall_tf = TH1F("hall_tf", "g4Tree; time [ns];",100,0,1e4) |
| |
| | truth_ana.henu_mu = TH1F("henu_mu","genTree;E_{#nu} [GeV]", 50,0,10) |
| |
| | truth_ana.henu_mubar = TH1F("henu_mubar","genTree;E_{#nu} [GeV]",50,0,10) |
| |
| | truth_ana.henu_e = TH1F("henu_e","genTree;E_{#nu} [GeV]", 50,0,10) |
| |
| | truth_ana.henu_ebar = TH1F("henu_ebar","genTree;E_{#nu} [GeV]", 50,0,10) |
| |
| | truth_ana.hx = TH1F("hx","g4Tree;x [cm]",100,-300,300) |
| |
| | truth_ana.hy = TH1F("hy","g4Tree;y [cm]",100,-450,120) |
| |
| | truth_ana.hz = TH1F("hz","g4Tree;z [cm]",100,1250,1750) |
| |
| dictionary | truth_ana.pdgcts = {13:0} |
| |
| dictionary | truth_ana.pdgcts_calo = {13:0} |
| |
| int | truth_ana.nnu = 0 |
| |
| int | truth_ana.nprim_mu = 0 |
| |
| int | truth_ana.nprim_p = 0 |
| |
| int | truth_ana.nstop_tpc = 0 |
| |
| int | truth_ana.nstop_calo = 0 |
| |
| int | truth_ana.nprim_chg_tpc = 0 |
| |
| int | truth_ana.nprim_chg_calo = 0 |
| |
| int | truth_ana.inu = -1 |
| |
| | truth_ana.theta = ROOT.TMath.RadToDeg()*g4.SimMomBegin(ig4).Theta() |
| |
| | truth_ana.pdgcode = g4.PDG(ig4) |
| |
| | truth_ana.partPDG = pdgDB.GetParticle(pdgcode) |
| |
| | truth_ana.enu = gen.NuP(igen).E() |
| |
| | truth_ana.hpdgtxt = TH1F("hpdgtxt","g4Tree: primaries;; counts per species",len(pdgcts),0,len(pdgcts)) |
| |
| | truth_ana.hpdgtxt_calo = TH1F("hpdgtxt_calo","g4Tree: primaries;; counts per species",len(pdgcts),0,len(pdgcts)) |
| |
| | truth_ana.cngen = TCanvas() |
| |
| | truth_ana.cnprim = TCanvas() |
| |
| | truth_ana.cg4nus = TCanvas() |
| |
| | truth_ana.cnureg = TCanvas() |
| |
| | truth_ana.cenu = TCanvas() |
| |
| | truth_ana.cprim_pdg = TCanvas() |
| |
| | truth_ana.cprim_mom = TCanvas() |
| |
| | truth_ana.cprim_mult = TCanvas() |
| |
| | truth_ana.leg3 = ROOT.TLegend(0.6,0.6,0.8,0.8) |
| |
| | truth_ana.cprim_stopmult = TCanvas() |
| |
| | truth_ana.leg4 = ROOT.TLegend(0.6,0.6,0.8,0.8) |
| |
| | truth_ana.cprim_chgmult = TCanvas() |
| |
| | truth_ana.cpdgtxt = TCanvas() |
| |
| | truth_ana.leg = ROOT.TLegend(0.6,0.2,0.8,0.4) |
| |
| | truth_ana.ctheta = TCanvas() |
| |
| | truth_ana.leg2 = ROOT.TLegend(0.6,0.6,0.8,0.8) |
| |
| | truth_ana.cx = TCanvas() |
| |
| | truth_ana.cy = TCanvas() |
| |
| | truth_ana.cz = TCanvas() |
| |
| | truth_ana.null = input("press <Enter> to close canvas and exit program.") |
| |
1.8.11
