1 | |
---|
2 | /** \class ExRootTreeReader |
---|
3 | * |
---|
4 | * Class simplifying access to ROOT tree branches |
---|
5 | * |
---|
6 | * $Date: 2008-06-04 13:57:57 $ |
---|
7 | * $Revision: 1.1 $ |
---|
8 | * |
---|
9 | * |
---|
10 | * \author P. Demin - UCL, Louvain-la-Neuve |
---|
11 | * |
---|
12 | */ |
---|
13 | |
---|
14 | #include "ExRootAnalysis/ExRootTreeReader.h" |
---|
15 | |
---|
16 | #include "TH2.h" |
---|
17 | #include "TStyle.h" |
---|
18 | #include "TCanvas.h" |
---|
19 | #include "TClonesArray.h" |
---|
20 | #include "TBranchElement.h" |
---|
21 | |
---|
22 | #include <iostream> |
---|
23 | |
---|
24 | using namespace std; |
---|
25 | |
---|
26 | //------------------------------------------------------------------------------ |
---|
27 | |
---|
28 | ExRootTreeReader::ExRootTreeReader(TTree *tree) : |
---|
29 | fChain(tree), fCurrentTree(-1) |
---|
30 | { |
---|
31 | } |
---|
32 | |
---|
33 | //------------------------------------------------------------------------------ |
---|
34 | |
---|
35 | ExRootTreeReader::~ExRootTreeReader() |
---|
36 | { |
---|
37 | TBranchMap::iterator itBranchMap; |
---|
38 | |
---|
39 | for(itBranchMap = fBranchMap.begin(); itBranchMap != fBranchMap.end(); ++itBranchMap) |
---|
40 | { |
---|
41 | delete itBranchMap->second.second; |
---|
42 | } |
---|
43 | } |
---|
44 | |
---|
45 | //------------------------------------------------------------------------------ |
---|
46 | |
---|
47 | Bool_t ExRootTreeReader::ReadEntry(Long64_t entry) |
---|
48 | { |
---|
49 | // Read contents of entry. |
---|
50 | if(!fChain) return kFALSE; |
---|
51 | |
---|
52 | Int_t treeEntry = fChain->LoadTree(entry); |
---|
53 | if(treeEntry < 0) return kFALSE; |
---|
54 | |
---|
55 | if(fChain->IsA() == TChain::Class()) |
---|
56 | { |
---|
57 | TChain *chain = static_cast<TChain*>(fChain); |
---|
58 | if(chain->GetTreeNumber() != fCurrentTree) |
---|
59 | { |
---|
60 | fCurrentTree = chain->GetTreeNumber(); |
---|
61 | Notify(); |
---|
62 | } |
---|
63 | } |
---|
64 | |
---|
65 | TBranchMap::iterator itBranchMap; |
---|
66 | TBranch *branch; |
---|
67 | |
---|
68 | for(itBranchMap = fBranchMap.begin(); itBranchMap != fBranchMap.end(); ++itBranchMap) |
---|
69 | { |
---|
70 | branch = itBranchMap->second.first; |
---|
71 | if(branch) |
---|
72 | { |
---|
73 | branch->GetEntry(treeEntry); |
---|
74 | } |
---|
75 | } |
---|
76 | |
---|
77 | return kTRUE; |
---|
78 | } |
---|
79 | |
---|
80 | //------------------------------------------------------------------------------ |
---|
81 | |
---|
82 | TClonesArray *ExRootTreeReader::UseBranch(const char *branchName) |
---|
83 | { |
---|
84 | TClonesArray *array = 0; |
---|
85 | |
---|
86 | TBranchMap::iterator itBranchMap = fBranchMap.find(branchName); |
---|
87 | |
---|
88 | if(itBranchMap != fBranchMap.end()) |
---|
89 | { |
---|
90 | cout << "** WARNING: branch '" << branchName << "' is already in use" << endl; |
---|
91 | array = itBranchMap->second.second; |
---|
92 | } |
---|
93 | else |
---|
94 | { |
---|
95 | TBranch *branch = fChain->GetBranch(branchName); |
---|
96 | if(branch) |
---|
97 | { |
---|
98 | if(branch->IsA() == TBranchElement::Class()) |
---|
99 | { |
---|
100 | TBranchElement *element = static_cast<TBranchElement*>(branch); |
---|
101 | const char *className = element->GetClonesName(); |
---|
102 | Int_t size = element->GetMaximum(); |
---|
103 | TClass *cl = gROOT->GetClass(className); |
---|
104 | if(cl) |
---|
105 | { |
---|
106 | array = new TClonesArray(cl, size); |
---|
107 | array->SetName(branchName); |
---|
108 | fBranchMap.insert(make_pair(branchName, make_pair(branch, array))); |
---|
109 | branch->SetAddress(&array); |
---|
110 | } |
---|
111 | } |
---|
112 | } |
---|
113 | } |
---|
114 | |
---|
115 | if(!array) |
---|
116 | { |
---|
117 | cout << "** WARNING: cannot access branch '" << branchName << "', return NULL pointer" << endl; |
---|
118 | } |
---|
119 | |
---|
120 | return array; |
---|
121 | } |
---|
122 | |
---|
123 | //------------------------------------------------------------------------------ |
---|
124 | |
---|
125 | Bool_t ExRootTreeReader::Notify() |
---|
126 | { |
---|
127 | // Called when loading a new file. |
---|
128 | // Get branch pointers. |
---|
129 | if(!fChain) return kFALSE; |
---|
130 | |
---|
131 | TBranchMap::iterator itBranchMap; |
---|
132 | TBranch *branch; |
---|
133 | |
---|
134 | for(itBranchMap = fBranchMap.begin(); itBranchMap != fBranchMap.end(); ++itBranchMap) |
---|
135 | { |
---|
136 | branch = fChain->GetBranch(itBranchMap->first); |
---|
137 | if(branch) |
---|
138 | { |
---|
139 | itBranchMap->second.first = branch; |
---|
140 | branch->SetAddress(&(itBranchMap->second.second)); |
---|
141 | } |
---|
142 | else |
---|
143 | { |
---|
144 | cout << "** WARNING: cannot get branch '" << itBranchMap->first << "'" << endl; |
---|
145 | } |
---|
146 | } |
---|
147 | return kTRUE; |
---|
148 | } |
---|
149 | |
---|
150 | //------------------------------------------------------------------------------ |
---|
151 | |
---|