24x = ROOT.RooRealVar(
"x",
"x", -5, 5)
25y = ROOT.RooRealVar(
"y",
"y", -5, 5)
26z = ROOT.RooRealVar(
"z",
"z", -5, 5)
30 "gx",
"gx", x, ROOT.RooFit.RooConst(0), ROOT.RooFit.RooConst(1))
32 "gy",
"gy", y, ROOT.RooFit.RooConst(0), ROOT.RooFit.RooConst(1))
34 "gz",
"gz", z, ROOT.RooFit.RooConst(0), ROOT.RooFit.RooConst(1))
35sig = ROOT.RooProdPdf(
"sig",
"sig", ROOT.RooArgList(gx, gy, gz))
38px = ROOT.RooPolynomial(
"px",
"px", x, ROOT.RooArgList(
39 ROOT.RooFit.RooConst(-0.1), ROOT.RooFit.RooConst(0.004)))
40py = ROOT.RooPolynomial(
"py",
"py", y, ROOT.RooArgList(
41 ROOT.RooFit.RooConst(0.1), ROOT.RooFit.RooConst(-0.004)))
42pz = ROOT.RooPolynomial(
"pz",
"pz", z)
43bkg = ROOT.RooProdPdf(
"bkg",
"bkg", ROOT.RooArgList(px, py, pz))
46fsig = ROOT.RooRealVar(
"fsig",
"signal fraction", 0.1, 0., 1.)
47model = ROOT.RooAddPdf(
"model",
"model", ROOT.RooArgList(sig, bkg), ROOT.RooArgList(fsig))
49data = model.generate(ROOT.RooArgSet(x, y, z), 20000)
55frame = x.frame(ROOT.RooFit.Title(
56 "Projection of 3D data and pdf on X"), ROOT.RooFit.Bins(40))
65sigyz = sig.createProjection(ROOT.RooArgSet(x))
66totyz = model.createProjection(ROOT.RooArgSet(x))
69llratio_func = ROOT.RooFormulaVar(
70 "llratio",
"log10(@0)-log10(@1)", ROOT.RooArgList(sigyz, totyz))
76data.addColumn(llratio_func)
79dataSel = data.reduce(ROOT.RooFit.Cut(
"llratio>0.7"))
82frame2 = x.frame(ROOT.RooFit.Title(
83 "Same projection on X with LLratio(y,z)>0.7"), ROOT.RooFit.Bins(40))
92mcprojData = model.generate(ROOT.RooArgSet(x, y, z), 10000)
96mcprojData.addColumn(llratio_func)
97mcprojDataSel = mcprojData.reduce(ROOT.RooFit.Cut(
"llratio>0.7"))
101model.plotOn(frame2, ROOT.RooFit.ProjWData(mcprojDataSel))
103c = ROOT.TCanvas(
"rf316_llratioplot",
"rf316_llratioplot", 800, 400)
106ROOT.gPad.SetLeftMargin(0.15)
107frame.GetYaxis().SetTitleOffset(1.4)
110ROOT.gPad.SetLeftMargin(0.15)
111frame2.GetYaxis().SetTitleOffset(1.4)
113c.SaveAs(
"rf316_llratioplot.png")