7 |
* Default constructor |
* Default constructor |
8 |
*/ |
*/ |
9 |
PamTrack::PamTrack(){ |
PamTrack::PamTrack(){ |
10 |
trk_track = this->TrkTrack::GetTrkTrack(); |
trk_track = 0; |
11 |
calo_track = this->CaloTrkVar::GetCaloTrkVar(); |
calo_track = 0; |
12 |
tof_track = this->ToFTrkVar::GetToFTrkVar(); |
tof_track = 0; |
13 |
|
candeleteobj = 0; |
14 |
}; |
}; |
15 |
//-------------------------------------- |
//-------------------------------------- |
16 |
// |
// |
20 |
* Constructor |
* Constructor |
21 |
*/ |
*/ |
22 |
PamTrack::PamTrack(TrkTrack* t, CaloTrkVar* c, ToFTrkVar* o){ |
PamTrack::PamTrack(TrkTrack* t, CaloTrkVar* c, ToFTrkVar* o){ |
23 |
trk_track = this->TrkTrack::GetTrkTrack(); |
|
24 |
calo_track = this->CaloTrkVar::GetCaloTrkVar(); |
trk_track = 0; |
25 |
tof_track = this->ToFTrkVar::GetToFTrkVar(); |
calo_track = 0; |
26 |
if(t) *trk_track = *t; |
tof_track = 0; |
27 |
if(c) *calo_track = *c; |
// if(t)trk_track = new TrkTrack(*t); |
28 |
if(o) *tof_track = *o; |
// if(c)calo_track = new CaloTrkVar(*c); |
29 |
|
// if(o)tof_track = new ToFTrkVar(*o); |
30 |
|
if(t)trk_track = t; |
31 |
|
if(c)calo_track = c; |
32 |
|
if(o)tof_track = o; |
33 |
|
candeleteobj = 0; |
34 |
|
|
35 |
}; |
}; |
36 |
|
PamTrack::PamTrack(const PamTrack& track){ |
37 |
|
|
38 |
|
TrkTrack *t = track.trk_track; |
39 |
|
CaloTrkVar *c = track.calo_track; |
40 |
|
ToFTrkVar *o = track.tof_track; |
41 |
|
|
42 |
|
trk_track = 0; |
43 |
|
calo_track = 0; |
44 |
|
tof_track = 0; |
45 |
|
if(t)trk_track = new TrkTrack(*t); |
46 |
|
if(c)calo_track = new CaloTrkVar(*c); |
47 |
|
if(o)tof_track = new ToFTrkVar(*o); |
48 |
|
candeleteobj = 1; |
49 |
|
|
50 |
|
} |
51 |
|
void PamTrack::Clear(){ |
52 |
|
|
53 |
|
// cout << "PamTrack::Clear() "<<candeleteobj<<endl; |
54 |
|
if(candeleteobj){ |
55 |
|
if(trk_track) trk_track->TrkTrack::Clear(); |
56 |
|
if(calo_track) calo_track->CaloTrkVar::Clear();//??? |
57 |
|
if(tof_track) tof_track->ToFTrkVar::Clear();//??? |
58 |
|
}else{ |
59 |
|
trk_track = 0; |
60 |
|
calo_track = 0; |
61 |
|
tof_track = 0; |
62 |
|
} |
63 |
|
} |
64 |
|
void PamTrack::Delete(){ |
65 |
|
// cout << "PamTrack::Delete() "<<candeleteobj<<endl; |
66 |
|
if(candeleteobj){ |
67 |
|
if(trk_track) { |
68 |
|
trk_track->TrkTrack::Clear(); |
69 |
|
delete trk_track; |
70 |
|
} |
71 |
|
if(calo_track){ |
72 |
|
calo_track->CaloTrkVar::Clear();//??? |
73 |
|
delete calo_track; |
74 |
|
} |
75 |
|
if(tof_track){ |
76 |
|
tof_track->ToFTrkVar::Clear();//??? |
77 |
|
delete tof_track; |
78 |
|
} |
79 |
|
}else{ |
80 |
|
Clear(); |
81 |
|
} |
82 |
|
} |
83 |
//-------------------------------------- |
//-------------------------------------- |
84 |
// |
// |
85 |
// |
// |
89 |
*/ |
*/ |
90 |
PamLevel2::PamLevel2(){ |
PamLevel2::PamLevel2(){ |
91 |
|
|
92 |
trk2_obj = TrkLevel2::GetTrkLevel2(); |
// trk2_obj = TrkLevel2::GetTrkLevel2(); |
93 |
trk1_obj = TrkLevel1::GetTrkLevel1(); |
// trk1_obj = TrkLevel1::GetTrkLevel1(); |
94 |
trkh_obj = TrkHough::GetTrkHough(); |
// trkh_obj = TrkHough::GetTrkHough(); |
95 |
calo1_obj = CaloLevel1::GetCaloLevel1(); |
// calo1_obj = CaloLevel1::GetCaloLevel1(); |
96 |
calo2_obj = CaloLevel2::GetCaloLevel2(); |
// calo2_obj = CaloLevel2::GetCaloLevel2(); |
97 |
tof_obj = ToFLevel2::GetToFLevel2(); |
// tof_obj = ToFLevel2::GetToFLevel2(); |
98 |
trig_obj = TrigLevel2::GetTrigLevel2(); |
// trig_obj = TrigLevel2::GetTrigLevel2(); |
99 |
s4_obj = S4Level2::GetS4Level2(); |
// s4_obj = S4Level2::GetS4Level2(); |
100 |
nd_obj = NDLevel2::GetNDLevel2(); |
// nd_obj = NDLevel2::GetNDLevel2(); |
101 |
ac_obj = AcLevel2::GetAcLevel2(); |
// ac_obj = AcLevel2::GetAcLevel2(); |
102 |
orb_obj = OrbitalInfo::GetOrbitalInfo(); |
// orb_obj = OrbitalInfo::GetOrbitalInfo(); |
103 |
|
trk2_obj = 0; |
104 |
run_obj = new GL_RUN(); |
trk1_obj = 0; |
105 |
|
trkh_obj = 0; |
106 |
// trk1_clone = 0; |
calo1_obj = 0; |
107 |
// trk2_clone = 0; |
calo2_obj = 0; |
108 |
// trkh_clone = 0; |
tof_obj = 0; |
109 |
// calo1_clone = 0; |
trig_obj = 0; |
110 |
// calo2_clone = 0; |
s4_obj = 0; |
111 |
// tof_clone = 0; |
nd_obj = 0; |
112 |
// trig_clone = 0; |
ac_obj = 0; |
113 |
// s4_clone = 0; |
orb_obj = 0; |
114 |
// nd_clone = 0; |
|
115 |
// ac_clone = 0; |
run_obj = 0;//new GL_RUN(); |
116 |
// orb_clone = 0; |
|
|
|
|
|
// T_clone = 0; |
|
|
// C_clone = 0; |
|
|
// O_clone = 0; |
|
|
// R_clone = 0; |
|
|
// S_clone = 0; |
|
|
// N_clone = 0; |
|
|
// A_clone = 0; |
|
|
// B_clone = 0; |
|
117 |
|
|
118 |
for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL; |
for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL; |
119 |
|
|
120 |
sorted_tracks = new TRefArray(); |
// sorted_tracks = 0;//new TRefArray(); |
121 |
|
|
122 |
CAL1 = true; |
CAL1 = true; |
123 |
CAL2 = true; |
CAL2 = true; |
131 |
AC = true; |
AC = true; |
132 |
ORB = true; |
ORB = true; |
133 |
|
|
134 |
|
RUN = true; |
135 |
|
|
136 |
|
tsorted=0; |
137 |
|
timage=0; |
138 |
|
|
139 |
}; |
}; |
140 |
/** |
/** |
141 |
* Destructor |
* Delete the event (act as Dtor) |
142 |
*/ |
*/ |
143 |
PamLevel2::~PamLevel2(){ |
void PamLevel2::Delete(){ |
144 |
delete run_obj; |
|
145 |
Clear(); |
if(run_obj)delete run_obj; |
146 |
|
|
147 |
|
// cout << "void PamLevel2::Clear()"<<endl; |
148 |
|
if(trk1_obj) delete trk1_obj; |
149 |
|
if(trk2_obj) delete trk2_obj; |
150 |
|
if(trkh_obj) delete trkh_obj; |
151 |
|
if(calo1_obj)delete calo1_obj; |
152 |
|
if(calo2_obj)delete calo2_obj; |
153 |
|
if(tof_obj) delete tof_obj; |
154 |
|
if(trig_obj) delete trig_obj; |
155 |
|
if(s4_obj) delete s4_obj; |
156 |
|
if(nd_obj) delete nd_obj; |
157 |
|
if(ac_obj) delete ac_obj; |
158 |
|
if(orb_obj) delete orb_obj; |
159 |
|
|
160 |
|
if(tsorted){ |
161 |
|
tsorted->Delete(); |
162 |
|
delete tsorted; |
163 |
|
} |
164 |
|
if(timage){ |
165 |
|
timage->Delete(); |
166 |
|
delete timage; |
167 |
|
} |
168 |
}; |
}; |
169 |
|
|
170 |
/** |
/** |
171 |
* Clear the event |
* Clear the event (NB! does not deallocate objects) |
172 |
*/ |
*/ |
173 |
void PamLevel2::Clear(){ |
void PamLevel2::Clear(){ |
174 |
|
|
175 |
TrkLevel1::Clear(); |
// cout << "void PamLevel2::Clear()"<<endl; |
176 |
|
if(run_obj)delete run_obj; |
177 |
TrkLevel2::Clear(); |
|
178 |
CaloLevel2::Clear(); |
if(trk1_obj) trk1_obj->Clear(); |
179 |
ToFLevel2::Clear(); |
if(trk2_obj) trk2_obj->Clear(); |
180 |
TrigLevel2::Clear(); |
if(trkh_obj) trkh_obj->Clear(); |
181 |
S4Level2::Clear(); |
if(calo1_obj)calo1_obj->Clear(); |
182 |
NDLevel2::Clear(); |
if(calo2_obj)calo2_obj->Clear(); |
183 |
AcLevel2::Clear(); |
if(tof_obj) tof_obj->Clear(); |
184 |
OrbitalInfo::Clear(); |
if(trig_obj) trig_obj->Clear(); |
185 |
|
if(s4_obj) s4_obj->Clear(); |
186 |
// if(trk1_clone)trk1_clone->Clear(); |
if(nd_obj) nd_obj->Clear(); |
187 |
// if(trk2_clone)trk2_clone->Clear(); |
if(ac_obj) ac_obj->Clear(); |
188 |
// if(trkh_clone)trkh_clone->Clear(); |
if(orb_obj) orb_obj->Clear(); |
|
// if(calo1_clone)calo1_clone->Clear(); |
|
|
// if(calo2_clone)calo2_clone->Clear(); |
|
|
// if(tof_clone)tof_clone->Clear(); |
|
|
// if(trig_clone)trig_clone->Clear(); |
|
|
// if(s4_clone)s4_clone->Clear(); |
|
|
// if(nd_clone)nd_clone->Clear(); |
|
|
// if(ac_clone)ac_clone->Clear(); |
|
|
// if(orb_clone)orb_clone->Clear(); |
|
189 |
|
|
190 |
sorted_tracks->Delete(); // clean the reference array |
// if(sorted_tracks)sorted_tracks->Clear(); |
191 |
|
// sorted_tracks.Clear(); |
192 |
|
|
193 |
|
if(tsorted){ |
194 |
|
tsorted->Delete(); |
195 |
|
} |
196 |
|
if(timage){ |
197 |
|
timage->Delete(); |
198 |
|
} |
199 |
|
|
200 |
}; |
}; |
201 |
|
|
207 |
void *PamLevel2::GetPointerTo(const char* c ){ |
void *PamLevel2::GetPointerTo(const char* c ){ |
208 |
|
|
209 |
TString objname = c; |
TString objname = c; |
210 |
if(!objname.CompareTo("TrkLevel1"))return &trk1_obj; |
|
211 |
if(!objname.CompareTo("TrkLevel2"))return &trk2_obj; |
if(!objname.CompareTo("TrkLevel1")) { |
212 |
if(!objname.CompareTo("TrkHough"))return &trkh_obj; |
if(!trk1_obj){ |
213 |
if(!objname.CompareTo("CaloLevel1"))return &calo1_obj; |
trk1_obj = new TrkLevel1(); |
214 |
if(!objname.CompareTo("CaloLevel2"))return &calo2_obj; |
trk1_obj->Set(); |
215 |
if(!objname.CompareTo("ToFLevel2"))return &tof_obj; |
} |
216 |
if(!objname.CompareTo("TrigLevel2"))return &trig_obj; |
return &trk1_obj; |
217 |
if(!objname.CompareTo("S4Level2"))return &s4_obj; |
}; |
218 |
if(!objname.CompareTo("NDLevel2"))return &nd_obj; |
if(!objname.CompareTo("TrkLevel2")) { |
219 |
if(!objname.CompareTo("AcLevel2"))return &ac_obj; |
if(!trk2_obj){ |
220 |
if(!objname.CompareTo("OrbitalInfo"))return &orb_obj; |
trk2_obj = new TrkLevel2(); |
221 |
|
trk2_obj->Set(); |
222 |
|
} |
223 |
|
return &trk2_obj; |
224 |
|
}; |
225 |
|
if(!objname.CompareTo("TrkHough")) { |
226 |
|
if(!trkh_obj) trkh_obj = new TrkHough(); |
227 |
|
return &trkh_obj; |
228 |
|
}; |
229 |
|
if(!objname.CompareTo("CaloLevel1")) { |
230 |
|
if(!calo1_obj) calo1_obj = new CaloLevel1(); |
231 |
|
return &calo1_obj; |
232 |
|
}; |
233 |
|
if(!objname.CompareTo("CaloLevel2")) { |
234 |
|
if(!calo2_obj){ |
235 |
|
calo2_obj = new CaloLevel2(); |
236 |
|
calo2_obj->Set(); |
237 |
|
}; |
238 |
|
return &calo2_obj; |
239 |
|
}; |
240 |
|
if(!objname.CompareTo("ToFLevel2")) { |
241 |
|
if(!tof_obj){ |
242 |
|
tof_obj = new ToFLevel2(); |
243 |
|
tof_obj->Set(); |
244 |
|
} |
245 |
|
return &tof_obj; |
246 |
|
}; |
247 |
|
if(!objname.CompareTo("TrigLevel2")) { |
248 |
|
if(!trig_obj) trig_obj = new TrigLevel2(); |
249 |
|
return &trig_obj; |
250 |
|
}; |
251 |
|
if(!objname.CompareTo("S4Level2")) { |
252 |
|
if(!s4_obj) s4_obj = new S4Level2(); |
253 |
|
return &s4_obj; |
254 |
|
}; |
255 |
|
if(!objname.CompareTo("NDLevel2")) { |
256 |
|
if(!nd_obj) nd_obj = new NDLevel2(); |
257 |
|
return &nd_obj; |
258 |
|
}; |
259 |
|
if(!objname.CompareTo("AcLevel2")) { |
260 |
|
if(!ac_obj) ac_obj = new AcLevel2(); |
261 |
|
return &ac_obj; |
262 |
|
}; |
263 |
|
if(!objname.CompareTo("OrbitalInfo")){ |
264 |
|
if(!orb_obj) orb_obj = new OrbitalInfo(); |
265 |
|
return &orb_obj; |
266 |
|
}; |
267 |
|
|
268 |
if(!objname.CompareTo("RunInfo"))return &run_obj; |
if(!objname.CompareTo("RunInfo"))return &run_obj; |
269 |
return NULL; |
return NULL; |
270 |
}; |
}; |
277 |
* (If seqno = -1 retrieves the self-trigger calorimeter track) |
* (If seqno = -1 retrieves the self-trigger calorimeter track) |
278 |
*/ |
*/ |
279 |
CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){ |
CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){ |
280 |
|
|
281 |
if( CaloLevel2::ntrk()==0 ){ |
if( !calo2_obj )return 0; |
282 |
|
|
283 |
|
if( calo2_obj->CaloLevel2::ntrk()==0 ){ |
284 |
cout << "PamLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no Calorimeter tracks are stored"<<endl; |
cout << "PamLevel2::GetCaloStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no Calorimeter tracks are stored"<<endl; |
285 |
return NULL; |
return NULL; |
286 |
}; |
}; |
289 |
Int_t it_calo=0; |
Int_t it_calo=0; |
290 |
|
|
291 |
do{ |
do{ |
292 |
c = CaloLevel2::GetCaloTrkVar(it_calo); |
c = calo2_obj->CaloLevel2::GetCaloTrkVar(it_calo); |
293 |
it_calo++; |
it_calo++; |
294 |
} while( c && seqno != c->trkseqno && it_calo < CaloLevel2::ntrk()); |
} while( c && seqno != c->trkseqno && it_calo < calo2_obj->CaloLevel2::ntrk()); |
295 |
|
|
296 |
if(!c || seqno != c->trkseqno){ |
if(!c || seqno != c->trkseqno){ |
297 |
c = 0; |
c = 0; |
310 |
*/ |
*/ |
311 |
ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){ |
ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){ |
312 |
|
|
313 |
if( ToFLevel2::ntrk()==0 ){ |
if( !tof_obj )return 0; |
314 |
|
|
315 |
|
if( tof_obj->ToFLevel2::ntrk()==0 ){ |
316 |
cout << "PamLevel2::GetToFStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no ToF tracks are stored"<<endl; |
cout << "PamLevel2::GetToFStoredTrack(int) : requested tracker SeqNo "<< seqno <<" but no ToF tracks are stored"<<endl; |
317 |
return NULL; |
return NULL; |
318 |
}; |
}; |
321 |
Int_t it_tof=0; |
Int_t it_tof=0; |
322 |
|
|
323 |
do{ |
do{ |
324 |
c = ToFLevel2::GetToFTrkVar(it_tof); |
c = tof_obj->ToFLevel2::GetToFTrkVar(it_tof); |
325 |
it_tof++; |
it_tof++; |
326 |
} while( c && seqno != c->trkseqno && it_tof < ToFLevel2::ntrk()); |
} while( c && seqno != c->trkseqno && it_tof < tof_obj->ToFLevel2::ntrk()); |
327 |
|
|
328 |
if(!c || seqno != c->trkseqno){ |
if(!c || seqno != c->trkseqno){ |
329 |
c = 0; |
c = 0; |
341 |
* Give the pamela track associated to a tracker track, retrieving related calorimeter and tof track information. |
* Give the pamela track associated to a tracker track, retrieving related calorimeter and tof track information. |
342 |
*/ |
*/ |
343 |
PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){ |
PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){ |
344 |
|
|
345 |
|
cout <<"PamLevel2::GetPamTrackAlong(TrkTrack* t) **obsolete** "<<endl; |
346 |
|
cout <<"(if you use it, remember to delete the PamTrack object)"<<endl; |
347 |
|
|
348 |
CaloTrkVar *c = 0; |
CaloTrkVar *c = 0; |
349 |
ToFTrkVar *o = 0; |
ToFTrkVar *o = 0; |
350 |
|
|
369 |
|
|
370 |
PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){ |
PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){ |
371 |
|
|
372 |
|
cout <<"PamLevel2::GetStoredTrack(Int_t itrk) **to-be-updated** "<<endl; |
373 |
|
cout <<"for the moment, better use separately the methods: TrkLevel2::GetStoredTrack(seqno) CaloLevel2::GetCaloTrkVar(Int_t notrack) ToFLevel2::GetToFTrkVar(Int_t notrack)"<<endl; |
374 |
|
cout <<"(if you use it, remember to delete the PamTrack object)"<<endl; |
375 |
PamTrack *track = 0; |
PamTrack *track = 0; |
376 |
|
|
377 |
if( itrk >=0 && itrk < TrkLevel2::ntrk() ){ |
if( itrk >=0 && itrk < trk2_obj->TrkLevel2::ntrk() ){ |
378 |
|
|
379 |
TrkTrack *t = TrkLevel2::GetStoredTrack(itrk); |
TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(itrk); |
380 |
track = GetPamTrackAlong(t); |
track = GetPamTrackAlong(t); |
381 |
|
|
382 |
|
|
383 |
|
|
384 |
}else{ |
}else{ |
385 |
cout << "PamLevel2::GetStoredTrack(int) : tracker track SeqNo "<< itrk <<" does not exist (GetNTracks() = "<<TrkLevel2::GetNTracks()<<")"<<endl; |
cout << "PamLevel2::GetStoredTrack(int) : tracker track SeqNo "<< itrk <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl; |
386 |
}; |
}; |
387 |
|
|
388 |
return track; |
return track; |
405 |
*/ |
*/ |
406 |
void PamLevel2::SortTracks(TString how){ |
void PamLevel2::SortTracks(TString how){ |
407 |
|
|
408 |
|
// cout <<" PamLevel2::SortTracks(TString how) "<<endl; |
409 |
|
if( !trk2_obj ){ |
410 |
|
cout << "void PamLevel2::SortTracks(TString how): TrkLevel2 not loaded !!!"; |
411 |
|
return; |
412 |
|
}; |
413 |
// cout << "call SortTracs() "<<endl; |
// cout << "call SortTracs() "<<endl; |
414 |
//Save current Object count |
//Save current Object count |
415 |
Int_t ObjectNumber = TProcessID::GetObjectCount(); |
Int_t ObjectNumber = TProcessID::GetObjectCount(); |
416 |
|
|
417 |
// cout << "ObjectNumber "<<ObjectNumber <<endl; |
// cout << "ObjectNumber "<<ObjectNumber <<endl; |
418 |
|
|
419 |
sorted_tracks->Delete(); //temporaneo??? |
// if(!sorted_tracks)sorted_tracks = new TRefArray(); |
420 |
|
// sorted_tracks->Clear(); |
421 |
|
// sorted_tracks.Clear(); |
422 |
|
|
423 |
|
if(!tsorted)tsorted = new TClonesArray("PamTrack",trk2_obj->GetNTracks()); |
424 |
|
tsorted->Delete(); |
425 |
|
TClonesArray &ttsorted = *tsorted; |
426 |
|
if(!timage)timage = new TClonesArray("PamTrack",trk2_obj->GetNTracks()); |
427 |
|
timage->Delete(); |
428 |
|
TClonesArray &ttimage = *timage; |
429 |
|
|
430 |
|
|
431 |
// loop over the tracks sorted by the tracker |
// loop over the tracks sorted by the tracker |
432 |
Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase); |
Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase); |
433 |
Bool_t use_CAL = how.Contains("CAL", TString::kIgnoreCase); |
Bool_t use_CAL = how.Contains("CAL", TString::kIgnoreCase); |
443 |
|
|
444 |
// cout << "use_CAL "<<use_CAL<<" use_TOF "<<use_TOF<<" use_TRK "<<use_TRK <<endl; |
// cout << "use_CAL "<<use_CAL<<" use_TOF "<<use_TOF<<" use_TRK "<<use_TRK <<endl; |
445 |
|
|
446 |
for(Int_t i=0; i < TrkLevel2::GetNTracks(); i++){ |
for(Int_t i=0; i < trk2_obj->TrkLevel2::GetNTracks(); i++){ |
447 |
|
|
448 |
TrkTrack *ts = 0; |
TrkTrack *ts = 0; |
449 |
|
CaloTrkVar *cs = 0; |
450 |
|
ToFTrkVar *os = 0; |
451 |
|
|
452 |
// get tracker tracks |
// get tracker tracks |
453 |
TrkTrack *tp = TrkLevel2::GetTrack(i); //tracker |
TrkTrack *tp = trk2_obj->TrkLevel2::GetTrack(i); //tracker |
454 |
CaloTrkVar *cp = 0; |
CaloTrkVar *cp = GetCaloStoredTrack(tp->GetSeqNo()); |
455 |
ToFTrkVar *op = 0; |
ToFTrkVar *op = GetToFStoredTrack(tp->GetSeqNo()); |
456 |
|
|
457 |
|
TrkTrack *ti = 0; //tracker (image) |
458 |
|
CaloTrkVar *ci = 0; |
459 |
|
ToFTrkVar *oi = 0; |
460 |
// cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl; |
// cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl; |
461 |
// if track has an image, check image selection |
// if track has an image, check image selection |
462 |
if(tp->HasImage()){ |
if(tp->HasImage()){ |
463 |
|
|
464 |
TrkTrack *ti = TrkLevel2::GetTrackImage(i); //tracker (image) |
ti = trk2_obj->TrkLevel2::GetTrackImage(i); //tracker (image) |
465 |
CaloTrkVar *ci = 0; |
ci = GetCaloStoredTrack(ti->GetSeqNo()); |
466 |
ToFTrkVar *oi = 0; |
oi = GetToFStoredTrack(ti->GetSeqNo()); |
467 |
|
|
468 |
// cout << "its image "<<i << " "<<hex<< ti <<dec<< endl; |
// cout << "its image "<<i << " "<<hex<< ti <<dec<< endl; |
469 |
|
|
476 |
// ------------------------ |
// ------------------------ |
477 |
// check the Y spatial residual on the first calorimeter plane |
// check the Y spatial residual on the first calorimeter plane |
478 |
// (cut on calorimeter variables from Emiliano) |
// (cut on calorimeter variables from Emiliano) |
479 |
|
if( use_CAL && !calo2_obj ){ |
480 |
|
cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but CaloLevel2 not loaded !!!"; |
481 |
|
return; |
482 |
|
}; |
483 |
if( |
if( |
484 |
use_CAL && |
use_CAL && |
485 |
npcfit[1] > 15 && //no. of fit planes on Y view |
calo2_obj->npcfit[1] > 15 && //no. of fit planes on Y view |
486 |
varcfit[1] < 1000. && //fit variance on Y view |
calo2_obj->varcfit[1] < 1000. && //fit variance on Y view |
487 |
true){ |
true){ |
488 |
|
|
|
cp = GetCaloStoredTrack(tp->GetSeqNo()); |
|
|
ci = GetCaloStoredTrack(ti->GetSeqNo()); |
|
489 |
|
|
490 |
Float_t resy_p = cp->tbar[0][1] - cbar[0][1]; if(resy_p < 0)resy_p= - resy_p; |
Float_t resy_p = cp->tbar[0][1] - calo2_obj->cbar[0][1]; if(resy_p < 0)resy_p= - resy_p; |
491 |
Float_t resy_i = ci->tbar[0][1] - cbar[0][1]; if(resy_i < 0)resy_i= - resy_i; |
Float_t resy_i = ci->tbar[0][1] - calo2_obj->cbar[0][1]; if(resy_i < 0)resy_i= - resy_i; |
492 |
|
|
493 |
if(resy_p <= resy_i) tp_score++; |
if(resy_p <= resy_i) tp_score++; |
494 |
else ti_score++; |
else ti_score++; |
501 |
// ------------------------ |
// ------------------------ |
502 |
// check the number of hit pmts along the track |
// check the number of hit pmts along the track |
503 |
// on S12 S21 and S32, where paddles are parallel to Y axis |
// on S12 S21 and S32, where paddles are parallel to Y axis |
504 |
|
if( use_TOF && !tof_obj ){ |
505 |
|
cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but ToFLevel2 not loaded !!!"; |
506 |
|
return; |
507 |
|
}; |
508 |
if( use_TOF ){ |
if( use_TOF ){ |
509 |
|
|
510 |
Int_t nphit_p =0; |
Int_t nphit_p =0; |
511 |
Int_t nphit_i =0; |
Int_t nphit_i =0; |
512 |
|
|
|
op = GetToFStoredTrack(tp->GetSeqNo()); |
|
|
oi = GetToFStoredTrack(ti->GetSeqNo()); |
|
513 |
|
|
514 |
/* cout << "track: npmtadc "<< op->npmtadc << endl; |
/* cout << "track: npmtadc "<< op->npmtadc << endl; |
515 |
cout << "track: npmttdc "<< op->npmttdc << endl; |
cout << "track: npmttdc "<< op->npmttdc << endl; |
517 |
cout << "image: npmttdc "<< oi->npmttdc << endl;*/ |
cout << "image: npmttdc "<< oi->npmttdc << endl;*/ |
518 |
|
|
519 |
for (Int_t ih=0; ih < op->npmtadc; ih++){ |
for (Int_t ih=0; ih < op->npmtadc; ih++){ |
520 |
Int_t pl = GetPlaneIndex( (op->pmtadc).At(ih) ); |
Int_t pl = tof_obj->GetPlaneIndex( (op->pmtadc).At(ih) ); |
521 |
if(pl == 1 || pl == 2 || pl == 5)nphit_p++; |
if(pl == 1 || pl == 2 || pl == 5)nphit_p++; |
522 |
}; |
}; |
523 |
|
|
524 |
for (Int_t ih=0; ih < oi->npmtadc; ih++){ |
for (Int_t ih=0; ih < oi->npmtadc; ih++){ |
525 |
Int_t pl = GetPlaneIndex( (oi->pmtadc).At(ih) ); |
Int_t pl = tof_obj->GetPlaneIndex( (oi->pmtadc).At(ih) ); |
526 |
if(pl == 1 || pl == 2 || pl == 5)nphit_i++; |
if(pl == 1 || pl == 2 || pl == 5)nphit_i++; |
527 |
}; |
}; |
528 |
|
|
552 |
// ------------------------ |
// ------------------------ |
553 |
// the winner is.... |
// the winner is.... |
554 |
// ------------------------ |
// ------------------------ |
555 |
if (tp_score > ti_score) ts = tp;//the track sorted by the tracker!! |
if (tp_score > ti_score) { |
556 |
else if (tp_score < ti_score) ts = ti;//its image!! |
// ts = tp;//the track sorted by the tracker!! |
557 |
else { |
// cs = cp; |
558 |
ts = tp; |
// os = op; |
559 |
|
}else if (tp_score < ti_score) { |
560 |
|
ts = ti;//its image!! |
561 |
|
cs = ci; |
562 |
|
os = oi; |
563 |
|
|
564 |
|
ti = tp;//its image!! |
565 |
|
ci = cp; |
566 |
|
oi = op; |
567 |
|
|
568 |
|
tp = ts;//its image!! |
569 |
|
cp = cs; |
570 |
|
op = os; |
571 |
|
|
572 |
|
|
573 |
|
}else { |
574 |
|
// ts = tp; |
575 |
|
// cs = cp; |
576 |
|
// os = op; |
577 |
// cout << "Warning - track image ambiguity not solved" << endl; |
// cout << "Warning - track image ambiguity not solved" << endl; |
578 |
// cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl; |
// cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl; |
579 |
}; |
}; |
580 |
|
|
581 |
}else{ |
}else{ |
582 |
ts = tp; |
// ts = tp; |
583 |
|
// cs = cp; |
584 |
|
// os = op; |
585 |
}; |
}; |
586 |
|
|
587 |
// cout <<" SortTracks() "<<i<<" -- "<<ts<<endl; |
// cout <<" SortTracks() "<<i<<" -- "<<ts<<endl; |
588 |
sorted_tracks->Add(ts);//save the track in the sorted array |
// sorted_tracks->Add(ts);//save the track in the sorted array |
589 |
|
// sorted_tracks.Add(ts);//save the track in the sorted array |
590 |
|
// sorted_tracks.Add(tp);//save the track in the sorted array |
591 |
// cout << "SortTracks:: sorted_tracks->Add(it) "<<i<<" "<<ts<<endl; |
// cout << "SortTracks:: sorted_tracks->Add(it) "<<i<<" "<<ts<<endl; |
592 |
|
// cout<<"o "<<tp<<endl; |
593 |
|
// cout<<"o "<<cp<<endl; |
594 |
|
// cout<<"o "<<op<<endl; |
595 |
|
new(ttsorted[i]) PamTrack(tp,cp,op); |
596 |
|
new(ttimage[i]) PamTrack(ti,ci,oi); |
597 |
}; |
}; |
598 |
|
|
599 |
|
if( tsorted->GetEntries() != trk2_obj->GetNTracks() ){ |
600 |
|
cout << "void PamLevel2::SortTracks(TString how): tsorted->GetEntries() "<<tsorted->GetEntries()<<" != trk2_obj->GetNTracks() = "<<trk2_obj->GetNTracks() <<endl; |
601 |
|
tsorted->Delete(); tsorted=0; |
602 |
|
timage->Delete(); timage=0; |
603 |
|
} |
604 |
|
|
605 |
//Restore Object count |
//Restore Object count |
606 |
//To save space in the table keeping track of all referenced objects |
//To save space in the table keeping track of all referenced objects |
607 |
//we assume that our events do not address each other. We reset the |
//We reset the object count to what it was at the beginning of the event. |
|
//object count to what it was at the beginning of the event. |
|
608 |
TProcessID::SetObjectCount(ObjectNumber); |
TProcessID::SetObjectCount(ObjectNumber); |
609 |
|
|
610 |
}; |
}; |
616 |
* This method overrides TrkLevel2::GetTracks(), where sorting is done by decreasing number of fit points and increasing chi^2. |
* This method overrides TrkLevel2::GetTracks(), where sorting is done by decreasing number of fit points and increasing chi^2. |
617 |
* PamLevel2::GetTracks() keeps the same track order given by TrkLevel2::GetTracks(), but checks image selection by using calorimeter and ToF tracking information. |
* PamLevel2::GetTracks() keeps the same track order given by TrkLevel2::GetTracks(), but checks image selection by using calorimeter and ToF tracking information. |
618 |
*/ |
*/ |
619 |
TRefArray *PamLevel2::GetTracks(){ |
// TRefArray *PamLevel2::GetTracks(){ |
620 |
|
|
621 |
|
// // *-*-*-*-*-*-*-*-*-*-*-*-* |
622 |
|
// SortTracks("+CAL+TOF"); |
623 |
|
// // *-*-*-*-*-*-*-*-*-*-*-*-* |
624 |
|
|
625 |
|
// // return sorted_tracks; |
626 |
|
// return &sorted_tracks; |
627 |
|
|
628 |
|
// }; |
629 |
|
TClonesArray *PamLevel2::GetTracks(){ |
630 |
|
|
631 |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
632 |
SortTracks("+CAL+TOF"); |
SortTracks("+CAL+TOF"); |
633 |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
634 |
|
|
635 |
return sorted_tracks; |
return tsorted; |
636 |
|
|
637 |
}; |
}; |
638 |
//-------------------------------------- |
//-------------------------------------- |
639 |
// |
// |
645 |
* @param it Track number, ranging from 0 to GetNTracks(). |
* @param it Track number, ranging from 0 to GetNTracks(). |
646 |
*/ |
*/ |
647 |
PamTrack *PamLevel2::GetTrack(int it){ |
PamTrack *PamLevel2::GetTrack(int it){ |
648 |
|
|
649 |
|
// if(!trk2_obj) return 0; |
650 |
|
|
651 |
|
// // *-*-*-*-*-*-*-*-*-*-*-*-* |
652 |
|
// SortTracks("+CAL+TOF"); |
653 |
|
// // *-*-*-*-*-*-*-*-*-*-*-*-* |
654 |
|
// // if(!sorted_tracks)return 0; |
655 |
|
// if(sorted_tracks.IsEmpty())return 0; |
656 |
|
|
657 |
|
// PamTrack *track = 0; |
658 |
|
|
659 |
|
// if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() && it<sorted_tracks.GetEntriesFast() ){ |
660 |
|
// TrkTrack *t = (TrkTrack*)sorted_tracks.At(it); |
661 |
|
// track = GetPamTrackAlong(t); |
662 |
|
// }else{ |
663 |
|
// cout << "PamLevel2::GetTrack(int) : tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl; |
664 |
|
// }; |
665 |
|
|
666 |
|
// return track; |
667 |
|
|
668 |
|
// cout << "PamLevel2::GetTrack(int it) "<<endl; |
669 |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
670 |
SortTracks("+CAL+TOF"); |
SortTracks("+CAL+TOF"); |
671 |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
672 |
|
if(!tsorted)return 0; |
673 |
|
if(!tsorted->GetEntries())return 0; |
674 |
|
|
675 |
PamTrack *track = 0; |
PamTrack *track = 0; |
676 |
|
|
677 |
if( it >=0 && it < TrkLevel2::GetNTracks() && it<sorted_tracks->GetEntriesFast() ){ |
if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks()){ |
678 |
TrkTrack *t = (TrkTrack*)sorted_tracks->At(it); |
// TrkTrack *t = (TrkTrack*)sorted_tracks.At(it); |
679 |
track = GetPamTrackAlong(t); |
// track = GetPamTrackAlong(t); |
680 |
|
TClonesArray &t = *(tsorted); |
681 |
|
track = (PamTrack*)t[it]; |
682 |
|
|
683 |
}else{ |
}else{ |
684 |
cout << "PamLevel2::GetTrack(int) : tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<TrkLevel2::GetNTracks()<<")"<<endl; |
cout << "PamLevel2::GetTrack(int) : tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl; |
685 |
}; |
}; |
686 |
|
|
687 |
return track; |
return track; |
698 |
*/ |
*/ |
699 |
PamTrack *PamLevel2::GetTrackImage(int it){ |
PamTrack *PamLevel2::GetTrackImage(int it){ |
700 |
|
|
701 |
|
// if(!trk2_obj) return 0; |
702 |
|
|
703 |
|
// // *-*-*-*-*-*-*-*-*-*-*-*-* |
704 |
|
// SortTracks("+CAL+TOF"); |
705 |
|
// // *-*-*-*-*-*-*-*-*-*-*-*-* |
706 |
|
// // if(!sorted_tracks)return 0; |
707 |
|
// if(sorted_tracks.IsEmpty())return 0; |
708 |
|
|
709 |
|
// PamTrack *image = 0; |
710 |
|
|
711 |
|
// if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){ |
712 |
|
// TrkTrack *temp = (TrkTrack*)sorted_tracks.At(it); |
713 |
|
// if( temp->HasImage() ){ |
714 |
|
// TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(temp->GetImageSeqNo()); |
715 |
|
// image = GetPamTrackAlong(t); |
716 |
|
// }else{ |
717 |
|
// cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl; |
718 |
|
// }; |
719 |
|
// }else{ |
720 |
|
// cout << "PamLevel2::GetTrackImage(int) : Tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl; |
721 |
|
// }; |
722 |
|
|
723 |
|
// return image; |
724 |
|
|
725 |
|
|
726 |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
727 |
SortTracks("+CAL+TOF"); |
SortTracks("+CAL+TOF"); |
728 |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
// *-*-*-*-*-*-*-*-*-*-*-*-* |
729 |
|
if(!timage)return 0; |
730 |
|
if(!timage->GetEntries())return 0; |
731 |
|
|
732 |
PamTrack *image = 0; |
PamTrack *image = 0; |
733 |
|
|
734 |
if( it >=0 && it < TrkLevel2::GetNTracks() ){ |
if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){ |
735 |
TrkTrack *temp = (TrkTrack*)sorted_tracks->At(it); |
TClonesArray &t = *(tsorted); |
736 |
if( temp->HasImage() ){ |
PamTrack *temp = (PamTrack*)t[it]; |
737 |
TrkTrack *t = TrkLevel2::GetStoredTrack(temp->GetImageSeqNo()); |
if( temp->GetTrkTrack()->HasImage() ){ |
738 |
image = GetPamTrackAlong(t); |
TClonesArray &t = *(timage); |
739 |
|
image = (PamTrack*)t[it]; |
740 |
}else{ |
}else{ |
741 |
cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl; |
// cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl; |
742 |
}; |
}; |
743 |
}else{ |
}else{ |
744 |
cout << "PamLevel2::GetTrackImage(int) : Tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<TrkLevel2::GetNTracks()<<")"<<endl; |
cout << "PamLevel2::GetTrackImage(int) : Tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl; |
745 |
}; |
}; |
746 |
|
|
747 |
return image; |
return image; |
757 |
* @param detlist String to select trees to be included |
* @param detlist String to select trees to be included |
758 |
* @return Pointer to a TTree |
* @return Pointer to a TTree |
759 |
*/ |
*/ |
760 |
TTree *PamLevel2::GetPamTree(TFile *f, TString detlist="+ALL"){ |
TTree *PamLevel2::GetPamTree(TFile *f, TString detlist ){ |
|
|
|
|
// cout << "WARNING!!! -- obsolete method -- \n"; |
|
|
// cout << "better use TChain *PamLevel2::GetPamTree(TList*, TString) \n"; |
|
761 |
|
|
762 |
|
// if( !detlist.IsNull() ) SetWhichTrees(detlist); |
763 |
|
// else GetWhichTrees(f); |
764 |
SetWhichTrees(detlist); |
SetWhichTrees(detlist); |
|
TTree *Trout =0; |
|
765 |
|
|
766 |
|
TTree *Trout =0; |
767 |
|
|
768 |
TString fname = f->GetName(); |
TString fname = f->GetName(); |
769 |
if(!CheckLevel2File(fname))return NULL; |
if(!CheckLevel2File(fname))return NULL; |
770 |
|
|
771 |
|
UInt_t *found=0; |
772 |
|
|
773 |
|
cout<< "GetPamTree(TFile*,TString): detector list --> "; |
774 |
|
if(TRK1)cout<<"TRK1 "; |
775 |
|
if(TRK2)cout<<"TRK2 "; |
776 |
|
if(TRKh)cout<<"TRKH "; |
777 |
|
if(CAL1)cout<<"CAL1 "; |
778 |
|
if(CAL2)cout<<"CAL2 "; |
779 |
|
if(TOF)cout<<"TOF "; |
780 |
|
if(TRG)cout<<"TRG "; |
781 |
|
if(AC)cout<<"AC "; |
782 |
|
if(ND)cout<<"ND "; |
783 |
|
if(S4)cout<<"S4 "; |
784 |
|
if(ORB)cout<<"ORB "; |
785 |
|
cout << endl; |
786 |
|
|
787 |
|
f->cd(); |
788 |
|
|
789 |
// Tracker |
// Tracker |
790 |
TTree *T = (TTree*)f->Get("Tracker"); |
TTree *T = (TTree*)f->Get("Tracker"); |
791 |
if(T && (TRK2||TRK1||TRKh)) { |
if(T && (TRK2||TRK1||TRKh)) { |
792 |
if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2")); |
if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2")); |
793 |
|
else T->SetBranchStatus("TrkLevel2",0,found); |
794 |
if(TRK2)cout << "Tracker : set branch address TrkLevel2"<<endl; |
if(TRK2)cout << "Tracker : set branch address TrkLevel2"<<endl; |
795 |
if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1")); |
if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1")); |
796 |
|
else T->SetBranchStatus("TrkLevel1",0,found); |
797 |
if(TRK1)cout << "Tracker : set branch address TrkLevel1"<<endl; |
if(TRK1)cout << "Tracker : set branch address TrkLevel1"<<endl; |
798 |
if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough")); |
if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough")); |
799 |
|
else T->SetBranchStatus("TrkHough",0,found); |
800 |
if(TRKh)cout << "Tracker : set branch address TrkHough"<<endl; |
if(TRKh)cout << "Tracker : set branch address TrkHough"<<endl; |
801 |
if(!Trout)Trout=T; |
if(!Trout)Trout=T; |
802 |
else Trout->AddFriend(T); |
else Trout->AddFriend(T); |
807 |
TTree *C = (TTree*)f->Get("Calorimeter"); |
TTree *C = (TTree*)f->Get("Calorimeter"); |
808 |
if(C && (CAL2||CAL1)) { |
if(C && (CAL2||CAL1)) { |
809 |
if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2")); |
if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2")); |
810 |
|
else C->SetBranchStatus("CaloLevel2",0,found); |
811 |
if(CAL2)cout << "Calorimeter : set branch address CaloLevel2"<<endl; |
if(CAL2)cout << "Calorimeter : set branch address CaloLevel2"<<endl; |
812 |
if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1")); |
if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1")); |
813 |
|
else C->SetBranchStatus("CaloLevel1",0,found); |
814 |
if(CAL1)cout << "Calorimeter : set branch address CaloLevel1"<<endl; |
if(CAL1)cout << "Calorimeter : set branch address CaloLevel1"<<endl; |
815 |
if(!Trout)Trout=C; |
if(!Trout)Trout=C; |
816 |
else Trout->AddFriend(C); |
else Trout->AddFriend(C); |
817 |
}else{ |
}else{ |
818 |
cout << "Calorimeter : missing tree"<<endl; |
cout << "Calorimeter : missing tree"<<endl; |
819 |
}; |
}; |
820 |
|
|
821 |
// ToF |
// ToF |
822 |
TTree *O = (TTree*)f->Get("ToF"); |
TTree *O = (TTree*)f->Get("ToF"); |
823 |
if(O && TOF) { |
if(O && TOF) { |
869 |
cout << "Anticounter : missing tree"<<endl; |
cout << "Anticounter : missing tree"<<endl; |
870 |
}; |
}; |
871 |
// Orbital Info |
// Orbital Info |
872 |
TTree *B = (TTree*)f->Get("OrbitalInfo"); |
TTree *B = (TTree*)f->Get("OrbitalInfo"); |
873 |
if(B && ORB) { |
if(B && ORB) { |
874 |
B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo")); |
B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo")); |
875 |
cout << "OrbitalInfo : set branch address OrbitalInfo"<<endl; |
cout << "OrbitalInfo : set branch address OrbitalInfo"<<endl; |
876 |
if(!Trout)Trout=O; |
if(!Trout)Trout=O; |
877 |
else Trout->AddFriend(B); |
else Trout->AddFriend(B); |
878 |
}else{ |
}else{ |
879 |
cout << "OrbitalInfo : missing tree"<<endl; |
cout << "OrbitalInfo : missing tree"<<endl; |
880 |
}; |
}; |
881 |
|
|
882 |
cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl; |
cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl; |
883 |
|
|
884 |
return Trout; |
return Trout; |
885 |
|
|
886 |
} |
} |
887 |
//-------------------------------------- |
//-------------------------------------- |
911 |
// if no input file list is given: |
// if no input file list is given: |
912 |
if ( flisttxt != "" ){ |
if ( flisttxt != "" ){ |
913 |
|
|
914 |
if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) )return 0; |
if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){ |
915 |
|
cout <<"File "<<flisttxt<<" not found"<<endl; |
916 |
|
return 0; |
917 |
|
} |
918 |
// flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt)); |
// flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt)); |
919 |
flisttxt = fullpath; |
flisttxt = fullpath; |
920 |
|
|
921 |
if( !gSystem->ChangeDirectory(ddir) )return 0; |
if( !gSystem->ChangeDirectory(ddir) ){ |
922 |
|
cout << "Cannot change directory : "<<ddir<<endl; |
923 |
|
return 0; |
924 |
|
} |
925 |
|
|
926 |
cout <<"Input file list : " << flisttxt <<endl; |
cout <<"Input file list : " << flisttxt <<endl; |
927 |
ifstream in; |
ifstream in; |
933 |
if (!in.good()) break; |
if (!in.good()) break; |
934 |
line++; |
line++; |
935 |
// cout <<"(1) " << file << endl; |
// cout <<"(1) " << file << endl; |
|
if(file.Contains("#"))file = file(0,file.First("#")); |
|
|
// cout <<"(2) " << file << endl; |
|
936 |
if(file.IsNull()){ |
if(file.IsNull()){ |
937 |
cout << "-- list interrupted at line "<<line <<endl; |
cout << "-- list interrupted at line "<<line <<endl; |
938 |
break; |
break; |
939 |
} |
} |
940 |
|
if(file.Contains("#"))file = file(0,file.First("#")); |
941 |
|
// cout <<"(2) " << file << endl; |
942 |
if( gSystem->IsFileInIncludePath(file,&fullpath) ){ |
if( gSystem->IsFileInIncludePath(file,&fullpath) ){ |
943 |
|
|
944 |
contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list |
contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list |
945 |
|
|
946 |
}else{ |
}else{ |
947 |
if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl; |
// if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl; |
948 |
}; |
}; |
949 |
}; |
}; |
950 |
in.close(); |
in.close(); |
994 |
* @param detlist String to select trees to be included |
* @param detlist String to select trees to be included |
995 |
* @return Pointer to a TChain |
* @return Pointer to a TChain |
996 |
*/ |
*/ |
997 |
TChain *PamLevel2::GetPamTree(TList *fl, TString detlist="+ALL"){ |
TChain *PamLevel2::GetPamTree(TList *fl, TString detlist ){ |
|
|
|
|
// TChain *Tout=0; |
|
|
// if(Tout)Tout->Delete(); |
|
|
// Tout = NULL; |
|
998 |
|
|
999 |
TChain *Tout =0; |
TChain *Trout =0; |
1000 |
|
|
1001 |
|
// if( !detlist.IsNull() )SetWhichTrees(detlist); |
1002 |
SetWhichTrees(detlist); |
SetWhichTrees(detlist); |
1003 |
|
|
1004 |
|
cout<< "GetPamTree(TList*,TString): input detector list --> "; |
1005 |
|
if(TRK1)cout<<"TRK1 "; |
1006 |
|
if(TRK2)cout<<"TRK2 "; |
1007 |
|
if(TRKh)cout<<"TRKH "; |
1008 |
|
if(CAL1)cout<<"CAL1 "; |
1009 |
|
if(CAL2)cout<<"CAL2 "; |
1010 |
|
if(TOF)cout<<"TOF "; |
1011 |
|
if(TRG)cout<<"TRG "; |
1012 |
|
if(AC)cout<<"AC "; |
1013 |
|
if(ND)cout<<"ND "; |
1014 |
|
if(S4)cout<<"S4 "; |
1015 |
|
if(ORB)cout<<"ORB "; |
1016 |
|
cout << endl; |
1017 |
|
|
1018 |
TChain *T = 0; |
TChain *T = 0; |
1019 |
TChain *C = 0; |
TChain *C = 0; |
1020 |
TChain *O = 0; |
TChain *O = 0; |
1025 |
TChain *B = 0; |
TChain *B = 0; |
1026 |
|
|
1027 |
if(TRK2||TRK1||TRKh) T = new TChain("Tracker"); |
if(TRK2||TRK1||TRKh) T = new TChain("Tracker"); |
1028 |
if(CAL2||CAL1) C = new TChain("Calorimeter"); |
if(CAL2||CAL1) C = new TChain("Calorimeter"); |
1029 |
if(TOF) O = new TChain("ToF"); |
if(TOF) O = new TChain("ToF"); |
1030 |
if(TRG) R = new TChain("Trigger"); |
if(TRG) R = new TChain("Trigger"); |
1031 |
if(S4) S = new TChain("S4"); |
if(S4) S = new TChain("S4"); |
1032 |
if(ND) N = new TChain("NeutronD"); |
if(ND) N = new TChain("NeutronD"); |
1033 |
if(AC) A = new TChain("Anticounter"); |
if(AC) A = new TChain("Anticounter"); |
1034 |
if(ORB) B = new TChain("OrbitalInfo"); |
if(ORB) B = new TChain("OrbitalInfo"); |
1035 |
|
|
1036 |
// loop over files and create chains |
// loop over files and create chains |
1037 |
TIter next(fl); |
TIter next(fl); |
1038 |
TSystemFile *questo = 0; |
TSystemFile *questo = 0; |
1039 |
while ( (questo = (TSystemFile*) next()) ) { |
while ( (questo = (TSystemFile*) next()) ) { |
1040 |
TString name = questo->GetName(); |
TString name = questo->GetName(); |
1041 |
// cout << "File: "<< name << endl; |
cout << "File: "<< name << endl; |
1042 |
if( CheckLevel2File(name) ){ |
if( CheckLevel2File(name) ){ |
1043 |
if(TRK2||TRK1||TRKh) T->Add(name); |
if(TRK2||TRK1||TRKh) T->Add(name); |
1044 |
if(CAL1||CAL2) C->Add(name); |
if(CAL1||CAL2) C->Add(name); |
1045 |
if(TOF) O->Add(name); |
if(TOF) O->Add(name); |
1046 |
if(TRG) R->Add(name); |
if(TRG) R->Add(name); |
1047 |
if(S4) S->Add(name); |
if(S4) S->Add(name); |
1048 |
if(ND) N->Add(name); |
if(ND) N->Add(name); |
1049 |
if(AC) A->Add(name); |
if(AC) A->Add(name); |
1050 |
if(ORB) B->Add(name); |
if(ORB) B->Add(name); |
1051 |
}; |
}; |
1052 |
} |
} |
1053 |
|
|
1054 |
cout << "done chain \n"; |
cout << "done chain \n"; |
1055 |
|
|
1056 |
|
UInt_t *found=0; |
1057 |
|
// Tracker |
1058 |
|
if(T && (TRK2||TRK1||TRKh)) { |
1059 |
|
if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2")); |
1060 |
|
else T->SetBranchStatus("TrkLevel2",0,found); |
1061 |
|
if(TRK2)cout << "Tracker : set branch address TrkLevel2"<<endl; |
1062 |
|
if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1")); |
1063 |
|
else T->SetBranchStatus("TrkLevel1",0,found); |
1064 |
|
if(TRK1)cout << "Tracker : set branch address TrkLevel1"<<endl; |
1065 |
|
if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough")); |
1066 |
|
else T->SetBranchStatus("TrkHough",0,found); |
1067 |
|
if(TRKh)cout << "Tracker : set branch address TrkHough"<<endl; |
1068 |
|
if(!Trout)Trout=T; |
1069 |
|
else Trout->AddFriend("Tracker"); |
1070 |
|
}else{ |
1071 |
|
cout << "Tracker : missing tree"<<endl; |
1072 |
|
}; |
1073 |
// Calorimeter |
// Calorimeter |
1074 |
if(CAL2||CAL1) { |
if(C && (CAL2||CAL1)) { |
1075 |
if(!Tout)Tout=C; |
if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2")); |
1076 |
else Tout->AddFriend("Calorimeter"); |
else C->SetBranchStatus("CaloLevel2",0,found); |
1077 |
|
if(CAL2)cout << "Calorimeter : set branch address CaloLevel2"<<endl; |
1078 |
|
if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1")); |
1079 |
|
else C->SetBranchStatus("CaloLevel1",0,found); |
1080 |
|
if(CAL1)cout << "Calorimeter : set branch address CaloLevel1"<<endl; |
1081 |
|
if(!Trout)Trout=C; |
1082 |
|
else Trout->AddFriend("Calorimeter"); |
1083 |
|
}else{ |
1084 |
|
cout << "Calorimeter : missing tree"<<endl; |
1085 |
}; |
}; |
1086 |
// ToF |
// ToF |
1087 |
if(TOF) { |
if(O && TOF) { |
1088 |
if(!Tout)Tout=O; |
O->SetBranchAddress("ToFLevel2", GetPointerTo("ToFLevel2")); |
1089 |
else Tout->AddFriend("ToF"); |
cout << "ToF : set branch address ToFLevel2"<<endl; |
1090 |
|
if(!Trout)Trout=O; |
1091 |
|
else Trout->AddFriend("ToF"); |
1092 |
|
}else{ |
1093 |
|
cout << "ToF : missing tree"<<endl; |
1094 |
}; |
}; |
|
// Tracker |
|
|
if(TRK2||TRK1||TRKh){ |
|
|
if(!Tout)Tout=T; |
|
|
else Tout->AddFriend("Tracker"); |
|
|
} |
|
1095 |
// Trigger |
// Trigger |
1096 |
if(TRG) { |
if(R && TRG) { |
1097 |
if(!Tout)Tout=R; |
R->SetBranchAddress("TrigLevel2", GetPointerTo("TrigLevel2")); |
1098 |
else Tout->AddFriend("Trigger"); |
cout << "Trigger : set branch address TrigLevel2"<<endl; |
1099 |
|
if(!Trout)Trout=O; |
1100 |
|
else Trout->AddFriend("Trigger"); |
1101 |
|
}else{ |
1102 |
|
cout << "Trigger : missing tree"<<endl; |
1103 |
}; |
}; |
1104 |
// S4 |
// S4 |
1105 |
if(S4) { |
if(S && S4) { |
1106 |
if(!Tout)Tout=S; |
S->SetBranchAddress("S4Level2", GetPointerTo("S4Level2")); |
1107 |
else Tout->AddFriend("S4"); |
cout << "S4 : set branch address S4Level2"<<endl; |
1108 |
|
if(!Trout)Trout=O; |
1109 |
|
else Trout->AddFriend("S4"); |
1110 |
|
}else{ |
1111 |
|
cout << "S4 : missing tree"<<endl; |
1112 |
}; |
}; |
1113 |
// Neutron Detector |
// Neutron Detector |
1114 |
if(ND) { |
if(N && ND) { |
1115 |
if(!Tout)Tout=N; |
N->SetBranchAddress("NDLevel2", GetPointerTo("NDLevel2")); |
1116 |
else Tout->AddFriend("NeutronD"); |
cout << "NeutronD : set branch address NDLevel2"<<endl; |
1117 |
|
if(!Trout)Trout=O; |
1118 |
|
else Trout->AddFriend("NeutronD"); |
1119 |
|
}else{ |
1120 |
|
cout << "NeutronD : missing tree"<<endl; |
1121 |
}; |
}; |
1122 |
// Anticounters |
// Anticounters |
1123 |
if(AC) { |
if(A && AC) { |
1124 |
if(!Tout)Tout=A; |
A->SetBranchAddress("AcLevel2", GetPointerTo("AcLevel2")); |
1125 |
else Tout->AddFriend("Anticounter"); |
cout << "Anticounter : set branch address AcLevel2"<<endl; |
1126 |
|
if(!Trout)Trout=O; |
1127 |
|
else Trout->AddFriend("Anticounter"); |
1128 |
|
}else{ |
1129 |
|
cout << "Anticounter : missing tree"<<endl; |
1130 |
}; |
}; |
1131 |
// OrbitalInfo |
// Orbital Info |
1132 |
if(ORB) { |
if(B && ORB) { |
1133 |
if(!Tout)Tout=B; |
B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo")); |
1134 |
else Tout->AddFriend("OrbitalInfo"); |
cout << "OrbitalInfo : set branch address OrbitalInfo"<<endl; |
1135 |
|
if(!Trout)Trout=O; |
1136 |
|
else Trout->AddFriend("OrbitalInfo"); |
1137 |
|
}else{ |
1138 |
|
cout << "OrbitalInfo : missing tree"<<endl; |
1139 |
}; |
}; |
1140 |
|
|
1141 |
|
// cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl; |
1142 |
|
|
1143 |
|
|
1144 |
|
// // Tracker |
1145 |
|
// if(TRK2||TRK1||TRKh){ |
1146 |
|
// if(!Tout)Tout=T; |
1147 |
|
// else Tout->AddFriend("Tracker"); |
1148 |
|
// } |
1149 |
|
// // Calorimeter |
1150 |
|
// if(CAL2||CAL1) { |
1151 |
|
// if(!Tout)Tout=C; |
1152 |
|
// else Tout->AddFriend("Calorimeter"); |
1153 |
|
// }; |
1154 |
|
// // ToF |
1155 |
|
// if(TOF) { |
1156 |
|
// if(!Tout)Tout=O; |
1157 |
|
// else Tout->AddFriend("ToF"); |
1158 |
|
// }; |
1159 |
|
// // Trigger |
1160 |
|
// if(TRG) { |
1161 |
|
// if(!Tout)Tout=R; |
1162 |
|
// else Tout->AddFriend("Trigger"); |
1163 |
|
// }; |
1164 |
|
// // S4 |
1165 |
|
// if(S4) { |
1166 |
|
// if(!Tout)Tout=S; |
1167 |
|
// else Tout->AddFriend("S4"); |
1168 |
|
// }; |
1169 |
|
// // Neutron Detector |
1170 |
|
// if(ND) { |
1171 |
|
// if(!Tout)Tout=N; |
1172 |
|
// else Tout->AddFriend("NeutronD"); |
1173 |
|
// }; |
1174 |
|
// // Anticounters |
1175 |
|
// if(AC) { |
1176 |
|
// if(!Tout)Tout=A; |
1177 |
|
// else Tout->AddFriend("Anticounter"); |
1178 |
|
// }; |
1179 |
|
// // OrbitalInfo |
1180 |
|
// if(ORB) { |
1181 |
|
// if(!Tout)Tout=B; |
1182 |
|
// else Tout->AddFriend("OrbitalInfo"); |
1183 |
|
// }; |
1184 |
|
|
1185 |
// cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl; |
// cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl; |
1186 |
|
|
1187 |
// if( Tout->GetEntries() )PamLevel2::SetBranchAddress(); |
// if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout ); |
|
if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout ); |
|
1188 |
|
|
1189 |
return Tout; |
return Trout; |
1190 |
} |
} |
1191 |
|
|
1192 |
|
|
1288 |
|
|
1289 |
// Tracker |
// Tracker |
1290 |
if(TRK2) { |
if(TRK2) { |
1291 |
t->SetBranchAddress("TrkLevel2", this->GetPointerTo("TrkLevel2")); |
t->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2")); |
1292 |
cout << "Tracker : set branch address TrkLevel2"<<endl; |
cout << "Tracker : set branch address TrkLevel2"<<endl; |
1293 |
}; |
}; |
1294 |
if(TRK1) { |
if(TRK1) { |
1343 |
|
|
1344 |
} |
} |
1345 |
|
|
|
// void* PamLevel2::GetPointerToTrk() {return &trk2_obj;}; |
|
|
// void* PamLevel2::GetPointerToTrk(int s){ |
|
|
// switch (s) { |
|
|
// case 1: return &trk1_obj; |
|
|
// case 2: return &trk2_obj; |
|
|
// case 3: return &trkh_obj; |
|
|
// default: return 0; |
|
|
// }; |
|
|
// }; |
|
|
// void* PamLevel2::GetPointerToCalo() {return &calo2_obj; }; |
|
|
// void* GetPointerToCalo(int s){ |
|
|
// switch (s) { |
|
|
// case 1: return &calo1_obj; |
|
|
// case 2: return &calo2_obj; |
|
|
// default: return 0; |
|
|
// }; |
|
|
// }; |
|
|
// void* GetPointerToToF() {return &tof_obj; }; |
|
|
// void* GetPointerTo("TrigLevel2") {return &trig_obj; }; |
|
|
// void* GetPointerTo("S4Level2") {return &s4_obj; }; |
|
|
// void* GetPointerTo("NDLevel2") {return &nd_obj; }; |
|
|
// void* GetPointerTo("AcLevel2") {return &ac_obj; }; |
|
|
// void* GetPointerTo("OrbitalInfo") {return &orb_obj; }; |
|
|
// void* GetPointerTo("RunInfo") {return &run_obj; }; |
|
1346 |
|
|
1347 |
//-------------------------------------- |
//-------------------------------------- |
1348 |
// |
// |
1404 |
*/ |
*/ |
1405 |
void PamLevel2::SetWhichTrees(TString detlist){ |
void PamLevel2::SetWhichTrees(TString detlist){ |
1406 |
|
|
1407 |
if(detlist.Contains("+ALL", TString::kIgnoreCase)){ |
if(detlist.IsNull() || detlist.Contains("+ALL", TString::kIgnoreCase)){ |
1408 |
CAL1 = true; |
CAL1 = true; |
1409 |
CAL2 = true; |
CAL2 = true; |
1410 |
TRK2 = true; |
TRK2 = true; |
1469 |
TRKh=false; |
TRKh=false; |
1470 |
} |
} |
1471 |
// ------------------------------------------------------------------------- |
// ------------------------------------------------------------------------- |
|
|
|
|
|
|
|
|
|
|
// if( detlist.Contains("-TRK", TString::kIgnoreCase) )TRK2 = false; |
|
|
// else if( detlist.Contains("+TRK", TString::kIgnoreCase) )TRK2 = true; |
|
|
|
|
|
// if( detlist.Contains("-TRK1", TString::kIgnoreCase) )TRK1 = false; |
|
|
// else if( detlist.Contains("+TRK1", TString::kIgnoreCase) )TRK1 = true; |
|
|
|
|
|
// if( detlist.Contains("-TRKh", TString::kIgnoreCase) )TRKh = false; |
|
|
// else if( detlist.Contains("+TRKh", TString::kIgnoreCase) )TRKh = true; |
|
1472 |
|
|
1473 |
if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false; |
if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false; |
1474 |
else if( detlist.Contains("+TRG", TString::kIgnoreCase) )TRG = true; |
else if( detlist.Contains("+TRG", TString::kIgnoreCase) )TRG = true; |
1487 |
|
|
1488 |
if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false; |
if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false; |
1489 |
else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true; |
else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true; |
1490 |
|
|
1491 |
|
// cout<< "Set detector list --> "; |
1492 |
|
// if(TRK1)cout<<"TRK1 "; |
1493 |
|
// if(TRK2)cout<<"TRK2 "; |
1494 |
|
// if(TRKh)cout<<"TRKH "; |
1495 |
|
// if(CAL1)cout<<"CAL1 "; |
1496 |
|
// if(CAL2)cout<<"CAL2 "; |
1497 |
|
// if(TOF)cout<<"TOF "; |
1498 |
|
// if(TRG)cout<<"TRG "; |
1499 |
|
// if(AC)cout<<"AC "; |
1500 |
|
// if(ND)cout<<"ND "; |
1501 |
|
// if(S4)cout<<"S4 "; |
1502 |
|
// if(ORB)cout<<"ORB "; |
1503 |
|
// cout << endl; |
1504 |
|
|
1505 |
}; |
}; |
1506 |
|
|
1507 |
|
|
1508 |
|
/** |
1509 |
|
* Set tree/branch detector flags from the content of a tree |
1510 |
|
*/ |
1511 |
|
void PamLevel2::GetWhichTrees(TFile* f){ |
1512 |
|
|
1513 |
|
// ----------- |
1514 |
|
// reset flags |
1515 |
|
// ----------- |
1516 |
|
CAL1 = false; |
1517 |
|
CAL2 = false; |
1518 |
|
TRK2 = false; |
1519 |
|
TRK1 = false; |
1520 |
|
TRKh = false; |
1521 |
|
TRG = false; |
1522 |
|
TOF = false; |
1523 |
|
S4 = false; |
1524 |
|
ND = false; |
1525 |
|
AC = false; |
1526 |
|
ORB = false; |
1527 |
|
|
1528 |
|
RUN = false; |
1529 |
|
|
1530 |
|
cout << "Checking file: "<<f->GetName()<<endl; |
1531 |
|
if( !f || f->IsZombie() ){ |
1532 |
|
cout << "File: "<< f->GetName() <<" Non valid root file"<< endl; |
1533 |
|
return; |
1534 |
|
} |
1535 |
|
|
1536 |
|
TList *lk = f->GetListOfKeys(); |
1537 |
|
TIter next(lk); |
1538 |
|
TKey *key =0; |
1539 |
|
|
1540 |
|
Int_t nev = 0; |
1541 |
|
|
1542 |
|
while( (key = (TKey*)next()) ){ |
1543 |
|
|
1544 |
|
if( !strcmp(key->GetName(),"Run" ) )RUN = true; |
1545 |
|
|
1546 |
|
//========================================================= |
1547 |
|
if( !strcmp(key->GetName(),"Trigger" ) ){ |
1548 |
|
TRG = true; |
1549 |
|
Int_t nevt = ((TTree*)f->Get("Trigger"))->GetEntries(); |
1550 |
|
if( nev && nevt!=nev){ |
1551 |
|
cout << "File: "<< f->GetName() <<" Trigger tree has "<<nevt<<" events instead of "<<nev<< endl; |
1552 |
|
TRG = false; |
1553 |
|
}else nev=nevt; |
1554 |
|
} |
1555 |
|
//========================================================= |
1556 |
|
if( !strcmp(key->GetName(),"ToF" ) ){ |
1557 |
|
TOF = true; |
1558 |
|
Int_t nevt = ((TTree*)f->Get("ToF"))->GetEntries(); |
1559 |
|
if( nev && nevt!=nev){ |
1560 |
|
cout << "File: "<< f->GetName() <<" ToF tree has "<<nevt<<" events instead of "<<nev<< endl; |
1561 |
|
TOF = false; |
1562 |
|
}else nev=nevt; |
1563 |
|
} |
1564 |
|
//========================================================= |
1565 |
|
if( !strcmp(key->GetName(),"S4" ) ){ |
1566 |
|
S4 = true; |
1567 |
|
Int_t nevt = ((TTree*)f->Get("S4"))->GetEntries(); |
1568 |
|
if( nev && nevt!=nev){ |
1569 |
|
cout << "File: "<< f->GetName() <<" S4 tree has "<<nevt<<" events instead of "<<nev<< endl; |
1570 |
|
S4 = false; |
1571 |
|
}else nev=nevt; |
1572 |
|
} |
1573 |
|
//========================================================= |
1574 |
|
|
1575 |
|
if( !strcmp(key->GetName(),"NeutronD" ) ){ |
1576 |
|
ND = true; |
1577 |
|
Int_t nevt = ((TTree*)f->Get("NeutronD"))->GetEntries(); |
1578 |
|
if( nev && nevt!=nev){ |
1579 |
|
cout << "File: "<< f->GetName() <<"NeutronD tree has "<<nevt<<" events instead of "<<nev<< endl; |
1580 |
|
ND =false; |
1581 |
|
}else nev=nevt; |
1582 |
|
} |
1583 |
|
//========================================================= |
1584 |
|
if( !strcmp(key->GetName(),"Anticounter") ){ |
1585 |
|
AC = true; |
1586 |
|
Int_t nevt = ((TTree*)f->Get("Anticounter"))->GetEntries(); |
1587 |
|
if( nev && nevt!=nev){ |
1588 |
|
cout << "File: "<< f->GetName() <<" Anticounter tree has "<<nevt<<" events instead of "<<nev<< endl; |
1589 |
|
AC =false; |
1590 |
|
}else nev=nevt; |
1591 |
|
} |
1592 |
|
//========================================================= |
1593 |
|
if( !strcmp(key->GetName(),"OrbitalInfo") ){ |
1594 |
|
ORB = true; |
1595 |
|
Int_t nevt = ((TTree*)f->Get("OrbitalInfo"))->GetEntries(); |
1596 |
|
if( nev && nevt!=nev){ |
1597 |
|
cout << "File: "<< f->GetName() <<" OrbitalInfo tree has "<<nevt<<" events instead of "<<nev<< endl; |
1598 |
|
ORB = false; |
1599 |
|
}else nev=nevt; |
1600 |
|
} |
1601 |
|
//========================================================= |
1602 |
|
if( !strcmp(key->GetName(),"Tracker" ) ){ |
1603 |
|
TTree *T = (TTree*)f->Get("Tracker"); |
1604 |
|
for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){ |
1605 |
|
TString name = T->GetListOfBranches()->At(i)->GetName(); |
1606 |
|
if( !name.CompareTo("TrkLevel1") )TRK1=true; |
1607 |
|
if( !name.CompareTo("TrkLevel2") )TRK2=true; |
1608 |
|
if( !name.CompareTo("TrkHough") )TRKh=true; |
1609 |
|
}; |
1610 |
|
Int_t nevt = T->GetEntries(); |
1611 |
|
if( nev && nevt!=nev){ |
1612 |
|
cout << "File: "<< f->GetName() <<" Tracker tree has "<<nevt<<" events instead of "<<nev<< endl; |
1613 |
|
TRK1 = false; |
1614 |
|
TRK2 = false; |
1615 |
|
TRKh = false; |
1616 |
|
}else nev=nevt; |
1617 |
|
T->Delete(); |
1618 |
|
}; |
1619 |
|
//========================================================= |
1620 |
|
if( !strcmp(key->GetName(),"Calorimeter" ) ){ |
1621 |
|
TTree *T = (TTree*)f->Get("Calorimeter"); |
1622 |
|
for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){ |
1623 |
|
TString name = T->GetListOfBranches()->At(i)->GetName(); |
1624 |
|
if( !name.CompareTo("CaloLevel1") )CAL1=true; |
1625 |
|
if( !name.CompareTo("CaloLevel2") )CAL2=true; |
1626 |
|
}; |
1627 |
|
Int_t nevt = T->GetEntries(); |
1628 |
|
if( nev && nevt!=nev){ |
1629 |
|
cout << "File: "<< f->GetName() <<" Calorimeter tree has "<<nevt<<" events instead of "<<nev<< endl; |
1630 |
|
CAL1 = false; |
1631 |
|
CAL2 = false; |
1632 |
|
}else nev=nevt; |
1633 |
|
T->Delete(); |
1634 |
|
}; |
1635 |
|
|
1636 |
|
}; |
1637 |
|
|
1638 |
|
delete lk; |
1639 |
|
|
1640 |
|
// cout<< "Get detector list from input file --> "; |
1641 |
|
// if(TRK1)cout<<"TRK1 "; |
1642 |
|
// if(TRK2)cout<<"TRK2 "; |
1643 |
|
// if(TRKh)cout<<"TRKH "; |
1644 |
|
// if(CAL1)cout<<"CAL1 "; |
1645 |
|
// if(CAL2)cout<<"CAL2 "; |
1646 |
|
// if(TOF)cout<<"TOF "; |
1647 |
|
// if(TRG)cout<<"TRG "; |
1648 |
|
// if(AC)cout<<"AC "; |
1649 |
|
// if(ND)cout<<"ND "; |
1650 |
|
// if(S4)cout<<"S4 "; |
1651 |
|
// if(ORB)cout<<"ORB "; |
1652 |
|
// cout << endl; |
1653 |
|
|
1654 |
|
return ; |
1655 |
|
|
1656 |
|
}; |
1657 |
|
|
1658 |
|
|
1659 |
//-------------------------------------- |
//-------------------------------------- |
1660 |
// |
// |
1661 |
// |
// |
1698 |
|
|
1699 |
while( (key = (TKey*)next()) ){ |
while( (key = (TKey*)next()) ){ |
1700 |
|
|
1701 |
// cout << key->GetName() << endl; |
// cout << key->GetName() << ""<<key->GetClassName()<<endl; |
1702 |
// cout << " Get tree: " << f->Get(key->GetName())<<endl; |
// cout << " Get tree: " << f->Get(key->GetName())<<endl; |
1703 |
// nev_previous = nev; |
// nev_previous = nev; |
1704 |
// cout << " n.entries "<< nev <<endl; |
// cout << " n.entries "<< nev <<endl; |
1705 |
// if( key->GetClassName()=="TTree" && nev_previous && nev != nev_previous ){ |
// if( key->GetClassName()=="TTree" && nev_previous && nev != nev_previous ){ |
1823 |
|
|
1824 |
}; |
}; |
1825 |
|
|
1826 |
|
// cout<< "CheckLevel2File(TString): detector list --> "; |
1827 |
|
// if(TRK1__ok)cout<<"TRK1 "; |
1828 |
|
// if(TRK2__ok)cout<<"TRK2 "; |
1829 |
|
// if(TRKh__ok)cout<<"TRKH "; |
1830 |
|
// if(CAL1__ok)cout<<"CAL1 "; |
1831 |
|
// if(CAL2__ok)cout<<"CAL2 "; |
1832 |
|
// if(TOF__ok)cout<<"TOF "; |
1833 |
|
// if(TRG__ok)cout<<"TRG "; |
1834 |
|
// if(AC__ok)cout<<"AC "; |
1835 |
|
// if(ND__ok)cout<<"ND "; |
1836 |
|
// if(S4__ok)cout<<"S4 "; |
1837 |
|
// if(ORB__ok)cout<<"ORB "; |
1838 |
|
// cout << endl; |
1839 |
|
|
1840 |
|
|
1841 |
|
if(TRK2 && TRK1__ok)TRK1=1; |
1842 |
|
// ---------------------------------------------------------------------------- |
1843 |
|
// NOTA |
1844 |
|
// se c'e` il level1, lo devo necessarimente leggere. |
1845 |
|
// infatti (non ho capito perche`) i cluster vengono letti e allocati in memoria |
1846 |
|
// comunque, ma non vengono disallocati da PamLevel2::Clear() |
1847 |
|
// ---------------------------------------------------------------------------- |
1848 |
|
|
1849 |
|
|
1850 |
if(!RUN__ok) { |
if(!RUN__ok) { |
1889 |
return false; |
return false; |
1890 |
}; |
}; |
1891 |
|
|
1892 |
if(ND && !ND__ok){ |
if(ND && !ND__ok){ |
1893 |
cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl; |
cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl; |
1894 |
return false; |
return false; |
1895 |
}; |
}; |
1896 |
if(TRG && !TRG__ok){ |
if(TRG && !TRG__ok){ |
1897 |
cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl; |
cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl; |
1898 |
return false; |
return false; |
1899 |
}; |
}; |
1900 |
|
|
|
|
|
|
|
|
|
|
|
1901 |
|
|
1902 |
lk->Delete(); |
// lk->Delete(); |
1903 |
|
// delete lk; |
1904 |
f->Close(); |
f->Close(); |
1905 |
|
|
1906 |
|
// cout<< "CheckLevel2File(TString): detector list --> "; |
1907 |
|
// if(TRK1)cout<<"TRK1 "; |
1908 |
|
// if(TRK2)cout<<"TRK2 "; |
1909 |
|
// if(TRKh)cout<<"TRKH "; |
1910 |
|
// if(CAL1)cout<<"CAL1 "; |
1911 |
|
// if(CAL2)cout<<"CAL2 "; |
1912 |
|
// if(TOF)cout<<"TOF "; |
1913 |
|
// if(TRG)cout<<"TRG "; |
1914 |
|
// if(AC)cout<<"AC "; |
1915 |
|
// if(ND)cout<<"ND "; |
1916 |
|
// if(S4)cout<<"S4 "; |
1917 |
|
// if(ORB)cout<<"ORB "; |
1918 |
|
// cout << endl; |
1919 |
|
|
1920 |
return true; |
return true; |
1921 |
|
|
1925 |
/** |
/** |
1926 |
* Create clone-trees |
* Create clone-trees |
1927 |
*/ |
*/ |
1928 |
void PamLevel2::CreateCloneTrees( TChain *fChain ){ |
void PamLevel2::CreateCloneTrees0( TChain *fChain, TFile *ofile ){ |
1929 |
|
|
1930 |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
1931 |
cout << "Create clones of PAMELA trees "<<endl; |
cout << "Create clones of PAMELA trees "<<endl; |
1941 |
TList *li = fChain->GetListOfFriends(); |
TList *li = fChain->GetListOfFriends(); |
1942 |
TIter next(li); |
TIter next(li); |
1943 |
TFriendElement* T_friend=0; |
TFriendElement* T_friend=0; |
1944 |
|
ofile->cd(); |
1945 |
while( (T_friend = (TFriendElement*)next()) ){ |
while( (T_friend = (TFriendElement*)next()) ){ |
1946 |
// cout<<T_friend->IsA()->GetName()<<" "<<T_friend->GetName()<<hex << T_friend->GetTree() << dec<<endl; |
// cout<<T_friend->IsA()->GetName()<<" "<<T_friend->GetName()<<hex << T_friend->GetTree() << dec<<endl; |
1947 |
// cout<<T_friend->GetTree()->GetName()<< endl; |
// cout<<T_friend->GetTree()->GetName()<< endl; |
1948 |
tree_clone[i] = T_friend->GetTree()->CloneTree(0); |
tree_clone[i] = T_friend->GetTree()->CloneTree(0); |
1949 |
|
tree_clone[i]->SetAutoSave(1000000); |
1950 |
name = tree_clone[i]->GetName(); |
name = tree_clone[i]->GetName(); |
1951 |
name.Append("_clone"); |
name.Append("_clone"); |
1952 |
// tree_clone[i]->SetName(name.Data()); |
// tree_clone[i]->SetName(name.Data()); |
1953 |
cout << tree_clone[i]->GetName() << endl; |
cout << tree_clone[i]->GetName() << endl; |
1954 |
i++; |
i++; |
1955 |
} |
} |
1956 |
|
|
1957 |
|
delete li; |
1958 |
|
|
1959 |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
1960 |
|
|
1961 |
} |
} |
1962 |
|
|
1963 |
// /** |
/** |
1964 |
// * Create clone-trees |
* Create clone-trees |
1965 |
// */ |
*/ |
1966 |
// void PamLevel2::CreateCloneTrees(){ |
void PamLevel2::CreateCloneTrees(TFile *ofile){ |
|
|
|
|
// cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
|
|
// cout << "Create new PAMELA trees "<<endl; |
|
|
|
|
|
// if(TRK1||TRK2||TRKh){ |
|
|
// T_clone = new TTree("Tracker_clone","PAMELA tracker level2 data "); |
|
|
// if(TRK1) { |
|
|
// trk1_clone = new TrkLevel1(); |
|
|
// T_clone->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1")); |
|
|
// T_clone->BranchRef(); |
|
|
// cout << "Tracker : branch TrkLevel1"<<endl; |
|
|
// }; |
|
|
// if(TRK2) { |
|
|
// trk2_clone = new TrkLevel2(); |
|
|
// T_clone->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2")); |
|
|
// cout << "Tracker : branch TrkLevel2"<<endl; |
|
|
// }; |
|
|
// if(TRKh) { |
|
|
// trkh_clone = new TrkHough(); |
|
|
// T_clone->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough")); |
|
|
// cout << "Tracker : branch TrkHough"<<endl; |
|
|
// }; |
|
|
// cout << hex << T_clone << dec << endl; |
|
|
// } |
|
|
|
|
|
// // Calorimeter |
|
|
// if(CAL1||CAL2){ |
|
|
// C_clone = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data "); |
|
|
// if(CAL1) { |
|
|
// calo1_clone = new CaloLevel1(); |
|
|
// C_clone->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1")); |
|
|
// cout << "Calorimeter : branch CaloLevel1"<<endl; |
|
|
// }; |
|
|
// if(CAL2) { |
|
|
// calo2_clone = new CaloLevel2(); |
|
|
// C_clone->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2")); |
|
|
// cout << "Calorimeter : branch CaloLevel2"<<endl; |
|
|
// }; |
|
|
// } |
|
1967 |
|
|
1968 |
// // ToF |
ofile->cd(); |
|
// if(TOF) { |
|
|
// O_clone = new TTree("ToF_clone","PAMELA ToF level2 data "); |
|
|
// tof_clone = new ToFLevel2(); |
|
|
// O_clone->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2")); |
|
|
// cout << "ToF : branch ToFLevel2"<<endl; |
|
|
// }; |
|
|
// // Trigger |
|
|
// if(TRG) { |
|
|
// R_clone = new TTree("Trigger_clone","PAMELA trigger level2 data "); |
|
|
// trig_clone = new TrigLevel2(); |
|
|
// R_clone->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2")); |
|
|
// cout << "Trigger : branch TrigLevel2"<<endl; |
|
|
// }; |
|
|
// // S4 |
|
|
// if(S4) { |
|
|
// S_clone = new TTree("S4_clone","PAMELA S4 level2 data "); |
|
|
// s4_clone = new S4Level2(); |
|
|
// S_clone->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2")); |
|
|
// cout << "S4 : branch S4Level2"<<endl; |
|
|
// }; |
|
|
// // Neutron Detector |
|
|
// if(ND) { |
|
|
// N_clone = new TTree("NeutronD_clone","PAMELA neutron detector level2 data "); |
|
|
// nd_clone = new NDLevel2(); |
|
|
// N_clone->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2")); |
|
|
// cout << "NeutronD : branch NDLevel2"<<endl; |
|
|
// }; |
|
|
// // Anticounters |
|
|
// if(AC) { |
|
|
// A_clone = new TTree("Anticounter_clone","PAMELA anticounter detector level2 data "); |
|
|
// ac_clone = new AcLevel2(); |
|
|
// A_clone->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2")); |
|
|
// cout << "Anticounter : branch AcLevel2"<<endl; |
|
|
// }; |
|
|
// // OrbitalInfo |
|
|
// if(ORB) { |
|
|
// B_clone = new TTree("OrbitalInfo_clone","PAMELA oribital info "); |
|
|
// orb_clone = new OrbitalInfo(); |
|
|
// B_clone->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo")); |
|
|
// cout << "OrbitalInfo : branch OrbitalInfo"<<endl; |
|
|
// }; |
|
|
// cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
|
1969 |
|
|
1970 |
// } |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
1971 |
|
cout << "Create new PAMELA trees "<<endl; |
1972 |
|
|
1973 |
|
Int_t i=0; |
1974 |
|
|
1975 |
// /** |
if(TRK1||TRK2||TRKh){ |
1976 |
// * Create clone-trees |
tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data "); |
1977 |
// */ |
if(TRK1) { |
1978 |
// void PamLevel2::CreateCloneTrees(){ |
tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1")); |
1979 |
|
tree_clone[i]->BranchRef(); |
1980 |
|
cout << "Tracker : branch TrkLevel1"<<endl; |
1981 |
|
}; |
1982 |
|
if(TRK2) { |
1983 |
|
tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2")); |
1984 |
|
cout << "Tracker : branch TrkLevel2"<<endl; |
1985 |
|
}; |
1986 |
|
if(TRKh) { |
1987 |
|
tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough")); |
1988 |
|
cout << "Tracker : branch TrkHough"<<endl; |
1989 |
|
}; |
1990 |
|
i++; |
1991 |
|
} |
1992 |
|
|
1993 |
// cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
// Calorimeter |
1994 |
// cout << "Create new PAMELA trees "<<endl; |
if(CAL1||CAL2){ |
1995 |
|
tree_clone[i] = new TTree("Calorimeter","PAMELA calorimeter level2 data "); |
1996 |
// if(TRK1||TRK2||TRKh){ |
if(CAL1) { |
1997 |
// T_clone = new TTree("Tracker_clone","PAMELA tracker level2 data "); |
tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1")); |
1998 |
// if(TRK1) { |
cout << "Calorimeter : branch CaloLevel1"<<endl; |
1999 |
// trk1_clone = new TrkLevel1(); |
}; |
2000 |
// T_clone->Branch("TrkLevel1","TrkLevel1", &trk1_clone); |
if(CAL2) { |
2001 |
// T_clone->BranchRef(); |
tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2")); |
2002 |
// cout << "Tracker : branch TrkLevel1"<<endl; |
cout << "Calorimeter : branch CaloLevel2"<<endl; |
2003 |
// }; |
}; |
2004 |
// if(TRK2) { |
i++; |
2005 |
// trk2_clone = new TrkLevel2(); |
} |
|
// T_clone->Branch("TrkLevel2", "TrkLevel2",&trk2_clone); |
|
|
// cout << "Tracker : branch TrkLevel2"<<endl; |
|
|
// }; |
|
|
// if(TRKh) { |
|
|
// trkh_clone = new TrkHough(); |
|
|
// T_clone->Branch("TrkHough","TrkHough", &trkh_clone); |
|
|
// cout << "Tracker : branch TrkHough"<<endl; |
|
|
// }; |
|
|
// cout << hex << T_clone << dec << endl; |
|
|
// } |
|
2006 |
|
|
2007 |
// // Calorimeter |
// ToF |
2008 |
// if(CAL1||CAL2){ |
if(TOF) { |
2009 |
// C_clone = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data "); |
tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data "); |
2010 |
// if(CAL1) { |
tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2")); |
2011 |
// calo1_clone = new CaloLevel1(); |
cout << "ToF : branch ToFLevel2"<<endl; |
2012 |
// C_clone->Branch("CaloLevel1", "CaloLevel1", &calo1_clone); |
i++; |
2013 |
// cout << "Calorimeter : branch CaloLevel1"<<endl; |
}; |
2014 |
// }; |
// Trigger |
2015 |
// if(CAL2) { |
if(TRG) { |
2016 |
// calo2_clone = new CaloLevel2(); |
tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data "); |
2017 |
// C_clone->Branch("CaloLevel2","CaloLevel2", &calo2_clone); |
tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2")); |
2018 |
// cout << "Calorimeter : branch CaloLevel2"<<endl; |
cout << "Trigger : branch TrigLevel2"<<endl; |
2019 |
// }; |
i++; |
2020 |
// } |
}; |
2021 |
|
// S4 |
2022 |
|
if(S4) { |
2023 |
|
tree_clone[i] = new TTree("S4","PAMELA S4 level2 data "); |
2024 |
|
tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2")); |
2025 |
|
cout << "S4 : branch S4Level2"<<endl; |
2026 |
|
i++; |
2027 |
|
}; |
2028 |
|
// Neutron Detector |
2029 |
|
if(ND) { |
2030 |
|
tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data "); |
2031 |
|
tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2")); |
2032 |
|
cout << "NeutronD : branch NDLevel2"<<endl; |
2033 |
|
i++; |
2034 |
|
}; |
2035 |
|
// Anticounters |
2036 |
|
if(AC) { |
2037 |
|
tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data "); |
2038 |
|
tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2")); |
2039 |
|
cout << "Anticounter : branch AcLevel2"<<endl; |
2040 |
|
i++; |
2041 |
|
}; |
2042 |
|
// OrbitalInfo |
2043 |
|
if(ORB) { |
2044 |
|
tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info "); |
2045 |
|
tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo")); |
2046 |
|
cout << "OrbitalInfo : branch OrbitalInfo"<<endl; |
2047 |
|
i++; |
2048 |
|
}; |
2049 |
|
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
2050 |
|
|
2051 |
// // ToF |
} |
|
// if(TOF) { |
|
|
// O_clone = new TTree("ToF_clone","PAMELA ToF level2 data "); |
|
|
// tof_clone = new ToFLevel2(); |
|
|
// O_clone->Branch("ToFLevel2","ToFLevel2", &tof_clone); |
|
|
// cout << "ToF : branch ToFLevel2"<<endl; |
|
|
// }; |
|
|
// // Trigger |
|
|
// if(TRG) { |
|
|
// R_clone = new TTree("Trigger_clone","PAMELA trigger level2 data "); |
|
|
// trig_clone = new TrigLevel2(); |
|
|
// R_clone->Branch("TrigLevel2","TrigLevel2", &trig_clone); |
|
|
// cout << "Trigger : branch TrigLevel2"<<endl; |
|
|
// }; |
|
|
// // S4 |
|
|
// if(S4) { |
|
|
// S_clone = new TTree("S4_clone","PAMELA S4 level2 data "); |
|
|
// s4_clone = new S4Level2(); |
|
|
// S_clone->Branch("S4Level2","S4Level2", &s4_clone); |
|
|
// cout << "S4 : branch S4Level2"<<endl; |
|
|
// }; |
|
|
// // Neutron Detector |
|
|
// if(ND) { |
|
|
// N_clone = new TTree("NeutronD_clone","PAMELA neutron detector level2 data "); |
|
|
// nd_clone = new NDLevel2(); |
|
|
// N_clone->Branch("NDLevel2","NDLevel2", &nd_clone); |
|
|
// cout << "NeutronD : branch NDLevel2"<<endl; |
|
|
// }; |
|
|
// // Anticounters |
|
|
// if(AC) { |
|
|
// A_clone = new TTree("Anticounter_clone","PAMELA anticounter detector level2 data "); |
|
|
// ac_clone = new AcLevel2(); |
|
|
// A_clone->Branch("AcLevel2","AcLevel2", &ac_clone); |
|
|
// cout << "Anticounter : branch AcLevel2"<<endl; |
|
|
// }; |
|
|
// // OrbitalInfo |
|
|
// if(ORB) { |
|
|
// B_clone = new TTree("OrbitalInfo_clone","PAMELA oribital info "); |
|
|
// orb_clone = new OrbitalInfo(); |
|
|
// B_clone->Branch("OrbitalInfo","OrbitalInfo", &orb_clone); |
|
|
// cout << "OrbitalInfo : branch OrbitalInfo"<<endl; |
|
|
// }; |
|
|
// cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl; |
|
2052 |
|
|
|
// } |
|
2053 |
|
|
2054 |
/** |
/** |
2055 |
* Fill tree (created with CreateCloneTrees) |
* Fill tree (created with CreateCloneTrees) |
2058 |
//void PamLevel2::FillNewPamTree(TTree *T){ |
//void PamLevel2::FillNewPamTree(TTree *T){ |
2059 |
void PamLevel2::FillCloneTrees(){ |
void PamLevel2::FillCloneTrees(){ |
2060 |
|
|
|
// // -------------------------------------- |
|
|
// // copy the event |
|
|
// // -------------------------------------- |
|
|
// if(trk1_clone) *trk1_clone = *trk1_obj; |
|
|
// if(trk2_clone){ |
|
|
// trk2_clone->Clear(); |
|
|
// trk2_obj->Copy(*trk2_clone); |
|
|
// // *trk2_clone = *trk2_obj; |
|
|
// } |
|
|
// if(trkh_clone) *trkh_clone = *trkh_obj; |
|
|
// if(calo1_clone){ |
|
|
// // *calo1_clone = *calo1_obj; |
|
|
// calo1_clone->Clear(); |
|
|
// calo1_obj->Copy(*calo1_clone); |
|
|
// } |
|
|
// if(calo2_clone){ |
|
|
// // *calo2_clone = *calo2_obj; |
|
|
// calo2_clone->Clear(); |
|
|
// calo2_obj->Copy(*calo2_clone); |
|
|
// } |
|
|
// if(tof_clone) *tof_clone = *tof_obj; |
|
|
// if(trig_clone) *trig_clone = *trig_obj; |
|
|
// if(s4_clone) *s4_clone = *s4_obj; |
|
|
// if(nd_clone) *nd_clone = *nd_obj; |
|
|
// if(ac_clone) *ac_clone = *ac_obj; |
|
|
// if(orb_clone) *orb_clone = *orb_obj; |
|
|
|
|
|
// if(T_clone)T_clone->Fill(); |
|
|
// if(C_clone)C_clone->Fill(); |
|
|
// if(O_clone)O_clone->Fill(); |
|
|
// if(R_clone)R_clone->Fill(); |
|
|
// if(S_clone)S_clone->Fill(); |
|
|
// if(N_clone)N_clone->Fill(); |
|
|
// if(A_clone)A_clone->Fill(); |
|
|
// if(O_clone)O_clone->Fill(); |
|
|
|
|
2061 |
for(Int_t i=0; i<8; i++){ |
for(Int_t i=0; i<8; i++){ |
2062 |
if(tree_clone[i])tree_clone[i]->Fill(); |
if(tree_clone[i])tree_clone[i]->Fill(); |
2063 |
} |
} |