21 x = ROOT.RooRealVar(
"x",
"x", 0, 10)
25 mean = ROOT.RooRealVar(
"mean",
"mean of gaussians", 5)
26 sigma1 = ROOT.RooRealVar(
"sigma1",
"width of gaussians", 0.5)
27 sigma2 = ROOT.RooRealVar(
"sigma2",
"width of gaussians", 1)
29 sig1 = ROOT.RooGaussian(
"sig1",
"Signal component 1", x, mean, sigma1)
30 sig2 = ROOT.RooGaussian(
"sig2",
"Signal component 2", x, mean, sigma2)
33 a0 = ROOT.RooRealVar(
"a0",
"a0", 0.5, 0., 1.)
34 a1 = ROOT.RooRealVar(
"a1",
"a1", -0.2, 0., 1.)
35 bkg = ROOT.RooChebychev(
"bkg",
"Background", x, ROOT.RooArgList(a0, a1))
43 sig1frac = ROOT.RooRealVar(
44 "sig1frac",
"fraction of component 1 in signal", 0.8, 0., 1.)
45 sig = ROOT.RooAddPdf(
"sig",
"Signal", ROOT.RooArgList(
46 sig1, sig2), ROOT.RooArgList(sig1frac))
52 bkgfrac = ROOT.RooRealVar(
"bkgfrac",
"fraction of background", 0.5, 0., 1.)
53 model = ROOT.RooAddPdf(
54 "model",
"g1+g2+a", ROOT.RooArgList(bkg, sig), ROOT.RooArgList(bkgfrac))
60 data = model.generate(ROOT.RooArgSet(x), 1000)
66 xframe = x.frame(ROOT.RooFit.Title(
67 "Example of composite pdf=(sig1+sig2)+bkg"))
72 ras_bkg = ROOT.RooArgSet(bkg)
73 model.plotOn(xframe, ROOT.RooFit.Components(ras_bkg),
74 ROOT.RooFit.LineStyle(ROOT.kDashed))
77 ras_bkg_sig2 = ROOT.RooArgSet(bkg, sig2)
78 model.plotOn(xframe, ROOT.RooFit.Components(ras_bkg_sig2),
79 ROOT.RooFit.LineStyle(ROOT.kDotted))
91 model2 = ROOT.RooAddPdf(
112 model2.plotOn(xframe, ROOT.RooFit.LineColor(ROOT.kRed),
113 ROOT.RooFit.LineStyle(ROOT.kDashed))
116 ROOT.RooFit.Components(ras_bkg_sig2),
117 ROOT.RooFit.LineColor(
119 ROOT.RooFit.LineStyle(
124 c = ROOT.TCanvas(
"rf201_composite",
"rf201_composite", 600, 600)
125 ROOT.gPad.SetLeftMargin(0.15)
126 xframe.GetYaxis().SetTitleOffset(1.4)
129 c.SaveAs(
"rf201_composite.png")