24ROOT.RooAbsPdf.defaultIntegratorConfig().setEpsRel(1e-8)
25ROOT.RooAbsPdf.defaultIntegratorConfig().setEpsAbs(1e-8)
28x = ROOT.RooRealVar(
"x",
"x", -5, 5)
29y = ROOT.RooRealVar(
"y",
"y", -2, 2)
32a0 = ROOT.RooRealVar(
"a0",
"a0", 0)
33a1 = ROOT.RooRealVar(
"a1",
"a1", -1.5, -3, 1)
34fy = ROOT.RooPolyVar(
"fy",
"fy", y, ROOT.RooArgList(a0, a1))
37sigmax = ROOT.RooRealVar(
"sigmax",
"width of gaussian", 0.5)
38gaussx = ROOT.RooGaussian(
39 "gaussx",
"Gaussian in x with shifting mean in y", x, fy, sigmax)
42gaussy = ROOT.RooGaussian(
43 "gaussy",
"Gaussian in y", y, ROOT.RooFit.RooConst(0), ROOT.RooFit.RooConst(2))
46model = ROOT.RooProdPdf(
"model",
"gaussx(x|y)*gaussy(y)", ROOT.RooArgSet(
47 gaussy), ROOT.RooFit.Conditional(ROOT.RooArgSet(gaussx), ROOT.RooArgSet(x)))
53modelx = model.createProjection(ROOT.RooArgSet(y))
59data = modelx.generateBinned(ROOT.RooArgSet(x), 1000)
62modelx.fitTo(data, ROOT.RooFit.Verbose())
70hh = model.createHistogram(
"x,y")
71hh.SetLineColor(ROOT.kBlue)
73c = ROOT.TCanvas(
"rf315_projectpdf",
"rf315_projectpdf", 800, 400)
76ROOT.gPad.SetLeftMargin(0.15)
77frame.GetYaxis().SetTitleOffset(1.4)
80ROOT.gPad.SetLeftMargin(0.20)
81hh.GetZaxis().SetTitleOffset(2.5)
83c.SaveAs(
"rf315_projectpdf.png")