/[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.19 by pam-fi, Thu Jan 18 11:56:12 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 656  TList*  PamLevel2::GetListOfLevel2Files( Line 955  TList*  PamLevel2::GetListOfLevel2Files(
955              if (!in.good()) break;              if (!in.good()) break;
956              line++;              line++;
957  //          cout <<"(1) " << file << endl;  //          cout <<"(1) " << file << endl;
             if(file.Contains("#"))file = file(0,file.First("#"));  
 //          cout <<"(2) " << file << endl;  
958              if(file.IsNull()){              if(file.IsNull()){
959                  cout << "-- list interrupted at line "<<line <<endl;                  cout << "-- list interrupted at line "<<line <<endl;
960                  break;                  break;
961              }              }
962              if( gSystem->IsFileInIncludePath(file,&fullpath) ){              if(file.Contains("#"))file = file(0,file.First("#"));
963    //          cout <<"(2) " << file << endl;
964    //          if( gSystem->IsFileInIncludePath(file,&fullpath) ){
965    //          if( (fullpath = gSystem->FindFile(ddir,file)) ){
966                if( file.EndsWith(".root") ){
967                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(file));
968                  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
969                    delete fullpath;
970              }else{              }
971                  if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;  //          }else{
972              };  //              if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;
973    //          };
974          };          };
975          in.close();          in.close();
976                    
# Line 690  TList*  PamLevel2::GetListOfLevel2Files( Line 992  TList*  PamLevel2::GetListOfLevel2Files(
992          while ( (questo = (TSystemFile*) next()) ) {          while ( (questo = (TSystemFile*) next()) ) {
993              TString name =  questo-> GetName();              TString name =  questo-> GetName();
994              if( name.EndsWith(".root") ){              if( name.EndsWith(".root") ){
995                  char *fullpath;  //              const char *fullpath = gSystem->FindFile(ddir,name);
996                  gSystem->IsFileInIncludePath(name,&fullpath);  //              char *fullpath;
997    //              gSystem->IsFileInIncludePath(name,&fullpath);
998                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(name));
999                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));
1000                    delete fullpath;
1001              };              };
1002          }          }
1003          delete temp;          delete temp;
# Line 702  TList*  PamLevel2::GetListOfLevel2Files( Line 1007  TList*  PamLevel2::GetListOfLevel2Files(
1007      gSystem->ChangeDirectory(wdir); // back to the working directory      gSystem->ChangeDirectory(wdir); // back to the working directory
1008  //      cout << endl << "Selected files:" << endl;  //      cout << endl << "Selected files:" << endl;
1009  //      contents->Print();  //      contents->Print();
1010      cout << contents->GetEntries()<<" files selected\n";      cout << contents->GetEntries()<<" files \n";
1011  //      cout << endl;  //      cout << endl;
1012  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;
1013      return contents;      return contents;
# Line 717  TList*  PamLevel2::GetListOfLevel2Files( Line 1022  TList*  PamLevel2::GetListOfLevel2Files(
1022   * @param detlist String to select trees to be included   * @param detlist String to select trees to be included
1023   * @return Pointer to a TChain   * @return Pointer to a TChain
1024   */   */
1025  TChain *PamLevel2::GetPamTree(TList *fl, TString detlist="+ALL"){  TChain *PamLevel2::GetPamTree(TList *fl, TString detlist ){
   
 //      TChain *Tout=0;  
 //    if(Tout)Tout->Delete();  
 //    Tout = NULL;  
1026            
1027      TChain *Tout =0;      TChain *Trout =0;
1028    
1029    //    if( !detlist.IsNull() )SetWhichTrees(detlist);
1030      SetWhichTrees(detlist);      SetWhichTrees(detlist);
1031            
1032        cout<< "GetPamTree(TList*,TString): input detector list --> ";
1033        if(TRK1)cout<<"TRK1 ";
1034        if(TRK2)cout<<"TRK2 ";
1035        if(TRKh)cout<<"TRKH ";
1036        if(CAL1)cout<<"CAL1 ";
1037        if(CAL2)cout<<"CAL2 ";
1038        if(TOF)cout<<"TOF ";
1039        if(TRG)cout<<"TRG ";
1040        if(AC)cout<<"AC ";
1041        if(ND)cout<<"ND ";
1042        if(S4)cout<<"S4 ";
1043        if(ORB)cout<<"ORB ";
1044        cout << endl;
1045    
1046      TChain *T = 0;            TChain *T = 0;      
1047      TChain *C = 0;      TChain *C = 0;
1048      TChain *O = 0;      TChain *O = 0;
# Line 737  TChain *PamLevel2::GetPamTree(TList *fl, Line 1053  TChain *PamLevel2::GetPamTree(TList *fl,
1053      TChain *B = 0;      TChain *B = 0;
1054            
1055      if(TRK2||TRK1||TRKh) T = new TChain("Tracker");          if(TRK2||TRK1||TRKh) T = new TChain("Tracker");    
1056      if(CAL2||CAL1) C = new TChain("Calorimeter");      if(CAL2||CAL1)       C = new TChain("Calorimeter");
1057      if(TOF) O = new TChain("ToF");      if(TOF)              O = new TChain("ToF");
1058      if(TRG) R = new TChain("Trigger");      if(TRG)              R = new TChain("Trigger");
1059      if(S4)  S = new TChain("S4");      if(S4)               S = new TChain("S4");
1060      if(ND)  N = new TChain("NeutronD");      if(ND)               N = new TChain("NeutronD");
1061      if(AC)  A = new TChain("Anticounter");      if(AC)               A = new TChain("Anticounter");
1062      if(ORB) B = new TChain("OrbitalInfo");      if(ORB)              B = new TChain("OrbitalInfo");
1063            
1064      // loop over files and create chains              // loop over files and create chains        
1065      TIter next(fl);      TIter next(fl);
1066      TSystemFile *questo = 0;      TSystemFile *questo = 0;
1067      while ( (questo = (TSystemFile*) next()) ) {      while ( (questo = (TSystemFile*) next()) ) {
1068          TString name =  questo->GetName();          TString name =  questo->GetName();
1069  //              cout << "File: "<< name << endl;          cout << "File: "<< name << endl;
1070          if( CheckLevel2File(name) ){          if( CheckLevel2File(name) ){
1071              if(TRK2||TRK1||TRKh) T->Add(name);              if(TRK2||TRK1||TRKh) T->Add(name);
1072              if(CAL1||CAL2) C->Add(name);              if(CAL1||CAL2)       C->Add(name);
1073              if(TOF) O->Add(name);              if(TOF)              O->Add(name);
1074              if(TRG) R->Add(name);              if(TRG)              R->Add(name);
1075              if(S4)  S->Add(name);              if(S4)               S->Add(name);
1076              if(ND)  N->Add(name);              if(ND)               N->Add(name);
1077              if(AC)  A->Add(name);              if(AC)               A->Add(name);
1078              if(ORB) B->Add(name);              if(ORB)              B->Add(name);
1079          };          };
1080      }      }
1081            
1082      cout << "done chain \n";      cout << "done chain \n";
1083        
1084        UInt_t *found=0;
1085    // Tracker
1086        if(T && (TRK2||TRK1||TRKh)) {
1087            if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1088            else    T->SetBranchStatus("TrkLevel2",0,found);
1089            if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
1090            if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
1091            else    T->SetBranchStatus("TrkLevel1",0,found);
1092            if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
1093            if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
1094            else    T->SetBranchStatus("TrkHough",0,found);
1095            if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
1096            if(!Trout)Trout=T;
1097            else Trout->AddFriend("Tracker");
1098        }else{
1099            cout << "Tracker      : missing tree"<<endl;
1100        };
1101      // Calorimeter      // Calorimeter
1102      if(CAL2||CAL1) {                  if(C && (CAL2||CAL1)) {
1103          if(!Tout)Tout=C;          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
1104          else Tout->AddFriend("Calorimeter");          else    C->SetBranchStatus("CaloLevel2",0,found);
1105            if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
1106            if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
1107            else    C->SetBranchStatus("CaloLevel1",0,found);
1108            if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
1109            if(!Trout)Trout=C;
1110            else Trout->AddFriend("Calorimeter");
1111        }else{
1112            cout << "Calorimeter  : missing tree"<<endl;
1113      };      };
1114      // ToF          // ToF    
1115      if(TOF) {      if(O && TOF) {
1116          if(!Tout)Tout=O;          O->SetBranchAddress("ToFLevel2", GetPointerTo("ToFLevel2"));
1117          else Tout->AddFriend("ToF");          cout << "ToF          : set branch address ToFLevel2"<<endl;
1118            if(!Trout)Trout=O;
1119            else Trout->AddFriend("ToF");
1120        }else{
1121            cout << "ToF         : missing tree"<<endl;
1122      };      };
     // Tracker  
     if(TRK2||TRK1||TRKh){  
         if(!Tout)Tout=T;  
         else Tout->AddFriend("Tracker");  
     }  
1123      // Trigger      // Trigger
1124      if(TRG) {      if(R && TRG) {
1125          if(!Tout)Tout=R;          R->SetBranchAddress("TrigLevel2", GetPointerTo("TrigLevel2"));
1126          else Tout->AddFriend("Trigger");          cout << "Trigger      : set branch address TrigLevel2"<<endl;
1127            if(!Trout)Trout=O;
1128            else Trout->AddFriend("Trigger");
1129        }else{
1130            cout << "Trigger      : missing tree"<<endl;
1131      };      };
1132      // S4      // S4
1133      if(S4) {      if(S && S4) {
1134          if(!Tout)Tout=S;          S->SetBranchAddress("S4Level2", GetPointerTo("S4Level2"));
1135          else Tout->AddFriend("S4");          cout << "S4           : set branch address S4Level2"<<endl;
1136            if(!Trout)Trout=O;
1137            else Trout->AddFriend("S4");
1138        }else{
1139            cout << "S4           : missing tree"<<endl;
1140      };      };
1141      // Neutron Detector      // Neutron Detector
1142      if(ND) {      if(N && ND) {
1143          if(!Tout)Tout=N;          N->SetBranchAddress("NDLevel2", GetPointerTo("NDLevel2"));
1144          else Tout->AddFriend("NeutronD");          cout << "NeutronD     : set branch address NDLevel2"<<endl;
1145            if(!Trout)Trout=O;
1146            else Trout->AddFriend("NeutronD");
1147        }else{
1148            cout << "NeutronD     : missing tree"<<endl;
1149      };      };
1150      // Anticounters      // Anticounters
1151      if(AC) {      if(A && AC) {
1152          if(!Tout)Tout=A;          A->SetBranchAddress("AcLevel2", GetPointerTo("AcLevel2"));
1153          else Tout->AddFriend("Anticounter");          cout << "Anticounter  : set branch address AcLevel2"<<endl;
1154            if(!Trout)Trout=O;
1155            else Trout->AddFriend("Anticounter");
1156        }else{
1157            cout << "Anticounter  : missing tree"<<endl;
1158      };      };
1159      // OrbitalInfo      // Orbital Info
1160      if(ORB) {      if(B && ORB) {
1161          if(!Tout)Tout=B;          B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1162          else Tout->AddFriend("OrbitalInfo");          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1163            if(!Trout)Trout=O;
1164            else Trout->AddFriend("OrbitalInfo");
1165        }else{
1166            cout << "OrbitalInfo  : missing tree"<<endl;
1167      };      };
1168            
1169    //    cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
1170    
1171        
1172    //     // Tracker
1173    //     if(TRK2||TRK1||TRKh){
1174    //      if(!Tout)Tout=T;
1175    //      else Tout->AddFriend("Tracker");
1176    //     }
1177    //     // Calorimeter
1178    //     if(CAL2||CAL1) {        
1179    //      if(!Tout)Tout=C;
1180    //      else Tout->AddFriend("Calorimeter");
1181    //     };
1182    //     // ToF    
1183    //     if(TOF) {
1184    //      if(!Tout)Tout=O;
1185    //      else Tout->AddFriend("ToF");
1186    //     };
1187    //     // Trigger
1188    //     if(TRG) {
1189    //      if(!Tout)Tout=R;
1190    //      else Tout->AddFriend("Trigger");
1191    //     };
1192    //     // S4
1193    //     if(S4) {
1194    //      if(!Tout)Tout=S;
1195    //      else Tout->AddFriend("S4");
1196    //     };
1197    //     // Neutron Detector
1198    //     if(ND) {
1199    //      if(!Tout)Tout=N;
1200    //      else Tout->AddFriend("NeutronD");
1201    //     };
1202    //     // Anticounters
1203    //     if(AC) {
1204    //      if(!Tout)Tout=A;
1205    //      else Tout->AddFriend("Anticounter");
1206    //     };
1207    //     // OrbitalInfo
1208    //     if(ORB) {
1209    //      if(!Tout)Tout=B;
1210    //      else Tout->AddFriend("OrbitalInfo");
1211    //     };
1212        
1213  //    cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl;  //    cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl;
1214            
1215  //    if( Tout->GetEntries() )PamLevel2::SetBranchAddress();  //    if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout );
     if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout );  
1216            
1217      return Tout;      return Trout;
1218  }  }
1219    
1220    
# Line 913  void PamLevel2::SetBranchAddress(TChain Line 1316  void PamLevel2::SetBranchAddress(TChain
1316    
1317      // Tracker      // Tracker
1318       if(TRK2) {       if(TRK2) {
1319          t->SetBranchAddress("TrkLevel2", this->GetPointerTo("TrkLevel2"));          t->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1320          cout << "Tracker      : set branch address TrkLevel2"<<endl;          cout << "Tracker      : set branch address TrkLevel2"<<endl;
1321      };      };
1322      if(TRK1) {      if(TRK1) {
# Line 968  void PamLevel2::SetBranchAddress(TChain Line 1371  void PamLevel2::SetBranchAddress(TChain
1371            
1372  }  }
1373    
 // 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;  };  
1374    
1375  //--------------------------------------  //--------------------------------------
1376  //  //
# Line 1019  TChain *PamLevel2::GetRunTree(TList *fl) Line 1398  TChain *PamLevel2::GetRunTree(TList *fl)
1398            
1399      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1400      cout << "Run         : set branch address RunInfo"<<endl;      cout << "Run         : set branch address RunInfo"<<endl;
1401        R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1402        cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1403            
1404      return R;      return R;
1405            
# Line 1039  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1420  TTree *PamLevel2::GetRunTree(TFile *f){
1420            
1421      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1422      cout << "Run         : set branch address RunInfo"<<endl;      cout << "Run         : set branch address RunInfo"<<endl;
1423        R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1424        cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1425            
1426      return R;      return R;
1427            
1428  }  }
1429    /**
1430     * Update the runinfo informations (to be used to have Run infos event by event basis)
1431     * @param run Pointer to the chain/tree which contains run infos
1432     * @return true if a new run has been read, false if it is still the same run
1433     */
1434    Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev){
1435      //
1436      // check if we have already called once GetEntry, if not call it
1437      //
1438      if ( run->GetEntries() <= 0 ) return(false);
1439      //
1440      Int_t oldrun = irun;
1441      if ( irun < 0 ){
1442        irun = 0;
1443        run->GetEntry(0);
1444        runfirstentry = 0ULL;
1445        runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS) - 1ULL;
1446      };
1447      while ( iev > (runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS-1)) && irun < run->GetEntries() ){
1448        //    printf(" iev %llu %u %llu \n",iev,this->GetRunInfo()->NEVENTS,(ULong64_t)(runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS)));
1449        irun++;
1450        run->GetEntry(irun);
1451        runfirstentry = runlastentry+1ULL;
1452        runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS);
1453      };
1454      //
1455      if ( irun == oldrun || irun >= run->GetEntries() ) return(false);
1456      //
1457      //  printf(" iev %llu irun %i nevents %u 1st %llu last %llu \n",iev,irun,this->GetRunInfo()->NEVENTS,(ULong64_t)runfirstentry,(ULong64_t)runlastentry);
1458      //
1459      return(true);    
1460      //
1461    };
1462    /**
1463     * Update the runinfo informations (to be used to have Run infos event by event basis)
1464     * @param run Pointer to the chain/tree which contains run infos
1465     * @return true if a new run has been read, false if it is still the same run
1466     */
1467    Bool_t PamLevel2::UpdateRunInfo(TTree *run, ULong64_t iev){
1468      return(UpdateRunInfo((TChain*)run,iev));
1469    };
1470    
1471  //--------------------------------------  //--------------------------------------
1472  //  //
1473  //  //
# Line 1053  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1478  TTree *PamLevel2::GetRunTree(TFile *f){
1478  */  */
1479  void PamLevel2::SetWhichTrees(TString detlist){  void PamLevel2::SetWhichTrees(TString detlist){
1480                    
1481      if(detlist.Contains("+ALL", TString::kIgnoreCase)){      if(detlist.IsNull() || detlist.Contains("+ALL", TString::kIgnoreCase)){
1482          CAL1 = true;          CAL1 = true;
1483          CAL2 = true;          CAL2 = true;
1484          TRK2 = true;          TRK2 = true;
# Line 1118  void PamLevel2::SetWhichTrees(TString de Line 1543  void PamLevel2::SetWhichTrees(TString de
1543          TRKh=false;          TRKh=false;
1544      }      }
1545  //  -------------------------------------------------------------------------  //  -------------------------------------------------------------------------
   
   
       
 //     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;  
1546            
1547      if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false;      if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false;
1548      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 1561  void PamLevel2::SetWhichTrees(TString de
1561            
1562      if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false;      if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false;
1563      else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true;      else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true;
1564    
1565    //     cout<< "Set detector list --> ";
1566    //     if(TRK1)cout<<"TRK1 ";
1567    //     if(TRK2)cout<<"TRK2 ";
1568    //     if(TRKh)cout<<"TRKH ";
1569    //     if(CAL1)cout<<"CAL1 ";
1570    //     if(CAL2)cout<<"CAL2 ";
1571    //     if(TOF)cout<<"TOF ";
1572    //     if(TRG)cout<<"TRG ";
1573    //     if(AC)cout<<"AC ";
1574    //     if(ND)cout<<"ND ";
1575    //     if(S4)cout<<"S4 ";
1576    //     if(ORB)cout<<"ORB ";
1577    //     cout << endl;
1578            
1579  };  };
1580    
1581    
1582    /**
1583     * Set tree/branch detector flags from the content of a tree
1584     */
1585    void  PamLevel2::GetWhichTrees(TFile* f){
1586            
1587        // -----------
1588        // reset flags
1589        // -----------
1590        CAL1   = false;    
1591        CAL2   = false;    
1592        TRK2   = false;    
1593        TRK1   = false;    
1594        TRKh   = false;    
1595        TRG    = false;    
1596        TOF    = false;    
1597        S4     = false;    
1598        ND     = false;    
1599        AC     = false;    
1600        ORB    = false;    
1601        
1602        RUN    = false;
1603            
1604        cout << "Checking file: "<<f->GetName()<<endl;
1605        if( !f || f->IsZombie() ){
1606            cout << "File: "<< f->GetName() <<" Non valid root file"<< endl;
1607            return;
1608        }
1609    
1610        TList *lk = f->GetListOfKeys();
1611        TIter next(lk);
1612        TKey *key =0;
1613    
1614        Int_t nev = 0;
1615    
1616        while( (key = (TKey*)next()) ){
1617            
1618            if( !strcmp(key->GetName(),"Run"        ) )RUN = true;
1619    
1620            //=========================================================    
1621            if( !strcmp(key->GetName(),"Trigger"    ) ){
1622                TRG = true;
1623                Int_t nevt = ((TTree*)f->Get("Trigger"))->GetEntries();
1624                if( nev && nevt!=nev){
1625                    cout << "File: "<< f->GetName() <<" Trigger tree has "<<nevt<<" events instead of "<<nev<< endl;
1626                    TRG = false;
1627                }else nev=nevt;
1628            }
1629            //=========================================================    
1630            if( !strcmp(key->GetName(),"ToF"        ) ){
1631                TOF = true;
1632                Int_t nevt = ((TTree*)f->Get("ToF"))->GetEntries();
1633                if( nev && nevt!=nev){
1634                    cout << "File: "<< f->GetName() <<"     ToF tree has "<<nevt<<" events instead of "<<nev<< endl;
1635                    TOF = false;
1636                }else nev=nevt;
1637            }
1638            //=========================================================  
1639            if( !strcmp(key->GetName(),"S4"         ) ){
1640                S4 = true;
1641                Int_t nevt = ((TTree*)f->Get("S4"))->GetEntries();
1642                if( nev && nevt!=nev){
1643                    cout << "File: "<< f->GetName() <<"      S4 tree has "<<nevt<<" events instead of "<<nev<< endl;
1644                    S4 = false;
1645                }else nev=nevt;
1646            }
1647            //=========================================================  
1648    
1649            if( !strcmp(key->GetName(),"NeutronD"   ) ){
1650                ND = true;
1651                Int_t nevt = ((TTree*)f->Get("NeutronD"))->GetEntries();
1652                if( nev && nevt!=nev){
1653                    cout << "File: "<< f->GetName() <<"NeutronD tree has "<<nevt<<" events instead of "<<nev<< endl;
1654                    ND =false;
1655                }else nev=nevt;
1656            }      
1657            //=========================================================  
1658            if( !strcmp(key->GetName(),"Anticounter") ){
1659                AC = true;
1660                Int_t nevt = ((TTree*)f->Get("Anticounter"))->GetEntries();
1661                if( nev && nevt!=nev){
1662                    cout << "File: "<< f->GetName() <<" Anticounter tree has "<<nevt<<" events instead of "<<nev<< endl;
1663                    AC =false;
1664                }else nev=nevt;
1665            }
1666            //=========================================================  
1667            if( !strcmp(key->GetName(),"OrbitalInfo") ){
1668                ORB = true;
1669                Int_t nevt = ((TTree*)f->Get("OrbitalInfo"))->GetEntries();
1670                if( nev && nevt!=nev){
1671                    cout << "File: "<< f->GetName() <<" OrbitalInfo tree has "<<nevt<<" events instead of "<<nev<< endl;
1672                    ORB = false;
1673                }else nev=nevt;
1674            }
1675            //=========================================================  
1676            if( !strcmp(key->GetName(),"Tracker"    ) ){
1677                TTree *T = (TTree*)f->Get("Tracker");
1678                for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){
1679                    TString name = T->GetListOfBranches()->At(i)->GetName();
1680                    if( !name.CompareTo("TrkLevel1") )TRK1=true;
1681                    if( !name.CompareTo("TrkLevel2") )TRK2=true;
1682                    if( !name.CompareTo("TrkHough") )TRKh=true;
1683                };      
1684                Int_t nevt = T->GetEntries();
1685                if( nev && nevt!=nev){
1686                    cout << "File: "<< f->GetName() <<" Tracker tree has "<<nevt<<" events instead of "<<nev<< endl;
1687                    TRK1 = false;
1688                    TRK2 = false;
1689                    TRKh = false;
1690                }else nev=nevt;
1691                T->Delete();
1692            };
1693            //=========================================================  
1694            if( !strcmp(key->GetName(),"Calorimeter"    ) ){
1695                TTree *T = (TTree*)f->Get("Calorimeter");
1696                for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){
1697                    TString name = T->GetListOfBranches()->At(i)->GetName();
1698                    if( !name.CompareTo("CaloLevel1") )CAL1=true;
1699                    if( !name.CompareTo("CaloLevel2") )CAL2=true;
1700                };    
1701                Int_t nevt = T->GetEntries();
1702                if( nev && nevt!=nev){
1703                    cout << "File: "<< f->GetName() <<"  Calorimeter tree has "<<nevt<<" events instead of "<<nev<< endl;
1704                    CAL1 = false;
1705                    CAL2 = false;
1706                }else nev=nevt;
1707                T->Delete();
1708            };      
1709    
1710        };
1711        
1712        delete lk;
1713    
1714    //     cout<< "Get detector list from input file --> ";
1715    //     if(TRK1)cout<<"TRK1 ";
1716    //     if(TRK2)cout<<"TRK2 ";
1717    //     if(TRKh)cout<<"TRKH ";
1718    //     if(CAL1)cout<<"CAL1 ";
1719    //     if(CAL2)cout<<"CAL2 ";
1720    //     if(TOF)cout<<"TOF ";
1721    //     if(TRG)cout<<"TRG ";
1722    //     if(AC)cout<<"AC ";
1723    //     if(ND)cout<<"ND ";
1724    //     if(S4)cout<<"S4 ";
1725    //     if(ORB)cout<<"ORB ";
1726    //     cout << endl;
1727          
1728        return ;
1729            
1730    };
1731    
1732    
1733  //--------------------------------------  //--------------------------------------
1734  //  //
1735  //  //
# Line 1191  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1772  Bool_t  PamLevel2::CheckLevel2File(TStri
1772    
1773      while( (key = (TKey*)next()) ){      while( (key = (TKey*)next()) ){
1774                    
1775  //      cout << key->GetName() << endl;  //      cout << key->GetName() << ""<<key->GetClassName()<<endl;
1776  //      cout << " Get tree: " << f->Get(key->GetName())<<endl;  //              cout << " Get tree: " << f->Get(key->GetName())<<endl;
1777  //      nev_previous = nev;  //      nev_previous = nev;
1778  //      cout << " n.entries  "<< nev <<endl;  //      cout << " n.entries  "<< nev <<endl;
1779  //      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 1897  Bool_t  PamLevel2::CheckLevel2File(TStri
1897    
1898      };      };
1899            
1900    //     cout<< "CheckLevel2File(TString): detector list --> ";
1901    //     if(TRK1__ok)cout<<"TRK1 ";
1902    //     if(TRK2__ok)cout<<"TRK2 ";
1903    //     if(TRKh__ok)cout<<"TRKH ";
1904    //     if(CAL1__ok)cout<<"CAL1 ";
1905    //     if(CAL2__ok)cout<<"CAL2 ";
1906    //     if(TOF__ok)cout<<"TOF ";
1907    //     if(TRG__ok)cout<<"TRG ";
1908    //     if(AC__ok)cout<<"AC ";
1909    //     if(ND__ok)cout<<"ND ";
1910    //     if(S4__ok)cout<<"S4 ";
1911    //     if(ORB__ok)cout<<"ORB ";
1912    //     cout << endl;
1913    
1914    
1915        if(TRK2 && TRK1__ok)TRK1=1;
1916    // ----------------------------------------------------------------------------
1917    // NOTA
1918    // se c'e` il level1, lo devo necessarimente leggere.
1919    // infatti (non ho capito perche`) i cluster vengono letti e allocati in memoria
1920    // comunque, ma non vengono disallocati da PamLevel2::Clear()
1921    // ----------------------------------------------------------------------------
1922    
1923    
1924      if(!RUN__ok) {      if(!RUN__ok) {
# Line 1360  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1963  Bool_t  PamLevel2::CheckLevel2File(TStri
1963          return false;          return false;
1964      };      };
1965    
1966          if(ND && !ND__ok){      if(ND && !ND__ok){
1967              cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl;          cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl;
1968              return false;          return false;
1969          };      };
1970          if(TRG && !TRG__ok){      if(TRG && !TRG__ok){
1971              cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl;          cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl;
1972              return false;          return false;
1973          };      };
1974        
   
   
   
1975    
1976      lk->Delete();  //    lk->Delete();
1977    //    delete lk;
1978      f->Close();      f->Close();
1979    
1980    //     cout<< "CheckLevel2File(TString): detector list --> ";
1981    //     if(TRK1)cout<<"TRK1 ";
1982    //     if(TRK2)cout<<"TRK2 ";
1983    //     if(TRKh)cout<<"TRKH ";
1984    //     if(CAL1)cout<<"CAL1 ";
1985    //     if(CAL2)cout<<"CAL2 ";
1986    //     if(TOF)cout<<"TOF ";
1987    //     if(TRG)cout<<"TRG ";
1988    //     if(AC)cout<<"AC ";
1989    //     if(ND)cout<<"ND ";
1990    //     if(S4)cout<<"S4 ";
1991    //     if(ORB)cout<<"ORB ";
1992    //     cout << endl;
1993                
1994      return true;      return true;
1995                    
# Line 1384  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1999  Bool_t  PamLevel2::CheckLevel2File(TStri
1999  /**  /**
2000   * Create clone-trees   * Create clone-trees
2001   */   */
2002  void PamLevel2::CreateCloneTrees( TChain *fChain ){  void PamLevel2::CreateCloneTrees0( TChain *fChain, TFile *ofile ){
2003    
2004      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2005      cout << "Create clones of PAMELA trees "<<endl;      cout << "Create clones of PAMELA trees "<<endl;
# Line 1400  void PamLevel2::CreateCloneTrees( TChain Line 2015  void PamLevel2::CreateCloneTrees( TChain
2015      TList *li = fChain->GetListOfFriends();      TList *li = fChain->GetListOfFriends();
2016      TIter next(li);      TIter next(li);
2017      TFriendElement* T_friend=0;      TFriendElement* T_friend=0;
2018        ofile->cd();
2019      while( (T_friend = (TFriendElement*)next()) ){      while( (T_friend = (TFriendElement*)next()) ){
2020  //      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;
2021  //      cout<<T_friend->GetTree()->GetName()<< endl;  //      cout<<T_friend->GetTree()->GetName()<< endl;
2022          tree_clone[i] = T_friend->GetTree()->CloneTree(0);                tree_clone[i] = T_friend->GetTree()->CloneTree(0);
2023            tree_clone[i]->SetAutoSave(1000000);
2024          name =  tree_clone[i]->GetName();          name =  tree_clone[i]->GetName();
2025          name.Append("_clone");          name.Append("_clone");
2026  //      tree_clone[i]->SetName(name.Data());  //      tree_clone[i]->SetName(name.Data());
2027          cout << tree_clone[i]->GetName() << endl;          cout << tree_clone[i]->GetName() << endl;
2028          i++;          i++;
2029      }      }
2030        
2031        delete li;
2032    
2033      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2034    
2035  }  }
2036    
2037  // /**  /**
2038  //  * Create clone-trees   * Create clone-trees
2039  //  */   */
2040  // 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;  
 //     }  
2041    
2042  //     // Calorimeter      ofile->cd();
 //     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;  
 //      };  
 //     }      
2043    
2044  //     // ToF          cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2045  //     if(TOF) {      cout << "Create new PAMELA trees "<<endl;
2046  //      O_clone = new TTree("ToF_clone","PAMELA ToF level2 data ");              
2047  //      tof_clone = new ToFLevel2();      Int_t i=0;
 //      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;  
   
 // }  
2048    
2049  // /**      if(TRK1||TRK2||TRKh){
2050  //  * Create clone-trees          tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data ");    
2051  //  */          if(TRK1) {
2052  // void PamLevel2::CreateCloneTrees(){              tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1"));
2053                tree_clone[i]->BranchRef();
2054                cout << "Tracker      : branch TrkLevel1"<<endl;
2055            };
2056            if(TRK2) {
2057                tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2"));
2058                cout << "Tracker      : branch TrkLevel2"<<endl;        
2059            };
2060            if(TRKh) {
2061                tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough"));
2062                cout << "Tracker      : branch TrkHough"<<endl;
2063            };
2064            i++;
2065        }
2066    
2067  //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      // Calorimeter
2068  //     cout << "Create new PAMELA trees "<<endl;      if(CAL1||CAL2){
2069                    tree_clone[i] = new TTree("Calorimeter","PAMELA calorimeter level2 data ");    
2070  //     if(TRK1||TRK2||TRKh){          if(CAL1) {
2071  //      T_clone = new TTree("Tracker_clone","PAMELA tracker level2 data ");                  tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1"));
2072  //      if(TRK1) {              cout << "Calorimeter  : branch CaloLevel1"<<endl;
2073  //          trk1_clone = new TrkLevel1();          };
2074  //          T_clone->Branch("TrkLevel1","TrkLevel1", &trk1_clone);          if(CAL2) {
2075  //          T_clone->BranchRef();              tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2"));
2076  //          cout << "Tracker      : branch TrkLevel1"<<endl;              cout << "Calorimeter  : branch CaloLevel2"<<endl;
2077  //      };          };
2078  //      if(TRK2) {          i++;
2079  //          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;  
 //     }  
2080    
2081  //     // Calorimeter      // ToF    
2082  //     if(CAL1||CAL2){      if(TOF) {
2083  //      C_clone = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data ");              tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data ");    
2084  //      if(CAL1) {          tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2"));
2085  //          calo1_clone = new CaloLevel1();          cout << "ToF          : branch ToFLevel2"<<endl;
2086  //          C_clone->Branch("CaloLevel1", "CaloLevel1", &calo1_clone);          i++;
2087  //          cout << "Calorimeter  : branch CaloLevel1"<<endl;      };
2088  //      };      // Trigger
2089  //      if(CAL2) {      if(TRG) {
2090  //          calo2_clone = new CaloLevel2();          tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data ");    
2091  //          C_clone->Branch("CaloLevel2","CaloLevel2", &calo2_clone);          tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2"));
2092  //          cout << "Calorimeter  : branch CaloLevel2"<<endl;          cout << "Trigger      : branch TrigLevel2"<<endl;
2093  //      };          i++;
2094  //     }          };
2095        // S4
2096        if(S4) {
2097            tree_clone[i] = new TTree("S4","PAMELA S4 level2 data ");      
2098            tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2"));
2099            cout << "S4           : branch S4Level2"<<endl;
2100            i++;
2101        };
2102        // Neutron Detector
2103        if(ND) {
2104            tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data ");  
2105            tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2"));
2106            cout << "NeutronD     : branch NDLevel2"<<endl;
2107            i++;
2108        };
2109        // Anticounters
2110        if(AC) {
2111            tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data ");    
2112            tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2"));
2113            cout << "Anticounter  : branch AcLevel2"<<endl;
2114            i++;
2115        };
2116        // OrbitalInfo
2117        if(ORB) {
2118            tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info  ");      
2119            tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo"));
2120            cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;
2121            i++;
2122        };
2123        cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2124    
2125  //     // 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;  
2126    
 // }  
2127    
2128  /**  /**
2129   * Fill tree (created with CreateCloneTrees)   * Fill tree (created with CreateCloneTrees)
# Line 1602  void PamLevel2::CreateCloneTrees( TChain Line 2132  void PamLevel2::CreateCloneTrees( TChain
2132  //void PamLevel2::FillNewPamTree(TTree *T){  //void PamLevel2::FillNewPamTree(TTree *T){
2133  void PamLevel2::FillCloneTrees(){  void PamLevel2::FillCloneTrees(){
2134                    
 // //  --------------------------------------  
 // //  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();  
   
2135      for(Int_t i=0; i<8; i++){      for(Int_t i=0; i<8; i++){
2136          if(tree_clone[i])tree_clone[i]->Fill();          if(tree_clone[i])tree_clone[i]->Fill();
2137      }      }
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.23