18 x = ROOT.RooRealVar(
"x",
"x", -10, 10)
19 y = ROOT.RooRealVar(
"y",
"y", -10, 10)
22 d = ROOT.RooDataSet(
"d",
"d", coord)
24 for i
in range(10000):
25 tmpy = ROOT.gRandom.Gaus(0, 10)
26 tmpx = ROOT.gRandom.Gaus(0.5 * tmpy, 1)
27 if (
abs(tmpy) < 10)
and (
abs(tmpx) < 10):
39x = ROOT.RooRealVar(
"x",
"x", -10, 10)
40y = ROOT.RooRealVar(
"y",
"y", -10, 10)
43a0 = ROOT.RooRealVar(
"a0",
"a0", -0.5, -5, 5)
44a1 = ROOT.RooRealVar(
"a1",
"a1", -0.5, -1, 1)
45fy = ROOT.RooPolyVar(
"fy",
"fy", y, [a0, a1])
48sigma = ROOT.RooRealVar(
"sigma",
"width of gaussian", 0.5, 0.1, 2.0)
49model = ROOT.RooGaussian(
"model",
"Gaussian with shifting mean", x, fy, sigma)
52expDataXY = makeFakeDataXY()
58expDataY = expDataXY.reduce({y})
62data = model.generate({x}, ProtoData=expDataY)
68model.fitTo(expDataXY, ConditionalObservables={y})
76expDataXY.plotOn(xframe)
77model.plotOn(xframe, ProjWData=expDataY)
81binnedDataY = expDataY.binnedClone()
82model.plotOn(xframe, ProjWData=binnedDataY, LineColor=
"c", LineStyle=
":")
85(expDataY.get().find(
"y")).setBins(5)
86binnedDataY2 = expDataY.binnedClone()
87model.plotOn(xframe, ProjWData=binnedDataY2, LineColor=
"r")
90c = ROOT.TCanvas(
"rf303_conditional",
"rf303_conditional", 600, 460)
91ROOT.gPad.SetLeftMargin(0.15)
92xframe.GetYaxis().SetTitleOffset(1.2)
95c.SaveAs(
"rf303_conditional.png")
RVec< PromoteType< T > > abs(const RVec< T > &v)