import ROOT
form1 = ROOT.TFormula("form1", "abs(sin(x)/x)")
rangeMin = 0.0
rangeMax = 10.0
sqroot = ROOT.TF1("sqroot", "x*gaus(0) + [3]*form1", rangeMin, rangeMax)
sqroot.SetLineColor(4)
sqroot.SetLineWidth(6)
sqroot.SetParameters(10.0, 4.0, 1.0, 20.0)
nBins = 200
h1d = ROOT.TH1D("h1d", "Test random numbers", nBins, rangeMin, rangeMax)
h1d.FillRandom("sqroot", 10000)
topX = 200
topY = 10
width = 700
height = 900
c1 = ROOT.TCanvas("c1", "The FillRandom example", topX, topY, width, height)
pad1 = ROOT.TPad("pad1", "The pad with the function", 0.05, 0.50, 0.95, 0.95)
pad2 = ROOT.TPad("pad2", "The pad with the histogram", 0.05, 0.05, 0.95, 0.45)
pad1.Draw()
pad2.Draw()
pad1.cd()
pad1.SetGridx()
pad1.SetGridy()
pad1.GetFrame().SetBorderMode(-1)
sqroot.Draw()
lfunction = ROOT.TPaveLabel(5, 39, 9.8, 46, "The sqroot function")
lfunction.Draw()
c1.Update()
pad2.cd()
pad2.GetFrame().SetBorderMode(-1)
h1d.SetFillColor(45)
h1d.Draw()
c1.Update()