/[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.9 by pam-fi, Tue Dec 5 12:13:39 2006 UTC revision 1.14 by pam-fi, Tue Jan 16 10:20:21 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;      
   
 //     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;
# Line 87  PamLevel2::PamLevel2(){ Line 131  PamLevel2::PamLevel2(){
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    
# Line 136  void PamLevel2::Clear(){ Line 207  void PamLevel2::Clear(){
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  };  };
# Line 159  void *PamLevel2::GetPointerTo(const char Line 277  void *PamLevel2::GetPointerTo(const char
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       };       };
# Line 169  void *PamLevel2::GetPointerTo(const char Line 289  void *PamLevel2::GetPointerTo(const char
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;
# Line 190  void *PamLevel2::GetPointerTo(const char Line 310  void *PamLevel2::GetPointerTo(const char
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       };       };
# Line 199  void *PamLevel2::GetPointerTo(const char Line 321  void *PamLevel2::GetPointerTo(const char
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;
# Line 219  void *PamLevel2::GetPointerTo(const char Line 341  void *PamLevel2::GetPointerTo(const char
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            
# Line 244  void *PamLevel2::GetPointerTo(const char Line 369  void *PamLevel2::GetPointerTo(const char
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;
# Line 275  PamTrack* PamLevel2::GetStoredTrack(Int_ Line 405  PamTrack* PamLevel2::GetStoredTrack(Int_
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);
# Line 298  void PamLevel2::SortTracks(TString how){ Line 443  void PamLevel2::SortTracks(TString how){
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    
# Line 326  void PamLevel2::SortTracks(TString how){ Line 476  void PamLevel2::SortTracks(TString how){
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++;
# Line 349  void PamLevel2::SortTracks(TString how){ Line 501  void PamLevel2::SortTracks(TString how){
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;
# Line 363  void PamLevel2::SortTracks(TString how){ Line 517  void PamLevel2::SortTracks(TString how){
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                                    
# Line 398  void PamLevel2::SortTracks(TString how){ Line 552  void PamLevel2::SortTracks(TString how){
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  };  };
# Line 430  void PamLevel2::SortTracks(TString how){ Line 616  void PamLevel2::SortTracks(TString how){
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   //   //
# Line 448  TRefArray *PamLevel2::GetTracks(){ Line 645  TRefArray *PamLevel2::GetTracks(){
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;
# Line 476  PamTrack *PamLevel2::GetTrack(int it){ Line 698  PamTrack *PamLevel2::GetTrack(int it){
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;
# Line 507  PamTrack *PamLevel2::GetTrackImage(int i Line 757  PamTrack *PamLevel2::GetTrackImage(int i
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);
# Line 538  TTree *PamLevel2::GetPamTree(TFile *f, T Line 807  TTree *PamLevel2::GetPamTree(TFile *f, 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) {
# Line 597  TTree *PamLevel2::GetPamTree(TFile *f, T Line 869  TTree *PamLevel2::GetPamTree(TFile *f, T
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  //--------------------------------------  //--------------------------------------
# Line 639  TList*  PamLevel2::GetListOfLevel2Files( Line 911  TList*  PamLevel2::GetListOfLevel2Files(
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;
# Line 656  TList*  PamLevel2::GetListOfLevel2Files( Line 933  TList*  PamLevel2::GetListOfLevel2Files(
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();
# Line 717  TList*  PamLevel2::GetListOfLevel2Files( Line 994  TList*  PamLevel2::GetListOfLevel2Files(
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;
# Line 737  TChain *PamLevel2::GetPamTree(TList *fl, Line 1025  TChain *PamLevel2::GetPamTree(TList *fl,
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    
# Line 913  void PamLevel2::SetBranchAddress(TChain Line 1288  void PamLevel2::SetBranchAddress(TChain
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) {
# Line 968  void PamLevel2::SetBranchAddress(TChain Line 1343  void PamLevel2::SetBranchAddress(TChain
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  //  //
# Line 1053  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1404  TTree *PamLevel2::GetRunTree(TFile *f){
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;
# Line 1118  void PamLevel2::SetWhichTrees(TString de Line 1469  void PamLevel2::SetWhichTrees(TString de
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;
# Line 1147  void PamLevel2::SetWhichTrees(TString de Line 1487  void PamLevel2::SetWhichTrees(TString de
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  //  //
# Line 1191  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1698  Bool_t  PamLevel2::CheckLevel2File(TStri
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 ){
# Line 1316  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1823  Bool_t  PamLevel2::CheckLevel2File(TStri
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) {
# Line 1360  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1889  Bool_t  PamLevel2::CheckLevel2File(TStri
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                    
# Line 1384  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1925  Bool_t  PamLevel2::CheckLevel2File(TStri
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;
# Line 1400  void PamLevel2::CreateCloneTrees( TChain Line 1941  void PamLevel2::CreateCloneTrees( TChain
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)
# Line 1602  void PamLevel2::CreateCloneTrees( TChain Line 2058  void PamLevel2::CreateCloneTrees( TChain
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      }      }
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23