/[PAMELA software]/PamelaLevel2/inc/PamLevel2.h
ViewVC logotype

Diff of /PamelaLevel2/inc/PamLevel2.h

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

revision 1.8 by pam-fi, Mon Dec 11 17:49:03 2006 UTC revision 1.53 by emocchiutti, Fri Oct 17 19:55:10 2014 UTC
# Line 9  Line 9 
9    
10  #include <TRefArray.h>  #include <TRefArray.h>
11    
12    #include <TrkStruct.h>
13    #include <TrkParams.h>
14    #include <TrkLevel0.h>
15  #include <TrkLevel1.h>  #include <TrkLevel1.h>
16  #include <TrkLevel2.h>  #include <TrkLevel2.h>
17  #include <TrkHough.h>  #include <TrkHough.h>
18    #include <CaloLevel0.h>
19  #include <CaloLevel1.h>  #include <CaloLevel1.h>
20  #include <CaloLevel2.h>  #include <CaloLevel2.h>
21  #include <ToFLevel2.h>  #include <ToFLevel2.h>
# Line 20  Line 24 
24  #include <NDLevel2.h>  #include <NDLevel2.h>
25  #include <AcLevel2.h>  #include <AcLevel2.h>
26  #include <OrbitalInfo.h>  #include <OrbitalInfo.h>
27  #include <OrbitalInfo.h>  #include <InclinationInfo.h> // EMILIANO, now this library is included in libDarthVader.so
28                             // and the header installed in $PAM_INC. If you get the
29                             // error "Cannot find InclinationInfo.h" update your DV installation
30    #include <OrientationInfo.h> // EMILIANO, now this library is included in libDarthVader.so
31                             // and the header installed in $PAM_INC. If you get the
32                             // error "Cannot find InclinationInfo.h" update your DV installation
33  #include <RunInfo.h>  #include <RunInfo.h>
34  #include <GLTables.h>  #include <GLTables.h>
35    #include <ExtTrack.h>
36    
37    
38  #include <TTree.h>  #include <TTree.h>
39  #include <TFriendElement.h>  #include <TFriendElement.h>
40  #include <TChain.h>  #include <TChain.h>
41    #include <TChainElement.h>
42  #include <TFile.h>  #include <TFile.h>
43    #include <TTimeStamp.h>
44  #include <TList.h>  #include <TList.h>
45  #include <TKey.h>  #include <TKey.h>
46  #include <TSystemFile.h>  #include <TSystemFile.h>
47  #include <TSystemDirectory.h>  #include <TSystemDirectory.h>
48    #include <TSQLServer.h>
49    
50    #include <stdlib.h>
51  #include <iostream>  #include <iostream>
52    #include <iomanip>
53    using namespace std;
54    
55    #define NCLONES 10
56    
57  using namespace std;  using namespace std;
58    
59    
60    /////////////////////////////////////////////////////////////////////////////
61    /////////////////////////////////////////////////////////////////////////////
62    /////////////////////////////////////////////////////////////////////////////
63    /////////////////////////////////////////////////////////////////////////////
64    /**
65     * \brief Class to describe GP output (to be moved to PamelaDigitizer)
66     *
67     */
68    
69    class GPamela : public TObject {
70    
71     private:
72    
73     public:
74    
75        // Declaration of leave types
76        Int_t           Irun;
77        Int_t           Ievnt;
78        Int_t           Ipa;
79        Float_t         X0;
80        Float_t         Y0;
81        Float_t         Z0;
82        Float_t         Theta;
83        Float_t         Phi;
84        Float_t         P0;
85        Int_t           Nthtof;
86        UChar_t         Ipltof[100];   //[Nthtof]
87        UChar_t         Ipaddle[100];   //[Nthtof]
88        UChar_t         Ipartof[100];   //[Nthtof]
89        Float_t         Xintof[100];   //[Nthtof]
90        Float_t         Yintof[100];   //[Nthtof]
91        Float_t         Zintof[100];   //[Nthtof]
92        Float_t         Xouttof[100];   //[Nthtof]
93        Float_t         Youttof[100];   //[Nthtof]
94        Float_t         Zouttof[100];   //[Nthtof]
95        Float_t         Ereltof[100];   //[Nthtof]
96        Float_t         Timetof[100];   //[Nthtof]
97        Float_t         Pathtof[100];   //[Nthtof]
98        Float_t         P0tof[100];   //[Nthtof]
99        Int_t           Nthcat;
100        UChar_t         Iparcat[50];   //[Nthcat]
101        UChar_t         Icat[50];   //[Nthcat]
102        Float_t         Xincat[50];   //[Nthcat]
103        Float_t         Yincat[50];   //[Nthcat]
104        Float_t         Zincat[50];   //[Nthcat]
105        Float_t         Xoutcat[50];   //[Nthcat]
106        Float_t         Youtcat[50];   //[Nthcat]
107        Float_t         Zoutcat[50];   //[Nthcat]
108        Float_t         Erelcat[50];   //[Nthcat]
109        Float_t         Timecat[50];   //[Nthcat]
110        Float_t         Pathcat[50];   //[Nthcat]
111        Float_t         P0cat[50];   //[Nthcat]
112        Int_t           Nthcas;
113        UChar_t         Iparcas[50];   //[Nthcas]
114        UChar_t         Icas[50];   //[Nthcas]
115        Float_t         Xincas[50];   //[Nthcas]
116        Float_t         Yincas[50];   //[Nthcas]
117        Float_t         Zincas[50];   //[Nthcas]
118        Float_t         Xoutcas[50];   //[Nthcas]
119        Float_t         Youtcas[50];   //[Nthcas]
120        Float_t         Zoutcas[50];   //[Nthcas]
121        Float_t         Erelcas[50];   //[Nthcas]
122        Float_t         Timecas[50];   //[Nthcas]
123        Float_t         Pathcas[50];   //[Nthcas]
124        Float_t         P0cas[50];   //[Nthcas]
125        Int_t           Nthspe;
126        UChar_t         Iparspe[100];   //[Nthspe]
127        UChar_t         Itrpb[100];   //[Nthspe]
128        UChar_t         Itrsl[100];   //[Nthspe]
129        UChar_t         Itspa[100];   //[Nthspe]
130        Float_t         Xinspe[100];   //[Nthspe]
131        Float_t         Yinspe[100];   //[Nthspe]
132        Float_t         Zinspe[100];   //[Nthspe]
133        Float_t         Xoutspe[100];   //[Nthspe]
134        Float_t         Youtspe[100];   //[Nthspe]
135        Float_t         Zoutspe[100];   //[Nthspe]
136        Float_t         Xavspe[100];   //[Nthspe]
137        Float_t         Yavspe[100];   //[Nthspe]
138        Float_t         Zavspe[100];   //[Nthspe]
139        Float_t         Erelspe[100];   //[Nthspe]
140        Float_t         Pathspe[100];   //[Nthspe]
141        Float_t         P0spe[100];   //[Nthspe]
142        UChar_t         Nxmult[100];   //[Nthspe]
143        UChar_t         Nymult[100];   //[Nthspe]
144        Int_t           Nstrpx;
145        UChar_t         Npstripx[1000];   //[Nstrpx]
146        UChar_t         Ntstripx[1000];   //[Nstrpx]
147        UShort_t        Istripx[1000];   //[Nstrpx]
148        //UInt_t        Istripx[1000];   //[Nstrpx]
149        Float_t         Qstripx[1000];   //[Nstrpx]
150        Float_t         Xstripx[1000];   //[Nstrpx]
151        Int_t           Nstrpy;
152        UChar_t         Npstripy[1000];   //[Nstrpy]
153        UChar_t         Ntstripy[1000];   //[Nstrpy]
154        UShort_t        Istripy[1000];   //[Nstrpy]
155        //UInt_t        Istripy[1000];   //[Nstrpy]
156        Float_t         Qstripy[1000];   //[Nstrpy]
157        Float_t         Ystripy[1000];   //[Nstrpy]
158        Int_t           Nthcali;
159        UChar_t         Icaplane[4224];   //[Nthcali]
160        UChar_t         Icastrip[4224];   //[Nthcali]
161        Int_t           Icamod[4224];   //[Nthcali]
162        Float_t         Enestrip[4224];   //[Nthcali]
163        Int_t           Nthcal;
164        UChar_t         Icapl[1000];   //[Nthcal]
165        UChar_t         Icasi[1000];   //[Nthcal]
166        UChar_t         Icast[1000];   //[Nthcal]
167        Float_t         Xincal[1000];   //[Nthcal]
168        Float_t         Yincal[1000];   //[Nthcal]
169        Float_t         Zincal[1000];   //[Nthcal]
170        Float_t         Erelcal[1000];   //[Nthcal]
171        Int_t           Nthnd;
172        UChar_t         Itubend[200];   //[Nthnd]
173        UChar_t         Iparnd[200];   //[Nthnd]
174        Float_t         Xinnd[200];   //[Nthnd]
175        Float_t         Yinnd[200];   //[Nthnd]
176        Float_t         Zinnd[200];   //[Nthnd]
177        Float_t         Xoutnd[200];   //[Nthnd]
178        Float_t         Youtnd[200];   //[Nthnd]
179        Float_t         Zoutnd[200];   //[Nthnd]
180        Float_t         Erelnd[200];   //[Nthnd]
181        Float_t         Timend[200];   //[Nthnd]
182        Float_t         Pathnd[200];   //[Nthnd]
183        Float_t         P0nd[200];   //[Nthnd]
184        Int_t           Nthcard;
185        UChar_t         Iparcard[100];   //[Nthcard]
186        UChar_t         Icard[100];   //[Nthcard]
187        Float_t         Xincard[100];   //[Nthcard]
188        Float_t         Yincard[100];   //[Nthcard]
189        Float_t         Zincard[100];   //[Nthcard]
190        Float_t         Xoutcard[100];   //[Nthcard]
191        Float_t         Youtcard[100];   //[Nthcard]
192        Float_t         Zoutcard[100];   //[Nthcard]
193        Float_t         Erelcard[100];   //[Nthcard]
194        Float_t         Timecard[100];   //[Nthcard]
195        Float_t         Pathcard[100];   //[Nthcard]
196        Float_t         P0card[100];   //[Nthcard]
197    
198        GPamela(){ Clear(); };
199    //    ~GPamela(){ Delete(); };
200        ~GPamela(){ Clear(); };
201    
202        void Clear();
203        void Delete();
204        void SetBranchAddress( TChain* );
205        void SetBranchAddress( TTree* t){ SetBranchAddress( (TChain*)t );};
206    
207        GPamela*    GetGPamela(){return this;};
208    
209        ClassDef(GPamela,1);
210    
211    };
212    /////////////////////////////////////////////////////////////////////////////
213    /////////////////////////////////////////////////////////////////////////////
214    /////////////////////////////////////////////////////////////////////////////
215    /////////////////////////////////////////////////////////////////////////////
216    
217    
218    
219  /**  /**
220   * \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF).   * \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF).
221   *   *
222   * This class inherites from TrkTrack, CaloTrkVar ecc...   * This class inherites from TrkTrack, CaloTrkVar ecc...
223   *   *
224   */   */
225  class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {  //class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {
226    class PamTrack : public TObject {
227    
228   private:   private:
229    
230        ExtTrack*   trk_ext_track;
231      TrkTrack*   trk_track;      TrkTrack*   trk_track;
232      CaloTrkVar* calo_track;      CaloTrkVar* calo_track;
233      ToFTrkVar*  tof_track;      ToFTrkVar*  tof_track;
234        OrbitalInfoTrkVar*  orb_track;
235        Int_t pscore;
236        Int_t iscore;
237        bool candeleteobj;
238    
239    //    TClonesArray *Cluster;        ///< clusters (if stored)
240    
241   public:   public:
242    
243      PamTrack();      PamTrack();
244        PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
245        PamTrack(ExtTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
246        PamTrack(const PamTrack&);
247        ~PamTrack(){ Delete(); };
248    
249        void Clear(Option_t *option);
250        void Clear(){Clear("");};
251        void Delete();
252    
253        Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
254        Int_t GetScore(){return (pscore-iscore);};
255        Int_t GetPScore(){return (pscore);};
256        Int_t GetIScore(){return (iscore);};
257        void SetPScore(Int_t p){pscore=p;};
258        void SetIScore(Int_t i){iscore=i;};
259    
260      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );      ExtTrack*   GetExtTrack(){return trk_ext_track;};
261        TrkTrack*   GetTrkTrack(){return trk_track;};
     TrkTrack* GetTrkTrack(){return trk_track;};  
262      CaloTrkVar* GetCaloTrack(){return calo_track;};      CaloTrkVar* GetCaloTrack(){return calo_track;};
263      ToFTrkVar* GetToFTrack(){return tof_track;};      ToFTrkVar*  GetToFTrack(){return tof_track;};
264        OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};
265    
266      ClassDef(PamTrack,1);      ClassDef(PamTrack,3);
267  };  };
268    
269  /**  /**
270   * \brief Class to describe PAMELA LEVEL2 data.   * \brief Class to describe PAMELA LEVEL2 data.
271   *   *
272   * This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.   * This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.
273   * This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors.   * This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors.
274   */   */
275                    
276  class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {  //class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {
277    class PamLevel2 : public TObject {
278    
279   private:   private:
  public:  
280    
281        // ------------------------------
282        // variables to retrieve run info
283        // ------------------------------
284        Long64_t  irun;          //current run entry, relative to the TChain
285        Long64_t  irunt;         //current run entry, relative to the TTree in the TChain
286        Long64_t  runfirstentry; //first entry of run, relative to the TChain
287        Long64_t  runlastentry;  //last  entry of run, relative to the TChain
288        Long64_t  irunentry;     //position of current entry, relative to the run
289        Long64_t  totrunentry;
290        Long64_t  totrunentrymax;
291        Long64_t  totrunentrymin;
292        ULong64_t    prevabstime; // abs time of previous event
293        UInt_t    prevpktnum; // packet number of previous event
294        ULong64_t    abstime; // current abs time
295        UInt_t    pktnum; // current packet number
296        Bool_t isFragment;
297        Bool_t isSync; ///< the file is fully sinchronized between L2 and L0/DB data
298            //    Bool_t hasL0EE; ///< backward compatibility flag for preselected event list
299        Long64_t il0entry; ///< entry of level0 file in selection list
300        TChain*   run_tree;
301        TTree*    run_tree_clone;
302        TChain*   proc_tree;
303        TTree*    proc_tree_clone;
304        TChain*   sel_tree;
305        TTree*    sel_tree_clone;
306        int *irunoffset;
307        GL_TIMESYNC *gltsync; // Emiliano
308        Bool_t fUpdateRunInfo; // Emiliano
309        Bool_t fUseDBinRunInfo; // Emiliano
310    
311        Bool_t fDiscarded; ///< true if any file from a list is missing
312        TString customString; ///< String passed to ProcessingInfo
313    
314        // -------------------------------------------
315        // variables to store total live and dead time
316        // -------------------------------------------
317        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events
318    
319        // ---------------------------------
320        // variables to retrieve level0 info
321        // ---------------------------------
322        TFile*      l0_file;
323        TTree*      l0_tree;
324        Int_t       iroot;
325        Long64_t       prevshift; // shift needed to sync level2 data with the DB (nested+DV_skipped packets) - global position in level2 file EM
326        Long64_t       yprevshift; // shift needed to sync level2 data with level0 data (nested+DV_skipped packets) - run wise EM
327        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
328        // --------------------------------------------
329        // variables to store Level0 info (if required)
330        // --------------------------------------------
331        EventHeader* h0_obj;
332        TrkLevel0*   trk0_obj;
333        CaloLevel0*  calo0_obj;
334    
335        // ---------------------------
336        // variables to connect the DB
337        // ---------------------------
338        TSQLServer* dbc;
339        TString host;
340        TString user;
341        TString psw;
342    
343        // ---------------------------------------------
344        // variables to retrieve and clone level1-2 info
345        // ---------------------------------------------
346        TChain* pam_tree;
347        TTree*  pam_tree_clone[NCLONES];
348    
349        // --------------------------------------------
350        // variables to store Level1-2 info
351        // --------------------------------------------
352      TrkLevel1*   trk1_obj;      TrkLevel1*   trk1_obj;
353      TrkLevel2*   trk2_obj;      TrkLevel2*   trk2_obj;
354      TrkHough*    trkh_obj;      TrkHough*    trkh_obj;
   
355      CaloLevel1*  calo1_obj;      CaloLevel1*  calo1_obj;
356      CaloLevel2*  calo2_obj;      CaloLevel2*  calo2_obj;
357        ToFLevel2*   tof2_obj;
     ToFLevel2*   tof_obj;  
358      TrigLevel2*  trig_obj;      TrigLevel2*  trig_obj;
359      S4Level2*    s4_obj;      S4Level2*    s4_obj;
360      NDLevel2*    nd_obj;      NDLevel2*    nd_obj;
361      AcLevel2*    ac_obj;      AcLevel2*    ac_obj;
362      OrbitalInfo* orb_obj;      OrbitalInfo* orb2_obj;
           
363      GL_RUN*      run_obj;      GL_RUN*      run_obj;
364            SoftInfo*    soft_obj; // Emiliano
365  /*     TrkLevel1*   trk1_clone; */      ProcInfo*    proc_obj; //EM
366  /*     TrkLevel2*   trk2_clone; */  
367  /*     TrkHough*    trkh_clone; */  
368    
369  /*     CaloLevel1*  calo1_clone; */      Int_t extAlgFlag;
370  /*     CaloLevel2*  calo2_clone; */  
371        TClonesArray* trk_ext_obj;
372  /*     ToFLevel2*   tof_clone; */      TClonesArray* trk_ext_nuc_obj;
373  /*     TrigLevel2*  trig_clone; */      TClonesArray* trk_nuc_obj;
374  /*     S4Level2*    s4_clone; */  
375  /*     NDLevel2*    nd_clone; */  
376  /*     AcLevel2*    ac_clone; */      TClonesArray* calo_ext_obj;
377  /*     OrbitalInfo* orb_clone; */      TClonesArray* calo_ext_nuc_obj;
378        TClonesArray* calo_nuc_obj;
379  /*     TTree*       T_clone; //tracker */  
380  /*     TTree*       C_clone; //calo */      TClonesArray* tof_ext_obj;
381  /*     TTree*       O_clone; //tof */      TClonesArray* tof_ext_nuc_obj;
382  /*     TTree*       R_clone; //trigger */      TClonesArray* tof_nuc_obj;
383  /*     TTree*       S_clone; //s4 */  
384  /*     TTree*       A_clone; //anti */      TClonesArray* orb_ext_obj;
385  /*     TTree*       N_clone; //nd */      TClonesArray* orb_ext_nuc_obj;
386  /*     TTree*       B_clone; //orbital */      TClonesArray* orb_nuc_obj;
387    
388    
389        TrkLevel2*   trk2_nuc_obj;
390        CaloLevel2*  calo2_nuc_obj;
391        ToFLevel2*   tof2_nuc_obj;
392        OrbitalInfo* orb2_nuc_obj;
393    
394    
395    
396        // --------------------------------------------
397        // variables to handle GP
398        // --------------------------------------------
399        GPamela*     gp_obj;
400    
401        // -----------------------------------
402        // variables to sort the track
403        // -----------------------------------
404        TClonesArray *tsorted;
405        TClonesArray *timage;
406        TClonesArray *text; //array of pamtrack objects built from extended algorythm
407        TClonesArray *tsorted_nuc;
408        TClonesArray *timage_nuc;
409        TClonesArray *text_nuc; //array of pamtrack objects built from extended algorythm
410    
411        PamTrack    *track; //store the last required
412        TString howtosort;
413        Float_t sortthr;
414    
415        Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria
416        Long64_t lastsorted; // the last event whose tracks has been sorted
417    
     TTree*       tree_clone[8];  
418    
419      TRefArray*   sorted_tracks;      Bool_t issorted_new;     // true if the current event has already been sorted with the current sort criteria
420                Long64_t lastsorted_new; // the last event whose tracks has been sorted
421    
422        // --------------------------------------
423        // flags to load level0-1-2 trees/braches
424        // --------------------------------------
425      Bool_t CAL2; //level2 branch      Bool_t CAL2; //level2 branch
426      Bool_t CAL1; //level1 branch      Bool_t CAL1; //level1 branch
427        Bool_t CAL0; //yoda
428    
429      Bool_t TRK2; //level2 branch      Bool_t TRK2; //level2 branch
430      Bool_t TRK1; //level1 branch      Bool_t TRK1; //level1 branch
431      Bool_t TRKh; //hough branch      Bool_t TRKh; //hough branch
432        Bool_t TRK0; //yoda
433    
     Bool_t TRG;  
434      Bool_t TOF;      Bool_t TOF;
435        Bool_t TOF0;
436    
437        Bool_t TRG;
438      Bool_t S4;      Bool_t S4;
439      Bool_t ND;      Bool_t ND;
440      Bool_t AC;      Bool_t AC;
441      Bool_t ORB;      Bool_t ORB;
442            
443        Bool_t GP; //GP simulation tree
444    
445    
446        Bool_t EXT;//get tracks from extended tracking algorythm
447        Bool_t NUC;//get tracks from nuclei algorythm
448    
449        Bool_t PROC;//EM
450    
451        const char*  trkAlg;
452    
453        // -----------
454        // other flags
455        // -----------
456        Bool_t RUN;   // RunInfo tree
457        Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
458        Int_t  SELLI; // tree of selected events (there is a selection list)
459    
460        Bool_t DBG;
461    
462        void Initialize();
463    
464    
465     public:
466    
467        PamLevel2(TString ddir,TList *list,TString detlist);
468        PamLevel2(TString ddir,TString list,TString detlist);
469        PamLevel2(TString ddir,TString list);
470      PamLevel2();      PamLevel2();
471      ~PamLevel2();      ~PamLevel2(){Delete();};
472            
473      void Clear();      void Clear();
474        void Reset();
475        void Delete();
476    
477    
478  /*     TrkLevel1*   GetTrkLevel1()  { return trk1_obj;  }; */  //  =============================================
479  /*     TrkLevel2*   GetTrkLevel2()  { return trk2_obj;  }; */  //  methods needed to read/write Level0-1-2 trees
480  /*     TrkHough*    GetTrkHough()   { return trkh_obj;  }; */  //  =============================================
 /*     CaloLevel1*  GetCaloLevel1() { return calo1_obj; }; */  
 /*     CaloLevel2*  GetCaloLevel2() { return calo2_obj; }; */  
 /*     ToFLevel2*   GetToFLevel2()  { return tof_obj;  }; */  
 /*     TrigLevel2*  GetTrigLevel2() { return trig_obj; }; */  
 /*     S4Level2*    GetS4Level2()   { return s4_obj;   }; */  
 /*     NDLevel2*    GetNDLevel2()   { return nd_obj;   }; */  
 /*     AcLevel2*    GetAcLevel2()   { return ac_obj;   }; */  
 /*     OrbitalInfo* GetOrbitalInfo(){ return orb_obj;  }; */  
 /*     TrkLevel1*   GetTrkLevel1()  { return ; }; */  
 /*     TrkLevel2*   GetTrkLevel2()  { return TrkLevel2::GetTrkLevel2(); }; */  
 /*     TrkHough*    GetTrkHough()   { return TrkHough::GetTrkHough(); }; */  
 /*     CaloLevel1*  GetCaloLevel1() { return CaloLevel1::GetCaloLevel1(); }; */  
 /*     CaloLevel2*  GetCaloLevel2() { return CaloLevel2::GetCaloLevel2(); }; */  
 /*     ToFLevel2*   GetToFLevel2()  { return ToFLevel2::GetToFLevel2(); }; */  
 /*     TrigLevel2*  GetTrigLevel2() { return TrigLevel2::GetTrigLevel2(); }; */  
 /*     S4Level2*    GetS4Level2()   { return S4Level2::GetS4Level2(); }; */  
 /*     NDLevel2*    GetNDLevel2()   { return NDLevel2::GetNDLevel2(); }; */  
 /*     AcLevel2*    GetAcLevel2()   { return AcLevel2::GetAcLevel2(); }; */  
 /*     OrbitalInfo* GetOrbitalInfo(){ return OrbitalInfo::GetOrbitalInfo(); }; */  
     GL_RUN*      GetRunInfo()    { return run_obj; };  
           
 //  =========================================  
 //  methods needed to read/write Level2 trees  
 //  =========================================  
481      void    SetWhichTrees(TString);      void    SetWhichTrees(TString);
482      void    GetWhichTrees(TFile*);      void    GetWhichTrees(TFile*);
483    
# Line 176  class PamLevel2 : public CaloLevel2, pub Line 489  class PamLevel2 : public CaloLevel2, pub
489      TChain* GetPamTree(TList*,TString);      TChain* GetPamTree(TList*,TString);
490      TChain* GetPamTree(TString ddir,TString list,TString detlist){ return GetPamTree(GetListOfLevel2Files(ddir,list),detlist); };      TChain* GetPamTree(TString ddir,TString list,TString detlist){ return GetPamTree(GetListOfLevel2Files(ddir,list),detlist); };
491      TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };      TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
492        TChain* AddFriend(TString);
493    
494      TTree*  GetRunTree(TFile*);      TTree*  GetRunTree(TFile*);
495      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
496        TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
497        Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug
498        Bool_t  UpdateRunInfo(TTree *,Long64_t);
499        Bool_t  UpdateRunInfo(TChain *,Long64_t);
500    
501        void    SkipRunInfoUpdate(); // Emiliano
502        void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano
503        void    UseDBinRunInfo(){fUseDBinRunInfo = true;}; // [default] Emiliano
504        void    NoDBinRunInfo(){fUseDBinRunInfo = false;}; // EM
505        void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano
506    
507      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
508      void    CreateCloneTrees(TFile* file);      void    CreateCloneTrees(TFile* file);
509        void    CreateCloneTrees(){CreateCloneTrees(NULL);};
510      void    FillCloneTrees();      void    FillCloneTrees();
511      TTree  *GetCloneTree(TString);      TTree  *GetCloneTree(TString);
512      void    WriteCloneTrees();      void    WriteCloneTrees();
# Line 189  class PamLevel2 : public CaloLevel2, pub Line 514  class PamLevel2 : public CaloLevel2, pub
514      void    SetBranchAddress(TTree *);      void    SetBranchAddress(TTree *);
515      void    SetBranchAddress(TChain *);      void    SetBranchAddress(TChain *);
516    
517  //  =====================================  //    Int_t   GetEntry(Int_t);
518  //  methods needed to set brach addresses      Int_t   GetEntry(Long64_t);
519  //  =====================================      Long64_t GetReadEntry();
520  /*     void*       GetPointerToTrk(); */      Int_t   GetYodaEntry();
521  /*     void*       GetPointerToTrk(int s); */  
522  /*     void*       GetPointerToCalo(); */      void SetMaxShift(Int_t sh); // EMILIANO method to change the maxshift variable
523  /*     void*       GetPointerToCalo(int s); */  
524  /*     void*       GetPointerToToF(); */      Long64_t GetEntries(){return pam_tree->GetEntries();};
525  /*     void*       GetPointerToTrig(); */      TChain* GetPamTree(){return pam_tree;};
526  /*     void*       GetPointerToS4(); */      TChain* GetRunTree(){return run_tree;};
527  /*     void*       GetPointerToND(); */      TChain* GetProcessingInfoTree(){return proc_tree;};
528  /*     void*       GetPointerToAc(); */      TTree*  GetYodaTree();
529  /*     void*       GetPointerToOrb(); */  
530  /*     void*       GetPointerToRun(); */      TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things
531    
532    //  =============
533    //  DB connection
534    //  =============
535        void SetDBHost(TString str){host=str;};
536        void SetDBUser(TString str){user=str;};
537        void SetDBPsw(TString str){psw=str;};
538        TSQLServer* GetSQLServer(){return dbc;};
539        Bool_t SetDBConnection();
540    
541    //  ========================================================
542    //  methods needed to retrieve subdetector general variables
543    //  ========================================================
544      void*       GetPointerTo(const char*);      void*       GetPointerTo(const char*);
545        EventHeader  *GetEventHeader()  { return h0_obj; };
546    
547    /*     TrkLevel0    *GetTrkLevel0()   { return trk0_obj; }; */
548    /*     CaloLevel0   *GetCaloLevel0()  { return calo0_obj; }; */
549        TrkLevel0    *GetTrkLevel0();
550        CaloLevel0   *GetCaloLevel0();
551    
552    
553        TrkLevel1    *GetTrkLevel1()  { return trk1_obj; };
554        TrkLevel2    *GetTrkLevel2()  { return trk2_obj; };
555        TrkHough     *GetTrkHough()   { return trkh_obj; };
556        CaloLevel1   *GetCaloLevel1() { return calo1_obj;};
557        CaloLevel2   *GetCaloLevel2() { return calo2_obj;};
558        ToFLevel2    *GetToFLevel2()  { return tof2_obj;  };
559        TrigLevel2   *GetTrigLevel2() { return trig_obj; };
560        S4Level2     *GetS4Level2()   { return s4_obj;   };
561        NDLevel2     *GetNDLevel2()   { return nd_obj;   };
562        AcLevel2     *GetAcLevel2()   { return ac_obj;   };
563        OrbitalInfo  *GetOrbitalInfo(){ return orb2_obj;  };
564        GL_RUN       *GetRunInfo()    { return run_obj;  };
565        SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
566        ProcInfo     *GetProcInfo()   { return proc_obj; }; // Emiliano
567        GPamela      *GetGPamela()    { return gp_obj;   };
568    
569  //  ==============================================  //  ==============================================
570  //  methods to retrieve subdetectors stored tracks  //  methods to retrieve subdetectors stored tracks
571  //  ==============================================  //  ==============================================
572      TrkTrack*   GetTrkStoredTrack(int seqno){return TrkLevel2::GetStoredTrack(seqno); };      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
573      CaloTrkVar* GetCaloStoredTrack(int);      CaloTrkVar* GetCaloStoredTrack(int);
574      ToFTrkVar*  GetToFStoredTrack(int);      ToFTrkVar*  GetToFStoredTrack(int);
575        OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);
576      //      //
577      PamTrack*   GetPamTrackAlong(TrkTrack*);      PamTrack*   GetPamTrackAlong(TrkTrack*);
578  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
579  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
580            
581  //  =================================  //  =================================
582  //  methods to retrieve pamela tracks  //  methods to retrieve pamela tracks
583  //  =================================  //  =================================
584      PamTrack*     GetStoredTrack(Int_t);      PamTrack*     GetStoredTrack(Int_t);
585      //      //
586      void          SortTracks(TString);          //assign value to variable sorted_tracks      void          SortTracks(TString);          //assign value to variable sorted_tracks
587      TRefArray*    GetTracks();                  //return sorted_tracks      void          SortTracks();          //assign value to variable sorted_tracks
588      PamTrack*     GetTrack(int);  //    TRefArray*    GetTracks();                  //return sorted_tracks
589      PamTrack*     GetTrackImage(int);      TClonesArray* GetTracks();                  //return sorted_tracks
590                PamTrack*     GetTrackOld(int);
591        PamTrack*     GetTrackImageOld(int);
592        PamTrack*     GetTrack(int i){ return GetTrackOld(i); };
593        PamTrack*     GetTrackImage(int i){ return GetTrackImageOld(i); };
594       void          SetSortingMethod(TString how);
595        TString       GetSortingMethod(){return howtosort;};
596        void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
597        Float_t       GetNucleiSortingThreshold(){ return sortthr;};
598    
599    
600    
601    
602    //  =================================
603    //  methods to handle new tracking algorythms
604    //  =================================
605    /*     Int_t GetTrackingAlgorythmFlag(){return extAlgFlag;}; */
606    /*     void  SetTrackingAlgorythm(const char * alg); */
607    /*     const char* GetTrackingAlgorythm(); */
608    
609        Int_t  GetNTracks(const char*);
610        PamTrack*     GetTrack(int,const char*);
611        PamTrack*     GetTrackImage(int, const char*);
612    
613    
614    /*     Int_t  GetNTracks(){return GetNTracks(trkAlg); }; */
615    /*     PamTrack*     GetTrack(int i) { return GetTrack(i,trkAlg); }; */
616    /*     PamTrack*     GetTrackImage(int i){ return GetTrackImage(i,trkAlg); }; */
617    
618        void          SortTracksNew();          //assign value to variable sorted_tracks
619    
620    
621  //  =================================  //  =================================
622  //  other methods  //  other methods
623  //  =================================  //  =================================
624      Int_t ntrk(){ return TrkLevel2::ntrk(); };  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
625    
626    
627        Long64_t  GetRunID(){ return irun;};
628        Long64_t  GetRunEntryID(){ return irunentry; };
629    
630        // -----
631        // flags methods
632        // -----
633        Bool_t IsCAL2(){return CAL2;}; //level2 branch
634        Bool_t IsCAL1(){return CAL1;}; //level1 branch
635        Bool_t IsCAL0(){return CAL0;}; //yoda
636    
637        Bool_t IsTRK2(){return TRK2;}; //level2 branch
638        Bool_t IsTRK1(){return TRK1;}; //level1 branch
639        Bool_t IsTRKh(){return TRKh;}; //hough branch
640        Bool_t IsTRK0(){return TRK0;}; //yoda
641    
642        Bool_t IsTOF(){return TOF;};
643        Bool_t IsTOF0(){return TOF0;};
644    
645        Bool_t IsTRG(){return TRG;};
646        Bool_t IsS4(){return S4;};
647        Bool_t IsND(){return ND;};
648        Bool_t IsAC(){return AC;};
649        Bool_t IsORB(){return ORB;};
650    
651        Bool_t IsRUN(){return RUN;};
652    
653        Bool_t IsEXT(){return EXT;};
654        Bool_t IsNUC(){return NUC;};
655    
656        Bool_t IsGood(Bool_t strict=true);
657    
658        Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)
659        void  SetSELLI(Int_t selli){SELLI=selli;};
660    
661        Bool_t IsPROC(){return PROC;};//EM
662    
663        Bool_t HasDiscardedFiles(){return fDiscarded;}; ///< true if something went wrong and not all files were processed
664    
665        Bool_t IsGP(){return ISGP;};
666    
667        void   SetGP(Bool_t tt){ISGP=tt;};
668    
669        void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
670    
671      ClassDef(PamLevel2,1);      void SetCustomString(TString cus){customString=cus;};
672        TString GetCustomString(){return customString;};
673    
674        ClassDef(PamLevel2,4);
675  };  };
676    
677  //===============================================================================================  //===============================================================================================
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.53
  ViewVC Help
Powered by ViewVC 1.1.23