3 from ROOT
import TH1F, TH2F, TFile
7 (
"ndecay", 16, -0.5, 15.5),
8 (
"ntype", 31, -15.5, 15.5),
9 (
"ptype", 7000, -3500., 3500.),
10 (
"nimpwt",100, 0., 50.),
11 (
"vz", 1000, -10000.0, 90000.0)
14 (
"vx_vs_vy", 100, -500.0, 500., 100, -500.0, 500.0),
15 (
"pdpz_vs_pdpt",100, 0, 60., 100, 0, 1.0),
16 (
"pppz_vs_pppt",100, 0, 60., 100, 0, 1.0),
18 hists = hists1D + hists2D
31 wgt = (thisentry.dk2nu.decay.nimpwt)*(thisentry.dk2nu.nuray[self.
idxdet].wgt)
39 self.
thisgroup[5].
Fill(thisentry.dk2nu.decay.vx,thisentry.dk2nu.decay.vy,wgt)
40 pdpt =
sqrt( thisentry.dk2nu.decay.pdpx* thisentry.dk2nu.decay.pdpx +
41 thisentry.dk2nu.decay.pdpy * thisentry.dk2nu.decay.pdpy )
43 pppt = thisentry.dk2nu.decay.pppz *
sqrt( thisentry.dk2nu.decay.ppdxdz* thisentry.dk2nu.decay.ppdxdz +
44 thisentry.dk2nu.decay.ppdydz* thisentry.dk2nu.decay.ppdydz )
62 hist_group.append(TH1F(
"{}_{}".
format(self.
thistag, i[0]),
";{};".
format(i[0]), i[1], i[2], i[3]))
64 hist_group.append(TH2F(
"{}_{}".
format(self.
thistag, i[0]),
";{};".
format(i[0]), i[1], i[2], i[3], i[4], i[5], i[6]))
71 if detname ==
"random":
73 elif detname ==
"minosnd":
75 elif detname ==
"novand":
77 elif detname ==
"minosfd":
79 elif detname ==
"novafd":
82 print "Bad input, no detector found!" def __init__(self, root_file, det)
def save(self, root_file)
correl_yv Fill(-(dy[iP-1][iC-1]), hyv->GetBinContent(iP, iC))
def fill(self, thisentry)
def GetIdxDet(self, detname)
std::string format(const int32_t &value, const int &ndigits=8)