Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
hsum.C
Go to the documentation of this file.
1/// \file
2/// \ingroup tutorial_hist
3/// \notebook -js
4/// Histograms filled and drawn in a loop.
5/// Simple example illustrating how to use the C++ interpreter
6/// to fill histograms in a loop and show the graphics results
7///
8/// \macro_image
9/// \macro_code
10///
11/// \author Rene Brun
12
13void hsum()
14{
15 TCanvas *c1 = new TCanvas("c1","The HSUM example",200,10,600,400);
16 c1->SetGrid();
17
18 gBenchmark->Start("hsum");
19
20 // Create some histograms.
21 auto total = new TH1F("total","This is the total distribution",100,-4,4);
22 auto main = new TH1F("main","Main contributor",100,-4,4);
23 auto s1 = new TH1F("s1","This is the first signal",100,-4,4);
24 auto s2 = new TH1F("s2","This is the second signal",100,-4,4);
25 total->Sumw2(); // store the sum of squares of weights
26 total->SetMarkerStyle(21);
27 total->SetMarkerSize(0.7);
28 main->SetFillColor(16);
29 s1->SetFillColor(42);
30 s2->SetFillColor(46);
31 TSlider *slider = nullptr;
32
33 // Fill histograms randomly
35 const Int_t kUPDATE = 500;
36 Float_t xs1, xs2, xmain;
37 for ( Int_t i=0; i<10000; i++) {
38 xmain = gRandom->Gaus(-1,1.5);
39 xs1 = gRandom->Gaus(-0.5,0.5);
40 xs2 = gRandom->Landau(1,0.15);
41 main->Fill(xmain);
42 s1->Fill(xs1,0.3);
43 s2->Fill(xs2,0.2);
44 total->Fill(xmain);
45 total->Fill(xs1,0.3);
46 total->Fill(xs2,0.2);
47 if (i && (i%kUPDATE) == 0) {
48 if (i == kUPDATE) {
49 total->Draw("e1p");
50 main->Draw("same");
51 s1->Draw("same");
52 s2->Draw("same");
53 c1->Update();
54 slider = new TSlider("slider",
55 "test",4.2,0,4.6,total->GetMaximum(),38);
56 slider->SetFillColor(46);
57 }
58 if (slider) slider->SetRange(0., 1.*i/10000.);
59 c1->Modified();
60 c1->Update();
61 }
62 }
63 slider->SetRange(0., 1.);
64 c1->Modified();
65 gBenchmark->Show("hsum");
66}
int main()
Definition Prototype.cxx:12
#define s1(x)
Definition RSha256.hxx:91
int Int_t
Definition RtypesCore.h:45
float Float_t
Definition RtypesCore.h:57
R__EXTERN TBenchmark * gBenchmark
Definition TBenchmark.h:59
static unsigned int total
R__EXTERN TRandom * gRandom
Definition TRandom.h:62
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
Definition TAttFill.h:37
virtual void Start(const char *name)
Starts Benchmark with the specified name.
virtual void Show(const char *name)
Stops Benchmark name and Prints results.
The Canvas class.
Definition TCanvas.h:23
1-D histogram with a float per channel (see TH1 documentation)
Definition TH1.h:622
virtual Double_t Gaus(Double_t mean=0, Double_t sigma=1)
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigm...
Definition TRandom.cxx:275
virtual void SetSeed(ULong_t seed=0)
Set the random generator seed.
Definition TRandom.cxx:615
virtual Double_t Landau(Double_t mean=0, Double_t sigma=1)
Generate a random number following a Landau distribution with location parameter mu and scale paramet...
Definition TRandom.cxx:381
A specialized TPad including a TSliderBox object.
Definition TSlider.h:17
virtual void SetRange(Double_t xmin=0, Double_t xmax=1)
Set Slider range in [0,1].
Definition TSlider.cxx:183
return c1
Definition legend1.C:41
Definition hsum.py:1