1 |
|
// $Id$ |
2 |
|
|
3 |
|
#include "MitAna/DataUtil/interface/TreeWriter.h" |
4 |
– |
|
4 |
|
#include <Riostream.h> |
5 |
|
#include <TObject.h> |
6 |
|
#include <TSystem.h> |
7 |
|
#include <TProcessID.h> |
8 |
< |
|
8 |
> |
#include <TBranchRef.h> |
9 |
|
#include "MitAna/DataUtil/interface/Debug.h" |
10 |
|
|
11 |
|
using namespace mithep; |
27 |
|
fEvtObjNum(-1), |
28 |
|
fIsInit(kFALSE), |
29 |
|
fDoObjNumReset(doreset), |
30 |
+ |
fDoBranchRef(0), |
31 |
|
fFile(0), |
32 |
|
fTrees(0) |
33 |
|
{ |
56 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
57 |
|
// to object pointer for class name "cname" using a given buffer size and splitlevel. |
58 |
|
|
59 |
< |
MyTree *t = AddOrGetMyTree(GetName()); |
60 |
< |
t->Bronch(name, cname, obj, bsize, level); |
59 |
> |
MyTree *t = AddOrGetMyTree(GetName()); |
60 |
> |
TBranch *b = t->Bronch(name, cname, obj, bsize, level); |
61 |
> |
b->SetCompressionLevel(GetCompressLevel()); |
62 |
|
} |
63 |
|
|
64 |
|
|
65 |
< |
//__________________________________________________________________________________________________ |
65 |
> |
//-------------------------------------------------------------------------------------------------- |
66 |
|
void TreeWriter::AddBranch(const char *name, void *obj, Int_t bsize, Int_t level) |
67 |
|
{ |
68 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
71 |
|
AddBranch(name, CName(obj), obj, bsize, level); |
72 |
|
} |
73 |
|
|
74 |
< |
//__________________________________________________________________________________________________ |
74 |
> |
//-------------------------------------------------------------------------------------------------- |
75 |
|
void TreeWriter::AddBranch(const char *name, const char *cname, |
76 |
|
void *obj, Int_t bsize) |
77 |
|
{ |
78 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
79 |
|
// to object pointer for class name "cname" using a given buffer size and default splitlevel. |
80 |
|
|
81 |
< |
MyTree *t = AddOrGetMyTree(GetName()); |
82 |
< |
t->Bronch(name, cname, obj, bsize, fDefSL); |
81 |
> |
MyTree *t = AddOrGetMyTree(GetName()); |
82 |
> |
TBranch *b = t->Bronch(name, cname, obj, bsize, fDefSL); |
83 |
> |
b->SetCompressionLevel(GetCompressLevel()); |
84 |
|
} |
85 |
|
|
86 |
< |
//__________________________________________________________________________________________________ |
86 |
> |
//-------------------------------------------------------------------------------------------------- |
87 |
|
void TreeWriter::AddBranch(const char *name, void *obj, Int_t bsize) |
88 |
|
{ |
89 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
92 |
|
AddBranch(name, CName(obj), obj, bsize); |
93 |
|
} |
94 |
|
|
95 |
< |
//__________________________________________________________________________________________________ |
95 |
> |
//-------------------------------------------------------------------------------------------------- |
96 |
|
void TreeWriter::AddBranch(const char *name, const char *cname, |
97 |
|
void *obj) |
98 |
|
{ |
99 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
100 |
|
// to object pointer for class name "cname" using a default buffer size and splitlevel. |
101 |
|
|
102 |
< |
MyTree *t = AddOrGetMyTree(GetName()); |
103 |
< |
t->Bronch(name, cname, obj, fDefBrSize, fDefSL); |
102 |
> |
MyTree *t = AddOrGetMyTree(GetName()); |
103 |
> |
TBranch *b = t->Bronch(name, cname, obj, fDefBrSize, fDefSL); |
104 |
> |
b->SetCompressionLevel(GetCompressLevel()); |
105 |
|
} |
106 |
|
|
107 |
< |
//__________________________________________________________________________________________________ |
107 |
> |
//-------------------------------------------------------------------------------------------------- |
108 |
|
void TreeWriter::AddBranch(const char *name, void *obj) |
109 |
|
{ |
110 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
113 |
|
AddBranch(name, CName(obj), obj); |
114 |
|
} |
115 |
|
|
116 |
< |
//__________________________________________________________________________________________________ |
116 |
> |
//-------------------------------------------------------------------------------------------------- |
117 |
|
void TreeWriter::AddBranchToTree(const char *tname, const char *name, const char *cname, |
118 |
|
void *obj, Int_t bsize, Int_t level) |
119 |
|
{ |
120 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
121 |
|
// to object pointer for class name "cname" using a given buffer size and splitlevel. |
122 |
|
|
123 |
< |
MyTree *t = AddOrGetMyTree(tname); |
124 |
< |
t->Bronch(name, cname, obj, bsize, level); |
123 |
> |
MyTree *t = AddOrGetMyTree(tname); |
124 |
> |
TBranch *b = t->Bronch(name, cname, obj, bsize, level); |
125 |
> |
b->SetCompressionLevel(GetCompressLevel()); |
126 |
|
} |
127 |
|
|
128 |
< |
//__________________________________________________________________________________________________ |
128 |
> |
//-------------------------------------------------------------------------------------------------- |
129 |
|
void TreeWriter::AddBranchToTree(const char *tname, const char *name, void *obj, |
130 |
|
Int_t bsize, Int_t level) |
131 |
|
{ |
135 |
|
AddBranchToTree(tname, name, CName(obj), obj, bsize, level); |
136 |
|
} |
137 |
|
|
138 |
< |
//__________________________________________________________________________________________________ |
138 |
> |
//-------------------------------------------------------------------------------------------------- |
139 |
|
void TreeWriter::AddBranchToTree(const char *tname, const char *name, const char *cname, |
140 |
|
void *obj, Int_t bsize) |
141 |
|
{ |
142 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
143 |
|
// to object pointer for class name "cname" using a given buffer size and default splitlevel. |
144 |
|
|
145 |
< |
MyTree *t = AddOrGetMyTree(tname); |
146 |
< |
t->Bronch(name, cname, obj, bsize, fDefSL); |
145 |
> |
MyTree *t = AddOrGetMyTree(tname); |
146 |
> |
TBranch *b = t->Bronch(name, cname, obj, bsize, fDefSL); |
147 |
> |
b->SetCompressionLevel(GetCompressLevel()); |
148 |
|
} |
149 |
|
|
150 |
< |
//__________________________________________________________________________________________________ |
150 |
> |
//-------------------------------------------------------------------------------------------------- |
151 |
|
void TreeWriter::AddBranchToTree(const char *tname, const char *name, void *obj, |
152 |
|
Int_t bsize) |
153 |
|
{ |
157 |
|
AddBranchToTree(tname, name, CName(obj), obj, bsize); |
158 |
|
} |
159 |
|
|
160 |
< |
//__________________________________________________________________________________________________ |
160 |
> |
//-------------------------------------------------------------------------------------------------- |
161 |
|
void TreeWriter::AddBranchToTree(const char *tname, const char *name, const char *cname, |
162 |
|
void *obj) |
163 |
|
{ |
164 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
165 |
|
// to object pointer for class name "cname" using a default buffer size and splitlevel. |
166 |
|
|
167 |
< |
MyTree *t = AddOrGetMyTree(tname); |
168 |
< |
t->Bronch(name, cname, obj, fDefBrSize, fDefSL); |
167 |
> |
MyTree *t = AddOrGetMyTree(tname); |
168 |
> |
TBranch *b = t->Bronch(name, cname, obj, fDefBrSize, fDefSL); |
169 |
> |
b->SetCompressionLevel(GetCompressLevel()); |
170 |
|
} |
171 |
|
|
172 |
< |
//__________________________________________________________________________________________________ |
172 |
> |
//-------------------------------------------------------------------------------------------------- |
173 |
|
void TreeWriter::AddBranchToTree(const char *tname, const char *name, void *obj) |
174 |
|
{ |
175 |
|
// Add branch with name "name" into tree with name "tname" and set its address |
178 |
|
AddBranchToTree(tname, name, CName(obj), obj); |
179 |
|
} |
180 |
|
|
181 |
< |
//__________________________________________________________________________________________________ |
181 |
> |
//-------------------------------------------------------------------------------------------------- |
182 |
|
MyTree *TreeWriter::AddOrGetMyTree(const char *tn) |
183 |
|
{ |
184 |
|
// Add new tree if not present in array of trees or return |
185 |
|
// present tree. |
186 |
|
|
187 |
< |
MyTree *tree = dynamic_cast<MyTree*>(fTrees.FindObject(tn)); |
187 |
> |
MyTree *tree = GetMyTree(tn); //dynamic_cast<MyTree*>(fTrees.FindObject(tn)); |
188 |
|
if (tree) |
189 |
|
return tree; |
190 |
|
|
191 |
|
TDirectory::TContext context(fFile); |
192 |
|
tree = new MyTree(tn, tn); |
193 |
|
tree->SetDirectory(fFile); |
194 |
+ |
if (fDoBranchRef) |
195 |
+ |
tree->BranchRef(); |
196 |
|
fTrees.AddLast(tree); |
197 |
|
return tree; |
198 |
|
} |
199 |
|
|
200 |
< |
//__________________________________________________________________________________________________ |
200 |
> |
//-------------------------------------------------------------------------------------------------- |
201 |
|
Bool_t TreeWriter::BeginEvent(Bool_t doreset) |
202 |
|
{ |
203 |
|
// Prepare for the next event. If doreset or fDoObjNumReset is kTRUE |
214 |
|
return kTRUE; |
215 |
|
} |
216 |
|
|
217 |
< |
//__________________________________________________________________________________________________ |
217 |
> |
//-------------------------------------------------------------------------------------------------- |
218 |
|
void TreeWriter::CloseFile() |
219 |
|
{ |
220 |
|
// Write tree(s) and close file. |
230 |
|
for (Int_t i=0;i<fTrees.GetEntries();++i) { |
231 |
|
MyTree *mt = static_cast<MyTree*>(fTrees.At(i)); |
232 |
|
mt->Write(mt->GetName(),TObject::kOverwrite); |
233 |
< |
mt->Reset(); |
233 |
> |
|
234 |
> |
// backup and restore list of branch pointers from TRefTable (needed for autoloading) |
235 |
> |
if (mt->GetBranchRef()) { |
236 |
> |
TObjArray *parents = mt->GetBranchRef()->GetRefTable()->GetParents(); |
237 |
> |
TObjArray parentsBak(*parents); |
238 |
> |
mt->Reset(); |
239 |
> |
for (Int_t j=0; j<parentsBak.GetEntries(); ++j) |
240 |
> |
parents->Add(parentsBak.At(j)); |
241 |
> |
} |
242 |
> |
else |
243 |
> |
mt->Reset(); |
244 |
|
mt->SetDirectory(0); |
245 |
|
} |
246 |
|
|
252 |
|
fFileNumber++; |
253 |
|
} |
254 |
|
|
255 |
< |
//__________________________________________________________________________________________________ |
255 |
> |
//-------------------------------------------------------------------------------------------------- |
256 |
|
const char *TreeWriter::CName(void *obj) const |
257 |
|
{ |
258 |
|
// Dereference void* pointer into TObject* pointer |
264 |
|
return tobj->ClassName(); |
265 |
|
} |
266 |
|
|
267 |
< |
//__________________________________________________________________________________________________ |
267 |
> |
//-------------------------------------------------------------------------------------------------- |
268 |
|
Bool_t TreeWriter::EndEvent(Bool_t doreset) |
269 |
|
{ |
270 |
|
// Store the event in the tree. If doreset or fDoObjNumReset is kTRUE |
291 |
|
if (fEvtObjNum<0) { |
292 |
|
Error("EndEvent", "Object counter is zero. Did you call BeginEvent(kTRUE)?"); |
293 |
|
} else { |
294 |
< |
// Reset the TRef table. keep it from growing with each event (see doc) |
294 |
> |
// Reset the TRef table. Keep it from growing with each event (see doc) |
295 |
|
TProcessID::SetObjectCount(fEvtObjNum); |
296 |
|
} |
297 |
|
} |
299 |
|
return (r >= 0); |
300 |
|
} |
301 |
|
|
302 |
< |
//------------------------------------------------------------------------------------------------- |
302 |
> |
//-------------------------------------------------------------------------------------------------- |
303 |
|
Long64_t TreeWriter::GetEntries(const char *tn) const |
304 |
|
{ |
305 |
|
// Return entries of tree with given name. If no tree is given, return sum of entries |
323 |
|
return ret; |
324 |
|
} |
325 |
|
|
326 |
< |
//------------------------------------------------------------------------------------------------- |
326 |
> |
//-------------------------------------------------------------------------------------------------- |
327 |
|
MyTree *mithep::TreeWriter::GetMyTree(const char *tn) |
328 |
|
{ |
329 |
|
// Return MyTree with given name from array. |
343 |
|
return 0; |
344 |
|
} |
345 |
|
|
346 |
< |
//------------------------------------------------------------------------------------------------- |
346 |
> |
//-------------------------------------------------------------------------------------------------- |
347 |
|
const TTree *mithep::TreeWriter::GetTree(const char *tn) const |
348 |
|
{ |
349 |
|
// Return TTree with given name from array. |
363 |
|
return 0; |
364 |
|
} |
365 |
|
|
366 |
< |
//------------------------------------------------------------------------------------------------- |
366 |
> |
//-------------------------------------------------------------------------------------------------- |
367 |
|
TTree *mithep::TreeWriter::GetTree(const char *tn) |
368 |
|
{ |
369 |
|
// Return TTree with given name from array. |
383 |
|
return 0; |
384 |
|
} |
385 |
|
|
386 |
< |
//__________________________________________________________________________________________________ |
386 |
> |
//-------------------------------------------------------------------------------------------------- |
387 |
|
Bool_t TreeWriter::IsFull() const |
388 |
|
{ |
389 |
|
// Check if the maximum file size has been reached. |
399 |
|
return (GetFileSize() + avgSize + fkMinFreeSpace) > fMaxSize; |
400 |
|
} |
401 |
|
|
402 |
< |
//__________________________________________________________________________________________________ |
402 |
> |
//-------------------------------------------------------------------------------------------------- |
403 |
|
void TreeWriter::OpenFile() |
404 |
|
{ |
405 |
|
// Open the file and attach the tree. |
430 |
|
fIsInit = kTRUE; |
431 |
|
} |
432 |
|
|
433 |
< |
//__________________________________________________________________________________________________ |
433 |
> |
//-------------------------------------------------------------------------------------------------- |
434 |
|
void TreeWriter::Print(Option_t *option) const |
435 |
|
{ |
436 |
|
// Print the contents of the tree writer. |
452 |
|
<< (GetEntries() == 1 ? " event" : " events") << endl; |
453 |
|
} |
454 |
|
|
455 |
< |
//------------------------------------------------------------------------------------------------- |
455 |
> |
//-------------------------------------------------------------------------------------------------- |
456 |
|
void TreeWriter::SetAutoFill(const char *tn, Bool_t b) |
457 |
|
{ |
458 |
|
// Set auto-fill mode of tree with given name. |
467 |
|
mt->SetAutoFill(b); |
468 |
|
} |
469 |
|
|
470 |
< |
|
453 |
< |
//------------------------------------------------------------------------------------------------- |
470 |
> |
//-------------------------------------------------------------------------------------------------- |
471 |
|
void TreeWriter::SetMaxSize(Long64_t s) |
472 |
|
{ |
473 |
|
// Set maximum file size. Check if this exceeds the ROOT file size and if, |
483 |
|
fMaxSize=s; |
484 |
|
} |
485 |
|
|
486 |
< |
//__________________________________________________________________________________________________ |
486 |
> |
//-------------------------------------------------------------------------------------------------- |
487 |
|
void TreeWriter::StoreObject(const TObject *obj) |
488 |
|
{ |
489 |
|
// Store object next to tree in file. Used to store the |