/[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.26 by pam-fi, Fri Mar 16 20:27:55 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        
104      run_obj   = new GL_RUN();      h0_obj    = 0;
105            trk0_obj  = 0;
106  //     trk1_clone  = 0;  
107  //     trk2_clone  = 0;      trk2_obj  = 0;
108  //     trkh_clone  = 0;      trk1_obj  = 0;
109  //     calo1_clone = 0;      trkh_obj  = 0;
110  //     calo2_clone = 0;      calo1_obj = 0;
111  //     tof_clone   = 0;      calo2_obj = 0;
112  //     trig_clone  = 0;      tof_obj   = 0;
113  //     s4_clone    = 0;      trig_obj  = 0;
114  //     nd_clone    = 0;      s4_obj    = 0;
115  //     ac_clone    = 0;      nd_obj    = 0;
116  //     orb_clone   = 0;      ac_obj    = 0;
117        orb_obj   = 0;
118  //     T_clone  = 0;  
119  //     C_clone  = 0;      run_obj   = 0;//new GL_RUN();
120  //     O_clone  = 0;      soft_obj   = 0;// Emiliano
121  //     R_clone  = 0;      irun = -1;
122  //     S_clone  = 0;      runfirstentry = 0ULL;
123  //     N_clone  = 0;      runlastentry = 0ULL;
124  //     A_clone  = 0;  
125  //     B_clone  = 0;      l0_file = NULL;
126        l0_tree = NULL;
127        iroot   = -1;
128        dbc     = 0;
129        
130        irun = -1;
131        run_tree = NULL;
132        run_tree_clone = NULL;
133        sel_tree = NULL;
134        sel_tree_clone = NULL;
135        
136        irunentry = -1;
137        pam_tree = NULL;
138        for(Int_t i=0; i<NCLONES; i++ )pam_tree_clone[i]=NULL;
139    
140        host = "mysql://localhost/pamelaprod";
141        user = "anonymous";
142        psw = "";
143        const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
144        const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
145        const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
146        if ( !pamdbhost ) pamdbhost = "";
147        if ( !pamdbuser ) pamdbuser = "";
148        if ( !pamdbpsw ) pamdbpsw = "";
149        if ( strcmp(pamdbhost,"") ) host = pamdbhost;
150        if ( strcmp(pamdbuser,"") ) user = pamdbuser;
151        if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
152    
153      for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL;      
154    
155      sorted_tracks = new TRefArray();  //    sorted_tracks = 0;//new TRefArray();
156            
157        CAL0 = false;
158      CAL1 = true;      CAL1 = true;
159      CAL2 = true;      CAL2 = true;
160      TRK2 = true;      TRK2 = true;
161      TRK1 = false;      TRK1 = false;
162      TRKh = false;      TRKh = false;
163      TRG = true;      TRKh = false;
164      TOF = true;      TRG  = true;
165      S4  = true;      TOF  = true;
166      ND  = true;      TOF0 = false;
167      AC  = true;      S4   = true;
168      ORB = true;      ND   = true;
169        AC   = true;
170        ORB  = true;
171            
172        RUN  = true;
173    
174        SELLI = -1;
175    
176        tsorted=0;
177        timage=0;
178    
179  };  };
180  /**  /**
181   * Destructor   * Delete the event (act as Dtor)
182   */   */
183  PamLevel2::~PamLevel2(){                  void PamLevel2::Delete(){
184      delete run_obj;          
185      Clear();      if(run_obj)delete run_obj;
186        if(soft_obj)delete soft_obj; //Emiliano
187    
188    //    cout << "void PamLevel2::Clear()"<<endl;
189        if(h0_obj)   delete h0_obj;
190        if(trk0_obj) delete trk0_obj;      
191        if(trk1_obj) delete trk1_obj;      
192        if(trk2_obj) delete trk2_obj;
193        if(trkh_obj) delete trkh_obj;
194        if(calo1_obj)delete calo1_obj;
195        if(calo2_obj)delete calo2_obj;
196        if(tof_obj)  delete tof_obj;
197        if(trig_obj) delete trig_obj;
198        if(s4_obj)   delete s4_obj;
199        if(nd_obj)   delete nd_obj;
200        if(ac_obj)   delete ac_obj;
201        if(orb_obj)  delete orb_obj;
202        
203        if(tsorted){
204            tsorted->Delete();
205            delete tsorted;
206        }
207        if(timage){
208            timage->Delete();
209            delete timage;
210        }
211    
212        if(dbc){
213            dbc->Close();
214            delete dbc;
215        }
216        
217        if(l0_file)l0_file->Close();
218        if(pam_tree)pam_tree->Delete();;
219        if(run_tree)run_tree->Delete();;
220        if(sel_tree)sel_tree->Delete();;
221        for(Int_t i=0; i<NCLONES; i++ )if(pam_tree_clone[i])pam_tree_clone[i]->Delete();;
222        if(run_tree_clone)run_tree_clone->Delete();;
223        if(sel_tree_clone)sel_tree_clone->Delete();;
224        
225    
226  };  };
227    
228  /**  /**
229   * Clear the event   * Clear the event (NB! does not deallocate objects)
230   */   */
231  void PamLevel2::Clear(){  void PamLevel2::Clear(){
232                    
233      TrkLevel1::Clear();  //    cout << "void PamLevel2::Clear()"<<endl;
234            
235      TrkLevel2::Clear();  //
236      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
237      ToFLevel2::Clear();  // want to load them for each event even if they are the same...
238      TrigLevel2::Clear();  //
239      S4Level2::Clear();  //    if(run_obj)delete run_obj;
240      NDLevel2::Clear();  //    if(run_obj) run_obj->Clear();  // Emiliano: Do not deallocate run_obj here, it will give segmentation fault! call clear instead
241      AcLevel2::Clear();  //    if(soft_obj) soft_obj->Clear();
242      OrbitalInfo::Clear();  
243        if(h0_obj)   h0_obj->Clear();      
244  //     if(trk1_clone)trk1_clone->Clear();      if(trk0_obj) trk0_obj->Clear();    
245  //     if(trk2_clone)trk2_clone->Clear();      if(trk1_obj) trk1_obj->Clear();    
246  //     if(trkh_clone)trkh_clone->Clear();      if(trk2_obj) trk2_obj->Clear();
247  //     if(calo1_clone)calo1_clone->Clear();      if(trkh_obj) trkh_obj->Clear();
248  //     if(calo2_clone)calo2_clone->Clear();      if(calo1_obj)calo1_obj->Clear();
249  //     if(tof_clone)tof_clone->Clear();      if(calo2_obj)calo2_obj->Clear();
250  //     if(trig_clone)trig_clone->Clear();      if(tof_obj)  tof_obj->Clear();
251  //     if(s4_clone)s4_clone->Clear();      if(trig_obj) trig_obj->Clear();
252  //     if(nd_clone)nd_clone->Clear();      if(s4_obj)   s4_obj->Clear();
253  //     if(ac_clone)ac_clone->Clear();      if(nd_obj)   nd_obj->Clear();
254  //     if(orb_clone)orb_clone->Clear();      if(ac_obj)   ac_obj->Clear();
255        if(orb_obj)  orb_obj->Clear();
256            
257      sorted_tracks->Delete(); // clean the reference array  //    if(sorted_tracks)sorted_tracks->Clear();
258    //    sorted_tracks.Clear();
259    
260        if(tsorted){
261            tsorted->Delete();
262        }
263        if(timage){
264            timage->Delete();
265        }
266            
267  };  };
268    
# Line 136  void PamLevel2::Clear(){ Line 274  void PamLevel2::Clear(){
274  void *PamLevel2::GetPointerTo(const char* c ){  void *PamLevel2::GetPointerTo(const char* c ){
275    
276      TString objname = c;      TString objname = c;
277      if(!objname.CompareTo("TrkLevel1"))return &trk1_obj;  
278      if(!objname.CompareTo("TrkLevel2"))return &trk2_obj;      if(!objname.CompareTo("TrkLevel1"))  {
279      if(!objname.CompareTo("TrkHough"))return &trkh_obj;          if(!trk1_obj){
280      if(!objname.CompareTo("CaloLevel1"))return &calo1_obj;              trk1_obj  = new TrkLevel1();
281      if(!objname.CompareTo("CaloLevel2"))return &calo2_obj;              trk1_obj->Set();
282      if(!objname.CompareTo("ToFLevel2"))return &tof_obj;          }
283      if(!objname.CompareTo("TrigLevel2"))return &trig_obj;          return &trk1_obj;
284      if(!objname.CompareTo("S4Level2"))return &s4_obj;      };
285      if(!objname.CompareTo("NDLevel2"))return &nd_obj;      if(!objname.CompareTo("TrkLevel2"))  {
286      if(!objname.CompareTo("AcLevel2"))return &ac_obj;          if(!trk2_obj){  
287      if(!objname.CompareTo("OrbitalInfo"))return &orb_obj;              trk2_obj  = new TrkLevel2();
288                trk2_obj->Set();
289            }
290            return &trk2_obj;
291        };
292        if(!objname.CompareTo("TrkHough"))   {
293            if(!trkh_obj)  trkh_obj  = new TrkHough();
294            return &trkh_obj;
295        };
296        if(!objname.CompareTo("CaloLevel1")) {
297            if(!calo1_obj) calo1_obj = new CaloLevel1();
298            return &calo1_obj;
299        };
300        if(!objname.CompareTo("CaloLevel2")) {
301            if(!calo2_obj){
302                calo2_obj = new CaloLevel2();
303                calo2_obj->Set();
304            };
305            return &calo2_obj;
306        };
307        if(!objname.CompareTo("ToFLevel2"))  {
308            if(!tof_obj){
309                tof_obj   = new ToFLevel2();
310                tof_obj->Set();
311            }
312            return &tof_obj;
313        };
314        if(!objname.CompareTo("TrigLevel2")) {
315            if(!trig_obj)  trig_obj  = new TrigLevel2();
316            return &trig_obj;
317        };
318        if(!objname.CompareTo("S4Level2"))   {
319            if(!s4_obj)    s4_obj    = new S4Level2();
320            return &s4_obj;
321        };
322        if(!objname.CompareTo("NDLevel2"))   {
323            if(!nd_obj)    nd_obj    = new NDLevel2();
324            return &nd_obj;
325        };
326        if(!objname.CompareTo("AcLevel2"))   {
327            if(!ac_obj)    ac_obj    = new AcLevel2();
328            return &ac_obj;
329        };
330        if(!objname.CompareTo("OrbitalInfo")){
331            if(!orb_obj)   orb_obj   = new OrbitalInfo();
332            return &orb_obj;
333        };
334        
335      if(!objname.CompareTo("RunInfo"))return &run_obj;      if(!objname.CompareTo("RunInfo"))return &run_obj;
336    
337        if(!objname.CompareTo("SoftInfo"))return &soft_obj; // Emiliano
338    
339      return NULL;      return NULL;
340  };  };
341  //--------------------------------------  //--------------------------------------
# Line 159  void *PamLevel2::GetPointerTo(const char Line 347  void *PamLevel2::GetPointerTo(const char
347   * (If seqno = -1 retrieves the self-trigger calorimeter track)   * (If seqno = -1 retrieves the self-trigger calorimeter track)
348   */   */
349   CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){   CaloTrkVar *PamLevel2::GetCaloStoredTrack(int seqno){
350                
351       if( CaloLevel2::ntrk()==0 ){       if( !calo2_obj )return 0;
352    
353         if( calo2_obj->CaloLevel2::ntrk()==0 ){
354           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;
355           return NULL;           return NULL;
356       };       };
# Line 169  void *PamLevel2::GetPointerTo(const char Line 359  void *PamLevel2::GetPointerTo(const char
359       Int_t it_calo=0;       Int_t it_calo=0;
360            
361       do{       do{
362           c = CaloLevel2::GetCaloTrkVar(it_calo);           c = calo2_obj->CaloLevel2::GetCaloTrkVar(it_calo);
363           it_calo++;           it_calo++;
364       } while( c && seqno != c->trkseqno && it_calo < CaloLevel2::ntrk());             } while( c && seqno != c->trkseqno && it_calo < calo2_obj->CaloLevel2::ntrk());    
365            
366       if(!c || seqno != c->trkseqno){       if(!c || seqno != c->trkseqno){
367           c = 0;           c = 0;
# Line 190  void *PamLevel2::GetPointerTo(const char Line 380  void *PamLevel2::GetPointerTo(const char
380   */   */
381   ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){   ToFTrkVar *PamLevel2::GetToFStoredTrack(int seqno){
382                    
383       if( ToFLevel2::ntrk()==0 ){       if( !tof_obj )return 0;
384    
385         if( tof_obj->ToFLevel2::ntrk()==0 ){
386           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;
387           return NULL;           return NULL;
388       };       };
# Line 199  void *PamLevel2::GetPointerTo(const char Line 391  void *PamLevel2::GetPointerTo(const char
391       Int_t it_tof=0;       Int_t it_tof=0;
392            
393       do{       do{
394           c = ToFLevel2::GetToFTrkVar(it_tof);           c = tof_obj->ToFLevel2::GetToFTrkVar(it_tof);
395           it_tof++;           it_tof++;
396       } while( c && seqno != c->trkseqno && it_tof < ToFLevel2::ntrk());       } while( c && seqno != c->trkseqno && it_tof < tof_obj->ToFLevel2::ntrk());        
397            
398       if(!c || seqno != c->trkseqno){       if(!c || seqno != c->trkseqno){
399           c = 0;           c = 0;
# Line 219  void *PamLevel2::GetPointerTo(const char Line 411  void *PamLevel2::GetPointerTo(const char
411   * 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.
412   */   */
413   PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){   PamTrack* PamLevel2::GetPamTrackAlong(TrkTrack* t){
414                
415         cout <<"PamLevel2::GetPamTrackAlong(TrkTrack* t) **obsolete** "<<endl;
416         cout <<"(if you use it, remember to delete the PamTrack object)"<<endl;
417    
418       CaloTrkVar *c = 0;       CaloTrkVar *c = 0;
419       ToFTrkVar  *o = 0;       ToFTrkVar  *o = 0;
420            
# Line 244  void *PamLevel2::GetPointerTo(const char Line 439  void *PamLevel2::GetPointerTo(const char
439    
440  PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){  PamTrack* PamLevel2::GetStoredTrack(Int_t itrk){
441            
442        cout <<"PamLevel2::GetStoredTrack(Int_t itrk) **to-be-updated** "<<endl;
443        cout <<"for the moment, better use separately the methods: TrkLevel2::GetStoredTrack(seqno) CaloLevel2::GetCaloTrkVar(Int_t notrack) ToFLevel2::GetToFTrkVar(Int_t notrack)"<<endl;
444        cout <<"(if you use it, remember to delete the PamTrack object)"<<endl;
445      PamTrack *track = 0;      PamTrack *track = 0;
446            
447      if( itrk >=0 && itrk < TrkLevel2::ntrk() ){      if( itrk >=0 && itrk < trk2_obj->TrkLevel2::ntrk() ){
448                    
449          TrkTrack *t = TrkLevel2::GetStoredTrack(itrk);          TrkTrack *t = trk2_obj->TrkLevel2::GetStoredTrack(itrk);
450          track = GetPamTrackAlong(t);          track = GetPamTrackAlong(t);
451                    
452            
453    
454      }else{      }else{
455          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;
456      };      };
457            
458      return track;      return track;
# Line 275  PamTrack* PamLevel2::GetStoredTrack(Int_ Line 475  PamTrack* PamLevel2::GetStoredTrack(Int_
475   */   */
476  void PamLevel2::SortTracks(TString how){  void PamLevel2::SortTracks(TString how){
477    
478    //    cout <<" PamLevel2::SortTracks(TString how) "<<endl;
479        if( !trk2_obj ){
480            cout << "void PamLevel2::SortTracks(TString how):  TrkLevel2 not loaded !!!";
481            return;
482        };
483  //    cout << "call SortTracs() "<<endl;  //    cout << "call SortTracs() "<<endl;
484    //Save current Object count    //Save current Object count
485      Int_t ObjectNumber = TProcessID::GetObjectCount();      Int_t ObjectNumber = TProcessID::GetObjectCount();
486    
487  //    cout << "ObjectNumber  "<<ObjectNumber <<endl;  //    cout << "ObjectNumber  "<<ObjectNumber <<endl;
488        
489      sorted_tracks->Delete(); //temporaneo???  //     if(!sorted_tracks)sorted_tracks = new TRefArray();
490            //     sorted_tracks->Clear();
491    //    sorted_tracks.Clear();
492    
493        if(!tsorted)tsorted = new TClonesArray("PamTrack",trk2_obj->GetNTracks());
494        tsorted->Delete();
495        TClonesArray &ttsorted = *tsorted;
496        if(!timage)timage = new TClonesArray("PamTrack",trk2_obj->GetNTracks());
497        timage->Delete();
498        TClonesArray &ttimage = *timage;
499    
500    
501      // loop over the tracks sorted by the tracker      // loop over the tracks sorted by the tracker
502      Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase);      Bool_t use_TRK = how.Contains("TRK", TString::kIgnoreCase);
503      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 513  void PamLevel2::SortTracks(TString how){
513    
514      //   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;
515            
516      for(Int_t i=0; i < TrkLevel2::GetNTracks(); i++){      for(Int_t i=0; i < trk2_obj->TrkLevel2::GetNTracks(); i++){
517                    
518          TrkTrack *ts = 0;          TrkTrack   *ts = 0;
519            CaloTrkVar *cs = 0;
520            ToFTrkVar  *os = 0;
521                    
522          // get tracker tracks          // get tracker tracks
523          TrkTrack   *tp = TrkLevel2::GetTrack(i);                    //tracker          TrkTrack   *tp = trk2_obj->TrkLevel2::GetTrack(i);                    //tracker
524          CaloTrkVar *cp = 0;          CaloTrkVar *cp = GetCaloStoredTrack(tp->GetSeqNo());
525          ToFTrkVar  *op = 0;          ToFTrkVar  *op = GetToFStoredTrack(tp->GetSeqNo());
526    
527            TrkTrack   *ti = 0;              //tracker (image)
528            CaloTrkVar *ci = 0;
529            ToFTrkVar  *oi = 0;
530  //      cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl;  //      cout << "trk track n. "<<i << " "<<hex<< tp <<dec<< endl;
531          // if track has an image, check image selection          // if track has an image, check image selection
532          if(tp->HasImage()){          if(tp->HasImage()){
533                            
534              TrkTrack   *ti = TrkLevel2::GetTrackImage(i);              //tracker (image)              ti = trk2_obj->TrkLevel2::GetTrackImage(i);              //tracker (image)
535              CaloTrkVar *ci = 0;              ci = GetCaloStoredTrack(ti->GetSeqNo());
536              ToFTrkVar  *oi = 0;              oi = GetToFStoredTrack(ti->GetSeqNo());
537                            
538  //          cout << "its image "<<i << " "<<hex<< ti <<dec<< endl;  //          cout << "its image "<<i << " "<<hex<< ti <<dec<< endl;
539    
# Line 326  void PamLevel2::SortTracks(TString how){ Line 546  void PamLevel2::SortTracks(TString how){
546              // ------------------------              // ------------------------
547              // check the Y spatial residual on the first calorimeter plane              // check the Y spatial residual on the first calorimeter plane
548              // (cut on calorimeter variables from Emiliano)              // (cut on calorimeter variables from Emiliano)
549                if( use_CAL && !calo2_obj ){
550                    cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but CaloLevel2 not loaded !!!";
551                    return;
552                };
553              if(              if(
554                  use_CAL            &&                  use_CAL            &&
555                  npcfit[1] > 15     &&   //no. of fit planes on Y view                  calo2_obj->npcfit[1] > 15     &&   //no. of fit planes on Y view
556                  varcfit[1] < 1000. &&  //fit variance on Y view                  calo2_obj->varcfit[1] < 1000. &&  //fit variance on Y view
557                    cp && ci &&
558                  true){                  true){
559    
                 cp = GetCaloStoredTrack(tp->GetSeqNo());  
                 ci = GetCaloStoredTrack(ti->GetSeqNo());  
560                                    
561                  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;
562                  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;
563                                    
564                  if(resy_p <= resy_i) tp_score++;                  if(resy_p <= resy_i) tp_score++;
565                  else                 ti_score++;                  else                 ti_score++;
# Line 349  void PamLevel2::SortTracks(TString how){ Line 572  void PamLevel2::SortTracks(TString how){
572              // ------------------------                  // ------------------------    
573              // check the number of hit pmts along the track              // check the number of hit pmts along the track
574              // on S12 S21 and S32, where paddles are parallel to Y axis              // on S12 S21 and S32, where paddles are parallel to Y axis
575              if( use_TOF ){              if( use_TOF && !tof_obj ){
576                    cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but ToFLevel2 not loaded !!!";
577                    return;
578                };
579                if( use_TOF && op && oi ){
580                                    
581                  Int_t nphit_p =0;                  Int_t nphit_p =0;
582                  Int_t nphit_i =0;                  Int_t nphit_i =0;
583                                    
                 op = GetToFStoredTrack(tp->GetSeqNo());  
                 oi = GetToFStoredTrack(ti->GetSeqNo());  
584                                    
585  /*                              cout << "track: npmtadc "<< op->npmtadc << endl;  /*                              cout << "track: npmtadc "<< op->npmtadc << endl;
586                                  cout << "track: npmttdc "<< op->npmttdc << endl;                                  cout << "track: npmttdc "<< op->npmttdc << endl;
# Line 363  void PamLevel2::SortTracks(TString how){ Line 588  void PamLevel2::SortTracks(TString how){
588                                  cout << "image: npmttdc "<< oi->npmttdc << endl;*/                                  cout << "image: npmttdc "<< oi->npmttdc << endl;*/
589                                    
590                  for (Int_t ih=0; ih < op->npmtadc; ih++){                  for (Int_t ih=0; ih < op->npmtadc; ih++){
591                      Int_t pl = GetPlaneIndex( (op->pmtadc).At(ih) );                      Int_t pl = tof_obj->GetPlaneIndex( (op->pmtadc).At(ih) );
592                      if(pl == 1 || pl == 2 || pl == 5)nphit_p++;                      if(pl == 1 || pl == 2 || pl == 5)nphit_p++;
593                  };                  };
594                                    
595                  for (Int_t ih=0; ih < oi->npmtadc; ih++){                  for (Int_t ih=0; ih < oi->npmtadc; ih++){
596                      Int_t pl = GetPlaneIndex( (oi->pmtadc).At(ih) );                      Int_t pl = tof_obj->GetPlaneIndex( (oi->pmtadc).At(ih) );
597                      if(pl == 1 || pl == 2 || pl == 5)nphit_i++;                      if(pl == 1 || pl == 2 || pl == 5)nphit_i++;
598                  };                  };
599                                    
# Line 398  void PamLevel2::SortTracks(TString how){ Line 623  void PamLevel2::SortTracks(TString how){
623              // ------------------------              // ------------------------
624              // the winner is....              // the winner is....
625              // ------------------------                  // ------------------------    
626              if      (tp_score > ti_score) ts = tp;//the track sorted by the tracker!!              if      (tp_score > ti_score) {
627              else if (tp_score < ti_score) ts = ti;//its image!!  //              ts = tp;//the track sorted by the tracker!!
628              else {  //              cs = cp;
629                  ts = tp;  //              os = op;
630                }else if (tp_score < ti_score) {
631                    ts = ti;//its image!!
632                    cs = ci;
633                    os = oi;
634    
635                    ti = tp;//its image!!
636                    ci = cp;
637                    oi = op;
638    
639                    tp = ts;//its image!!
640                    cp = cs;
641                    op = os;
642    
643                    
644                }else {
645    //              ts = tp;
646    //              cs = cp;
647    //              os = op;
648  //                              cout << "Warning - track image ambiguity not solved" << endl;  //                              cout << "Warning - track image ambiguity not solved" << endl;
649  //                              cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl;  //                              cout << ts->GetNtot() << " "<< ts->chi2 << " " << npcfit[1] << " "<< nphit_p << endl;
650              };              };
651                            
652          }else{          }else{
653              ts = tp;  //          ts = tp;
654    //          cs = cp;
655    //          os = op;
656          };          };
657                    
658  //      cout <<" SortTracks() "<<i<<" -- "<<ts<<endl;  //      cout <<" SortTracks() "<<i<<" -- "<<ts<<endl;
659          sorted_tracks->Add(ts);//save the track in the sorted array  //      sorted_tracks->Add(ts);//save the track in the sorted array
660    //      sorted_tracks.Add(ts);//save the track in the sorted array
661    //      sorted_tracks.Add(tp);//save the track in the sorted array
662  //      cout << "SortTracks:: sorted_tracks->Add(it) "<<i<<" "<<ts<<endl;  //      cout << "SortTracks:: sorted_tracks->Add(it) "<<i<<" "<<ts<<endl;
663            //      cout<<"o "<<tp<<endl;
664    //      cout<<"o "<<cp<<endl;
665    //      cout<<"o "<<op<<endl;
666            new(ttsorted[i]) PamTrack(tp,cp,op);
667            new(ttimage[i])  PamTrack(ti,ci,oi);
668      };      };
669    
670        if( tsorted->GetEntries() != trk2_obj->GetNTracks() ){
671            cout << "void PamLevel2::SortTracks(TString how): tsorted->GetEntries() "<<tsorted->GetEntries()<<" != trk2_obj->GetNTracks() = "<<trk2_obj->GetNTracks() <<endl;
672            tsorted->Delete(); tsorted=0;
673            timage->Delete(); timage=0;    
674        }
675    
676      //Restore Object count      //Restore Object count
677      //To save space in the table keeping track of all referenced objects      //To save space in the table keeping track of all referenced objects
678      //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.  
679      TProcessID::SetObjectCount(ObjectNumber);      TProcessID::SetObjectCount(ObjectNumber);
680            
681  };  };
# Line 430  void PamLevel2::SortTracks(TString how){ Line 687  void PamLevel2::SortTracks(TString how){
687   * 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.
688   * 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.
689   */   */
690  TRefArray *PamLevel2::GetTracks(){  // TRefArray *PamLevel2::GetTracks(){
691    
692    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
693    //     SortTracks("+CAL+TOF");
694    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
695    
696    // //   return  sorted_tracks;
697    //     return &sorted_tracks;
698        
699    //  };
700    TClonesArray *PamLevel2::GetTracks(){
701    
702  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
703      SortTracks("+CAL+TOF");      SortTracks("+CAL+TOF");
704  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
705    
706      return sorted_tracks;      return tsorted;
707        
708   };   };
709  //--------------------------------------  //--------------------------------------
710   //   //
# Line 448  TRefArray *PamLevel2::GetTracks(){ Line 716  TRefArray *PamLevel2::GetTracks(){
716   * @param it Track number, ranging from 0 to GetNTracks().   * @param it Track number, ranging from 0 to GetNTracks().
717   */   */
718  PamTrack *PamLevel2::GetTrack(int it){  PamTrack *PamLevel2::GetTrack(int it){
719    
720    //    if(!trk2_obj) return 0;
721                    
722    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
723    //     SortTracks("+CAL+TOF");
724    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
725    // //    if(!sorted_tracks)return 0;
726    //     if(sorted_tracks.IsEmpty())return 0;
727    
728    //     PamTrack *track = 0;
729        
730    //     if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() && it<sorted_tracks.GetEntriesFast() ){    
731    //      TrkTrack   *t = (TrkTrack*)sorted_tracks.At(it);
732    //      track = GetPamTrackAlong(t);
733    //     }else{
734    //      cout << "PamLevel2::GetTrack(int) : tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl;
735    //     };
736        
737    //     return track;
738    
739    //    cout << "PamLevel2::GetTrack(int it) "<<endl;
740  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
741      SortTracks("+CAL+TOF");      SortTracks("+CAL+TOF");
742  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
743            if(!tsorted)return 0;
744        if(!tsorted->GetEntries())return 0;
745    
746      PamTrack *track = 0;      PamTrack *track = 0;
747            
748      if( it >=0 && it < TrkLevel2::GetNTracks() && it<sorted_tracks->GetEntriesFast() ){      if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks()){
749          TrkTrack   *t = (TrkTrack*)sorted_tracks->At(it);  //      TrkTrack   *t = (TrkTrack*)sorted_tracks.At(it);
750          track = GetPamTrackAlong(t);  //      track = GetPamTrackAlong(t);
751            TClonesArray &t = *(tsorted);
752            track = (PamTrack*)t[it];
753            
754      }else{      }else{
755          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;
756      };      };
757            
758      return track;      return track;
# Line 476  PamTrack *PamLevel2::GetTrack(int it){ Line 769  PamTrack *PamLevel2::GetTrack(int it){
769   */   */
770  PamTrack *PamLevel2::GetTrackImage(int it){  PamTrack *PamLevel2::GetTrackImage(int it){
771    
772    //     if(!trk2_obj) return 0;
773    
774    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
775    //     SortTracks("+CAL+TOF");
776    // //  *-*-*-*-*-*-*-*-*-*-*-*-*
777    // //    if(!sorted_tracks)return 0;
778    //     if(sorted_tracks.IsEmpty())return 0;
779        
780    //     PamTrack *image = 0;
781        
782    //     if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){
783    //      TrkTrack *temp = (TrkTrack*)sorted_tracks.At(it);
784    //      if( temp->HasImage() ){
785    //          TrkTrack   *t = trk2_obj->TrkLevel2::GetStoredTrack(temp->GetImageSeqNo());
786    //          image = GetPamTrackAlong(t);
787    //      }else{
788    //          cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl;
789    //      };
790    //     }else{
791    //      cout << "PamLevel2::GetTrackImage(int) : Tracker track SeqNo "<< it <<" does not exist (GetNTracks() = "<<trk2_obj->TrkLevel2::GetNTracks()<<")"<<endl;
792    //     };
793        
794    //     return image;
795    
796    
797  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
798      SortTracks("+CAL+TOF");      SortTracks("+CAL+TOF");
799  //  *-*-*-*-*-*-*-*-*-*-*-*-*  //  *-*-*-*-*-*-*-*-*-*-*-*-*
800        if(!timage)return 0;
801        if(!timage->GetEntries())return 0;
802            
803      PamTrack *image = 0;      PamTrack *image = 0;
804            
805      if( it >=0 && it < TrkLevel2::GetNTracks() ){      if( it >=0 && it < trk2_obj->TrkLevel2::GetNTracks() ){
806          TrkTrack *temp = (TrkTrack*)sorted_tracks->At(it);          TClonesArray &t = *(tsorted);
807          if( temp->HasImage() ){          PamTrack *temp = (PamTrack*)t[it];
808              TrkTrack   *t = TrkLevel2::GetStoredTrack(temp->GetImageSeqNo());          if( temp->GetTrkTrack()->HasImage() ){
809              image = GetPamTrackAlong(t);              TClonesArray &t = *(timage);
810                image = (PamTrack*)t[it];
811          }else{          }else{
812              cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl;  //          cout <<"PamLevel2::GetTrackImage(int) : Track SeqNo "<<it<<" does not have image"<<endl;
813          };          };
814      }else{      }else{
815          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;
816      };      };
817            
818      return image;      return image;
# Line 507  PamTrack *PamLevel2::GetTrackImage(int i Line 828  PamTrack *PamLevel2::GetTrackImage(int i
828   * @param detlist String to select trees to be included   * @param detlist String to select trees to be included
829   * @return Pointer to a TTree   * @return Pointer to a TTree
830   */   */
831  TTree *PamLevel2::GetPamTree(TFile *f, TString detlist="+ALL"){  TTree *PamLevel2::GetPamTree(TFile *f, TString detlist ){
832    
833  //      cout << "WARNING!!! -- obsolete method -- \n";  //     if( !detlist.IsNull() ) SetWhichTrees(detlist);
834  //      cout << "better use TChain *PamLevel2::GetPamTree(TList*, TString) \n";  //     else                    GetWhichTrees(f);
835        
836        cout << "TTree *PamLevel2::GetPamTree(TFile *f, TString detlist ) -- obsolte "<<endl;
837    
838      SetWhichTrees(detlist);      SetWhichTrees(detlist);
     TTree *Trout =0;  
839    
840        TTree *Trout =0;
841    
842      TString fname = f->GetName();      TString fname = f->GetName();
843      if(!CheckLevel2File(fname))return NULL;      if(!CheckLevel2File(fname))return NULL;
844    
845        //    UInt_t *found=0;    
846    
847        cout<< "GetPamTree(TFile*,TString): detector list --> ";
848        if(TRK1)cout<<"TRK1 ";
849        if(TRK2)cout<<"TRK2 ";
850        if(TRKh)cout<<"TRKH ";
851        if(CAL1)cout<<"CAL1 ";
852        if(CAL2)cout<<"CAL2 ";
853        if(TOF)cout<<"TOF ";
854        if(TRG)cout<<"TRG ";
855        if(AC)cout<<"AC ";
856        if(ND)cout<<"ND ";
857        if(S4)cout<<"S4 ";
858        if(ORB)cout<<"ORB ";
859        cout << endl;
860        if(SELLI)cout<<">>> Found selection-list <<<"<<endl;
861    
862        f->cd();
863    
864  // Tracker  // Tracker
865      TTree *T = (TTree*)f->Get("Tracker");      TTree *T = (TTree*)f->Get("Tracker");
866      if(T && (TRK2||TRK1||TRKh)) {      if(T && (TRK2||TRK1||TRKh)) {
867          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
868    //      else    T->SetBranchStatus("TrkLevel2",0,found);
869          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
870          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
871    //      else    T->SetBranchStatus("TrkLevel1",0,found);
872          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
873          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
874    //      else    T->SetBranchStatus("TrkHough",0,found);
875          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
876          if(!Trout)Trout=T;          if(!Trout)Trout=T;
877          else Trout->AddFriend(T);          else Trout->AddFriend(T);
# Line 538  TTree *PamLevel2::GetPamTree(TFile *f, T Line 882  TTree *PamLevel2::GetPamTree(TFile *f, T
882      TTree *C = (TTree*)f->Get("Calorimeter");      TTree *C = (TTree*)f->Get("Calorimeter");
883      if(C && (CAL2||CAL1)) {      if(C && (CAL2||CAL1)) {
884          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
885    //      else    C->SetBranchStatus("CaloLevel2",0,found);
886          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
887          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
888    //      else    C->SetBranchStatus("CaloLevel1",0,found);
889          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
890          if(!Trout)Trout=C;          if(!Trout)Trout=C;
891          else Trout->AddFriend(C);          else Trout->AddFriend(C);
892      }else{      }else{
893          cout << "Calorimeter  : missing tree"<<endl;          cout << "Calorimeter  : missing tree"<<endl;
894      };      };
895    
896      // ToF          // ToF    
897      TTree *O = (TTree*)f->Get("ToF");      TTree *O = (TTree*)f->Get("ToF");
898      if(O && TOF) {      if(O && TOF) {
# Line 597  TTree *PamLevel2::GetPamTree(TFile *f, T Line 944  TTree *PamLevel2::GetPamTree(TFile *f, T
944          cout << "Anticounter  : missing tree"<<endl;          cout << "Anticounter  : missing tree"<<endl;
945      };      };
946      // Orbital Info      // Orbital Info
947          TTree *B = (TTree*)f->Get("OrbitalInfo");      TTree *B = (TTree*)f->Get("OrbitalInfo");
948          if(B && ORB) {      if(B && ORB) {
949                  B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));          B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
950                  cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
951                  if(!Trout)Trout=O;          if(!Trout)Trout=O;
952                  else Trout->AddFriend(B);          else Trout->AddFriend(B);
953          }else{      }else{
954                  cout << "OrbitalInfo  : missing tree"<<endl;          cout << "OrbitalInfo  : missing tree"<<endl;
955          };      };
956    
957        TTree *L = (TTree*)f->Get("SelectionList");
958        if(L && SELLI==1) {
959    //      L->SetBranchAddress("RunEntry",&irun);
960    //      cout << "SelectionList: set branch address RunEntry"<<endl;
961    //      L->SetBranchAddress("EventEntry",&irunentry);
962    //      cout << "SelectionList: set branch address EventEntry"<<endl;
963    //      if(!Trout)Trout=O;
964    //      else Trout->AddFriend("SelectionList");
965            cout << " TTree *PamLevel2::GetPamTree(TFile, TString) >>> SelectionList not implemente!!!!"<<endl;
966            sel_tree = 0;
967        }else{
968            cout << "SelectionList  : missing tree"<<endl;
969        };
970            
971          cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;      cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
972                        
973          return Trout;      return Trout;
974            
975  }  }
976  //--------------------------------------  //--------------------------------------
# Line 634  TList*  PamLevel2::GetListOfLevel2Files( Line 995  TList*  PamLevel2::GetListOfLevel2Files(
995      TList *contents  = new TList; // create output list      TList *contents  = new TList; // create output list
996      contents->SetOwner();      contents->SetOwner();
997            
998      char *fullpath;  //    char *fullpath;
999    //    const char *fullpath;
1000            
1001      // if no input file list is given:        // if no input file list is given:  
1002      if ( flisttxt != "" ){      if ( flisttxt != "" ){
1003                    
1004          if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) )return 0;  //      if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){        
1005            //      if( !(fullpath = gSystem->FindFile("./",flisttxt)) ){
1006  //              flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));  //          cout <<"File "<<flisttxt<<" not found"<<endl;
1007          flisttxt = fullpath;  //          return 0;
1008            //      }      
1009          if( !gSystem->ChangeDirectory(ddir) )return 0;  //      flisttxt = fullpath;
1010    
1011            flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));
1012    
1013            if( !gSystem->ChangeDirectory(ddir) ){
1014                cout << "Cannot change directory : "<<ddir<<endl;
1015                return 0;
1016            }
1017                    
1018          cout <<"Input file list : " << flisttxt <<endl;          cout <<"Input file list : " << flisttxt <<endl;
1019          ifstream in;          ifstream in;
1020          in.open(flisttxt, ios::in); //open input file list          in.open(flisttxt, ios::in); //open input file list
1021            if(!in.good()){
1022                cout <<" ERROR opening the file "<<endl;
1023                gSystem->ChangeDirectory(wdir); // back to the working directory
1024                return 0;
1025            }      
1026          int line=0;          int line=0;
1027          while (1) {          while (1) {
1028              TString file;              TString file;
# Line 656  TList*  PamLevel2::GetListOfLevel2Files( Line 1030  TList*  PamLevel2::GetListOfLevel2Files(
1030              if (!in.good()) break;              if (!in.good()) break;
1031              line++;              line++;
1032  //          cout <<"(1) " << file << endl;  //          cout <<"(1) " << file << endl;
             if(file.Contains("#"))file = file(0,file.First("#"));  
 //          cout <<"(2) " << file << endl;  
1033              if(file.IsNull()){              if(file.IsNull()){
1034                  cout << "-- list interrupted at line "<<line <<endl;                  cout << "-- list interrupted at line "<<line <<endl;
1035                  break;                  break;
1036              }              }
1037              if( gSystem->IsFileInIncludePath(file,&fullpath) ){              if(file.Contains("#"))file = file(0,file.First("#"));
1038    //          cout <<"(2) " << file << endl;
1039    //          if( gSystem->IsFileInIncludePath(file,&fullpath) ){
1040    //          if( (fullpath = gSystem->FindFile(ddir,file)) ){
1041                if( file.EndsWith(".root") ){
1042                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(file));
1043                  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
1044                    delete fullpath;
1045              }else{              }
1046                  if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;  //          }else{
1047              };  //              if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;
1048    //          };
1049          };          };
1050          in.close();          in.close();
1051                    
# Line 690  TList*  PamLevel2::GetListOfLevel2Files( Line 1067  TList*  PamLevel2::GetListOfLevel2Files(
1067          while ( (questo = (TSystemFile*) next()) ) {          while ( (questo = (TSystemFile*) next()) ) {
1068              TString name =  questo-> GetName();              TString name =  questo-> GetName();
1069              if( name.EndsWith(".root") ){              if( name.EndsWith(".root") ){
1070                  char *fullpath;  //              const char *fullpath = gSystem->FindFile(ddir,name);
1071                  gSystem->IsFileInIncludePath(name,&fullpath);  //              char *fullpath;
1072    //              gSystem->IsFileInIncludePath(name,&fullpath);
1073                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(name));
1074                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));
1075                    delete fullpath;
1076              };              };
1077          }          }
1078          delete temp;          delete temp;
# Line 702  TList*  PamLevel2::GetListOfLevel2Files( Line 1082  TList*  PamLevel2::GetListOfLevel2Files(
1082      gSystem->ChangeDirectory(wdir); // back to the working directory      gSystem->ChangeDirectory(wdir); // back to the working directory
1083  //      cout << endl << "Selected files:" << endl;  //      cout << endl << "Selected files:" << endl;
1084  //      contents->Print();  //      contents->Print();
1085      cout << contents->GetEntries()<<" files selected\n";      cout << contents->GetEntries()<<" files \n";
1086  //      cout << endl;  //      cout << endl;
1087  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;
1088      return contents;      return contents;
# Line 717  TList*  PamLevel2::GetListOfLevel2Files( Line 1097  TList*  PamLevel2::GetListOfLevel2Files(
1097   * @param detlist String to select trees to be included   * @param detlist String to select trees to be included
1098   * @return Pointer to a TChain   * @return Pointer to a TChain
1099   */   */
1100  TChain *PamLevel2::GetPamTree(TList *fl, TString detlist="+ALL"){  TChain *PamLevel2::GetPamTree(TList *fl, TString detlist ){
   
 //      TChain *Tout=0;  
 //    if(Tout)Tout->Delete();  
 //    Tout = NULL;  
1101            
1102      TChain *Tout =0;      TChain *Trout =0;
1103    
1104    //    if( !detlist.IsNull() )SetWhichTrees(detlist);
1105      SetWhichTrees(detlist);      SetWhichTrees(detlist);
1106            
1107        cout<< "GetPamTree(TList*,TString): input detector list --> ";
1108        if(TRK1)cout<<"TRK1 ";
1109        if(TRK2)cout<<"TRK2 ";
1110        if(TRKh)cout<<"TRKH ";
1111        if(CAL1)cout<<"CAL1 ";
1112        if(CAL2)cout<<"CAL2 ";
1113        if(TOF)cout<<"TOF ";
1114        if(TRG)cout<<"TRG ";
1115        if(AC)cout<<"AC ";
1116        if(ND)cout<<"ND ";
1117        if(S4)cout<<"S4 ";
1118        if(ORB)cout<<"ORB ";
1119        cout << endl;
1120    
1121      TChain *T = 0;            TChain *T = 0;      
1122      TChain *C = 0;      TChain *C = 0;
1123      TChain *O = 0;      TChain *O = 0;
# Line 735  TChain *PamLevel2::GetPamTree(TList *fl, Line 1126  TChain *PamLevel2::GetPamTree(TList *fl,
1126      TChain *N = 0;      TChain *N = 0;
1127      TChain *A = 0;      TChain *A = 0;
1128      TChain *B = 0;      TChain *B = 0;
1129    
1130        TChain *L = 0;
1131            
1132      if(TRK2||TRK1||TRKh) T = new TChain("Tracker");          if(TRK2||TRK1||TRKh) T = new TChain("Tracker");    
1133      if(CAL2||CAL1) C = new TChain("Calorimeter");      if(CAL2||CAL1)       C = new TChain("Calorimeter");
1134      if(TOF) O = new TChain("ToF");      if(TOF)              O = new TChain("ToF");
1135      if(TRG) R = new TChain("Trigger");      if(TRG)              R = new TChain("Trigger");
1136      if(S4)  S = new TChain("S4");      if(S4)               S = new TChain("S4");
1137      if(ND)  N = new TChain("NeutronD");      if(ND)               N = new TChain("NeutronD");
1138      if(AC)  A = new TChain("Anticounter");      if(AC)               A = new TChain("Anticounter");
1139      if(ORB) B = new TChain("OrbitalInfo");      if(ORB)              B = new TChain("OrbitalInfo");
1140    
1141         L = new TChain("SelectionList");
1142            
1143      // loop over files and create chains              // loop over files and create chains        
1144      TIter next(fl);      TIter next(fl);
1145      TSystemFile *questo = 0;      TSystemFile *questo = 0;
1146      while ( (questo = (TSystemFile*) next()) ) {      while ( (questo = (TSystemFile*) next()) ) {
1147          TString name =  questo->GetName();          TString name =  questo->GetName();
1148  //              cout << "File: "<< name << endl;          cout << "File: "<< name << endl;
1149          if( CheckLevel2File(name) ){          if( CheckLevel2File(name) ){
1150              if(TRK2||TRK1||TRKh) T->Add(name);              if(TRK2||TRK1||TRKh) T->Add(name);
1151              if(CAL1||CAL2) C->Add(name);              if(CAL1||CAL2)       C->Add(name);
1152              if(TOF) O->Add(name);              if(TOF)              O->Add(name);
1153              if(TRG) R->Add(name);              if(TRG)              R->Add(name);
1154              if(S4)  S->Add(name);              if(S4)               S->Add(name);
1155              if(ND)  N->Add(name);              if(ND)               N->Add(name);
1156              if(AC)  A->Add(name);              if(AC)               A->Add(name);
1157              if(ORB) B->Add(name);              if(ORB)              B->Add(name);
1158                if(SELLI==1)         L->Add(name);
1159          };          };
1160      }      }
1161            
1162      cout << "done chain \n";      cout << "done chain \n";
1163        
1164    //    UInt_t *found=0;
1165    // Tracker
1166        if(T && (TRK2||TRK1||TRKh)) {
1167            if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1168    //      else    T->SetBranchStatus("TrkLevel2",0,found);
1169            if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
1170            if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
1171    //      else    T->SetBranchStatus("TrkLevel1",0,found);
1172            if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
1173            if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
1174    //      else    T->SetBranchStatus("TrkHough",0,found);
1175            if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
1176            if(!Trout)Trout=T;
1177            else Trout->AddFriend("Tracker");
1178        }else{
1179            cout << "Tracker      : missing tree"<<endl;
1180        };
1181      // Calorimeter      // Calorimeter
1182      if(CAL2||CAL1) {                  if(C && (CAL2||CAL1)) {
1183          if(!Tout)Tout=C;          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
1184          else Tout->AddFriend("Calorimeter");  //      else    C->SetBranchStatus("CaloLevel2",0,found);
1185            if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
1186            if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
1187    //      else    C->SetBranchStatus("CaloLevel1",0,found);
1188            if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
1189            if(!Trout)Trout=C;
1190            else Trout->AddFriend("Calorimeter");
1191        }else{
1192            cout << "Calorimeter  : missing tree"<<endl;
1193      };      };
1194      // ToF          // ToF    
1195      if(TOF) {      if(O && TOF) {
1196          if(!Tout)Tout=O;          O->SetBranchAddress("ToFLevel2", GetPointerTo("ToFLevel2"));
1197          else Tout->AddFriend("ToF");          cout << "ToF          : set branch address ToFLevel2"<<endl;
1198            if(!Trout)Trout=O;
1199            else Trout->AddFriend("ToF");
1200        }else{
1201            cout << "ToF         : missing tree"<<endl;
1202      };      };
     // Tracker  
     if(TRK2||TRK1||TRKh){  
         if(!Tout)Tout=T;  
         else Tout->AddFriend("Tracker");  
     }  
1203      // Trigger      // Trigger
1204      if(TRG) {      if(R && TRG) {
1205          if(!Tout)Tout=R;          R->SetBranchAddress("TrigLevel2", GetPointerTo("TrigLevel2"));
1206          else Tout->AddFriend("Trigger");          cout << "Trigger      : set branch address TrigLevel2"<<endl;
1207            if(!Trout)Trout=O;
1208            else Trout->AddFriend("Trigger");
1209        }else{
1210            cout << "Trigger      : missing tree"<<endl;
1211      };      };
1212      // S4      // S4
1213      if(S4) {      if(S && S4) {
1214          if(!Tout)Tout=S;          S->SetBranchAddress("S4Level2", GetPointerTo("S4Level2"));
1215          else Tout->AddFriend("S4");          cout << "S4           : set branch address S4Level2"<<endl;
1216            if(!Trout)Trout=O;
1217            else Trout->AddFriend("S4");
1218        }else{
1219            cout << "S4           : missing tree"<<endl;
1220      };      };
1221      // Neutron Detector      // Neutron Detector
1222      if(ND) {      if(N && ND) {
1223          if(!Tout)Tout=N;          N->SetBranchAddress("NDLevel2", GetPointerTo("NDLevel2"));
1224          else Tout->AddFriend("NeutronD");          cout << "NeutronD     : set branch address NDLevel2"<<endl;
1225            if(!Trout)Trout=O;
1226            else Trout->AddFriend("NeutronD");
1227        }else{
1228            cout << "NeutronD     : missing tree"<<endl;
1229      };      };
1230      // Anticounters      // Anticounters
1231      if(AC) {      if(A && AC) {
1232          if(!Tout)Tout=A;          A->SetBranchAddress("AcLevel2", GetPointerTo("AcLevel2"));
1233          else Tout->AddFriend("Anticounter");          cout << "Anticounter  : set branch address AcLevel2"<<endl;
1234            if(!Trout)Trout=O;
1235            else Trout->AddFriend("Anticounter");
1236        }else{
1237            cout << "Anticounter  : missing tree"<<endl;
1238      };      };
1239      // OrbitalInfo      // Orbital Info
1240      if(ORB) {      if(B && ORB) {
1241          if(!Tout)Tout=B;          B->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1242          else Tout->AddFriend("OrbitalInfo");          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1243            if(!Trout)Trout=O;
1244            else Trout->AddFriend("OrbitalInfo");
1245        }else{
1246            cout << "OrbitalInfo  : missing tree"<<endl;
1247        };
1248    
1249        // Selection List
1250        if(L && SELLI==1) {
1251            cout<<">>> Found selection-list <<<"<<endl;
1252            L->SetBranchAddress("RunEntry",&irun);
1253            cout << "SelectionList: set branch address RunEntry"<<endl;
1254            L->SetBranchAddress("EventEntry",&irunentry);
1255            cout << "SelectionList: set branch address EventEntry"<<endl;
1256            sel_tree = L;
1257    //      if(!Trout)Trout=O;
1258    //      else Trout->AddFriend("SelectionList");
1259        }else{
1260    //      cout << "SelectionList  : missing tree"<<endl;
1261            if(L)L->Delete();
1262      };      };
1263            
1264  //    cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl;  //    cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
1265        
1266  //    if( Tout->GetEntries() )PamLevel2::SetBranchAddress();      pam_tree = Trout;
1267      if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout );  
1268            return Trout;
     return Tout;  
1269  }  }
1270    
1271    
# Line 892  void PamLevel2::SetBranchAddress(TTree * Line 1346  void PamLevel2::SetBranchAddress(TTree *
1346          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1347          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1348      };      };
1349        // SelectionList
1350        if(SELLI==1) {
1351            t->SetBranchAddress("RunEntry", &irun);
1352            cout << "SelectionList: set branch address RunEntry"<<endl;
1353            t->SetBranchAddress("EventEntry", &irunentry);
1354            cout << "SelectionList: set branch address EventEntry"<<endl;
1355        };
1356            
1357  }  }
1358  /**  /**
# Line 913  void PamLevel2::SetBranchAddress(TChain Line 1374  void PamLevel2::SetBranchAddress(TChain
1374    
1375      // Tracker      // Tracker
1376       if(TRK2) {       if(TRK2) {
1377          t->SetBranchAddress("TrkLevel2", this->GetPointerTo("TrkLevel2"));          t->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1378          cout << "Tracker      : set branch address TrkLevel2"<<endl;          cout << "Tracker      : set branch address TrkLevel2"<<endl;
1379      };      };
1380      if(TRK1) {      if(TRK1) {
# Line 965  void PamLevel2::SetBranchAddress(TChain Line 1426  void PamLevel2::SetBranchAddress(TChain
1426          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1427          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1428      };      };
1429        // SelectionList
1430        if(SELLI==1) {
1431            t->SetBranchAddress("RunEntry", &irun);
1432            cout << "SelectionList: set branch address RunEntry"<<endl;
1433            t->SetBranchAddress("EventEntry", &irunentry);
1434            cout << "SelectionList: set branch address EventEntry"<<endl;
1435        };
1436            
1437  }  }
1438    
 // 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;  };  
1439    
1440  //--------------------------------------  //--------------------------------------
1441  //  //
# Line 1016  TChain *PamLevel2::GetRunTree(TList *fl) Line 1460  TChain *PamLevel2::GetRunTree(TList *fl)
1460              R->Add(name);              R->Add(name);
1461          };          };
1462      }      }
1463    
1464        if(R->GetNtrees()){
1465            
1466      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));          R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1467      cout << "Run         : set branch address RunInfo"<<endl;          cout << "Run         : set branch address RunInfo"<<endl;
1468                R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1469            cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1470        }else{
1471            delete R;
1472            R=0;
1473        }
1474    
1475        run_tree = R;
1476    
1477      return R;      return R;
1478            
1479  }  }
# Line 1034  TChain *PamLevel2::GetRunTree(TList *fl) Line 1488  TChain *PamLevel2::GetRunTree(TList *fl)
1488   */   */
1489  TTree *PamLevel2::GetRunTree(TFile *f){  TTree *PamLevel2::GetRunTree(TFile *f){
1490    
1491        cout << "TTree *PamLevel2::GetRunTree(TFile *f) -- obsolte "<<endl;
1492    
1493      TTree *R = (TTree*)f->Get("Run");      TTree *R = (TTree*)f->Get("Run");
1494            
1495      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));      if(R){
1496      cout << "Run         : set branch address RunInfo"<<endl;          R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1497                cout << "Run         : set branch address RunInfo"<<endl;
1498            R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1499            cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1500        }
1501      return R;      return R;
1502            
1503  }  }
1504    /**
1505     * Update the runinfo informations (to be used to have Run infos event by event basis)
1506     * @param run Pointer to the chain/tree which contains run infos
1507     * @return true if a new run has been read, false if it is still the same run
1508     */
1509    Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev){
1510      //
1511      // check if we have already called once GetEntry, if not call it
1512      //
1513      if ( run->GetEntries() <= 0 ) return(false);
1514      //
1515      
1516      Int_t oldrun = irun;
1517      // --------------------------------------
1518      // if it is a full file (not preselected)
1519      // --------------------------------------
1520      if(SELLI==0){
1521          if ( irun < 0 ){
1522              irun = 0;
1523              run->GetEntry(irun);
1524              runfirstentry = 0ULL;
1525              runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS) - 1ULL;
1526              return(true);
1527          };      
1528    
1529          if( !GetOrbitalInfo() ){
1530              cout << "Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev) -- ERROR -- missing OrbitalInfo "<<endl;
1531              return(false);
1532          }
1533          // modificato il controllo sull'aggiornamento del run, per evitare problemi
1534          // dovuti agli eventi annidati (NB! NEVENTS conta anche questi!!)
1535          while ( GetOrbitalInfo()->absTime > GetRunInfo()->RUNTRAILER_TIME && irun < run->GetEntries() ){
1536    //      while ( iev > (runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS-1)) && irun < run->GetEntries() ){
1537              //    printf(" iev %llu %u %llu \n",iev,this->GetRunInfo()->NEVENTS,(ULong64_t)(runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS)));
1538              irun++;
1539              run->GetEntry(irun);
1540              runfirstentry = runlastentry+1ULL;
1541              runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS);
1542          };
1543          
1544          //
1545          if
1546     ( irun == oldrun || irun >= run->GetEntries() ) return(false);
1547          //
1548          //  printf(" iev %llu irun %i nevents %u 1st %llu last %llu \n",iev,irun,this->GetRunInfo()->NEVENTS,(ULong64_t)runfirstentry,(ULong64_t)runlastentry);
1549          //
1550          return(true);    
1551      };
1552      // ----------------------------------------------------
1553      // if it is a preselected file (there is SelectionList)
1554      // NBNB - the event tree MUST be read first
1555      // ----------------------------------------------------
1556      if(SELLI==1){      
1557          sel_tree->GetEntry(iev);
1558          if(irun != oldrun){
1559              run->GetEntry(irun);
1560              return true;
1561          }
1562          return false;
1563      }
1564    
1565      return false;
1566      //
1567    };
1568    /**
1569     * Update the runinfo informations (to be used to have Run infos event by event basis)
1570     * @param run Pointer to the chain/tree which contains run infos
1571     * @return true if a new run has been read, false if it is still the same run
1572     */
1573    Bool_t PamLevel2::UpdateRunInfo(TTree *run, ULong64_t iev){
1574      return(UpdateRunInfo((TChain*)run,iev));
1575    };
1576    
1577  //--------------------------------------  //--------------------------------------
1578  //  //
1579  //  //
# Line 1053  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1584  TTree *PamLevel2::GetRunTree(TFile *f){
1584  */  */
1585  void PamLevel2::SetWhichTrees(TString detlist){  void PamLevel2::SetWhichTrees(TString detlist){
1586                    
1587      if(detlist.Contains("+ALL", TString::kIgnoreCase)){      if(detlist.IsNull() || detlist.Contains("+ALL", TString::kIgnoreCase)){
1588            CAL0 = false;
1589          CAL1 = true;          CAL1 = true;
1590          CAL2 = true;          CAL2 = true;
1591          TRK2 = true;          TRK2 = true;
1592          TRK1 = false;          TRK1 = false;
1593          TRKh = false;          TRKh = false;
1594            TRK0 = false;
1595          TRG = true;          TRG = true;
1596          TOF = true;          TOF = true;
1597            TOF0 = false;
1598          S4  = true;          S4  = true;
1599          ND  = true;          ND  = true;
1600          AC  = true;          AC  = true;
1601          ORB = true;          ORB = true;
1602      }else if( detlist.Contains("-ALL", TString::kIgnoreCase) ){      }else if( detlist.Contains("-ALL", TString::kIgnoreCase) ){
1603            CAL0 = false;
1604          CAL1 = false;          CAL1 = false;
1605          CAL2 = false;          CAL2 = false;
1606          TRK2 = false;          TRK2 = false;
1607          TRK1 = false;          TRK1 = false;
1608          TRKh = false;          TRKh = false;
1609            TRK0 = false;
1610          TRG = false;          TRG = false;
1611          TOF = false;          TOF = false;
1612            TOF0 = false;
1613          S4  = false;          S4  = false;
1614          ND  = false;          ND  = false;
1615          AC  = false;          AC  = false;
# Line 1084  void PamLevel2::SetWhichTrees(TString de Line 1621  void PamLevel2::SetWhichTrees(TString de
1621          if ( detlist.Contains("-CAL1", TString::kIgnoreCase) )CAL1=false;          if ( detlist.Contains("-CAL1", TString::kIgnoreCase) )CAL1=false;
1622          if ( detlist.Contains("+CAL1", TString::kIgnoreCase) )CAL1=true;          if ( detlist.Contains("+CAL1", TString::kIgnoreCase) )CAL1=true;
1623      };      };
1624    
1625        if( detlist.Contains("CAL0", TString::kIgnoreCase) ){
1626            if ( detlist.Contains("-CAL0", TString::kIgnoreCase) )CAL0=false;
1627            if ( detlist.Contains("+CAL0", TString::kIgnoreCase) )CAL0=true;
1628        };
1629                    
1630      if( detlist.Contains("CAL2", TString::kIgnoreCase)){      if( detlist.Contains("CAL2", TString::kIgnoreCase)){
1631          if ( detlist.Contains("-CAL2", TString::kIgnoreCase) )CAL2=false;          if ( detlist.Contains("-CAL2", TString::kIgnoreCase) )CAL2=false;
# Line 1096  void PamLevel2::SetWhichTrees(TString de Line 1638  void PamLevel2::SetWhichTrees(TString de
1638          CAL1=false;          CAL1=false;
1639      }      }
1640  //  -------------------------------------------------------------------------  //  -------------------------------------------------------------------------
1641        if( detlist.Contains("TRK0", TString::kIgnoreCase) ){
1642            if ( detlist.Contains("-TRK0", TString::kIgnoreCase) )TRK0=false;
1643            if ( detlist.Contains("+TRK0", TString::kIgnoreCase) )TRK0=true;
1644        };
1645    
1646      if( detlist.Contains("TRK1", TString::kIgnoreCase) ){      if( detlist.Contains("TRK1", TString::kIgnoreCase) ){
1647          if ( detlist.Contains("-TRK1", TString::kIgnoreCase) )TRK1=false;          if ( detlist.Contains("-TRK1", TString::kIgnoreCase) )TRK1=false;
1648          if ( detlist.Contains("+TRK1", TString::kIgnoreCase) )TRK1=true;          if ( detlist.Contains("+TRK1", TString::kIgnoreCase) )TRK1=true;
# Line 1118  void PamLevel2::SetWhichTrees(TString de Line 1665  void PamLevel2::SetWhichTrees(TString de
1665          TRKh=false;          TRKh=false;
1666      }      }
1667  //  -------------------------------------------------------------------------  //  -------------------------------------------------------------------------
   
   
       
 //     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;  
1668            
1669      if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false;      if( detlist.Contains("-TRG", TString::kIgnoreCase) )TRG = false;
1670      else if( detlist.Contains("+TRG", TString::kIgnoreCase) )TRG = true;      else if( detlist.Contains("+TRG", TString::kIgnoreCase) )TRG = true;
1671            
1672      if( detlist.Contains("-TOF", TString::kIgnoreCase) )TOF = false;      if( detlist.Contains("-TOF", TString::kIgnoreCase) )TOF = false;
1673      else if( detlist.Contains("+TOF", TString::kIgnoreCase) )TOF = true;      else if( detlist.Contains("+TOF", TString::kIgnoreCase) )TOF = true;
1674    
1675        if( detlist.Contains("-TOF0", TString::kIgnoreCase) )TOF0 = false;
1676        else if( detlist.Contains("+TOF0", TString::kIgnoreCase) )TOF0 = true;
1677            
1678      if( detlist.Contains("-S4",  TString::kIgnoreCase) )S4  = false;      if( detlist.Contains("-S4",  TString::kIgnoreCase) )S4  = false;
1679      else if( detlist.Contains("+S4",  TString::kIgnoreCase) )S4  = true;      else if( detlist.Contains("+S4",  TString::kIgnoreCase) )S4  = true;
# Line 1147  void PamLevel2::SetWhichTrees(TString de Line 1686  void PamLevel2::SetWhichTrees(TString de
1686            
1687      if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false;      if( detlist.Contains("-ORB", TString::kIgnoreCase) )ORB = false;
1688      else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true;      else if( detlist.Contains("+ORB", TString::kIgnoreCase) )ORB = true;
1689    
1690    //     cout<< "Set detector list --> ";
1691    //     if(TRK1)cout<<"TRK1 ";
1692    //     if(TRK2)cout<<"TRK2 ";
1693    //     if(TRKh)cout<<"TRKH ";
1694    //     if(CAL1)cout<<"CAL1 ";
1695    //     if(CAL2)cout<<"CAL2 ";
1696    //     if(TOF)cout<<"TOF ";
1697    //     if(TRG)cout<<"TRG ";
1698    //     if(AC)cout<<"AC ";
1699    //     if(ND)cout<<"ND ";
1700    //     if(S4)cout<<"S4 ";
1701    //     if(ORB)cout<<"ORB ";
1702    //     cout << endl;
1703            
1704  };  };
1705    
1706    
1707    /**
1708     * Set tree/branch detector flags from the content of a tree
1709     */
1710    void  PamLevel2::GetWhichTrees(TFile* f){
1711        
1712    
1713    
1714        cout << "void  PamLevel2::GetWhichTrees(TFile* f) --- WARNING!! --- ...potrebbe non funzionare "<<endl;
1715        // -----------
1716        // reset flags
1717        // -----------
1718        CAL1   = false;    
1719        CAL2   = false;    
1720        TRK2   = false;    
1721        TRK1   = false;    
1722        TRKh   = false;    
1723        TRG    = false;    
1724        TOF    = false;    
1725        S4     = false;    
1726        ND     = false;    
1727        AC     = false;    
1728        ORB    = false;    
1729        
1730        RUN    = false;
1731            
1732        cout << "Checking file: "<<f->GetName()<<endl;
1733        if( !f || f->IsZombie() ){
1734            cout << "File: "<< f->GetName() <<" Non valid root file"<< endl;
1735            return;
1736        }
1737    
1738        TList *lk = f->GetListOfKeys();
1739        TIter next(lk);
1740        TKey *key =0;
1741    
1742        Int_t nev = 0;
1743    
1744        while( (key = (TKey*)next()) ){
1745            
1746            if( !strcmp(key->GetName(),"Run"        ) )RUN = true;
1747    
1748            //=========================================================    
1749            if( !strcmp(key->GetName(),"Trigger"    ) ){
1750                TRG = true;
1751                Int_t nevt = ((TTree*)f->Get("Trigger"))->GetEntries();
1752                if( nev && nevt!=nev){
1753                    cout << "File: "<< f->GetName() <<" Trigger tree has "<<nevt<<" events instead of "<<nev<< endl;
1754                    TRG = false;
1755                }else nev=nevt;
1756            }
1757            //=========================================================    
1758            if( !strcmp(key->GetName(),"ToF"        ) ){
1759                TOF = true;
1760                Int_t nevt = ((TTree*)f->Get("ToF"))->GetEntries();
1761                if( nev && nevt!=nev){
1762                    cout << "File: "<< f->GetName() <<"     ToF tree has "<<nevt<<" events instead of "<<nev<< endl;
1763                    TOF = false;
1764                }else nev=nevt;
1765            }
1766            //=========================================================  
1767            if( !strcmp(key->GetName(),"S4"         ) ){
1768                S4 = true;
1769                Int_t nevt = ((TTree*)f->Get("S4"))->GetEntries();
1770                if( nev && nevt!=nev){
1771                    cout << "File: "<< f->GetName() <<"      S4 tree has "<<nevt<<" events instead of "<<nev<< endl;
1772                    S4 = false;
1773                }else nev=nevt;
1774            }
1775            //=========================================================  
1776    
1777            if( !strcmp(key->GetName(),"NeutronD"   ) ){
1778                ND = true;
1779                Int_t nevt = ((TTree*)f->Get("NeutronD"))->GetEntries();
1780                if( nev && nevt!=nev){
1781                    cout << "File: "<< f->GetName() <<"NeutronD tree has "<<nevt<<" events instead of "<<nev<< endl;
1782                    ND =false;
1783                }else nev=nevt;
1784            }      
1785            //=========================================================  
1786            if( !strcmp(key->GetName(),"Anticounter") ){
1787                AC = true;
1788                Int_t nevt = ((TTree*)f->Get("Anticounter"))->GetEntries();
1789                if( nev && nevt!=nev){
1790                    cout << "File: "<< f->GetName() <<" Anticounter tree has "<<nevt<<" events instead of "<<nev<< endl;
1791                    AC =false;
1792                }else nev=nevt;
1793            }
1794            //=========================================================  
1795            if( !strcmp(key->GetName(),"OrbitalInfo") ){
1796                ORB = true;
1797                Int_t nevt = ((TTree*)f->Get("OrbitalInfo"))->GetEntries();
1798                if( nev && nevt!=nev){
1799                    cout << "File: "<< f->GetName() <<" OrbitalInfo tree has "<<nevt<<" events instead of "<<nev<< endl;
1800                    ORB = false;
1801                }else nev=nevt;
1802            }
1803            //=========================================================  
1804            if( !strcmp(key->GetName(),"Tracker"    ) ){
1805                TTree *T = (TTree*)f->Get("Tracker");
1806                for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){
1807                    TString name = T->GetListOfBranches()->At(i)->GetName();
1808                    if( !name.CompareTo("TrkLevel1") )TRK1=true;
1809                    if( !name.CompareTo("TrkLevel2") )TRK2=true;
1810                    if( !name.CompareTo("TrkHough") )TRKh=true;
1811                };      
1812                Int_t nevt = T->GetEntries();
1813                if( nev && nevt!=nev){
1814                    cout << "File: "<< f->GetName() <<" Tracker tree has "<<nevt<<" events instead of "<<nev<< endl;
1815                    TRK1 = false;
1816                    TRK2 = false;
1817                    TRKh = false;
1818                }else nev=nevt;
1819                T->Delete();
1820            };
1821            //=========================================================  
1822            if( !strcmp(key->GetName(),"Calorimeter"    ) ){
1823                TTree *T = (TTree*)f->Get("Calorimeter");
1824                for(Int_t i=0; i<T->GetListOfBranches()->GetEntries(); i++){
1825                    TString name = T->GetListOfBranches()->At(i)->GetName();
1826                    if( !name.CompareTo("CaloLevel1") )CAL1=true;
1827                    if( !name.CompareTo("CaloLevel2") )CAL2=true;
1828                };    
1829                Int_t nevt = T->GetEntries();
1830                if( nev && nevt!=nev){
1831                    cout << "File: "<< f->GetName() <<"  Calorimeter tree has "<<nevt<<" events instead of "<<nev<< endl;
1832                    CAL1 = false;
1833                    CAL2 = false;
1834                }else nev=nevt;
1835                T->Delete();
1836            };      
1837    
1838        };
1839        
1840        delete lk;
1841    
1842    //     cout<< "Get detector list from input file --> ";
1843    //     if(TRK1)cout<<"TRK1 ";
1844    //     if(TRK2)cout<<"TRK2 ";
1845    //     if(TRKh)cout<<"TRKH ";
1846    //     if(CAL1)cout<<"CAL1 ";
1847    //     if(CAL2)cout<<"CAL2 ";
1848    //     if(TOF)cout<<"TOF ";
1849    //     if(TRG)cout<<"TRG ";
1850    //     if(AC)cout<<"AC ";
1851    //     if(ND)cout<<"ND ";
1852    //     if(S4)cout<<"S4 ";
1853    //     if(ORB)cout<<"ORB ";
1854    //     cout << endl;
1855          
1856        return ;
1857            
1858    };
1859    
1860    
1861  //--------------------------------------  //--------------------------------------
1862  //  //
1863  //  //
# Line 1174  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1883  Bool_t  PamLevel2::CheckLevel2File(TStri
1883            
1884      Bool_t RUN__ok    = false;      Bool_t RUN__ok    = false;
1885            
1886            Bool_t SELLI__ok  = false;
1887    
1888      cout << "Checking file: "<<name<<endl;      cout << "Checking file: "<<name<<endl;
1889      TFile *f = new TFile(name.Data());      TFile *f = new TFile(name.Data());
# Line 1192  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1901  Bool_t  PamLevel2::CheckLevel2File(TStri
1901      while( (key = (TKey*)next()) ){      while( (key = (TKey*)next()) ){
1902                    
1903  //      cout << key->GetName() << endl;  //      cout << key->GetName() << endl;
1904  //      cout << " Get tree: " << f->Get(key->GetName())<<endl;  //      cout << key->GetName() << ""<<key->GetClassName()<<endl;
1905    //              cout << " Get tree: " << f->Get(key->GetName())<<endl;
1906  //      nev_previous = nev;  //      nev_previous = nev;
1907  //      cout << " n.entries  "<< nev <<endl;  //      cout << " n.entries  "<< nev <<endl;
1908  //      if( key->GetClassName()=="TTree" && nev_previous && nev != nev_previous ){  //      if( key->GetClassName()=="TTree" && nev_previous && nev != nev_previous ){
# Line 1201  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1911  Bool_t  PamLevel2::CheckLevel2File(TStri
1911  //          return false;  //          return false;
1912  //      };  //      };
1913    
1914    
1915          if( !strcmp(key->GetName(),"Run"        ) )RUN__ok = true;                if( !strcmp(key->GetName(),"Run"        ) )RUN__ok = true;      
1916    
1917          //=========================================================              //=========================================================    
1918            if( !strcmp(key->GetName(),"SelectionList"    ) ){
1919                SELLI__ok = true;
1920                if(SELLI==1){
1921                    Int_t nevt = ((TTree*)f->Get("SelectionList"))->GetEntries();
1922                    if( nev && nevt!=nev){
1923                        cout << "File: "<< f->GetName() <<" discarded ---- SelectionList tree has "<<nevt<<" events instead of "<<nev<< endl;
1924                        return false;
1925                    }
1926                    nev=nevt;
1927                }
1928            }
1929    
1930            //=========================================================    
1931          if( !strcmp(key->GetName(),"Trigger"    ) ){          if( !strcmp(key->GetName(),"Trigger"    ) ){
1932              TRG__ok = true;              TRG__ok = true;
1933              if(TRG){              if(TRG){
# Line 1315  Bool_t  PamLevel2::CheckLevel2File(TStri Line 2039  Bool_t  PamLevel2::CheckLevel2File(TStri
2039          };                };      
2040    
2041      };      };
2042    
2043        if( SELLI==-1 )SELLI = (Int_t)SELLI__ok;
2044        if( SELLI==0 && SELLI__ok ){
2045            cout << "File: "<< f->GetName() <<" discarded ---- found SelectionList (it is not a full-event file)" << endl;
2046            return false;  
2047        }
2048        if( SELLI==1 && !SELLI__ok ){
2049            cout << "File: "<< f->GetName() <<" discarded ---- SelectionList missing" << endl;
2050            return false;
2051        }
2052            
2053    //    cout << "SELLI "<<SELLI<<endl;
2054    
2055    //     cout<< "CheckLevel2File(TString): detector list --> ";
2056    //     if(TRK1__ok)cout<<"TRK1 ";
2057    //     if(TRK2__ok)cout<<"TRK2 ";
2058    //     if(TRKh__ok)cout<<"TRKH ";
2059    //     if(CAL1__ok)cout<<"CAL1 ";
2060    //     if(CAL2__ok)cout<<"CAL2 ";
2061    //     if(TOF__ok)cout<<"TOF ";
2062    //     if(TRG__ok)cout<<"TRG ";
2063    //     if(AC__ok)cout<<"AC ";
2064    //     if(ND__ok)cout<<"ND ";
2065    //     if(S4__ok)cout<<"S4 ";
2066    //     if(ORB__ok)cout<<"ORB ";
2067    //     cout << endl;
2068    
2069    
2070        if(TRK2 && TRK1__ok)TRK1=1;
2071    // ----------------------------------------------------------------------------
2072    // NOTA
2073    // se c'e` il level1, lo devo necessarimente leggere.
2074    // infatti (non ho capito perche`) i cluster vengono letti e allocati in memoria
2075    // comunque, ma non vengono disallocati da PamLevel2::Clear()
2076    // ----------------------------------------------------------------------------
2077    
2078    
2079      if(!RUN__ok) {      if(!RUN__ok) {
# Line 1360  Bool_t  PamLevel2::CheckLevel2File(TStri Line 2118  Bool_t  PamLevel2::CheckLevel2File(TStri
2118          return false;          return false;
2119      };      };
2120    
2121          if(ND && !ND__ok){      if(ND && !ND__ok){
2122              cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl;          cout << "File: "<< f->GetName() <<" discarded ---- Missing ND tree"<< endl;
2123              return false;          return false;
2124          };      };
2125          if(TRG && !TRG__ok){      if(TRG && !TRG__ok){
2126              cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl;          cout << "File: "<< f->GetName() <<" discarded ---- Missing Trigger tree"<< endl;
2127              return false;          return false;
2128          };      };
2129        
   
   
   
2130    
2131      lk->Delete();  //    lk->Delete();
2132    //    delete lk;
2133      f->Close();      f->Close();
2134    
2135    //     cout<< "CheckLevel2File(TString): detector list --> ";
2136    //     if(TRK1)cout<<"TRK1 ";
2137    //     if(TRK2)cout<<"TRK2 ";
2138    //     if(TRKh)cout<<"TRKH ";
2139    //     if(CAL1)cout<<"CAL1 ";
2140    //     if(CAL2)cout<<"CAL2 ";
2141    //     if(TOF)cout<<"TOF ";
2142    //     if(TRG)cout<<"TRG ";
2143    //     if(AC)cout<<"AC ";
2144    //     if(ND)cout<<"ND ";
2145    //     if(S4)cout<<"S4 ";
2146    //     if(ORB)cout<<"ORB ";
2147    //     cout << endl;
2148                
2149      return true;      return true;
2150                    
# Line 1384  Bool_t  PamLevel2::CheckLevel2File(TStri Line 2154  Bool_t  PamLevel2::CheckLevel2File(TStri
2154  /**  /**
2155   * Create clone-trees   * Create clone-trees
2156   */   */
2157  void PamLevel2::CreateCloneTrees( TChain *fChain ){  void PamLevel2::CreateCloneTrees0( TChain *fChain, TFile *ofile ){
2158    
2159      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2160      cout << "Create clones of PAMELA trees "<<endl;      cout << "Create clones of PAMELA trees "<<endl;
2161            
2162      Int_t i=0;      Int_t i=0;
2163      tree_clone[i] = fChain->GetTree()->CloneTree(0);      pam_tree_clone[i] = fChain->GetTree()->CloneTree(0);
2164      TString name =  tree_clone[i]->GetName();      TString name =  pam_tree_clone[i]->GetName();
2165      name.Append("_clone");      name.Append("_clone");
2166  //    tree_clone[i]->SetName(name.Data());  //    pam_tree_clone[i]->SetName(name.Data());
2167      cout << tree_clone[i]->GetName() <<endl;      cout << pam_tree_clone[i]->GetName() <<endl;
2168      i++;      i++;
2169    
2170      TList *li = fChain->GetListOfFriends();      TList *li = fChain->GetListOfFriends();
2171      TIter next(li);      TIter next(li);
2172      TFriendElement* T_friend=0;      TFriendElement* T_friend=0;
2173        ofile->cd();
2174      while( (T_friend = (TFriendElement*)next()) ){      while( (T_friend = (TFriendElement*)next()) ){
2175  //      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;
2176  //      cout<<T_friend->GetTree()->GetName()<< endl;  //      cout<<T_friend->GetTree()->GetName()<< endl;
2177          tree_clone[i] = T_friend->GetTree()->CloneTree(0);                pam_tree_clone[i] = T_friend->GetTree()->CloneTree(0);
2178          name =  tree_clone[i]->GetName();          pam_tree_clone[i]->SetAutoSave(1000000);
2179            name =  pam_tree_clone[i]->GetName();
2180          name.Append("_clone");          name.Append("_clone");
2181  //      tree_clone[i]->SetName(name.Data());  //      pam_tree_clone[i]->SetName(name.Data());
2182          cout << tree_clone[i]->GetName() << endl;          cout << pam_tree_clone[i]->GetName() << endl;
2183          i++;          i++;
2184      }      }
2185        
2186        delete li;
2187    
2188      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2189    
2190  }  }
2191    
2192  // /**  /**
2193  //  * Create clone-trees   * Create clone-trees
2194  //  */   */
2195  // void PamLevel2::CreateCloneTrees(){  void PamLevel2::CreateCloneTrees(TFile *ofile){
2196    
2197  //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      ofile->cd();
2198  //     cout << "Create new PAMELA trees "<<endl;  
2199        cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2200        cout << "Create new PAMELA trees "<<endl;
2201                    
 //     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;  
 //     }  
2202    
2203  //     // Calorimeter      run_tree_clone = new TTree("Run","PAMELA Level2 data from the GL_RUN table ");
2204  //     if(CAL1||CAL2){      run_tree_clone->Branch("RunInfo","GL_RUN",GetPointerTo("RunInfo"));
2205  //      C_clone = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data ");          cout << "Run          : branch RunInfo"<<endl;
2206  //      if(CAL1) {      run_tree_clone->Branch("SoftInfo","SoftInfo",GetPointerTo("SoftInfo"));
2207  //          calo1_clone = new CaloLevel1();      cout << "Run          : branch SoftInfo"<<endl;
2208  //          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;  
 //      };  
 //     }      
2209    
 //     // ToF      
 //     if(TOF) {  
 //      O_clone = new TTree("ToF_clone","PAMELA ToF level2 data ");      
 //      tof_clone = new ToFLevel2();  
 //      O_clone->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2"));  
 //      cout << "ToF          : branch ToFLevel2"<<endl;  
 //     };  
 //     // Trigger  
 //     if(TRG) {  
 //      R_clone = new TTree("Trigger_clone","PAMELA trigger level2 data ");      
 //      trig_clone = new TrigLevel2();  
 //      R_clone->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2"));  
 //      cout << "Trigger      : branch TrigLevel2"<<endl;  
 //     };  
 //     // S4  
 //     if(S4) {  
 //      S_clone = new TTree("S4_clone","PAMELA S4 level2 data ");        
 //      s4_clone = new S4Level2();  
 //      S_clone->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2"));  
 //      cout << "S4           : branch S4Level2"<<endl;  
 //     };  
 //     // Neutron Detector  
 //     if(ND) {  
 //      N_clone = new TTree("NeutronD_clone","PAMELA neutron detector level2 data ");    
 //      nd_clone = new NDLevel2();  
 //      N_clone->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2"));  
 //      cout << "NeutronD     : branch NDLevel2"<<endl;  
 //     };  
 //     // Anticounters  
 //     if(AC) {  
 //      A_clone = new TTree("Anticounter_clone","PAMELA anticounter detector level2 data ");      
 //      ac_clone = new AcLevel2();  
 //      A_clone->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2"));  
 //      cout << "Anticounter  : branch AcLevel2"<<endl;  
 //     };  
 //     // OrbitalInfo  
 //     if(ORB) {  
 //      B_clone = new TTree("OrbitalInfo_clone","PAMELA oribital info  ");        
 //      orb_clone = new OrbitalInfo();  
 //      B_clone->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo"));  
 //      cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;  
 //     };  
 //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;  
2210    
2211  // }      sel_tree_clone = new TTree("SelectionList","List of selected events ");
2212        sel_tree_clone->Branch("RunEntry",&irun,"runentry/I");
2213        sel_tree_clone->Branch("EventEntry",&irunentry,"eventry/I");
2214        
2215    
2216  // /**      Int_t i=0;
2217  //  * Create clone-trees      if(TRK1||TRK2||TRKh){
2218  //  */          pam_tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data ");
2219  // void PamLevel2::CreateCloneTrees(){          if(TRK1) {
2220                pam_tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1"));
2221                pam_tree_clone[i]->BranchRef();
2222                cout << "Tracker      : branch TrkLevel1"<<endl;
2223                cout << "CreateCloneTrees " << GetTrkLevel1()<<endl;
2224            };
2225            if(TRK2) {
2226                pam_tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2"));
2227                cout << "Tracker      : branch TrkLevel2"<<endl;        
2228            };
2229            if(TRKh) {
2230                pam_tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough"));
2231                cout << "Tracker      : branch TrkHough"<<endl;
2232            };
2233            i++;
2234        }
2235    
2236  //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      // Calorimeter
2237  //     cout << "Create new PAMELA trees "<<endl;      if(CAL1||CAL2){
2238                    pam_tree_clone[i] = new TTree("Calorimeter","PAMELA calorimeter level2 data ");
2239  //     if(TRK1||TRK2||TRKh){          if(CAL1) {
2240  //      T_clone = new TTree("Tracker_clone","PAMELA tracker level2 data ");                  pam_tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1"));
2241  //      if(TRK1) {              cout << "Calorimeter  : branch CaloLevel1"<<endl;
2242  //          trk1_clone = new TrkLevel1();          };
2243  //          T_clone->Branch("TrkLevel1","TrkLevel1", &trk1_clone);          if(CAL2) {
2244  //          T_clone->BranchRef();              pam_tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2"));
2245  //          cout << "Tracker      : branch TrkLevel1"<<endl;              cout << "Calorimeter  : branch CaloLevel2"<<endl;
2246  //      };          };
2247  //      if(TRK2) {          i++;
2248  //          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;  
 //     }  
2249    
2250  //     // Calorimeter      // ToF    
2251  //     if(CAL1||CAL2){      if(TOF) {
2252  //      C_clone = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data ");              pam_tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data ");
2253  //      if(CAL1) {          pam_tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2"));
2254  //          calo1_clone = new CaloLevel1();          cout << "ToF          : branch ToFLevel2"<<endl;
2255  //          C_clone->Branch("CaloLevel1", "CaloLevel1", &calo1_clone);          i++;
2256  //          cout << "Calorimeter  : branch CaloLevel1"<<endl;      };
2257  //      };      // Trigger
2258  //      if(CAL2) {      if(TRG) {
2259  //          calo2_clone = new CaloLevel2();          pam_tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data ");
2260  //          C_clone->Branch("CaloLevel2","CaloLevel2", &calo2_clone);          pam_tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2"));
2261  //          cout << "Calorimeter  : branch CaloLevel2"<<endl;          cout << "Trigger      : branch TrigLevel2"<<endl;
2262  //      };          i++;
2263  //     }          };
2264        // S4
2265        if(S4) {
2266            pam_tree_clone[i] = new TTree("S4","PAMELA S4 level2 data ");  
2267            pam_tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2"));
2268            cout << "S4           : branch S4Level2"<<endl;
2269            i++;
2270        };
2271        // Neutron Detector
2272        if(ND) {
2273            pam_tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data ");      
2274            pam_tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2"));
2275            cout << "NeutronD     : branch NDLevel2"<<endl;
2276            i++;
2277        };
2278        // Anticounters
2279        if(AC) {
2280            pam_tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data ");        
2281            pam_tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2"));
2282            cout << "Anticounter  : branch AcLevel2"<<endl;
2283            i++;
2284        };
2285        // OrbitalInfo
2286        if(ORB) {
2287            pam_tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info  ");  
2288            pam_tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo"));
2289            cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;
2290            i++;
2291        };
2292        cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2293    
2294  //     // 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;  
2295    
 // }  
2296    
2297  /**  /**
2298   * Fill tree (created with CreateCloneTrees)   * Fill tree (created with CreateCloneTrees)
# Line 1601  void PamLevel2::CreateCloneTrees( TChain Line 2300  void PamLevel2::CreateCloneTrees( TChain
2300   */   */
2301  //void PamLevel2::FillNewPamTree(TTree *T){  //void PamLevel2::FillNewPamTree(TTree *T){
2302  void PamLevel2::FillCloneTrees(){  void PamLevel2::FillCloneTrees(){
2303                
2304  // //  --------------------------------------  //    cout << "PamLevel2::FillCloneTrees()" << irunentry << endl;
 // //  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();  
2305    
2306      for(Int_t i=0; i<8; i++){      for(Int_t i=0; i<NCLONES; i++){
2307          if(tree_clone[i])tree_clone[i]->Fill();          if(pam_tree_clone[i])pam_tree_clone[i]->Fill();
2308      }      }
2309        if(sel_tree_clone)sel_tree_clone->Fill();    
2310    
2311  }  }
2312    
2313    
2314  TTree* PamLevel2::GetCloneTree(TString name){  TTree* PamLevel2::GetCloneTree(TString name){
2315    
2316      for(Int_t i=0; i<8; i++){        for(Int_t i=0; i<NCLONES; i++){    
2317          if(tree_clone[i]){          if(pam_tree_clone[i]){
2318              TString na = tree_clone[i]->GetName();              TString na = pam_tree_clone[i]->GetName();
2319              if(!name.CompareTo(na))return tree_clone[i];              if(!name.CompareTo(na))return pam_tree_clone[i];
2320          };          };
2321      }      }
2322        if(run_tree_clone){
2323            TString na = run_tree_clone->GetName();
2324            if(!name.CompareTo(na))return run_tree_clone;
2325        }
2326        if(sel_tree_clone){
2327            TString na = sel_tree_clone->GetName();
2328            if(!name.CompareTo(na))return sel_tree_clone;
2329        }
2330      return NULL;      return NULL;
2331    
2332  }  }
2333  void PamLevel2::WriteCloneTrees(){  void PamLevel2::WriteCloneTrees(){
2334      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2335      cout << "Write clones of PAMELA trees "<<endl;      cout << "Write clones of PAMELA trees "<<endl;
2336      for(Int_t i=0; i<8; i++){        cout << run_tree_clone->GetName()<<endl;    
2337          if(tree_clone[i]){      run_tree_clone->Write();
2338              cout << tree_clone[i]->GetName()<<endl;      cout << sel_tree_clone->GetName()<<endl;    
2339              tree_clone[i]->Write();      sel_tree_clone->Write();
2340        for(Int_t i=0; i<NCLONES; i++){    
2341            if(pam_tree_clone[i]){
2342                cout << pam_tree_clone[i]->GetName()<<endl;
2343                pam_tree_clone[i]->Write();
2344          };          };
2345      }      }
2346      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
# Line 1669  void PamLevel2::WriteCloneTrees(){ Line 2348  void PamLevel2::WriteCloneTrees(){
2348  }  }
2349    
2350  /**  /**
2351   * Create a new (empty) Pamela trees   * Method to get level2-trees entry, the corresponding run entry and (if required) the level0 entry.
2352   */   */
2353  // TTree* PamLevel2::GetNewPamTree(){  Int_t PamLevel2::GetEntry(Int_t iee){
2354        
2355        if(!pam_tree){
2356            cout << " Int_t PamLevel2::GetEntry(Int_t) -- ERROR -- level2 trees not loeaded"<<endl;
2357            return 0;
2358        }
2359        if(!run_tree ){
2360            cout << " Int_t PamLevel2::GetEntry(Int_t) -- ERROR -- run tree not loeaded"<<endl;
2361            return 0;
2362        }
2363    
2364  //     if(tree_clone)return tree_clone;      Int_t ii=0;
2365        //-------------------------------
2366        ii = iee;
2367        if( !pam_tree->GetEntry(ii) ) return 0;
2368        //-------------------------------
2369        ii = iee;
2370        Bool_t UPDATED = UpdateRunInfo(run_tree,ii);
2371        if(SELLI==0)irunentry = iee-runfirstentry;
2372        if(UPDATED && run_tree_clone)run_tree_clone->Fill();
2373    
2374  //     TTree *Tout = 0;  //    cout << "PamLevel2::GetEntry("<<iee<<") "<<irun<<" "<<runfirstentry<<" "<<irunentry<<endl;
2375        
2376  //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      if( TRK0 || CAL0 || TOF0 )GetYodaEntry( );
2377  //     cout << "Create new PAMELA trees "<<endl;  
2378            return 1;
2379    
2380    }
2381    
2382    /**
2383     * Method to retrieve the level0 tree (YODA tree) that contains the current event.
2384     * Given the run ID (...), if needed it query the DB and load the proper file.
2385     * @return Pointer to the tree
2386     */
2387    
2388    
2389    TTree* PamLevel2::GetYodaTree( ){
2390        
2391    //    cout << "TTree* PamLevel2::GetYodaTree( )"<<endl;
2392        //===================================
2393        // check if iroot has changed
2394        //===================================
2395        if( irun<0 ){
2396            cout << "PamLevel2::GetYodaTree() -- ERROR "<<endl;
2397            cout << "In order to use this method you have to load the RunInfo tree "<<endl;
2398            cout << "with the method TChain* PamLevel2::GetRunTree(TString, TString)"<<endl;
2399    //      cout << " - read the event with PamLevel2::GetEntry(Int_t)"<<endl;
2400            return NULL;
2401        }    
2402        Int_t irootnew = run_obj->ID_ROOT_L0;
2403    //     cout << "iroot    "<<iroot<<endl;
2404    //     cout << "irootnew "<<irootnew<<endl;
2405    
2406        //===================================
2407        // load the level0 file
2408        // (if not already loaded)
2409        //===================================
2410        if( iroot != irootnew || !l0_tree){
2411            iroot = irootnew;
2412            //===================================
2413            // open the DB connection
2414            // (if not already opened)
2415            //===================================
2416            if(!dbc || (dbc && !dbc->IsConnected())){
2417                cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl;
2418                cout<<"Connecting to DB"<<endl;
2419                cout<<"HOST "<<host<<endl;
2420                cout<<"USER "<<user<<endl;
2421                cout<<"PSW  "<<psw<<endl;
2422                dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
2423                if( !dbc )return NULL;
2424                if( !dbc->IsConnected() )return NULL;    
2425    //          cout<<"...done"<<endl;
2426                cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl;
2427            }else{
2428    //          cout<<"DB already connected"<<endl;
2429            }
2430            GL_ROOT glroot = GL_ROOT();
2431            if( glroot.Query_GL_ROOT(iroot,dbc) ){
2432                cout << "TTree* PamLevel2::GetYodaTree( ) -- ERROR -- level0 file iroot = "<<iroot<< " does not exists"<<endl;
2433                return NULL;
2434            };
2435            TString filename = glroot.PATH + glroot.NAME;
2436            if(l0_file){
2437                l0_file->Close();
2438                l0_file->Delete();
2439            }
2440            cout << "Opening LEVEL0 file: "<< filename << endl;
2441            FileStat_t t;
2442            if( gSystem->GetPathInfo(filename.Data(),t) ){
2443                cout << " PamLevel2::GetYodaTree() -- ERROR opening file "<<endl;
2444                return NULL;
2445            }
2446            l0_file = new TFile(filename);
2447            if( !l0_file )return NULL;
2448            l0_tree = (TTree*)l0_file->Get("Physics");
2449            if(!h0_obj)h0_obj = new EventHeader();
2450            l0_tree->SetBranchAddress("Header" ,&h0_obj);
2451            //---------------------------------------------------
2452            // TRACKER:
2453            if(TRK0){
2454                if(!trk0_obj){
2455                    trk0_obj = new TrkLevel0();
2456                    trk0_obj->Set();
2457                };
2458                l0_tree->SetBranchAddress("Tracker" ,trk0_obj->GetPointerToTrackerEvent());
2459            }
2460            //---------------------------------------------------
2461            // CALORIMETER:
2462            if(CAL0){
2463                cout << "PamLevel2::GetYodaTree() --- level0 calorimeter not implemented "<<endl;
2464            }
2465            //---------------------------------------------------
2466            // TOF:
2467            if(TOF0){
2468                cout << "PamLevel2::GetYodaTree() --- level0 TOF not implemented "<<endl;
2469            }
2470    
2471        };
2472    
2473        if(!dbc || (dbc && !dbc->IsConnected())){
2474            cout << " TTree* PamLevel2::GetYodaTree( ) -- no DB connected... hai fatto qualche cazzata "<<endl;
2475        }
2476    
2477        TrkParams::SetCalib(run_obj,dbc);
2478    
2479    
2480    //    cout << l0_tree << endl;
2481            
2482  //     if(TRK1||TRK2||TRKh){      return l0_tree;
 //      TTree *T = new TTree("Tracker_clone","PAMELA tracker level2 data ");      
 //      if(TRK1) {  
 //          trk1_clone = new TrkLevel1();  
 //          T->Branch("TrkLevel1","TrkLevel1", &trk1_clone);  
 //          T->BranchRef();  
 //          cout << "Tracker      : branch TrkLevel1"<<endl;  
 //      };  
 //      if(TRK2) {  
 //          trk2_clone = new TrkLevel2();  
 //          T->Branch("TrkLevel2", "TrkLevel2",&trk2_clone);  
 //          cout << "Tracker      : branch TrkLevel2"<<endl;          
 //      };  
 //      if(TRKh) {  
 //          trkh_clone = new TrkHough();  
 //          T->Branch("TrkHough","TrkHough", &trkh_clone);  
 //          cout << "Tracker      : branch TrkHough"<<endl;  
 //      };  
 //      if(!Tout)Tout=T;  
 //      else Tout->AddFriend("Tracker_clone");  
 //     }  
   
 //     // Calorimeter  
 //     if(CAL1||CAL2){  
 //      TTree *C = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data ");      
 //      if(CAL1) {  
 //          calo1_clone = new CaloLevel1();  
 //          C->Branch("CaloLevel1", "CaloLevel1", &calo1_clone);  
 //          cout << "Calorimeter  : branch CaloLevel1"<<endl;  
 //      };  
 //      if(CAL2) {  
 //          calo2_clone = new CaloLevel2();  
 //          C->Branch("CaloLevel2","CaloLevel2", &calo2_clone);  
 //          cout << "Calorimeter  : branch CaloLevel2"<<endl;  
 //      };  
 //      if(!Tout)Tout=C;  
 //      else Tout->AddFriend("Calorimeter_clone");  
 //     }      
   
 //     // ToF      
 //     if(TOF) {  
 //      TTree *O = new TTree("ToF_clone","PAMELA ToF level2 data ");      
 //      tof_clone = new ToFLevel2();  
 //      O->Branch("ToFLevel2","ToFLevel2", &tof_clone);  
 //      cout << "ToF          : branch ToFLevel2"<<endl;  
 //      if(!Tout)Tout=O;  
 //      else Tout->AddFriend("ToF_clone");  
 //     };  
 //     // Trigger  
 //     if(TRG) {  
 //      TTree *R = new TTree("Trigger_clone","PAMELA trigger level2 data ");      
 //      trig_clone = new TrigLevel2();  
 //      R->Branch("TrigLevel2","TrigLevel2", &trig_clone);  
 //      cout << "Trigger      : branch TrigLevel2"<<endl;  
 //      if(!Tout)Tout=R;  
 //      else Tout->AddFriend("Trigger_clone");  
 //     };  
 //     // S4  
 //     if(S4) {  
 //      TTree *S = new TTree("S4_clone","PAMELA S4 level2 data ");        
 //      s4_clone = new S4Level2();  
 //      S->Branch("S4Level2","S4Level2", &s4_clone);  
 //      cout << "S4           : branch S4Level2"<<endl;  
 //      if(!Tout)Tout=S;  
 //      else Tout->AddFriend("S4_clone");  
 //     };  
 //     // Neutron Detector  
 //     if(ND) {  
 //      TTree *N = new TTree("NeutronD_clone","PAMELA neutron detector level2 data ");    
 //      nd_clone = new NDLevel2();  
 //      N->Branch("NDLevel2","NDLevel2", &nd_clone);  
 //      cout << "NeutronD     : branch NDLevel2"<<endl;  
 //      if(!Tout)Tout=N;  
 //      else Tout->AddFriend("NeutronD_clone");  
 //     };  
 //     // Anticounters  
 //     if(AC) {  
 //      TTree *A = new TTree("Anticounter_clone","PAMELA anticounter detector level2 data ");    
 //      ac_clone = new AcLevel2();  
 //      A->Branch("AcLevel2","AcLevel2", &ac_clone);  
 //      cout << "Anticounter  : branch AcLevel2"<<endl;  
 //      if(!Tout)Tout=A;  
 //      else Tout->AddFriend("Anticounter_clone");  
 //     };  
 //     // OrbitalInfo  
 //     if(ORB) {  
 //      TTree *B = new TTree("OrbitalInfo_clone","PAMELA oribital info  ");      
 //      orb_clone = new OrbitalInfo();  
 //      B->Branch("OrbitalInfo","OrbitalInfo", &orb_clone);  
 //      cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;  
 //      if(!Tout)Tout=B;  
 //      else Tout->AddFriend("OrbitalInfo_clone");  
 //     };  
 //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;  
2483    
2484  //     tree_clone = Tout;  }
 //     tree_clone->SetDirectory(0);  
2485    
2486  //     return Tout;  /**
2487  // }   * Method to retrieve the level0 tree (YODA tree) that contains the current event.
2488  // /**   */
2489  //  * Fill a tree (created with GetNewPamTree)  Int_t PamLevel2::GetYodaEntry(){
 //  *  
 //  */  
 // //void PamLevel2::FillNewPamTree(TTree *T){  
 // void PamLevel2::FillNewPamTree(){  
           
2490    
2491  //     if(trk1_clone)  *trk1_clone = *trk1_obj;  //    cout << "Int_t PamLevel2::GetYodaEntry()"<<endl;
2492  //     if(trk2_clone){      if(!GetYodaTree())return 0;
2493  //      trk2_clone->Clear();      
2494  //      trk2_obj->Copy(*trk2_clone);      //  ---------------------------------
2495  // //   *trk2_clone = *trk2_obj;      //  if file is NOT a preselected file
2496  //     }      //  ---------------------------------
2497  //     if(trkh_clone)  *trkh_clone = *trkh_obj;      Int_t quellagiusta = irunentry + run_obj->EV_FROM;
2498  //     if(calo1_clone){  
2499  // //   *calo1_clone = *calo1_obj;  //     cout << " irun "<< irun << " irunentry "<< irunentry<<" run_obj->EV_FROM "<<run_obj->EV_FROM <<" quella giusta "<<quellagiusta << endl;
2500  //      calo1_clone->Clear();  //     cout << " iroot "<<iroot<<" run_obj->ID_ROOT_L0 "<<run_obj->ID_ROOT_L0<<endl;
2501  //      calo1_obj->Copy(*calo1_clone);  //     cout << " time "<< GetOrbitalInfo()->absTime << endl;
2502  //     }  //     cout << " trk_calib_used "<<run_obj->TRK_CALIB_USED<< endl;
2503  //     if(calo2_clone){      
2504  // //   *calo2_clone = *calo2_obj;      if( !GetOrbitalInfo() ){
2505  //      calo2_clone->Clear();          cout << "Int_t PamLevel2::GetYodaEntry() -- ERROR -- missing OrbitalInfo "<<endl;
2506  //      calo2_obj->Copy(*calo2_clone);          return 0;
2507  //     }      }
 //     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;  
   
 //     TTree *T = tree_clone;  
   
 //     T->Fill(); //fill main tree  
 // //    cout<<T->IsA()->GetName()<<" "<<T->GetName()<<endl;  
 //     TList *li = T->GetListOfFriends();  
 //     TIter next(li);  
 //     TFriendElement* T_friend=0;  
 //     while( (T_friend = (TFriendElement*)next()) ){  
 // //   cout<<T_friend->IsA()->GetName()<<" "<<T_friend->GetName()<<hex << T_friend->GetTree() << dec<<endl;  
 //      T_friend->GetTree()->Fill();//fill friends  
 //     }  
2508    
2509  // }      // ---------------------------------------------------------------------
2510        // ATTENTION!!!
2511        // If data are not pre-processed with cleaner, the level0 tree may contain
2512        // spurious nested physics packets.
2513        // The GL_RUN variables EV_FROM, EV_TO and NEVENTS counts also these entries
2514        // while level2 tree DOES NOT!!
2515        // This means that "quellagiusta" in these cases is not correct.
2516        // In order to retrieve the correct level0 event, I implemented a check
2517        // of the OBT and pkt-number. In case of mismatch, the level0 entry number
2518        // is shift forward until when the packets match.
2519        // ---------------------------------------------------------------------
2520        Int_t answer = 0;
2521        Int_t shift =0;
2522        do{
2523            if(shift>0){
2524                cout << " level0 <--> level2 mismatch ( irun "<<irun<<" irunentry "<<irunentry<<" shift "<<shift<<" )"<<endl;
2525            }
2526            answer = l0_tree->GetEntry(quellagiusta+shift);
2527            shift++;
2528            if( !GetEventHeader() ){
2529                cout << "Int_t PamLevel2::GetYodaEntry() -- ERROR -- missing EventHeader "<<endl;
2530                return 0;
2531            }
2532    //      if( (quellagiusta+shift) == l0_tree->GetEntries()+1 )cout << ">>> end of level0 tree <<<"<<endl;
2533    //      cout << " GetOrbitalInfo()->OBT "<< GetOrbitalInfo()->OBT << endl;
2534    //      cout << " GetEventHeader()->GetPscuHeader()->GetOrbitalTime() "<< GetEventHeader()->GetPscuHeader()->GetOrbitalTime() << endl;
2535    //      cout << " GetOrbitalInfo()->pkt_num "<< GetOrbitalInfo()->pkt_num << endl;
2536    //      cout << " GetEventHeader()->GetPscuHeader()->GetCounter() "<< GetEventHeader()->GetPscuHeader()->GetCounter() << endl;
2537    
2538        }while( ( GetOrbitalInfo()->OBT != (UInt_t)(GetEventHeader()->GetPscuHeader()->GetOrbitalTime()) || GetOrbitalInfo()->pkt_num != (UInt_t)(GetEventHeader()->GetPscuHeader()->GetCounter())) && (quellagiusta+shift) < GetYodaTree()->GetEntries() );
2539    
2540    //    cout << "LA ENTRY GIUSTA E`: "<<quellagiusta<<" (spero...)"<<endl;
2541    //    return GetYodaTree()->GetEntry(quellagiusta);
2542        
2543        return answer;
2544        
2545    }

Legend:
Removed from v.1.9  
changed lines
  Added in v.1.26

  ViewVC Help
Powered by ViewVC 1.1.23