/[PAMELA software]/PamelaLevel2/src/PamLevel2.cpp
ViewVC logotype

Diff of /PamelaLevel2/src/PamLevel2.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.10 by pam-fi, Wed Dec 6 11:07:35 2006 UTC revision 1.18 by pam-fi, Thu Jan 18 08:51:16 2007 UTC
# Line 7  Line 7 
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  //  //
# Line 19  PamTrack::PamTrack(){ Line 20  PamTrack::PamTrack(){
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  //  //
# Line 36  PamTrack::PamTrack(TrkTrack* t, CaloTrkV Line 89  PamTrack::PamTrack(TrkTrack* t, CaloTrkV
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;      soft_obj   = 0;// Emiliano
117        irun = -1;
118  //     T_clone  = 0;      runfirstentry = 0ULL;
119  //     C_clone  = 0;      runlastentry = 0ULL;
 //     O_clone  = 0;  
 //     R_clone  = 0;  
 //     S_clone  = 0;  
 //     N_clone  = 0;  
 //     A_clone  = 0;  
 //     B_clone  = 0;  
120    
121      for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL;      for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL;
122    
123      sorted_tracks = new TRefArray();  //    sorted_tracks = 0;//new TRefArray();
124            
125      CAL1 = true;      CAL1 = true;
126      CAL2 = true;      CAL2 = true;
# Line 87  PamLevel2::PamLevel2(){ Line 134  PamLevel2::PamLevel2(){
134      AC  = true;      AC  = true;
135      ORB = true;      ORB = true;
136            
137        RUN = true;
138    
139        tsorted=0;
140        timage=0;
141    
142  };  };
143  /**  /**
144   * Destructor   * Delete the event (act as Dtor)
145   */   */
146  PamLevel2::~PamLevel2(){                  void PamLevel2::Delete(){
147      delete run_obj;          
148      Clear();      if(run_obj)delete run_obj;
149        if(soft_obj)delete soft_obj; //Emiliano
150    
151    //    cout << "void PamLevel2::Clear()"<<endl;
152        if(trk1_obj) delete trk1_obj;      
153        if(trk2_obj) delete trk2_obj;
154        if(trkh_obj) delete trkh_obj;
155        if(calo1_obj)delete calo1_obj;
156        if(calo2_obj)delete calo2_obj;
157        if(tof_obj)  delete tof_obj;
158        if(trig_obj) delete trig_obj;
159        if(s4_obj)   delete s4_obj;
160        if(nd_obj)   delete nd_obj;
161        if(ac_obj)   delete ac_obj;
162        if(orb_obj)  delete orb_obj;
163        
164        if(tsorted){
165            tsorted->Delete();
166            delete tsorted;
167        }
168        if(timage){
169            timage->Delete();
170            delete timage;
171        }
172  };  };
173    
174  /**  /**
175   * Clear the event   * Clear the event (NB! does not deallocate objects)
176   */   */
177  void PamLevel2::Clear(){  void PamLevel2::Clear(){
178                    
179      TrkLevel1::Clear();  //    cout << "void PamLevel2::Clear()"<<endl;
180            
181      TrkLevel2::Clear();  //
182      CaloLevel2::Clear();  // This method is called once for every entry but RunInfo and SoftInfo do not change until the next run so we cannot clear them here unless we don't
183      ToFLevel2::Clear();  // want to load them for each event even if they are the same...
184      TrigLevel2::Clear();  //
185      S4Level2::Clear();  //    if(run_obj)delete run_obj;
186      NDLevel2::Clear();  //    if(run_obj) run_obj->Clear();  // Emiliano: Do not deallocate run_obj here, it will give segmentation fault! call clear instead
187      AcLevel2::Clear();  //    if(soft_obj) soft_obj->Clear();
188      OrbitalInfo::Clear();  
189        if(trk1_obj) trk1_obj->Clear();    
190  //     if(trk1_clone)trk1_clone->Clear();      if(trk2_obj) trk2_obj->Clear();
191  //     if(trk2_clone)trk2_clone->Clear();      if(trkh_obj) trkh_obj->Clear();
192  //     if(trkh_clone)trkh_clone->Clear();      if(calo1_obj)calo1_obj->Clear();
193  //     if(calo1_clone)calo1_clone->Clear();      if(calo2_obj)calo2_obj->Clear();
194  //     if(calo2_clone)calo2_clone->Clear();      if(tof_obj)  tof_obj->Clear();
195  //     if(tof_clone)tof_clone->Clear();      if(trig_obj) trig_obj->Clear();
196  //     if(trig_clone)trig_clone->Clear();      if(s4_obj)   s4_obj->Clear();
197  //     if(s4_clone)s4_clone->Clear();      if(nd_obj)   nd_obj->Clear();
198  //     if(nd_clone)nd_clone->Clear();      if(ac_obj)   ac_obj->Clear();
199  //     if(ac_clone)ac_clone->Clear();      if(orb_obj)  orb_obj->Clear();
 //     if(orb_clone)orb_clone->Clear();  
200            
201      sorted_tracks->Delete(); // clean the reference array  //    if(sorted_tracks)sorted_tracks->Clear();
202    //    sorted_tracks.Clear();
203    
204        if(tsorted){
205            tsorted->Delete();
206        }
207        if(timage){
208            timage->Delete();
209        }
210            
211  };  };
212    
# Line 136  void PamLevel2::Clear(){ Line 218  void PamLevel2::Clear(){
218  void *PamLevel2::GetPointerTo(const char* c ){  void *PamLevel2::GetPointerTo(const char* c ){
219    
220      TString objname = c;      TString objname = c;
221      if(!objname.CompareTo("TrkLevel1"))return &trk1_obj;  
222      if(!objname.CompareTo("TrkLevel2"))return &trk2_obj;      if(!objname.CompareTo("TrkLevel1"))  {
223      if(!objname.CompareTo("TrkHough"))return &trkh_obj;          if(!trk1_obj){
224      if(!objname.CompareTo("CaloLevel1"))return &calo1_obj;              trk1_obj  = new TrkLevel1();
225      if(!objname.CompareTo("CaloLevel2"))return &calo2_obj;              trk1_obj->Set();
226      if(!objname.CompareTo("ToFLevel2"))return &tof_obj;          }
227      if(!objname.CompareTo("TrigLevel2"))return &trig_obj;          return &trk1_obj;
228      if(!objname.CompareTo("S4Level2"))return &s4_obj;      };
229      if(!objname.CompareTo("NDLevel2"))return &nd_obj;      if(!objname.CompareTo("TrkLevel2"))  {
230      if(!objname.CompareTo("AcLevel2"))return &ac_obj;          if(!trk2_obj){  
231      if(!objname.CompareTo("OrbitalInfo"))return &orb_obj;              trk2_obj  = new TrkLevel2();
232                trk2_obj->Set();
233            }
234            return &trk2_obj;
235        };
236        if(!objname.CompareTo("TrkHough"))   {
237            if(!trkh_obj)  trkh_obj  = new TrkHough();
238            return &trkh_obj;
239        };
240        if(!objname.CompareTo("CaloLevel1")) {
241            if(!calo1_obj) calo1_obj = new CaloLevel1();
242            return &calo1_obj;
243        };
244        if(!objname.CompareTo("CaloLevel2")) {
245            if(!calo2_obj){
246                calo2_obj = new CaloLevel2();
247                calo2_obj->Set();
248            };
249            return &calo2_obj;
250        };
251        if(!objname.CompareTo("ToFLevel2"))  {
252            if(!tof_obj){
253                tof_obj   = new ToFLevel2();
254                tof_obj->Set();
255            }
256            return &tof_obj;
257        };
258        if(!objname.CompareTo("TrigLevel2")) {
259            if(!trig_obj)  trig_obj  = new TrigLevel2();
260            return &trig_obj;
261        };
262        if(!objname.CompareTo("S4Level2"))   {
263            if(!s4_obj)    s4_obj    = new S4Level2();
264            return &s4_obj;
265        };
266        if(!objname.CompareTo("NDLevel2"))   {
267            if(!nd_obj)    nd_obj    = new NDLevel2();
268            return &nd_obj;
269        };
270        if(!objname.CompareTo("AcLevel2"))   {
271            if(!ac_obj)    ac_obj    = new AcLevel2();
272            return &ac_obj;
273        };
274        if(!objname.CompareTo("OrbitalInfo")){
275            if(!orb_obj)   orb_obj   = new OrbitalInfo();
276            return &orb_obj;
277        };
278        
279      if(!objname.CompareTo("RunInfo"))return &run_obj;      if(!objname.CompareTo("RunInfo"))return &run_obj;
280    
281        if(!objname.CompareTo("SoftInfo"))return &soft_obj; // Emiliano
282    
283      return NULL;      return NULL;
284  };  };
285  //--------------------------------------  //--------------------------------------
# Line 159  void *PamLevel2::GetPointerTo(const char Line 291  void *PamLevel2::GetPointerTo(const char
291   * (If seqno = -1 retrieves the self-trigger calorimeter track)   * (If seqno = -1 retrieves the self-trigger calorimeter track)
292   */   */
293   CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){   CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){
294                
295       if( CaloLevel2::ntrk()==0 ){       if( !calo2_obj )return 0;
296    
297         if( calo2_obj->CaloLevel2::ntrk()==0 ){
298           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;
299           return NULL;           return NULL;
300       };       };
# Line 169  void *PamLevel2::GetPointerTo(const char Line 303  void *PamLevel2::GetPointerTo(const char
303       Int_t it_calo=0;       Int_t it_calo=0;
304            
305       do{       do{
306           c = CaloLevel2::GetCaloTrkVar(it_calo);           c = calo2_obj->CaloLevel2::GetCaloTrkVar(it_calo);
307           it_calo++;           it_calo++;
308       } while( c && seqno != c->trkseqno && it_calo < CaloLevel2::ntrk());             } while( c && seqno != c->trkseqno && it_calo < calo2_obj->CaloLevel2::ntrk());    
309            
310       if(!c || seqno != c->trkseqno){       if(!c || seqno != c->trkseqno){
311           c = 0;           c = 0;
# Line 190  void *PamLevel2::GetPointerTo(const char Line 324  void *PamLevel2::GetPointerTo(const char
324   */   */
325   ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){   ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){
326                    
327       if( ToFLevel2::ntrk()==0 ){       if( !tof_obj )return 0;
328    
329         if( tof_obj->ToFLevel2::ntrk()==0 ){
330           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;
331           return NULL;           return NULL;
332       };       };
# Line 199  void *PamLevel2::GetPointerTo(const char Line 335  void *PamLevel2::GetPointerTo(const char
335       Int_t it_tof=0;       Int_t it_tof=0;
336            
337       do{       do{
338           c = ToFLevel2::GetToFTrkVar(it_tof);           c = tof_obj->ToFLevel2::GetToFTrkVar(it_tof);
339           it_tof++;           it_tof++;
340       } while( c && seqno != c->trkseqno && it_tof < ToFLevel2::ntrk());       } while( c && seqno != c->trkseqno && it_tof < tof_obj->ToFLevel2::ntrk());        
341            
342       if(!c || seqno != c->trkseqno){       if(!c || seqno != c->trkseqno){
343           c = 0;           c = 0;
# Line 219  void *PamLevel2::GetPointerTo(const char Line 355  void *PamLevel2::GetPointerTo(const char
355   * 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.
356   */   */
357   PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){   PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){
358                
359         cout <<"PamLevel2::GetPamTrackAlong(TrkTrack* t) **obsolete** "<<endl;
360         cout <<"(if you use it, remember to delete the PamTrack object)"<<endl;
361    
362       CaloTrkVar *c = 0;       CaloTrkVar *c = 0;
363       ToFTrkVar  *o = 0;       ToFTrkVar  *o = 0;
364            
# Line 244  void *PamLevel2::GetPointerTo(const char Line 383  void *PamLevel2::GetPointerTo(const char
383    
384  PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){  PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){
385            
386        cout <<"PamLevel2::GetStoredTrack(Int_t itrk) **to-be-updated** "<<endl;
387        cout <<"for the moment, better use separately the methods: TrkLevel2::GetStoredTrack(seqno) CaloLevel2::GetCaloTrkVar(Int_t notrack) ToFLevel2::GetToFTrkVar(Int_t notrack)"<<endl;
388        cout <<"(if you use it, remember to delete the PamTrack object)"<<endl;
389      PamTrack *track = 0;      PamTrack *track = 0;
390            
391      if( itrk >=0 && itrk < TrkLevel2::ntrk() ){      if( itrk >=0 && itrk < trk2_obj->TrkLevel2::ntrk() ){
392                    
393          TrkTrack *t = TrkLevel2::GetStoredTrack(itrk);          TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(itrk);
394          track = GetPamTrackAlong(t);          track = GetPamTrackAlong(t);
395                    
396            
397    
398      }else{      }else{
399          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;
400      };      };
401            
402      return track;      return track;
# Line 275  PamTrack* PamLevel2::GetStoredTrack(Int_ Line 419  PamTrack* PamLevel2::GetStoredTrack(Int_
419   */   */
420  void PamLevel2::SortTracks(TString how){  void PamLevel2::SortTracks(TString how){
421    
422    //    cout <<" PamLevel2::SortTracks(TString how) "<<endl;
423        if( !trk2_obj ){
424            cout << "void PamLevel2::SortTracks(TString how):  TrkLevel2 not loaded !!!";
425            return;
426        };
427  //    cout << "call SortTracs() "<<endl;  //    cout << "call SortTracs() "<<endl;
428    //Save current Object count    //Save current Object count
429      Int_t ObjectNumber = TProcessID::GetObjectCount();      Int_t ObjectNumber = TProcessID::GetObjectCount();
430    
431  //    cout << "ObjectNumber  "<<ObjectNumber <<endl;  //    cout << "ObjectNumber  "<<ObjectNumber <<endl;
432        
433      sorted_tracks->Delete(); //temporaneo???  //     if(!sorted_tracks)sorted_tracks = new TRefArray();
434            //     sorted_tracks->Clear();
435    //    sorted_tracks.Clear();
436    
437        if(!tsorted)tsorted = new TClonesArray("PamTrack",trk2_obj->GetNTracks());
438        tsorted->Delete();
439        TClonesArray &ttsorted = *tsorted;
440        if(!timage)timage = new TClonesArray("PamTrack",trk2_obj->GetNTracks());
441        timage->Delete();
442        TClonesArray &ttimage = *timage;
443    
444    
445      // loop over the tracks sorted by the tracker      // loop over the tracks sorted by the tracker
446      Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase);      Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase);
447      Bool_t use_CAL = how.Contains("CAL", TString::kIgnoreCase);      Bool_t use_CAL = how.Contains("CAL", TString::kIgnoreCase);
# Line 298  void PamLevel2::SortTracks(TString how){ Line 457  void PamLevel2::SortTracks(TString how){
457    
458      //   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;
459            
460      for(Int_t i=0; i < TrkLevel2::GetNTracks(); i++){      for(Int_t i=0; i < trk2_obj->TrkLevel2::GetNTracks(); i++){
461                    
462          TrkTrack *ts = 0;          TrkTrack   *ts = 0;
463            CaloTrkVar *cs = 0;
464            ToFTrkVar  *os = 0;
465                    
466          // get tracker tracks          // get tracker tracks
467          TrkTrack   *tp = TrkLevel2::GetTrack(i);                    //tracker          TrkTrack   *tp = trk2_obj->TrkLevel2::GetTrack(i);                    //tracker
468          CaloTrkVar *cp = 0;          CaloTrkVar *cp = GetCaloStoredTrack(tp->GetSeqNo());
469          ToFTrkVar  *op = 0;          ToFTrkVar  *op = GetToFStoredTrack(tp->GetSeqNo());
470    
471            TrkTrack   *ti = 0;              //tracker (image)
472            CaloTrkVar *ci = 0;
473            ToFTrkVar  *oi = 0;
474  //      cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl;  //      cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl;
475          // if track has an image, check image selection          // if track has an image, check image selection
476          if(tp->HasImage()){          if(tp->HasImage()){
477                            
478              TrkTrack   *ti = TrkLevel2::GetTrackImage(i);              //tracker (image)              ti = trk2_obj->TrkLevel2::GetTrackImage(i);              //tracker (image)
479              CaloTrkVar *ci = 0;              ci = GetCaloStoredTrack(ti->GetSeqNo());
480              ToFTrkVar  *oi = 0;              oi = GetToFStoredTrack(ti->GetSeqNo());
481                            
482  //          cout << "its image "<<i << " "<<hex<< ti <<dec<< endl;  //          cout << "its image "<<i << " "<<hex<< ti <<dec<< endl;
483    
# Line 326  void PamLevel2::SortTracks(TString how){ Line 490  void PamLevel2::SortTracks(TString how){
490              // ------------------------              // ------------------------
491              // check the Y spatial residual on the first calorimeter plane              // check the Y spatial residual on the first calorimeter plane
492              // (cut on calorimeter variables from Emiliano)              // (cut on calorimeter variables from Emiliano)
493                if( use_CAL && !calo2_obj ){
494                    cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but CaloLevel2 not loaded !!!";
495                    return;
496                };
497              if(              if(
498                  use_CAL            &&                  use_CAL            &&
499                  npcfit[1] > 15     &&   //no. of fit planes on Y view                  calo2_obj->npcfit[1] > 15     &&   //no. of fit planes on Y view
500                  varcfit[1] < 1000. &&  //fit variance on Y view                  calo2_obj->varcfit[1] < 1000. &&  //fit variance on Y view
501                  true){                  true){
502    
                 cp = GetCaloStoredTrack(tp->GetSeqNo());  
                 ci = GetCaloStoredTrack(ti->GetSeqNo());  
503                                    
504                  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;
505                  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;
506                                    
507                  if(resy_p <= resy_i) tp_score++;                  if(resy_p <= resy_i) tp_score++;
508                  else                 ti_score++;                  else                 ti_score++;
# Line 349  void PamLevel2::SortTracks(TString how){ Line 515  void PamLevel2::SortTracks(TString how){
515              // ------------------------                  // ------------------------    
516              // check the number of hit pmts along the track              // check the number of hit pmts along the track
517              // on S12 S21 and S32, where paddles are parallel to Y axis              // on S12 S21 and S32, where paddles are parallel to Y axis
518                if( use_TOF && !tof_obj ){
519                    cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but ToFLevel2 not loaded !!!";
520                    return;
521                };
522              if( use_TOF ){              if( use_TOF ){
523                                    
524                  Int_t nphit_p =0;                  Int_t nphit_p =0;
525                  Int_t nphit_i =0;                  Int_t nphit_i =0;
526                                    
                 op = GetToFStoredTrack(tp->GetSeqNo());  
                 oi = GetToFStoredTrack(ti->GetSeqNo());  
527                                    
528  /*                              cout << "track: npmtadc "<< op->npmtadc << endl;  /*                              cout << "track: npmtadc "<< op->npmtadc << endl;
529                                  cout << "track: npmttdc "<< op->npmttdc << endl;                                  cout << "track: npmttdc "<< op->npmttdc << endl;
# Line 363  void PamLevel2::SortTracks(TString how){ Line 531  void PamLevel2::SortTracks(TString how){
531                                  cout << "image: npmttdc "<< oi->npmttdc << endl;*/                                  cout << "image: npmttdc "<< oi->npmttdc << endl;*/
532                                    
533                  for (Int_t ih=0; ih < op->npmtadc; ih++){                  for (Int_t ih=0; ih < op->npmtadc; ih++){
534                      Int_t pl = GetPlaneIndex( (op->pmtadc).At(ih) );                      Int_t pl = tof_obj->GetPlaneIndex( (op->pmtadc).At(ih) );
535                      if(pl == 1 || pl == 2 || pl == 5)nphit_p++;                      if(pl == 1 || pl == 2 || pl == 5)nphit_p++;
536                  };                  };
537                                    
538                  for (Int_t ih=0; ih < oi->npmtadc; ih++){                  for (Int_t ih=0; ih < oi->npmtadc; ih++){
539                      Int_t pl = GetPlaneIndex( (oi->pmtadc).At(ih) );                      Int_t pl = tof_obj->GetPlaneIndex( (oi->pmtadc).At(ih) );
540                      if(pl == 1 || pl == 2 || pl == 5)nphit_i++;                      if(pl == 1 || pl == 2 || pl == 5)nphit_i++;
541                  };                  };
542                                    
# Line 398  void PamLevel2::SortTracks(TString how){ Line 566  void PamLevel2::SortTracks(TString how){
566              // ------------------------              // ------------------------
567              // the winner is....              // the winner is....
568              // ------------------------                  // ------------------------    
569              if      (tp_score > ti_score) ts = tp;//the track sorted by the tracker!!              if      (tp_score > ti_score) {
570              else if (tp_score < ti_score) ts = ti;//its image!!  //              ts = tp;//the track sorted by the tracker!!
571              else {  //              cs = cp;
572                  ts = tp;  //              os = op;
573                }else if (tp_score < ti_score) {
574                    ts = ti;//its image!!
575                    cs = ci;
576                    os = oi;
577    
578                    ti = tp;//its image!!
579                    ci = cp;
580                    oi = op;
581    
582                    tp = ts;//its image!!
583                    cp = cs;
584                    op = os;
585    
586                    
587                }else {
588    //              ts = tp;
589    //              cs = cp;
590    //              os = op;
591  //                              cout << "Warning - track image ambiguity not solved" << endl;  //                              cout << "Warning - track image ambiguity not solved" << endl;
592  //                              cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl;  //                              cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl;
593              };              };
594                            
595          }else{          }else{
596              ts = tp;  //          ts = tp;
597    //          cs = cp;
598    //          os = op;
599          };          };
600                    
601  //      cout <<" SortTracks() "<<i<<" -- "<<ts<<endl;  //      cout <<" SortTracks() "<<i<<" -- "<<ts<<endl;
602          sorted_tracks->Add(ts);//save the track in the sorted array  //      sorted_tracks->Add(ts);//save the track in the sorted array
603    //      sorted_tracks.Add(ts);//save the track in the sorted array
604    //      sorted_tracks.Add(tp);//save the track in the sorted array
605  //      cout << "SortTracks:: sorted_tracks->Add(it) "<<i<<" "<<ts<<endl;  //      cout << "SortTracks:: sorted_tracks->Add(it) "<<i<<" "<<ts<<endl;
606            //      cout<<"o "<<tp<<endl;
607    //      cout<<"o "<<cp<<endl;
608    //      cout<<"o "<<op<<endl;
609            new(ttsorted[i]) PamTrack(tp,cp,op);
610            new(ttimage[i])  PamTrack(ti,ci,oi);
611      };      };
612    
613        if( tsorted->GetEntries() != trk2_obj->GetNTracks() ){
614            cout << "void PamLevel2::SortTracks(TString how): tsorted->GetEntries() "<<tsorted->GetEntries()<<" != trk2_obj->GetNTracks() = "<<trk2_obj->GetNTracks() <<endl;
615            tsorted->Delete(); tsorted=0;
616            timage->Delete(); timage=0;    
617        }
618    
619      //Restore Object count      //Restore Object count
620      //To save space in the table keeping track of all referenced objects      //To save space in the table keeping track of all referenced objects
621      //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.  
622      TProcessID::SetObjectCount(ObjectNumber);      TProcessID::SetObjectCount(ObjectNumber);
623            
624  };  };
# Line 430  void PamLevel2::SortTracks(TString how){ Line 630  void PamLevel2::SortTracks(TString how){
630   * 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.
631   * 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.
632   */   */
633  TRefArray *PamLevel2::GetTracks(){  // TRefArray *PamLevel2::GetTracks(){
634    
635    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
636    //     SortTracks("+CAL+TOF");
637    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
638    
639    // //   return  sorted_tracks;
640    //     return &sorted_tracks;
641        
642    //  };
643    TClonesArray *PamLevel2::GetTracks(){
644    
645  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
646      SortTracks("+CAL+TOF");      SortTracks("+CAL+TOF");
647  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
648    
649      return sorted_tracks;      return tsorted;
650        
651   };   };
652  //--------------------------------------  //--------------------------------------
653   //   //
# Line 448  TRefArray *PamLevel2::GetTracks(){ Line 659  TRefArray *PamLevel2::GetTracks(){
659   * @param it Track number, ranging from 0 to GetNTracks().   * @param it Track number, ranging from 0 to GetNTracks().
660   */   */
661  PamTrack *PamLevel2::GetTrack(int it){  PamTrack *PamLevel2::GetTrack(int it){
662    
663    //    if(!trk2_obj) return 0;
664                    
665    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
666    //     SortTracks("+CAL+TOF");
667    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
668    // //    if(!sorted_tracks)return 0;
669    //     if(sorted_tracks.IsEmpty())return 0;
670    
671    //     PamTrack *track = 0;
672        
673    //     if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() && it<sorted_tracks.GetEntriesFast() ){    
674    //      TrkTrack   *t = (TrkTrack*)sorted_tracks.At(it);
675    //      track = GetPamTrackAlong(t);
676    //     }else{
677    //      cout << "PamLevel2::GetTrack(int) : tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl;
678    //     };
679        
680    //     return track;
681    
682    //    cout << "PamLevel2::GetTrack(int it) "<<endl;
683  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
684      SortTracks("+CAL+TOF");      SortTracks("+CAL+TOF");
685  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
686            if(!tsorted)return 0;
687        if(!tsorted->GetEntries())return 0;
688    
689      PamTrack *track = 0;      PamTrack *track = 0;
690            
691      if( it >=0 && it < TrkLevel2::GetNTracks() && it<sorted_tracks->GetEntriesFast() ){      if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks()){
692          TrkTrack   *t = (TrkTrack*)sorted_tracks->At(it);  //      TrkTrack   *t = (TrkTrack*)sorted_tracks.At(it);
693          track = GetPamTrackAlong(t);  //      track = GetPamTrackAlong(t);
694            TClonesArray &t = *(tsorted);
695            track = (PamTrack*)t[it];
696            
697      }else{      }else{
698          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;
699      };      };
700            
701      return track;      return track;
# Line 476  PamTrack *PamLevel2::GetTrack(int it){ Line 712  PamTrack *PamLevel2::GetTrack(int it){
712   */   */
713  PamTrack *PamLevel2::GetTrackImage(int it){  PamTrack *PamLevel2::GetTrackImage(int it){
714    
715    //     if(!trk2_obj) return 0;
716    
717    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
718    //     SortTracks("+CAL+TOF");
719    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
720    // //    if(!sorted_tracks)return 0;
721    //     if(sorted_tracks.IsEmpty())return 0;
722        
723    //     PamTrack *image = 0;
724        
725    //     if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){
726    //      TrkTrack *temp = (TrkTrack*)sorted_tracks.At(it);
727    //      if( temp->HasImage() ){
728    //          TrkTrack   *t = trk2_obj->TrkLevel2::GetStoredTrack(temp->GetImageSeqNo());
729    //          image = GetPamTrackAlong(t);
730    //      }else{
731    //          cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl;
732    //      };
733    //     }else{
734    //      cout << "PamLevel2::GetTrackImage(int) : Tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl;
735    //     };
736        
737    //     return image;
738    
739    
740  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
741      SortTracks("+CAL+TOF");      SortTracks("+CAL+TOF");
742  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
743        if(!timage)return 0;
744        if(!timage->GetEntries())return 0;
745            
746      PamTrack *image = 0;      PamTrack *image = 0;
747            
748      if( it >=0 && it < TrkLevel2::GetNTracks() ){      if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){
749          TrkTrack *temp = (TrkTrack*)sorted_tracks->At(it);          TClonesArray &t = *(tsorted);
750          if( temp->HasImage() ){          PamTrack *temp = (PamTrack*)t[it];
751              TrkTrack   *t = TrkLevel2::GetStoredTrack(temp->GetImageSeqNo());          if( temp->GetTrkTrack()->HasImage() ){
752              image = GetPamTrackAlong(t);              TClonesArray &t = *(timage);
753                image = (PamTrack*)t[it];
754          }else{          }else{
755              cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl;  //          cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl;
756          };          };
757      }else{      }else{
758          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;
759      };      };
760            
761      return image;      return image;
# Line 507  PamTrack *PamLevel2::GetTrackImage(int i Line 771  PamTrack *PamLevel2::GetTrackImage(int i
771   * @param detlist String to select trees to be included   * @param detlist String to select trees to be included
772   * @return Pointer to a TTree   * @return Pointer to a TTree
773   */   */
774  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";  
775    
776    //     if( !detlist.IsNull() ) SetWhichTrees(detlist);
777    //     else                    GetWhichTrees(f);
778      SetWhichTrees(detlist);      SetWhichTrees(detlist);
     TTree *Trout =0;  
779    
780        TTree *Trout =0;
781    
782      TString fname = f->GetName();      TString fname = f->GetName();
783      if(!CheckLevel2File(fname))return NULL;      if(!CheckLevel2File(fname))return NULL;
784    
785            UInt_t *found=0;    
786    
787        cout<< "GetPamTree(TFile*,TString): detector list --> ";
788        if(TRK1)cout<<"TRK1 ";
789        if(TRK2)cout<<"TRK2 ";
790        if(TRKh)cout<<"TRKH ";
791        if(CAL1)cout<<"CAL1 ";
792        if(CAL2)cout<<"CAL2 ";
793        if(TOF)cout<<"TOF ";
794        if(TRG)cout<<"TRG ";
795        if(AC)cout<<"AC ";
796        if(ND)cout<<"ND ";
797        if(S4)cout<<"S4 ";
798        if(ORB)cout<<"ORB ";
799        cout << endl;
800    
801        f->cd();
802    
803  // Tracker  // Tracker
804      TTree *T = (TTree*)f->Get("Tracker");      TTree *T = (TTree*)f->Get("Tracker");
805      if(T && (TRK2||TRK1||TRKh)) {      if(T && (TRK2||TRK1||TRKh)) {
806          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
807            else    T->SetBranchStatus("TrkLevel2",0,found);
808          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
809          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
810            else    T->SetBranchStatus("TrkLevel1",0,found);
811          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
812          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
813            else    T->SetBranchStatus("TrkHough",0,found);
814          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
815          if(!Trout)Trout=T;          if(!Trout)Trout=T;
816          else Trout->AddFriend(T);          else Trout->AddFriend(T);
# Line 538  TTree *PamLevel2::GetPamTree(TFile *f, T Line 821  TTree *PamLevel2::GetPamTree(TFile *f, T
821      TTree *C = (TTree*)f->Get("Calorimeter");      TTree *C = (TTree*)f->Get("Calorimeter");
822      if(C && (CAL2||CAL1)) {      if(C && (CAL2||CAL1)) {
823          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
824            else    C->SetBranchStatus("CaloLevel2",0,found);
825          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
826          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
827            else    C->SetBranchStatus("CaloLevel1",0,found);
828          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
829          if(!Trout)Trout=C;          if(!Trout)Trout=C;
830          else Trout->AddFriend(C);          else Trout->AddFriend(C);
831      }else{      }else{
832          cout << "Calorimeter  : missing tree"<<endl;          cout << "Calorimeter  : missing tree"<<endl;
833      };      };
834    
835      // ToF          // ToF    
836      TTree *O = (TTree*)f->Get("ToF");      TTree *O = (TTree*)f->Get("ToF");
837      if(O && TOF) {      if(O && TOF) {
# Line 597  TTree *PamLevel2::GetPamTree(TFile *f, T Line 883  TTree *PamLevel2::GetPamTree(TFile *f, T
883          cout << "Anticounter  : missing tree"<<endl;          cout << "Anticounter  : missing tree"<<endl;
884      };      };
885      // Orbital Info      // Orbital Info
886          TTree *B = (TTree*)f->Get("OrbitalInfo");      TTree *B = (TTree*)f->Get("OrbitalInfo");
887          if(B && ORB) {      if(B && ORB) {
888                  B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));          B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
889                  cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
890                  if(!Trout)Trout=O;          if(!Trout)Trout=O;
891                  else Trout->AddFriend(B);          else Trout->AddFriend(B);
892          }else{      }else{
893                  cout << "OrbitalInfo  : missing tree"<<endl;          cout << "OrbitalInfo  : missing tree"<<endl;
894          };      };
895            
896          cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;      cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
897                        
898          return Trout;      return Trout;
899            
900  }  }
901  //--------------------------------------  //--------------------------------------
# Line 634  TList*  PamLevel2::GetListOfLevel2Files( Line 920  TList*  PamLevel2::GetListOfLevel2Files(
920      TList *contents  = new TList; // create output list      TList *contents  = new TList; // create output list
921      contents->SetOwner();      contents->SetOwner();
922            
923      char *fullpath;  //    char *fullpath;
924    //    const char *fullpath;
925            
926      // if no input file list is given:        // if no input file list is given:  
927      if ( flisttxt != "" ){      if ( flisttxt != "" ){
928                    
929          if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) )return 0;  //      if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){        
930            //      if( !(fullpath = gSystem->FindFile("./",flisttxt)) ){
931  //              flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));  //          cout <<"File "<<flisttxt<<" not found"<<endl;
932          flisttxt = fullpath;  //          return 0;
933            //      }      
934          if( !gSystem->ChangeDirectory(ddir) )return 0;  //      flisttxt = fullpath;
935    
936            flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));
937    
938            if( !gSystem->ChangeDirectory(ddir) ){
939                cout << "Cannot change directory : "<<ddir<<endl;
940                return 0;
941            }
942                    
943          cout <<"Input file list : " << flisttxt <<endl;          cout <<"Input file list : " << flisttxt <<endl;
944          ifstream in;          ifstream in;
945          in.open(flisttxt, ios::in); //open input file list          in.open(flisttxt, ios::in); //open input file list
946            if(!in.good()){
947                cout <<" ERROR opening the file "<<endl;
948                gSystem->ChangeDirectory(wdir); // back to the working directory
949                return 0;
950            }      
951          int line=0;          int line=0;
952          while (1) {          while (1) {
953              TString file;              TString file;
# Line 662  TList*  PamLevel2::GetListOfLevel2Files( Line 961  TList*  PamLevel2::GetListOfLevel2Files(
961              }              }
962              if(file.Contains("#"))file = file(0,file.First("#"));              if(file.Contains("#"))file = file(0,file.First("#"));
963  //          cout <<"(2) " << file << endl;  //          cout <<"(2) " << file << endl;
964              if( gSystem->IsFileInIncludePath(file,&fullpath) ){  //          if( gSystem->IsFileInIncludePath(file,&fullpath) ){
965    //          if( (fullpath = gSystem->FindFile(ddir,file)) ){
966                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list              char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(file));
967                contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list
968              }else{              delete fullpath;
969    //          }else{
970  //              if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;  //              if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;
971              };  //          };
972          };          };
973          in.close();          in.close();
974                    
# Line 690  TList*  PamLevel2::GetListOfLevel2Files( Line 990  TList*  PamLevel2::GetListOfLevel2Files(
990          while ( (questo = (TSystemFile*) next()) ) {          while ( (questo = (TSystemFile*) next()) ) {
991              TString name =  questo-> GetName();              TString name =  questo-> GetName();
992              if( name.EndsWith(".root") ){              if( name.EndsWith(".root") ){
993                  char *fullpath;  //              const char *fullpath = gSystem->FindFile(ddir,name);
994                  gSystem->IsFileInIncludePath(name,&fullpath);  //              char *fullpath;
995    //              gSystem->IsFileInIncludePath(name,&fullpath);
996                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(name));
997                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));
998                    delete fullpath;
999              };              };
1000          }          }
1001          delete temp;          delete temp;
# Line 702  TList*  PamLevel2::GetListOfLevel2Files( Line 1005  TList*  PamLevel2::GetListOfLevel2Files(
1005      gSystem->ChangeDirectory(wdir); // back to the working directory      gSystem->ChangeDirectory(wdir); // back to the working directory
1006  //      cout << endl << "Selected files:" << endl;  //      cout << endl << "Selected files:" << endl;
1007  //      contents->Print();  //      contents->Print();
1008      cout << contents->GetEntries()<<" files selected\n";      cout << contents->GetEntries()<<" files \n";
1009  //      cout << endl;  //      cout << endl;
1010  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;
1011      return contents;      return contents;
# Line 717  TList*  PamLevel2::GetListOfLevel2Files( Line 1020  TList*  PamLevel2::GetListOfLevel2Files(
1020   * @param detlist String to select trees to be included   * @param detlist String to select trees to be included
1021   * @return Pointer to a TChain   * @return Pointer to a TChain
1022   */   */
1023  TChain *PamLevel2::GetPamTree(TList *fl, TString detlist="+ALL"){  TChain *PamLevel2::GetPamTree(TList *fl, TString detlist ){
   
 //      TChain *Tout=0;  
 //    if(Tout)Tout->Delete();  
 //    Tout = NULL;  
1024            
1025      TChain *Tout =0;      TChain *Trout =0;
1026    
1027    //    if( !detlist.IsNull() )SetWhichTrees(detlist);
1028      SetWhichTrees(detlist);      SetWhichTrees(detlist);
1029            
1030        cout<< "GetPamTree(TList*,TString): input detector list --> ";
1031        if(TRK1)cout<<"TRK1 ";
1032        if(TRK2)cout<<"TRK2 ";
1033        if(TRKh)cout<<"TRKH ";
1034        if(CAL1)cout<<"CAL1 ";
1035        if(CAL2)cout<<"CAL2 ";
1036        if(TOF)cout<<"TOF ";
1037        if(TRG)cout<<"TRG ";
1038        if(AC)cout<<"AC ";
1039        if(ND)cout<<"ND ";
1040        if(S4)cout<<"S4 ";
1041        if(ORB)cout<<"ORB ";
1042        cout << endl;
1043    
1044      TChain *T = 0;            TChain *T = 0;      
1045      TChain *C = 0;      TChain *C = 0;
1046      TChain *O = 0;      TChain *O = 0;
# Line 737  TChain *PamLevel2::GetPamTree(TList *fl, Line 1051  TChain *PamLevel2::GetPamTree(TList *fl,
1051      TChain *B = 0;      TChain *B = 0;
1052            
1053      if(TRK2||TRK1||TRKh) T = new TChain("Tracker");          if(TRK2||TRK1||TRKh) T = new TChain("Tracker");    
1054      if(CAL2||CAL1) C = new TChain("Calorimeter");      if(CAL2||CAL1)       C = new TChain("Calorimeter");
1055      if(TOF) O = new TChain("ToF");      if(TOF)              O = new TChain("ToF");
1056      if(TRG) R = new TChain("Trigger");      if(TRG)              R = new TChain("Trigger");
1057      if(S4)  S = new TChain("S4");      if(S4)               S = new TChain("S4");
1058      if(ND)  N = new TChain("NeutronD");      if(ND)               N = new TChain("NeutronD");
1059      if(AC)  A = new TChain("Anticounter");      if(AC)               A = new TChain("Anticounter");
1060      if(ORB) B = new TChain("OrbitalInfo");      if(ORB)              B = new TChain("OrbitalInfo");
1061            
1062      // loop over files and create chains              // loop over files and create chains        
1063      TIter next(fl);      TIter next(fl);
1064      TSystemFile *questo = 0;      TSystemFile *questo = 0;
1065      while ( (questo = (TSystemFile*) next()) ) {      while ( (questo = (TSystemFile*) next()) ) {
1066          TString name =  questo->GetName();          TString name =  questo->GetName();
1067  //              cout << "File: "<< name << endl;          cout << "File: "<< name << endl;
1068          if( CheckLevel2File(name) ){          if( CheckLevel2File(name) ){
1069              if(TRK2||TRK1||TRKh) T->Add(name);              if(TRK2||TRK1||TRKh) T->Add(name);
1070              if(CAL1||CAL2) C->Add(name);              if(CAL1||CAL2)       C->Add(name);
1071              if(TOF) O->Add(name);              if(TOF)              O->Add(name);
1072              if(TRG) R->Add(name);              if(TRG)              R->Add(name);
1073              if(S4)  S->Add(name);              if(S4)               S->Add(name);
1074              if(ND)  N->Add(name);              if(ND)               N->Add(name);
1075              if(AC)  A->Add(name);              if(AC)               A->Add(name);
1076              if(ORB) B->Add(name);              if(ORB)              B->Add(name);
1077          };          };
1078      }      }
1079            
1080      cout << "done chain \n";      cout << "done chain \n";
1081        
1082        UInt_t *found=0;
1083    // Tracker
1084        if(T && (TRK2||TRK1||TRKh)) {
1085            if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1086            else    T->SetBranchStatus("TrkLevel2",0,found);
1087            if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
1088            if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
1089            else    T->SetBranchStatus("TrkLevel1",0,found);
1090            if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
1091            if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
1092            else    T->SetBranchStatus("TrkHough",0,found);
1093            if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
1094            if(!Trout)Trout=T;
1095            else Trout->AddFriend("Tracker");
1096        }else{
1097            cout << "Tracker      : missing tree"<<endl;
1098        };
1099      // Calorimeter      // Calorimeter
1100      if(CAL2||CAL1) {                  if(C && (CAL2||CAL1)) {
1101          if(!Tout)Tout=C;          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
1102          else Tout->AddFriend("Calorimeter");          else    C->SetBranchStatus("CaloLevel2",0,found);
1103            if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
1104            if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
1105            else    C->SetBranchStatus("CaloLevel1",0,found);
1106            if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
1107            if(!Trout)Trout=C;
1108            else Trout->AddFriend("Calorimeter");
1109        }else{
1110            cout << "Calorimeter  : missing tree"<<endl;
1111      };      };
1112      // ToF          // ToF    
1113      if(TOF) {      if(O && TOF) {
1114          if(!Tout)Tout=O;          O->SetBranchAddress("ToFLevel2", GetPointerTo("ToFLevel2"));
1115          else Tout->AddFriend("ToF");          cout << "ToF          : set branch address ToFLevel2"<<endl;
1116            if(!Trout)Trout=O;
1117            else Trout->AddFriend("ToF");
1118        }else{
1119            cout << "ToF         : missing tree"<<endl;
1120      };      };
     // Tracker  
     if(TRK2||TRK1||TRKh){  
         if(!Tout)Tout=T;  
         else Tout->AddFriend("Tracker");  
     }  
1121      // Trigger      // Trigger
1122      if(TRG) {      if(R && TRG) {
1123          if(!Tout)Tout=R;          R->SetBranchAddress("TrigLevel2", GetPointerTo("TrigLevel2"));
1124          else Tout->AddFriend("Trigger");          cout << "Trigger      : set branch address TrigLevel2"<<endl;
1125            if(!Trout)Trout=O;
1126            else Trout->AddFriend("Trigger");
1127        }else{
1128            cout << "Trigger      : missing tree"<<endl;
1129      };      };
1130      // S4      // S4
1131      if(S4) {      if(S && S4) {
1132          if(!Tout)Tout=S;          S->SetBranchAddress("S4Level2", GetPointerTo("S4Level2"));
1133          else Tout->AddFriend("S4");          cout << "S4           : set branch address S4Level2"<<endl;
1134            if(!Trout)Trout=O;
1135            else Trout->AddFriend("S4");
1136        }else{
1137            cout << "S4           : missing tree"<<endl;
1138      };      };
1139      // Neutron Detector      // Neutron Detector
1140      if(ND) {      if(N && ND) {
1141          if(!Tout)Tout=N;          N->SetBranchAddress("NDLevel2", GetPointerTo("NDLevel2"));
1142          else Tout->AddFriend("NeutronD");          cout << "NeutronD     : set branch address NDLevel2"<<endl;
1143            if(!Trout)Trout=O;
1144            else Trout->AddFriend("NeutronD");
1145        }else{
1146            cout << "NeutronD     : missing tree"<<endl;
1147      };      };
1148      // Anticounters      // Anticounters
1149      if(AC) {      if(A && AC) {
1150          if(!Tout)Tout=A;          A->SetBranchAddress("AcLevel2", GetPointerTo("AcLevel2"));
1151          else Tout->AddFriend("Anticounter");          cout << "Anticounter  : set branch address AcLevel2"<<endl;
1152            if(!Trout)Trout=O;
1153            else Trout->AddFriend("Anticounter");
1154        }else{
1155            cout << "Anticounter  : missing tree"<<endl;
1156      };      };
1157      // OrbitalInfo      // Orbital Info
1158      if(ORB) {      if(B && ORB) {
1159          if(!Tout)Tout=B;          B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1160          else Tout->AddFriend("OrbitalInfo");          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1161            if(!Trout)Trout=O;
1162            else Trout->AddFriend("OrbitalInfo");
1163        }else{
1164            cout << "OrbitalInfo  : missing tree"<<endl;
1165      };      };
1166            
1167    //    cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
1168    
1169        
1170    //     // Tracker
1171    //     if(TRK2||TRK1||TRKh){
1172    //      if(!Tout)Tout=T;
1173    //      else Tout->AddFriend("Tracker");
1174    //     }
1175    //     // Calorimeter
1176    //     if(CAL2||CAL1) {        
1177    //      if(!Tout)Tout=C;
1178    //      else Tout->AddFriend("Calorimeter");
1179    //     };
1180    //     // ToF    
1181    //     if(TOF) {
1182    //      if(!Tout)Tout=O;
1183    //      else Tout->AddFriend("ToF");
1184    //     };
1185    //     // Trigger
1186    //     if(TRG) {
1187    //      if(!Tout)Tout=R;
1188    //      else Tout->AddFriend("Trigger");
1189    //     };
1190    //     // S4
1191    //     if(S4) {
1192    //      if(!Tout)Tout=S;
1193    //      else Tout->AddFriend("S4");
1194    //     };
1195    //     // Neutron Detector
1196    //     if(ND) {
1197    //      if(!Tout)Tout=N;
1198    //      else Tout->AddFriend("NeutronD");
1199    //     };
1200    //     // Anticounters
1201    //     if(AC) {
1202    //      if(!Tout)Tout=A;
1203    //      else Tout->AddFriend("Anticounter");
1204    //     };
1205    //     // OrbitalInfo
1206    //     if(ORB) {
1207    //      if(!Tout)Tout=B;
1208    //      else Tout->AddFriend("OrbitalInfo");
1209    //     };
1210        
1211  //    cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl;  //    cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl;
1212            
1213  //    if( Tout->GetEntries() )PamLevel2::SetBranchAddress();  //    if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout );
     if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout );  
1214            
1215      return Tout;      return Trout;
1216  }  }
1217    
1218    
# Line 913  void PamLevel2::SetBranchAddress(TChain Line 1314  void PamLevel2::SetBranchAddress(TChain
1314    
1315      // Tracker      // Tracker
1316       if(TRK2) {       if(TRK2) {
1317          t->SetBranchAddress("TrkLevel2", this->GetPointerTo("TrkLevel2"));          t->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1318          cout << "Tracker      : set branch address TrkLevel2"<<endl;          cout << "Tracker      : set branch address TrkLevel2"<<endl;
1319      };      };
1320      if(TRK1) {      if(TRK1) {
# Line 968  void PamLevel2::SetBranchAddress(TChain Line 1369  void PamLevel2::SetBranchAddress(TChain
1369            
1370  }  }
1371    
 // 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;  };  
1372    
1373  //--------------------------------------  //--------------------------------------
1374  //  //
# Line 1019  TChain *PamLevel2::GetRunTree(TList *fl) Line 1396  TChain *PamLevel2::GetRunTree(TList *fl)
1396            
1397      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1398      cout << "Run         : set branch address RunInfo"<<endl;      cout << "Run         : set branch address RunInfo"<<endl;
1399        R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1400        cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1401            
1402      return R;      return R;
1403            
# Line 1039  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1418  TTree *PamLevel2::GetRunTree(TFile *f){
1418            
1419      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1420      cout << "Run         : set branch address RunInfo"<<endl;      cout << "Run         : set branch address RunInfo"<<endl;
1421        R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1422        cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1423            
1424      return R;      return R;
1425            
1426  }  }
1427    /**
1428     * Update the runinfo informations (to be used to have Run infos event by event basis)
1429     * @param run Pointer to the chain/tree which contains run infos
1430     * @return true if a new run has been read, false if it is still the same run
1431     */
1432    Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev){
1433      //
1434      // check if we have already called once GetEntry, if not call it
1435      //
1436      if ( run->GetEntries() <= 0 ) return(false);
1437      //
1438      Int_t oldrun = irun;
1439      if ( irun < 0 ){
1440        irun = 0;
1441        run->GetEntry(0);
1442        runfirstentry = 0ULL;
1443        runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS) - 1ULL;
1444      };
1445      while ( iev > (runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS-1)) && irun < run->GetEntries() ){
1446        //    printf(" iev %llu %u %llu \n",iev,this->GetRunInfo()->NEVENTS,(ULong64_t)(runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS)));
1447        irun++;
1448        run->GetEntry(irun);
1449        runfirstentry = runlastentry+1ULL;
1450        runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS);
1451      };
1452      //
1453      if ( irun == oldrun || irun >= run->GetEntries() ) return(false);
1454      //
1455      //  printf(" iev %llu irun %i nevents %u 1st %llu last %llu \n",iev,irun,this->GetRunInfo()->NEVENTS,(ULong64_t)runfirstentry,(ULong64_t)runlastentry);
1456      //
1457      return(true);    
1458      //
1459    };
1460    /**
1461     * Update the runinfo informations (to be used to have Run infos event by event basis)
1462     * @param run Pointer to the chain/tree which contains run infos
1463     * @return true if a new run has been read, false if it is still the same run
1464     */
1465    Bool_t PamLevel2::UpdateRunInfo(TTree *run, ULong64_t iev){
1466      return(UpdateRunInfo((TChain*)run,iev));
1467    };
1468    
1469  //--------------------------------------  //--------------------------------------
1470  //  //
1471  //  //
# Line 1053  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1476  TTree *PamLevel2::GetRunTree(TFile *f){
1476  */  */
1477  void PamLevel2::SetWhichTrees(TString detlist){  void PamLevel2::SetWhichTrees(TString detlist){
1478                    
1479      if(detlist.Contains("+ALL", TString::kIgnoreCase)){      if(detlist.IsNull() || detlist.Contains("+ALL", TString::kIgnoreCase)){
1480          CAL1 = true;          CAL1 = true;
1481          CAL2 = true;          CAL2 = true;
1482          TRK2 = true;          TRK2 = true;
# Line 1118  void PamLevel2::SetWhichTrees(TString de Line 1541  void PamLevel2::SetWhichTrees(TString de
1541          TRKh=false;          TRKh=false;
1542      }      }
1543  //  -------------------------------------------------------------------------  //  -------------------------------------------------------------------------
   
   
       
 //     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;  
1544            
1545      if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false;      if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false;
1546      else if( detlist.Contains("+TRG", TString::kIgnoreCase) )TRG = true;      else if( detlist.Contains("+TRG", TString::kIgnoreCase) )TRG = true;
# Line 1147  void PamLevel2::SetWhichTrees(TString de Line 1559  void PamLevel2::SetWhichTrees(TString de
1559            
1560      if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false;      if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false;
1561      else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true;      else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true;
1562    
1563    //     cout<< "Set detector list --> ";
1564    //     if(TRK1)cout<<"TRK1 ";
1565    //     if(TRK2)cout<<"TRK2 ";
1566    //     if(TRKh)cout<<"TRKH ";
1567    //     if(CAL1)cout<<"CAL1 ";
1568    //     if(CAL2)cout<<"CAL2 ";
1569    //     if(TOF)cout<<"TOF ";
1570    //     if(TRG)cout<<"TRG ";
1571    //     if(AC)cout<<"AC ";
1572    //     if(ND)cout<<"ND ";
1573    //     if(S4)cout<<"S4 ";
1574    //     if(ORB)cout<<"ORB ";
1575    //     cout << endl;
1576        
1577    };
1578    
1579    
1580    /**
1581     * Set tree/branch detector flags from the content of a tree
1582     */
1583    void  PamLevel2::GetWhichTrees(TFile* f){
1584            
1585        // -----------
1586        // reset flags
1587        // -----------
1588        CAL1   = false;    
1589        CAL2   = false;    
1590        TRK2   = false;    
1591        TRK1   = false;    
1592        TRKh   = false;    
1593        TRG    = false;    
1594        TOF    = false;    
1595        S4     = false;    
1596        ND     = false;    
1597        AC     = false;    
1598        ORB    = false;    
1599        
1600        RUN    = false;
1601            
1602        cout << "Checking file: "<<f->GetName()<<endl;
1603        if( !f || f->IsZombie() ){
1604            cout << "File: "<< f->GetName() <<" Non valid root file"<< endl;
1605            return;
1606        }
1607    
1608        TList *lk = f->GetListOfKeys();
1609        TIter next(lk);
1610        TKey *key =0;
1611    
1612        Int_t nev = 0;
1613    
1614        while( (key = (TKey*)next()) ){
1615            
1616            if( !strcmp(key->GetName(),"Run"        ) )RUN = true;
1617    
1618            //=========================================================    
1619            if( !strcmp(key->GetName(),"Trigger"    ) ){
1620                TRG = true;
1621                Int_t nevt = ((TTree*)f->Get("Trigger"))->GetEntries();
1622                if( nev && nevt!=nev){
1623                    cout << "File: "<< f->GetName() <<" Trigger tree has "<<nevt<<" events instead of "<<nev<< endl;
1624                    TRG = false;
1625                }else nev=nevt;
1626            }
1627            //=========================================================    
1628            if( !strcmp(key->GetName(),"ToF"        ) ){
1629                TOF = true;
1630                Int_t nevt = ((TTree*)f->Get("ToF"))->GetEntries();
1631                if( nev && nevt!=nev){
1632                    cout << "File: "<< f->GetName() <<"     ToF tree has "<<nevt<<" events instead of "<<nev<< endl;
1633                    TOF = false;
1634                }else nev=nevt;
1635            }
1636            //=========================================================  
1637            if( !strcmp(key->GetName(),"S4"         ) ){
1638                S4 = true;
1639                Int_t nevt = ((TTree*)f->Get("S4"))->GetEntries();
1640                if( nev && nevt!=nev){
1641                    cout << "File: "<< f->GetName() <<"      S4 tree has "<<nevt<<" events instead of "<<nev<< endl;
1642                    S4 = false;
1643                }else nev=nevt;
1644            }
1645            //=========================================================  
1646    
1647            if( !strcmp(key->GetName(),"NeutronD"   ) ){
1648                ND = true;
1649                Int_t nevt = ((TTree*)f->Get("NeutronD"))->GetEntries();
1650                if( nev && nevt!=nev){
1651                    cout << "File: "<< f->GetName() <<"NeutronD tree has "<<nevt<<" events instead of "<<nev<< endl;
1652                    ND =false;
1653                }else nev=nevt;
1654            }      
1655            //=========================================================  
1656            if( !strcmp(key->GetName(),"Anticounter") ){
1657                AC = true;
1658                Int_t nevt = ((TTree*)f->Get("Anticounter"))->GetEntries();
1659                if( nev && nevt!=nev){
1660                    cout << "File: "<< f->GetName() <<" Anticounter tree has "<<nevt<<" events instead of "<<nev<< endl;
1661                    AC =false;
1662                }else nev=nevt;
1663            }
1664            //=========================================================  
1665            if( !strcmp(key->GetName(),"OrbitalInfo") ){
1666                ORB = true;
1667                Int_t nevt = ((TTree*)f->Get("OrbitalInfo"))->GetEntries();
1668                if( nev && nevt!=nev){
1669                    cout << "File: "<< f->GetName() <<" OrbitalInfo tree has "<<nevt<<" events instead of "<<nev<< endl;
1670                    ORB = false;
1671                }else nev=nevt;
1672            }
1673            //=========================================================  
1674            if( !strcmp(key->GetName(),"Tracker"    ) ){
1675                TTree *T = (TTree*)f->Get("Tracker");
1676                for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){
1677                    TString name = T->GetListOfBranches()->At(i)->GetName();
1678                    if( !name.CompareTo("TrkLevel1") )TRK1=true;
1679                    if( !name.CompareTo("TrkLevel2") )TRK2=true;
1680                    if( !name.CompareTo("TrkHough") )TRKh=true;
1681                };      
1682                Int_t nevt = T->GetEntries();
1683                if( nev && nevt!=nev){
1684                    cout << "File: "<< f->GetName() <<" Tracker tree has "<<nevt<<" events instead of "<<nev<< endl;
1685                    TRK1 = false;
1686                    TRK2 = false;
1687                    TRKh = false;
1688                }else nev=nevt;
1689                T->Delete();
1690            };
1691            //=========================================================  
1692            if( !strcmp(key->GetName(),"Calorimeter"    ) ){
1693                TTree *T = (TTree*)f->Get("Calorimeter");
1694                for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){
1695                    TString name = T->GetListOfBranches()->At(i)->GetName();
1696                    if( !name.CompareTo("CaloLevel1") )CAL1=true;
1697                    if( !name.CompareTo("CaloLevel2") )CAL2=true;
1698                };    
1699                Int_t nevt = T->GetEntries();
1700                if( nev && nevt!=nev){
1701                    cout << "File: "<< f->GetName() <<"  Calorimeter tree has "<<nevt<<" events instead of "<<nev<< endl;
1702                    CAL1 = false;
1703                    CAL2 = false;
1704                }else nev=nevt;
1705                T->Delete();
1706            };      
1707    
1708        };
1709            
1710        delete lk;
1711    
1712    //     cout<< "Get detector list from input file --> ";
1713    //     if(TRK1)cout<<"TRK1 ";
1714    //     if(TRK2)cout<<"TRK2 ";
1715    //     if(TRKh)cout<<"TRKH ";
1716    //     if(CAL1)cout<<"CAL1 ";
1717    //     if(CAL2)cout<<"CAL2 ";
1718    //     if(TOF)cout<<"TOF ";
1719    //     if(TRG)cout<<"TRG ";
1720    //     if(AC)cout<<"AC ";
1721    //     if(ND)cout<<"ND ";
1722    //     if(S4)cout<<"S4 ";
1723    //     if(ORB)cout<<"ORB ";
1724    //     cout << endl;
1725          
1726        return ;
1727            
1728  };  };
1729    
1730    
1731  //--------------------------------------  //--------------------------------------
1732  //  //
1733  //  //
# Line 1191  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1770  Bool_t  PamLevel2::CheckLevel2File(TStri
1770    
1771      while( (key = (TKey*)next()) ){      while( (key = (TKey*)next()) ){
1772                    
1773  //      cout << key->GetName() << endl;  //      cout << key->GetName() << ""<<key->GetClassName()<<endl;
1774  //      cout << " Get tree: " << f->Get(key->GetName())<<endl;  //              cout << " Get tree: " << f->Get(key->GetName())<<endl;
1775  //      nev_previous = nev;  //      nev_previous = nev;
1776  //      cout << " n.entries  "<< nev <<endl;  //      cout << " n.entries  "<< nev <<endl;
1777  //      if( key->GetClassName()=="TTree" && nev_previous && nev != nev_previous ){  //      if( key->GetClassName()=="TTree" && nev_previous && nev != nev_previous ){
# Line 1316  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1895  Bool_t  PamLevel2::CheckLevel2File(TStri
1895    
1896      };      };
1897            
1898    //     cout<< "CheckLevel2File(TString): detector list --> ";
1899    //     if(TRK1__ok)cout<<"TRK1 ";
1900    //     if(TRK2__ok)cout<<"TRK2 ";
1901    //     if(TRKh__ok)cout<<"TRKH ";
1902    //     if(CAL1__ok)cout<<"CAL1 ";
1903    //     if(CAL2__ok)cout<<"CAL2 ";
1904    //     if(TOF__ok)cout<<"TOF ";
1905    //     if(TRG__ok)cout<<"TRG ";
1906    //     if(AC__ok)cout<<"AC ";
1907    //     if(ND__ok)cout<<"ND ";
1908    //     if(S4__ok)cout<<"S4 ";
1909    //     if(ORB__ok)cout<<"ORB ";
1910    //     cout << endl;
1911    
1912    
1913        if(TRK2 && TRK1__ok)TRK1=1;
1914    // ----------------------------------------------------------------------------
1915    // NOTA
1916    // se c'e` il level1, lo devo necessarimente leggere.
1917    // infatti (non ho capito perche`) i cluster vengono letti e allocati in memoria
1918    // comunque, ma non vengono disallocati da PamLevel2::Clear()
1919    // ----------------------------------------------------------------------------
1920    
1921    
1922      if(!RUN__ok) {      if(!RUN__ok) {
# Line 1360  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1961  Bool_t  PamLevel2::CheckLevel2File(TStri
1961          return false;          return false;
1962      };      };
1963    
1964          if(ND && !ND__ok){      if(ND && !ND__ok){
1965              cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl;          cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl;
1966              return false;          return false;
1967          };      };
1968          if(TRG && !TRG__ok){      if(TRG && !TRG__ok){
1969              cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl;          cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl;
1970              return false;          return false;
1971          };      };
1972        
   
   
   
1973    
1974      lk->Delete();  //    lk->Delete();
1975    //    delete lk;
1976      f->Close();      f->Close();
1977    
1978    //     cout<< "CheckLevel2File(TString): detector list --> ";
1979    //     if(TRK1)cout<<"TRK1 ";
1980    //     if(TRK2)cout<<"TRK2 ";
1981    //     if(TRKh)cout<<"TRKH ";
1982    //     if(CAL1)cout<<"CAL1 ";
1983    //     if(CAL2)cout<<"CAL2 ";
1984    //     if(TOF)cout<<"TOF ";
1985    //     if(TRG)cout<<"TRG ";
1986    //     if(AC)cout<<"AC ";
1987    //     if(ND)cout<<"ND ";
1988    //     if(S4)cout<<"S4 ";
1989    //     if(ORB)cout<<"ORB ";
1990    //     cout << endl;
1991                
1992      return true;      return true;
1993                    
# Line 1384  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1997  Bool_t  PamLevel2::CheckLevel2File(TStri
1997  /**  /**
1998   * Create clone-trees   * Create clone-trees
1999   */   */
2000  void PamLevel2::CreateCloneTrees( TChain *fChain ){  void PamLevel2::CreateCloneTrees0( TChain *fChain, TFile *ofile ){
2001    
2002      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2003      cout << "Create clones of PAMELA trees "<<endl;      cout << "Create clones of PAMELA trees "<<endl;
# Line 1400  void PamLevel2::CreateCloneTrees( TChain Line 2013  void PamLevel2::CreateCloneTrees( TChain
2013      TList *li = fChain->GetListOfFriends();      TList *li = fChain->GetListOfFriends();
2014      TIter next(li);      TIter next(li);
2015      TFriendElement* T_friend=0;      TFriendElement* T_friend=0;
2016        ofile->cd();
2017      while( (T_friend = (TFriendElement*)next()) ){      while( (T_friend = (TFriendElement*)next()) ){
2018  //      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;
2019  //      cout<<T_friend->GetTree()->GetName()<< endl;  //      cout<<T_friend->GetTree()->GetName()<< endl;
2020          tree_clone[i] = T_friend->GetTree()->CloneTree(0);                tree_clone[i] = T_friend->GetTree()->CloneTree(0);
2021            tree_clone[i]->SetAutoSave(1000000);
2022          name =  tree_clone[i]->GetName();          name =  tree_clone[i]->GetName();
2023          name.Append("_clone");          name.Append("_clone");
2024  //      tree_clone[i]->SetName(name.Data());  //      tree_clone[i]->SetName(name.Data());
2025          cout << tree_clone[i]->GetName() << endl;          cout << tree_clone[i]->GetName() << endl;
2026          i++;          i++;
2027      }      }
2028        
2029        delete li;
2030    
2031      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2032    
2033  }  }
2034    
2035  // /**  /**
2036  //  * Create clone-trees   * Create clone-trees
2037  //  */   */
2038  // 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;  
 //      };  
 //     }      
2039    
2040  //     // 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;  
2041    
2042  // }      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2043        cout << "Create new PAMELA trees "<<endl;
2044            
2045        Int_t i=0;
2046    
2047  // /**      if(TRK1||TRK2||TRKh){
2048  //  * Create clone-trees          tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data ");    
2049  //  */          if(TRK1) {
2050  // void PamLevel2::CreateCloneTrees(){              tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1"));
2051                tree_clone[i]->BranchRef();
2052                cout << "Tracker      : branch TrkLevel1"<<endl;
2053            };
2054            if(TRK2) {
2055                tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2"));
2056                cout << "Tracker      : branch TrkLevel2"<<endl;        
2057            };
2058            if(TRKh) {
2059                tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough"));
2060                cout << "Tracker      : branch TrkHough"<<endl;
2061            };
2062            i++;
2063        }
2064    
2065  //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      // Calorimeter
2066  //     cout << "Create new PAMELA trees "<<endl;      if(CAL1||CAL2){
2067                    tree_clone[i] = new TTree("Calorimeter","PAMELA calorimeter level2 data ");    
2068  //     if(TRK1||TRK2||TRKh){          if(CAL1) {
2069  //      T_clone = new TTree("Tracker_clone","PAMELA tracker level2 data ");                  tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1"));
2070  //      if(TRK1) {              cout << "Calorimeter  : branch CaloLevel1"<<endl;
2071  //          trk1_clone = new TrkLevel1();          };
2072  //          T_clone->Branch("TrkLevel1","TrkLevel1", &trk1_clone);          if(CAL2) {
2073  //          T_clone->BranchRef();              tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2"));
2074  //          cout << "Tracker      : branch TrkLevel1"<<endl;              cout << "Calorimeter  : branch CaloLevel2"<<endl;
2075  //      };          };
2076  //      if(TRK2) {          i++;
2077  //          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;  
 //     }  
2078    
2079  //     // Calorimeter      // ToF    
2080  //     if(CAL1||CAL2){      if(TOF) {
2081  //      C_clone = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data ");              tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data ");    
2082  //      if(CAL1) {          tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2"));
2083  //          calo1_clone = new CaloLevel1();          cout << "ToF          : branch ToFLevel2"<<endl;
2084  //          C_clone->Branch("CaloLevel1", "CaloLevel1", &calo1_clone);          i++;
2085  //          cout << "Calorimeter  : branch CaloLevel1"<<endl;      };
2086  //      };      // Trigger
2087  //      if(CAL2) {      if(TRG) {
2088  //          calo2_clone = new CaloLevel2();          tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data ");    
2089  //          C_clone->Branch("CaloLevel2","CaloLevel2", &calo2_clone);          tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2"));
2090  //          cout << "Calorimeter  : branch CaloLevel2"<<endl;          cout << "Trigger      : branch TrigLevel2"<<endl;
2091  //      };          i++;
2092  //     }          };
2093        // S4
2094        if(S4) {
2095            tree_clone[i] = new TTree("S4","PAMELA S4 level2 data ");      
2096            tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2"));
2097            cout << "S4           : branch S4Level2"<<endl;
2098            i++;
2099        };
2100        // Neutron Detector
2101        if(ND) {
2102            tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data ");  
2103            tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2"));
2104            cout << "NeutronD     : branch NDLevel2"<<endl;
2105            i++;
2106        };
2107        // Anticounters
2108        if(AC) {
2109            tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data ");    
2110            tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2"));
2111            cout << "Anticounter  : branch AcLevel2"<<endl;
2112            i++;
2113        };
2114        // OrbitalInfo
2115        if(ORB) {
2116            tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info  ");      
2117            tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo"));
2118            cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;
2119            i++;
2120        };
2121        cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2122    
2123  //     // 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;  
2124    
 // }  
2125    
2126  /**  /**
2127   * Fill tree (created with CreateCloneTrees)   * Fill tree (created with CreateCloneTrees)
# Line 1602  void PamLevel2::CreateCloneTrees( TChain Line 2130  void PamLevel2::CreateCloneTrees( TChain
2130  //void PamLevel2::FillNewPamTree(TTree *T){  //void PamLevel2::FillNewPamTree(TTree *T){
2131  void PamLevel2::FillCloneTrees(){  void PamLevel2::FillCloneTrees(){
2132                    
 // //  --------------------------------------  
 // //  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();  
   
2133      for(Int_t i=0; i<8; i++){      for(Int_t i=0; i<8; i++){
2134          if(tree_clone[i])tree_clone[i]->Fill();          if(tree_clone[i])tree_clone[i]->Fill();
2135      }      }
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.23