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