Logo ROOT   6.12/07
Reference Guide
ratioplotOld.C
Go to the documentation of this file.
1 /// \file
2 /// \ingroup tutorial_hist
3 /// Example displaying two histograms and their ratio. This macro does not use the
4 /// class TRatioPlot. For ROOT version >= 6.08 TRatioPlot should be used. See
5 /// the other ratio plots examples in this folder.
6 ///
7 /// \macro_image
8 /// \macro_code
9 ///
10 /// \author Olivier Couet
11 
12 void ratioplotOld( ) {
13  // Define two gaussian histograms. Note the X and Y title are defined
14  // at booking time using the convention "Hist_title ; X_title ; Y_title"
15  TH1F *h1 = new TH1F("h1", "Two gaussian plots and their ratio;x title; h1 and h2 gaussian histograms", 100, -5, 5);
16  TH1F *h2 = new TH1F("h2", "h2", 100, -5, 5);
17  h1->FillRandom("gaus");
18  h2->FillRandom("gaus");
19 
20  // Define the Canvas
21  TCanvas *c = new TCanvas("c", "canvas", 800, 800);
22 
23  // Upper plot will be in pad1
24  TPad *pad1 = new TPad("pad1", "pad1", 0, 0.3, 1, 1.0);
25  pad1->SetBottomMargin(0); // Upper and lower plot are joined
26  pad1->SetGridx(); // Vertical grid
27  pad1->Draw(); // Draw the upper pad: pad1
28  pad1->cd(); // pad1 becomes the current pad
29  h1->SetStats(0); // No statistics on upper plot
30  h1->Draw(); // Draw h1
31  h2->Draw("same"); // Draw h2 on top of h1
32 
33  // Do not draw the Y axis label on the upper plot and redraw a small
34  // axis instead, in order to avoid the first label (0) to be clipped.
35  h1->GetYaxis()->SetLabelSize(0.);
36  TGaxis *axis = new TGaxis( -5, 20, -5, 220, 20,220,510,"");
37  axis->SetLabelFont(43); // Absolute font size in pixel (precision 3)
38  axis->SetLabelSize(15);
39  axis->Draw();
40 
41  // lower plot will be in pad
42  c->cd(); // Go back to the main canvas before defining pad2
43  TPad *pad2 = new TPad("pad2", "pad2", 0, 0.05, 1, 0.3);
44  pad2->SetTopMargin(0);
45  pad2->SetBottomMargin(0.2);
46  pad2->SetGridx(); // vertical grid
47  pad2->Draw();
48  pad2->cd(); // pad2 becomes the current pad
49 
50  // Define the ratio plot
51  TH1F *h3 = (TH1F*)h1->Clone("h3");
52  h3->SetLineColor(kBlack);
53  h3->SetMinimum(0.8); // Define Y ..
54  h3->SetMaximum(1.35); // .. range
55  h3->Sumw2();
56  h3->SetStats(0); // No statistics on lower plot
57  h3->Divide(h2);
58  h3->SetMarkerStyle(21);
59  h3->Draw("ep"); // Draw the ratio plot
60 
61  // h1 settings
62  h1->SetLineColor(kBlue+1);
63  h1->SetLineWidth(2);
64 
65  // Y axis h1 plot settings
66  h1->GetYaxis()->SetTitleSize(20);
67  h1->GetYaxis()->SetTitleFont(43);
68  h1->GetYaxis()->SetTitleOffset(1.55);
69 
70  // h2 settings
71  h2->SetLineColor(kRed);
72  h2->SetLineWidth(2);
73 
74  // Ratio plot (h3) settings
75  h3->SetTitle(""); // Remove the ratio title
76 
77  // Y axis ratio plot settings
78  h3->GetYaxis()->SetTitle("ratio h1/h2 ");
79  h3->GetYaxis()->SetNdivisions(505);
80  h3->GetYaxis()->SetTitleSize(20);
81  h3->GetYaxis()->SetTitleFont(43);
82  h3->GetYaxis()->SetTitleOffset(1.55);
83  h3->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
84  h3->GetYaxis()->SetLabelSize(15);
85 
86  // X axis ratio plot settings
87  h3->GetXaxis()->SetTitleSize(20);
88  h3->GetXaxis()->SetTitleFont(43);
89  h3->GetXaxis()->SetTitleOffset(4.);
90  h3->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
91  h3->GetXaxis()->SetLabelSize(15);
92 }
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title Offset is a correction factor with respect to the "s...
Definition: TAttAxis.cxx:294
virtual void SetGridx(Int_t value=1)
Definition: TPad.h:328
virtual void SetLineWidth(Width_t lwidth)
Set the line width.
Definition: TAttLine.h:43
virtual void SetMaximum(Double_t maximum=-1111)
Definition: TH1.h:390
Definition: Rtypes.h:59
THist< 1, float, THistStatContent, THistStatUncertainty > TH1F
Definition: THist.hxx:285
Definition: Rtypes.h:58
virtual void SetNdivisions(Int_t n=510, Bool_t optim=kTRUE)
Set the number of divisions for this axis.
Definition: TAttAxis.cxx:229
TVirtualPad * cd(Int_t subpadnumber=0)
Set current canvas & pad.
Definition: TCanvas.cxx:688
virtual void SetMinimum(Double_t minimum=-1111)
Definition: TH1.h:391
1-D histogram with a float per channel (see TH1 documentation)}
Definition: TH1.h:567
virtual void SetTitleFont(Style_t font=62)
Set the title font.
Definition: TAttAxis.cxx:322
virtual void Draw(Option_t *option="")
Default Draw method for all objects.
Definition: TObject.cxx:195
virtual void SetTopMargin(Float_t topmargin)
Set Pad top margin in fraction of the pad height.
Definition: TAttPad.cxx:130
virtual void SetLabelFont(Style_t font=62)
Set labels&#39; font.
Definition: TAttAxis.cxx:183
TVirtualPad * cd(Int_t subpadnumber=0)
Set Current pad.
Definition: TPad.cxx:585
TH1F * h1
Definition: legend1.C:5
virtual void Draw(Option_t *option="")
Draw Pad in Current pad (re-parent pad if necessary).
Definition: TPad.cxx:1272
void SetLabelSize(Float_t labelsize)
Definition: TGaxis.h:108
virtual Bool_t Divide(TF1 *f1, Double_t c1=1)
Performs the operation: this = this/(c1*f1) if errors are defined (see TH1::Sumw2), errors are also recalculated.
Definition: TH1.cxx:2724
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition: TAttLine.h:40
virtual void SetBottomMargin(Float_t bottommargin)
Set Pad bottom margin in fraction of the pad height.
Definition: TAttPad.cxx:100
virtual void FillRandom(const char *fname, Int_t ntimes=5000)
Fill histogram following distribution in function fname.
Definition: TH1.cxx:3414
virtual void Draw(Option_t *option="")
Draw this histogram with options.
Definition: TH1.cxx:2969
The most important graphics class in the ROOT system.
Definition: TPad.h:29
The axis painter class.
Definition: TGaxis.h:24
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.
Definition: TAttMarker.h:40
TAxis * GetYaxis()
Definition: TH1.h:316
virtual void SetLabelSize(Float_t size=0.04)
Set size of axis labels The size is expressed in per cent of the pad width.
Definition: TAttAxis.cxx:204
virtual void SetTitleSize(Float_t size=0.04)
Set size of axis title The size is expressed in per cent of the pad width.
Definition: TAttAxis.cxx:304
The Canvas class.
Definition: TCanvas.h:31
void SetLabelFont(Int_t labelfont)
Definition: TGaxis.h:106
virtual void Sumw2(Bool_t flag=kTRUE)
Create structure to store sum of squares of weights.
Definition: TH1.cxx:8276
TObject * Clone(const char *newname=0) const
Make a complete copy of the underlying object.
Definition: TH1.cxx:2662
Definition: Rtypes.h:59
virtual void SetTitle(const char *title)
See GetStatOverflows for more information.
Definition: TH1.cxx:6154
virtual void SetTitle(const char *title="")
Set the title of the TNamed.
Definition: TNamed.cxx:164
virtual void SetStats(Bool_t stats=kTRUE)
Set statistics option on/off.
Definition: TH1.cxx:8247
TAxis * GetXaxis()
Get the behaviour adopted by the object about the statoverflows. See EStatOverflows for more informat...
Definition: TH1.h:315