[5] | 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 | |
---|