import ROOT

ROOT.gROOT.SetStyle("Plain")

x = ROOT.RooRealVar("x", "x", -10, 10)

mu1 = ROOT.RooRealVar("mu1", "average 1", 0, -1, 1)
sigma1 = ROOT.RooRealVar ("sigma1", "sigma 1", 1, 0, 5)
gauss1 = ROOT.RooGaussian("gauss1","gaussian 1 PDF", x, mu1, sigma1)

mu2 = ROOT.RooRealVar("mu2", "average 2", 0.5, -1, 1)
sigma2 = ROOT.RooRealVar("sigma2", "sigma 2", 0.5, 0, 5)
gauss2 = ROOT.RooGaussian("gauss2","gaussian 2 PDF", x, mu2, sigma2)

mu3 = ROOT.RooRealVar("mu3", "average 2", -0.5, -1, 1)
sigma3 = ROOT.RooRealVar("sigma3", "sigma 2", 1.5, 0, 5)
gauss3 = ROOT.RooGaussian("gauss3","gaussian 2 PDF", x, mu3, sigma3)

xFrame = x.frame()
gauss1.plotOn(xFrame, ROOT.RooFit.LineColor(ROOT.kRed))
gauss2.plotOn(xFrame, ROOT.RooFit.LineColor(ROOT.kBlue))
gauss3.plotOn(xFrame, ROOT.RooFit.LineColor(ROOT.kMagenta))
c = ROOT.TCanvas()
xFrame.Draw()
c.SaveAs("01-gaussians.pdf")

