13from array
import array
20 hh = ROOT.TH1D(
"hh",
"hh", 25, -10, 10)
22 hh.Fill(ROOT.gRandom.Gaus(0, 3))
30 tree = ROOT.TTree(
"tree",
"tree")
33 tree.Branch(
"x", px,
"x/D")
34 tree.Branch(
"y", py,
"y/D")
36 px[0] = ROOT.gRandom.Gaus(0, 3)
37 py[0] = ROOT.gRandom.Uniform() * 30 - 15
51x = ROOT.RooRealVar(
"x",
"x", -10, 10)
55dh = ROOT.RooDataHist(
"dh",
"dh", [x], Import=hh)
60frame = x.frame(Title=
"Imported ROOT.TH1 with Poisson error bars")
64mean = ROOT.RooRealVar(
"mean",
"mean", 0, -10, 10)
65sigma = ROOT.RooRealVar(
"sigma",
"sigma", 3, 0.1, 10)
66gauss = ROOT.RooGaussian(
"gauss",
"gauss", x, mean, sigma)
76frame2 = x.frame(Title=
"Imported ROOT.TH1 with internal errors")
77dh.plotOn(frame2, DataError=
"SumW2")
94y = ROOT.RooRealVar(
"y",
"y", -10, 10)
104ds = ROOT.RooDataSet(
"ds",
"ds", {x, y}, Import=tree)
110def write_dataset(ds, filename):
112 outstream = ROOT.std.ofstream(filename)
118write_dataset(ds,
"rf102_testData.txt")
122print(
"\n-----------------------\nReading data from ASCII")
123dataReadBack = ROOT.RooDataSet.read(
124 "rf102_testData.txt",
129dataReadBack.Print(
"V")
131print(
"\nOriginal data, line 20:")
134print(
"\nRead-back data, line 20:")
135dataReadBack.get(20).
Print(
"V")
144frame3 = y.frame(Title=
"Unbinned data shown in default frame binning")
148frame4 = y.frame(Title=
"Unbinned data shown with custom binning")
149ds.plotOn(frame4, Binning=20)
151frame5 = y.frame(Title=
"Unbinned data read back from ASCII file")
152ds.plotOn(frame5, Binning=20)
153dataReadBack.plotOn(frame5, Binning=20, MarkerColor=
"r", MarkerStyle=5)
156c = ROOT.TCanvas(
"rf102_dataimport",
"rf102_dataimport", 800, 800)
159ROOT.gPad.SetLeftMargin(0.15)
160frame.GetYaxis().SetTitleOffset(1.4)
163ROOT.gPad.SetLeftMargin(0.15)
164frame2.GetYaxis().SetTitleOffset(1.4)
167ROOT.gPad.SetLeftMargin(0.15)
168frame3.GetYaxis().SetTitleOffset(1.4)
171ROOT.gPad.SetLeftMargin(0.15)
172frame4.GetYaxis().SetTitleOffset(1.4)
175ROOT.gPad.SetLeftMargin(0.15)
176frame4.GetYaxis().SetTitleOffset(1.4)
179c.SaveAs(
"rf102_dataimport.png")
void Print(std::ostream &os, const OptionType &opt)