Sophie

Sophie

distrib > Fedora > 18 > x86_64 > by-pkgid > 8c86774a3e53d77cc119f53a2b94a57a > files > 1309

root-tutorial-5.34.14-2.fc18.noarch.rpm

# ruby-root testsuite
# port of the original $ROOT/hstack.C tutorial
# (20/01/2004)  --elathan  <elathan@phys.uoa.gr>
#
# original header:
# Example of stacked histograms: class THStack
# Author: Rene Brun

hs = THStack.new("hs","test stacked histograms")

#create three 1-d histograms
h1 = TH1F.new("h1","test hstack",100,-4,4)
    h1.FillRandom("gaus",20000)
    h1.SetFillColor(kRed)
    h1.SetMarkerStyle(21)
    h1.SetMarkerColor(kRed)

hs.Add(h1)

h2 = TH1F.new("h2","test hstack",100,-4,4)
    h2.FillRandom("gaus",15000)
    h2.SetFillColor(kBlue)
    h2.SetMarkerStyle(21)
    h2.SetMarkerColor(kBlue)

hs.Add(h2)

h3 = TH1F.new("h3","test hstack",100,-4,4)
    h3.FillRandom("gaus",10000)
    h3.SetFillColor(kGreen)
    h3.SetMarkerStyle(21)
    h3.SetMarkerColor(kGreen)

hs.Add(h3)
   
c1 = TCanvas.new("c1","stacked hists",10,10,1000,800)
    c1.SetFillColor(41)
    c1.Divide(2,2)
    # in top left pad, draw the stack with defaults
    c1.cd(1)
    
hs.Draw

    # in top right pad, draw the stack in non-stack mode and errors option
    c1.cd(2)
    gPad.SetGrid

hs.Draw("nostack,e1p")
   
    #in bottom left, draw in stack mode with "lego1" option
    c1.cd(3)
    gPad.SetFrameFillColor(17)
    gPad.SetTheta(3.77)
    gPad.SetPhi(2.9)

hs.Draw("lego1")

    c1.cd(4)
    #create two 2-D histograms and draw them in stack mode
    gPad.SetFrameFillColor(17)
   
a = THStack.new("a","test legos")
f1 = TF2.new("f1","xygaus + xygaus(5) + xylandau(10)",-4,4,-4,4)
    params = [130,-1.4,1.8,1.5,1, 150,2,0.5,-2,0.5, 3600,-2,0.7,-3,0.3]
    f1.SetParameters(params)
   
h2a = TH2F.new("h2a","h2a",20,-4,4,20,-4,4)
    h2a.SetFillColor(38)
    h2a.FillRandom("f1",4000)
   
    f2 = TF2.new("f2","xygaus + xygaus(5)",-4,4,-4,4)
    params = [100,-1.4,1.9,1.1,2, 80,2,0.7,-2,0.5]
    f2.SetParameters(params)
   
h2b = TH2F.new("h2b","h2b",20,-4,4,20,-4,4)
    h2b.SetFillColor(46)
    h2b.FillRandom("f2",3000)
    
a.Add(h2a)
a.Add(h2b)
a.Draw

gApplication.Run