Logo ROOT  
Reference Guide
ntuple1.C
Go to the documentation of this file.
1 /// \file
2 /// \ingroup tutorial_tree
3 /// \notebook
4 /// Simple tree analysis.
5 ///
6 /// \macro_image
7 /// \macro_code
8 ///
9 /// \author Rene Brun
10 
11 #include "TCanvas.h"
12 #include "TFile.h"
13 #include "TNtuple.h"
14 #include "TProfile.h"
15 #include "TBenchmark.h"
16 #include "TStyle.h"
17 #include "TPaveText.h"
18 #include "TFrame.h"
19 #include "TF1.h"
20 #include "TROOT.h"
21 #include "TSystem.h"
22 #include "TInterpreter.h"
23 
24 void ntuple1() {
25 
26  //just in case this script is executed multiple times
27  delete gROOT->GetListOfFiles()->FindObject("hsimple.root");
28  delete gROOT->GetListOfCanvases()->FindObject("c1");
29 
30  gBenchmark->Start("ntuple1");
31  //
32  // Connect ROOT histogram/ntuple demonstration file
33  // generated by example $ROOTSYS/tutorials/hsimple.C.
34  TString dir = gROOT->GetTutorialDir();
35  dir.Append("/hsimple.C");
36  dir.ReplaceAll("/./","/");
37  if (gBenchmark->GetBench("hsimple") < 0) gInterpreter->LoadMacro(dir.Data());
38  TFile *f1 = (TFile*)gROOT->ProcessLineFast("hsimple(1)");
39  if (!f1) return;
40  //
41  // Create a canvas, with 4 pads
42  //
43  TCanvas *c1 = new TCanvas("c1","The Ntuple canvas",200,10,700,780);
44  TPad *pad1 = new TPad("pad1","This is pad1",0.02,0.52,0.48,0.98,21);
45  TPad *pad2 = new TPad("pad2","This is pad2",0.52,0.52,0.98,0.98,21);
46  TPad *pad3 = new TPad("pad3","This is pad3",0.02,0.02,0.48,0.48,21);
47  TPad *pad4 = new TPad("pad4","This is pad4",0.52,0.02,0.98,0.48,1);
48  pad1->Draw();
49  pad2->Draw();
50  pad3->Draw();
51  pad4->Draw();
52  //
53  // Change default style for the statistics box
54  gStyle->SetStatW(0.30);
55  gStyle->SetStatH(0.20);
56  //
57  // Display a function of one ntuple column imposing a condition
58  // on another column.
59  pad1->cd();
60  pad1->SetGrid();
61  pad1->SetLogy();
62  TNtuple *ntuple = (TNtuple*)f1->Get("ntuple");
63  ntuple->SetLineColor(1);
64  ntuple->SetFillStyle(1001);
65  ntuple->SetFillColor(45);
66  ntuple->Draw("3*px+2","px**2+py**2>1");
67  ntuple->SetFillColor(38);
68  ntuple->Draw("2*px+2","pz>2","same");
69  ntuple->SetFillColor(5);
70  ntuple->Draw("1.3*px+2","(px^2+py^2>4) && py>0","same");
71  pad1->RedrawAxis();
72  //
73  // Display the profile of two columns
74  // The profile histogram produced is saved in the current directory with
75  // the name hprofs
76  pad2->cd();
77  pad2->SetGrid();
78  ntuple->Draw("pz:px>>hprofs","","goffprofs");
79  TProfile *hprofs = (TProfile*)gDirectory->Get("hprofs");
80  hprofs->SetMarkerColor(5);
81  hprofs->SetMarkerSize(0.7);
82  hprofs->SetMarkerStyle(21);
83  hprofs->Fit("pol2");
84  // Get pointer to fitted function and modify its attributes
85  TF1 *fpol2 = hprofs->GetFunction("pol2");
86  fpol2->SetLineWidth(4);
87  fpol2->SetLineColor(2);
88  //
89  // Display a scatter plot of two columns with a selection.
90  // Superimpose the result of another cut with a different marker color
91  pad3->cd();
92  pad3->GetFrame()->SetBorderSize(8);
93  ntuple->SetMarkerColor(1);
94  ntuple->Draw("py:px","pz>1");
95  ntuple->SetMarkerColor(2);
96  ntuple->Draw("py:px","pz<1","same");
97  //
98  // Display a 3-D scatter plot of 3 columns. Superimpose a different selection.
99  pad4->cd();
100  ntuple->Draw("pz:py:px","(pz<10 && pz>6)+(pz<4 && pz>3)");
101  ntuple->SetMarkerColor(4);
102  ntuple->Draw("pz:py:px","pz<6 && pz>4","same");
103  ntuple->SetMarkerColor(5);
104  ntuple->Draw("pz:py:px","pz<4 && pz>3","same");
105  TPaveText *l4 = new TPaveText(-0.9,0.5,0.9,0.95);
106  l4->SetFillColor(42);
107  l4->SetTextAlign(12);
108  l4->AddText("You can interactively rotate this view in 2 ways:");
109  l4->AddText(" - With the RotateCube in clicking in this pad");
110  l4->AddText(" - Selecting View with x3d in the View menu");
111  l4->Draw();
112  //
113  c1->cd();
114  c1->Update();
115  gStyle->SetStatColor(19);
116  gBenchmark->Show("ntuple1");
117 }
TAttMarker::SetMarkerSize
virtual void SetMarkerSize(Size_t msize=1)
Set the marker size.
Definition: TAttMarker.h:41
TStyle::SetStatW
void SetStatW(Float_t w=0.19)
Definition: TStyle.h:382
TH1::GetFunction
virtual TF1 * GetFunction(const char *name) const
Return pointer to function with name.
Definition: TH1.cxx:8726
TAttFill::SetFillStyle
virtual void SetFillStyle(Style_t fstyle)
Set the fill area style.
Definition: TAttFill.h:39
TNtuple
A simple TTree restricted to a list of float variables only.
Definition: TNtuple.h:28
TPad
The most important graphics class in the ROOT system.
Definition: TPad.h:26
TBenchmark::Start
virtual void Start(const char *name)
Starts Benchmark with the specified name.
Definition: TBenchmark.cxx:172
TString::Data
const char * Data() const
Definition: TString.h:369
gInterpreter
#define gInterpreter
Definition: TInterpreter.h:560
TStyle.h
TStyle::SetStatColor
void SetStatColor(Color_t color=19)
Definition: TStyle.h:373
TAttLine::SetLineColor
virtual void SetLineColor(Color_t lcolor)
Set the line color.
Definition: TAttLine.h:40
TCanvas.h
TProfile.h
TString
Basic string class.
Definition: TString.h:136
TFile.h
TString::ReplaceAll
TString & ReplaceAll(const TString &s1, const TString &s2)
Definition: TString.h:692
TPaveText.h
TNtuple.h
TROOT.h
TPaveText::Draw
virtual void Draw(Option_t *option="")
Draw this pavetext with its current attributes.
Definition: TPaveText.cxx:234
TBenchmark::GetBench
Int_t GetBench(const char *name) const
Returns index of Benchmark name.
Definition: TBenchmark.cxx:106
TPad::SetLogy
virtual void SetLogy(Int_t value=1)
Set Lin/Log scale for Y.
Definition: TPad.cxx:5913
gStyle
R__EXTERN TStyle * gStyle
Definition: TStyle.h:412
TSystem.h
TAttMarker::SetMarkerColor
virtual void SetMarkerColor(Color_t mcolor=1)
Set the marker color.
Definition: TAttMarker.h:38
gBenchmark
R__EXTERN TBenchmark * gBenchmark
Definition: TBenchmark.h:59
TString::Append
TString & Append(const char *cs)
Definition: TString.h:564
gDirectory
#define gDirectory
Definition: TDirectory.h:236
TTree::Draw
virtual void Draw(Option_t *opt)
Default Draw method for all objects.
Definition: TTree.h:427
TBenchmark.h
TWbox::SetBorderSize
virtual void SetBorderSize(Short_t bordersize)
Definition: TWbox.h:52
TStyle::SetStatH
void SetStatH(Float_t h=0.1)
Definition: TStyle.h:383
TFile
A ROOT file is a suite of consecutive data records (TKey instances) with a well defined format.
Definition: TFile.h:54
TPaveText::AddText
virtual TText * AddText(Double_t x1, Double_t y1, const char *label)
Add a new Text line to this pavetext at given coordinates.
Definition: TPaveText.cxx:183
TProfile
Profile Histogram.
Definition: TProfile.h:32
TAttLine::SetLineWidth
virtual void SetLineWidth(Width_t lwidth)
Set the line width.
Definition: TAttLine.h:43
TAttText::SetTextAlign
virtual void SetTextAlign(Short_t align=11)
Set the text alignment.
Definition: TAttText.h:41
f1
TF1 * f1
Definition: legend1.C:11
TF1.h
TCanvas
The Canvas class.
Definition: TCanvas.h:23
TInterpreter.h
TPad::cd
TVirtualPad * cd(Int_t subpadnumber=0)
Set Current pad.
Definition: TPad.cxx:602
TAttFill::SetFillColor
virtual void SetFillColor(Color_t fcolor)
Set the fill area color.
Definition: TAttFill.h:37
TFrame.h
TPad::GetFrame
TFrame * GetFrame()
Get frame.
Definition: TPad.cxx:2855
TPad::SetGrid
virtual void SetGrid(Int_t valuex=1, Int_t valuey=1)
Definition: TPad.h:327
ntuple1
Definition: ntuple1.py:1
TPaveText
A Pave (see TPave) with text, lines or/and boxes inside.
Definition: TPaveText.h:21
TH1::Fit
virtual TFitResultPtr Fit(const char *formula, Option_t *option="", Option_t *goption="", Double_t xmin=0, Double_t xmax=0)
Fit histogram with function fname.
Definition: TH1.cxx:3869
TAttMarker::SetMarkerStyle
virtual void SetMarkerStyle(Style_t mstyle=1)
Set the marker style.
Definition: TAttMarker.h:40
TF1
1-Dim function class
Definition: TF1.h:213
TBenchmark::Show
virtual void Show(const char *name)
Stops Benchmark name and Prints results.
Definition: TBenchmark.cxx:155
TPad::RedrawAxis
virtual void RedrawAxis(Option_t *option="")
Redraw the frame axis.
Definition: TPad.cxx:5284
gROOT
#define gROOT
Definition: TROOT.h:406
c1
return c1
Definition: legend1.C:41
TPad::Draw
virtual void Draw(Option_t *option="")
Draw Pad in Current pad (re-parent pad if necessary).
Definition: TPad.cxx:1293