| 1 | //
|
|---|
| 2 | // This script is a first version to read, with ROOT, a ELYSE.root
|
|---|
| 3 | // file having the "TTree within TTree" of the ELYSE_batch
|
|---|
| 4 | // progrem (see Analysis.cxx).
|
|---|
| 5 | // But the below way to read the data has clearly a HUGE
|
|---|
| 6 | // memory consumption problem that someone can see by using
|
|---|
| 7 | // the "top" program on a UNIX.
|
|---|
| 8 | //
|
|---|
| 9 | // Usage :
|
|---|
| 10 | // OS> <setup ROOT>
|
|---|
| 11 | // ( OS> root batch_dummy.C )
|
|---|
| 12 | // OS> root analyis_leaks.C
|
|---|
| 13 | //
|
|---|
| 14 | // G.Barrand
|
|---|
| 15 | //
|
|---|
| 16 |
|
|---|
| 17 | void analysis_leaks() {
|
|---|
| 18 |
|
|---|
| 19 | TFile* f = new TFile("ELYSE.root");
|
|---|
| 20 | TTree* tEvent = (TTree*)f->Get("Event");
|
|---|
| 21 |
|
|---|
| 22 | //declare Non subTuple part
|
|---|
| 23 | Int_t eventId,inputEvtId, interMode, vtxVol;
|
|---|
| 24 | tEvent->SetBranchAddress("eventId",&eventId);
|
|---|
| 25 | tEvent->SetBranchAddress("inputEvtId",&inputEvtId);
|
|---|
| 26 | tEvent->SetBranchAddress("interMode",&interMode);
|
|---|
| 27 | tEvent->SetBranchAddress("vtxVol",&vtxVol);
|
|---|
| 28 |
|
|---|
| 29 | Int_t nPart,leptonIndex, protonIndex;
|
|---|
| 30 | tEvent->SetBranchAddress("nPart",&nPart);
|
|---|
| 31 | tEvent->SetBranchAddress("leptonIndex",&leptonIndex);
|
|---|
| 32 | tEvent->SetBranchAddress("protonIndex",&protonIndex);
|
|---|
| 33 |
|
|---|
| 34 | Int_t nHits;
|
|---|
| 35 | tEvent->SetBranchAddress("nHits",&nHits);
|
|---|
| 36 |
|
|---|
| 37 | Int_t nDigits;
|
|---|
| 38 | Double_t sumPE;
|
|---|
| 39 | tEvent->SetBranchAddress("nDigits",&nDigits);
|
|---|
| 40 | tEvent->SetBranchAddress("sumPE",&sumPE);
|
|---|
| 41 |
|
|---|
| 42 | TTree* Event_vtxPos = new TTree();
|
|---|
| 43 | tEvent->SetBranchAddress("vtxPos",&Event_vtxPos);
|
|---|
| 44 | Double_t vx,vy,vz;
|
|---|
| 45 |
|
|---|
| 46 | //The Track Event:subTuple part
|
|---|
| 47 | TTree* Event_track = new TTree();
|
|---|
| 48 | tEvent->SetBranchAddress("track",&Event_track);
|
|---|
| 49 | Int_t pId,parent;
|
|---|
| 50 | Float_t timeStart;
|
|---|
| 51 | Double_t mass, pTot, ETot;
|
|---|
| 52 | Int_t startVol, stopVol;
|
|---|
| 53 |
|
|---|
| 54 | //The Direction Track:subTuple
|
|---|
| 55 | TTree* Track_direction = new TTree();
|
|---|
| 56 | Double_t dx,dy,dz;
|
|---|
| 57 |
|
|---|
| 58 | //The Momentum Track:subTuple
|
|---|
| 59 | TTree* Track_momentum = new TTree();
|
|---|
| 60 | Double_t px,py,pz;
|
|---|
| 61 |
|
|---|
| 62 | //The StartPos Track:subTuple
|
|---|
| 63 | TTree* Track_startPos = new TTree();
|
|---|
| 64 | Double_t start_x,start_y,start_z;
|
|---|
| 65 |
|
|---|
| 66 | //The StopPos Track:subTuple
|
|---|
| 67 | TTree* Track_stopPos = new TTree();
|
|---|
| 68 | Double_t stop_x,stop_y,stop_z;
|
|---|
| 69 |
|
|---|
| 70 | //The Hit Event:subTuple
|
|---|
| 71 | TTree* Event_hit = new TTree();
|
|---|
| 72 | tEvent->SetBranchAddress("hit",&Event_hit);
|
|---|
| 73 | Int_t tubeId_hit; //JEC 16/1/06
|
|---|
| 74 | Int_t totalPE;
|
|---|
| 75 |
|
|---|
| 76 | //The Pe Hit::subTuple
|
|---|
| 77 | TTree* Hit_pe = new TTree();
|
|---|
| 78 | Float_t hit_time;
|
|---|
| 79 | Float_t trk_length; //NV 13/6/06
|
|---|
| 80 |
|
|---|
| 81 | //The Digit Event:subTuple
|
|---|
| 82 | TTree* Event_digit = new TTree();
|
|---|
| 83 | tEvent->SetBranchAddress("digit",&Event_digit);
|
|---|
| 84 | Int_t tubeId;
|
|---|
| 85 | Double_t digit_pe, digit_time;
|
|---|
| 86 |
|
|---|
| 87 | Int_t nEvent = (Int_t)tEvent->GetEntries();
|
|---|
| 88 | std::cout << " nEvents = " << nEvent << endl;
|
|---|
| 89 |
|
|---|
| 90 | //Int_t nEvent = 10000;
|
|---|
| 91 | for (Int_t i=0; i<nEvent; ++i) {
|
|---|
| 92 |
|
|---|
| 93 | /*
|
|---|
| 94 | Event_vtxPos->Reset();
|
|---|
| 95 | Event_track->Reset();
|
|---|
| 96 | Track_direction->Reset();
|
|---|
| 97 | Track_momentum->Reset();
|
|---|
| 98 | Track_startPos->Reset();
|
|---|
| 99 | Track_stopPos->Reset();
|
|---|
| 100 | Event_hit->Reset();
|
|---|
| 101 | Hit_pe->Reset();
|
|---|
| 102 | Event_digit->Reset();
|
|---|
| 103 | //tEvent->ResetBranchAddresses();
|
|---|
| 104 | */
|
|---|
| 105 |
|
|---|
| 106 | tEvent->GetEntry(i);
|
|---|
| 107 | std::cout << ">>>>>>>>>>>>> Event{" << i << "}: "
|
|---|
| 108 | << " evt Id " << eventId
|
|---|
| 109 | // << " evt Input Id " << inputEvtId
|
|---|
| 110 | // << "\n interaction mode " << interMode
|
|---|
| 111 | // << " start in volume " << vtxVol << "\n"
|
|---|
| 112 | // <<" #tracks: " << nPart
|
|---|
| 113 | // <<" #hits: " << nHits
|
|---|
| 114 | // <<" #digits: " << nDigits
|
|---|
| 115 | << std::endl;
|
|---|
| 116 |
|
|---|
| 117 | Int_t nTracks = (Int_t)Event_track->GetEntries();
|
|---|
| 118 | Int_t nTubeHits = (Int_t)Event_hit->GetEntries();
|
|---|
| 119 | Int_t nTubeDigits = (Int_t)Event_digit->GetEntries();
|
|---|
| 120 | //std::cout << "Verif: nTracks = " << nTracks
|
|---|
| 121 | // << " nTube Hits = " << nTubeHits
|
|---|
| 122 | // << " nTube Digits = " << nTubeDigits
|
|---|
| 123 | // << std::endl;
|
|---|
| 124 |
|
|---|
| 125 | Event_vtxPos->SetBranchAddress("x",&vx);
|
|---|
| 126 | Event_vtxPos->SetBranchAddress("y",&vy);
|
|---|
| 127 | Event_vtxPos->SetBranchAddress("z",&vz);
|
|---|
| 128 |
|
|---|
| 129 | // Have a brand new overwritten track TTree ; we have
|
|---|
| 130 | // to rebind its user variables :
|
|---|
| 131 | Event_track->SetBranchAddress("pId",&pId);
|
|---|
| 132 | Event_track->SetBranchAddress("parent",&parent);
|
|---|
| 133 | Event_track->SetBranchAddress("timeStart",&timeStart);
|
|---|
| 134 |
|
|---|
| 135 | Event_track->SetBranchAddress("mass",&mass);
|
|---|
| 136 | Event_track->SetBranchAddress("pTot",&pTot);
|
|---|
| 137 | Event_track->SetBranchAddress("ETot",&ETot);
|
|---|
| 138 |
|
|---|
| 139 | Event_track->SetBranchAddress("startVol",&startVol);
|
|---|
| 140 | Event_track->SetBranchAddress("stopVol",&stopVol);
|
|---|
| 141 |
|
|---|
| 142 | Event_track->SetBranchAddress("direction",&Track_direction);
|
|---|
| 143 | Event_track->SetBranchAddress("momentum",&Track_momentum);
|
|---|
| 144 | Event_track->SetBranchAddress("startPos",&Track_startPos);
|
|---|
| 145 | Event_track->SetBranchAddress("stopPos",&Track_stopPos);
|
|---|
| 146 |
|
|---|
| 147 | // Have a brand new overwritten hit TTree ; we have
|
|---|
| 148 | // to rebind its user variables :
|
|---|
| 149 | Event_hit->SetBranchAddress("tubeId",&tubeId_hit); //JEC 16/1/06
|
|---|
| 150 | Event_hit->SetBranchAddress("totalPE",&totalPE);
|
|---|
| 151 | Event_hit->SetBranchAddress("pe",&Hit_pe);
|
|---|
| 152 |
|
|---|
| 153 | // Have a brand new overwritten digit TTree ; we have
|
|---|
| 154 | // to rebind its user variables :
|
|---|
| 155 | Event_digit->SetBranchAddress("tubeId",&tubeId);
|
|---|
| 156 | Event_digit->SetBranchAddress("pe",&digit_pe);
|
|---|
| 157 | Event_digit->SetBranchAddress("time",&digit_time);
|
|---|
| 158 |
|
|---|
| 159 |
|
|---|
| 160 | for (Int_t j=0; j<nTracks; ++j) {
|
|---|
| 161 | //for (Int_t j=0; j<0; ++j) {
|
|---|
| 162 | Event_track->GetEntry(j);
|
|---|
| 163 |
|
|---|
| 164 | Track_direction->SetBranchAddress("dx",&dx);
|
|---|
| 165 | Track_direction->SetBranchAddress("dy",&dy);
|
|---|
| 166 | Track_direction->SetBranchAddress("dz",&dz);
|
|---|
| 167 | Track_direction->GetEntry(0);
|
|---|
| 168 |
|
|---|
| 169 | Track_momentum->SetBranchAddress("px",&px);
|
|---|
| 170 | Track_momentum->SetBranchAddress("py",&py);
|
|---|
| 171 | Track_momentum->SetBranchAddress("pz",&pz);
|
|---|
| 172 | Track_momentum->GetEntry(0);
|
|---|
| 173 |
|
|---|
| 174 | Track_startPos->SetBranchAddress("x",&start_x);
|
|---|
| 175 | Track_startPos->SetBranchAddress("y",&start_y);
|
|---|
| 176 | Track_startPos->SetBranchAddress("z",&start_z);
|
|---|
| 177 | Track_startPos->GetEntry(0);
|
|---|
| 178 |
|
|---|
| 179 | Track_stopPos->SetBranchAddress("x",&stop_x);
|
|---|
| 180 | Track_stopPos->SetBranchAddress("y",&stop_y);
|
|---|
| 181 | Track_stopPos->SetBranchAddress("z",&stop_z);
|
|---|
| 182 | Track_stopPos->GetEntry(0);
|
|---|
| 183 |
|
|---|
| 184 | std::cout << "----> Tk{"<<j<<"}: "
|
|---|
| 185 | << " pId " << pId
|
|---|
| 186 | << " parent " << parent
|
|---|
| 187 | << " creation time " << timeStart
|
|---|
| 188 | << " Volumes " << startVol << " " << stopVol << "\n"
|
|---|
| 189 | << " Start Pos (" << start_x << "," << start_y << "," << start_z << ")\n"
|
|---|
| 190 | << " Stop Pos (" << stop_x << "," << stop_y << "," << stop_z << ")\n"
|
|---|
| 191 | << " dx,dy,dz " << dx << " " << dy << " " << dz << "\n"
|
|---|
| 192 | << " m " << mass
|
|---|
| 193 | << " ETot " << ETot
|
|---|
| 194 | << " pTot " << pTot
|
|---|
| 195 | << " px,py,pz " << px << " " << py << " " << pz << "\n"
|
|---|
| 196 | << std::endl;
|
|---|
| 197 | }//loop on Tracks
|
|---|
| 198 |
|
|---|
| 199 | //--------
|
|---|
| 200 | // The Hits
|
|---|
| 201 | //--------
|
|---|
| 202 |
|
|---|
| 203 |
|
|---|
| 204 | for (Int_t k=0; k<nTubeHits; ++k) {
|
|---|
| 205 |
|
|---|
| 206 | Event_hit->GetEntry(k);
|
|---|
| 207 |
|
|---|
| 208 | Hit_pe->SetBranchAddress("time",&hit_time);
|
|---|
| 209 | //Hit_pe->SetBranchAddress("length",&trk_length);
|
|---|
| 210 | //JEC 16/1/06 add the tubeId_hit info
|
|---|
| 211 | std::cout << "----> Hit{"<<k<<"}: tube[" << tubeId_hit << "] total #PE " << totalPE << std::endl;
|
|---|
| 212 | for (Int_t ki=0; ki<Hit_pe->GetEntries(); ++ki) {
|
|---|
| 213 | Hit_pe->GetEntry(ki);
|
|---|
| 214 | std::cout << "<" << hit_time << ">" << "<" << trk_length << ">" ;
|
|---|
| 215 | }
|
|---|
| 216 | std::cout << std::endl;
|
|---|
| 217 | }//Loop on Hits
|
|---|
| 218 |
|
|---|
| 219 | //--------
|
|---|
| 220 | // The Digits
|
|---|
| 221 | //--------
|
|---|
| 222 | //for (Int_t l=0; l<nTubeDigits; ++l)
|
|---|
| 223 | for (Int_t l=0; l<0; ++l) {
|
|---|
| 224 | Event_digit->GetEntry(l);
|
|---|
| 225 |
|
|---|
| 226 | std::cout << "----> Digit{"<<l<<"}: "
|
|---|
| 227 | << "tube[" << tubeId << "] = "
|
|---|
| 228 | << " pe: " << digit_pe
|
|---|
| 229 | << " time: " << digit_time
|
|---|
| 230 | << std::endl;
|
|---|
| 231 |
|
|---|
| 232 | }//Loop on Digits
|
|---|
| 233 | }//loop on event
|
|---|
| 234 | }//read_event
|
|---|
