/[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.32 by pam-fi, Tue Nov 20 18:16:29 2007 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>
 #include <OrbitalInfo.h>  
27  #include <RunInfo.h>  #include <RunInfo.h>
28  #include <GLTables.h>  #include <GLTables.h>
29    
30    
31  #include <TTree.h>  #include <TTree.h>
32    #include <TFriendElement.h>
33  #include <TChain.h>  #include <TChain.h>
34    #include <TChainElement.h>
35  #include <TFile.h>  #include <TFile.h>
36  #include <TList.h>  #include <TList.h>
37  #include <TKey.h>  #include <TKey.h>
38  #include <TSystemFile.h>  #include <TSystemFile.h>
39  #include <TSystemDirectory.h>  #include <TSystemDirectory.h>
40    #include <TSQLServer.h>
41    
42    #include <stdlib.h>
43  #include <iostream>  #include <iostream>
44    #include <iomanip>
45    using namespace std;
46    
47    #define NCLONES 10
48    
49  using namespace std;  using namespace std;
50    
51    
52    /////////////////////////////////////////////////////////////////////////////
53    /////////////////////////////////////////////////////////////////////////////
54    /////////////////////////////////////////////////////////////////////////////
55    /////////////////////////////////////////////////////////////////////////////
56    /**
57     * \brief Class to describe GP output (to be moved to PamelaDigitizer)
58     *
59     */
60    
61    class GPamela : public TObject {
62    
63     private:
64    
65     public:
66    
67        // Declaration of leave types
68        Int_t           Irun;
69        Int_t           Ievnt;
70        Int_t           Ipa;
71        Float_t         X0;
72        Float_t         Y0;
73        Float_t         Z0;
74        Float_t         Theta;
75        Float_t         Phi;
76        Float_t         P0;
77        Int_t           Nthtof;
78        UChar_t         Ipltof[100];   //[Nthtof]
79        UChar_t         Ipaddle[100];   //[Nthtof]
80        UChar_t         Ipartof[100];   //[Nthtof]
81        Float_t         Xintof[100];   //[Nthtof]
82        Float_t         Yintof[100];   //[Nthtof]
83        Float_t         Zintof[100];   //[Nthtof]
84        Float_t         Xouttof[100];   //[Nthtof]
85        Float_t         Youttof[100];   //[Nthtof]
86        Float_t         Zouttof[100];   //[Nthtof]
87        Float_t         Ereltof[100];   //[Nthtof]
88        Float_t         Timetof[100];   //[Nthtof]
89        Float_t         Pathtof[100];   //[Nthtof]
90        Float_t         P0tof[100];   //[Nthtof]
91        Int_t           Nthcat;
92        UChar_t         Iparcat[50];   //[Nthcat]
93        UChar_t         Icat[50];   //[Nthcat]
94        Float_t         Xincat[50];   //[Nthcat]
95        Float_t         Yincat[50];   //[Nthcat]
96        Float_t         Zincat[50];   //[Nthcat]
97        Float_t         Xoutcat[50];   //[Nthcat]
98        Float_t         Youtcat[50];   //[Nthcat]
99        Float_t         Zoutcat[50];   //[Nthcat]
100        Float_t         Erelcat[50];   //[Nthcat]
101        Float_t         Timecat[50];   //[Nthcat]
102        Float_t         Pathcat[50];   //[Nthcat]
103        Float_t         P0cat[50];   //[Nthcat]
104        Int_t           Nthcas;
105        UChar_t         Iparcas[50];   //[Nthcas]
106        UChar_t         Icas[50];   //[Nthcas]
107        Float_t         Xincas[50];   //[Nthcas]
108        Float_t         Yincas[50];   //[Nthcas]
109        Float_t         Zincas[50];   //[Nthcas]
110        Float_t         Xoutcas[50];   //[Nthcas]
111        Float_t         Youtcas[50];   //[Nthcas]
112        Float_t         Zoutcas[50];   //[Nthcas]
113        Float_t         Erelcas[50];   //[Nthcas]
114        Float_t         Timecas[50];   //[Nthcas]
115        Float_t         Pathcas[50];   //[Nthcas]
116        Float_t         P0cas[50];   //[Nthcas]
117        Int_t           Nthspe;
118        UChar_t         Iparspe[100];   //[Nthspe]
119        UChar_t         Itrpb[100];   //[Nthspe]
120        UChar_t         Itrsl[100];   //[Nthspe]
121        UChar_t         Itspa[100];   //[Nthspe]
122        Float_t         Xinspe[100];   //[Nthspe]
123        Float_t         Yinspe[100];   //[Nthspe]
124        Float_t         Zinspe[100];   //[Nthspe]
125        Float_t         Xoutspe[100];   //[Nthspe]
126        Float_t         Youtspe[100];   //[Nthspe]
127        Float_t         Zoutspe[100];   //[Nthspe]
128        Float_t         Xavspe[100];   //[Nthspe]
129        Float_t         Yavspe[100];   //[Nthspe]
130        Float_t         Zavspe[100];   //[Nthspe]
131        Float_t         Erelspe[100];   //[Nthspe]
132        Float_t         Pathspe[100];   //[Nthspe]
133        Float_t         P0spe[100];   //[Nthspe]
134        UChar_t         Nxmult[100];   //[Nthspe]
135        UChar_t         Nymult[100];   //[Nthspe]
136        Int_t           Nstrpx;
137        UChar_t         Npstripx[1000];   //[Nstrpx]
138        UChar_t         Ntstripx[1000];   //[Nstrpx]
139        UShort_t        Istripx[1000];   //[Nstrpx]
140        //UInt_t        Istripx[1000];   //[Nstrpx]
141        Float_t         Qstripx[1000];   //[Nstrpx]
142        Float_t         Xstripx[1000];   //[Nstrpx]
143        Int_t           Nstrpy;
144        UChar_t         Npstripy[1000];   //[Nstrpy]
145        UChar_t         Ntstripy[1000];   //[Nstrpy]
146        UShort_t        Istripy[1000];   //[Nstrpy]
147        //UInt_t        Istripy[1000];   //[Nstrpy]
148        Float_t         Qstripy[1000];   //[Nstrpy]
149        Float_t         Ystripy[1000];   //[Nstrpy]
150        Int_t           Nthcali;
151        UChar_t         Icaplane[4224];   //[Nthcali]
152        UChar_t         Icastrip[4224];   //[Nthcali]
153        Int_t           Icamod[4224];   //[Nthcali]
154        Float_t         Enestrip[4224];   //[Nthcali]
155        Int_t           Nthcal;
156        UChar_t         Icapl[1000];   //[Nthcal]
157        UChar_t         Icasi[1000];   //[Nthcal]
158        UChar_t         Icast[1000];   //[Nthcal]
159        Float_t         Xincal[1000];   //[Nthcal]
160        Float_t         Yincal[1000];   //[Nthcal]
161        Float_t         Zincal[1000];   //[Nthcal]
162        Float_t         Erelcal[1000];   //[Nthcal]
163        Int_t           Nthnd;
164        UChar_t         Itubend[200];   //[Nthnd]
165        UChar_t         Iparnd[200];   //[Nthnd]
166        Float_t         Xinnd[200];   //[Nthnd]
167        Float_t         Yinnd[200];   //[Nthnd]
168        Float_t         Zinnd[200];   //[Nthnd]
169        Float_t         Xoutnd[200];   //[Nthnd]
170        Float_t         Youtnd[200];   //[Nthnd]
171        Float_t         Zoutnd[200];   //[Nthnd]
172        Float_t         Erelnd[200];   //[Nthnd]
173        Float_t         Timend[200];   //[Nthnd]
174        Float_t         Pathnd[200];   //[Nthnd]
175        Float_t         P0nd[200];   //[Nthnd]
176        Int_t           Nthcard;
177        UChar_t         Iparcard[100];   //[Nthcard]
178        UChar_t         Icard[100];   //[Nthcard]
179        Float_t         Xincard[100];   //[Nthcard]
180        Float_t         Yincard[100];   //[Nthcard]
181        Float_t         Zincard[100];   //[Nthcard]
182        Float_t         Xoutcard[100];   //[Nthcard]
183        Float_t         Youtcard[100];   //[Nthcard]
184        Float_t         Zoutcard[100];   //[Nthcard]
185        Float_t         Erelcard[100];   //[Nthcard]
186        Float_t         Timecard[100];   //[Nthcard]
187        Float_t         Pathcard[100];   //[Nthcard]
188        Float_t         P0card[100];   //[Nthcard]
189    
190        GPamela(){ Clear(); };
191    //    ~GPamela(){ Delete(); };
192        ~GPamela(){ Clear(); };
193    
194        void Clear();
195        void Delete();
196        void SetBranchAddress( TChain* );
197        void SetBranchAddress( TTree* t){ SetBranchAddress( (TChain*)t );};
198    
199        GPamela*    GetGPamela(){return this;};
200    
201        ClassDef(GPamela,1);
202    
203    };
204    /////////////////////////////////////////////////////////////////////////////
205    /////////////////////////////////////////////////////////////////////////////
206    /////////////////////////////////////////////////////////////////////////////
207    /////////////////////////////////////////////////////////////////////////////
208    
209    
210    
211  /**  /**
212   * \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).
213   *   *
214   * This class inherites from TrkTrack, CaloTrkVar ecc...   * This class inherites from TrkTrack, CaloTrkVar ecc...
215   *   *
216   */   */
217  class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {  //class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {
218    class PamTrack : public TObject {
219    
220   private:   private:
221    
222      TrkTrack*   trk_track;      TrkTrack*   trk_track;
223      CaloTrkVar* calo_track;      CaloTrkVar* calo_track;
224      ToFTrkVar*  tof_track;      ToFTrkVar*  tof_track;
225        Int_t pscore;
226        Int_t iscore;
227        bool candeleteobj;
228    
229    //    TClonesArray *Cluster;        ///< clusters (if stored)
230    
231   public:   public:
232    
233      PamTrack();      PamTrack();
   
234      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );
235        PamTrack(const PamTrack&);
236        ~PamTrack(){ Delete(); };
237    
238        void Clear();
239        void Clear(Option_t *option){Clear();};
240        void Delete();
241    
242        Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
243        Int_t GetScore(){return (pscore-iscore);};
244        Int_t GetPScore(){return (pscore);};
245        Int_t GetIScore(){return (iscore);};
246        void SetPScore(Int_t p){pscore=p;};
247        void SetIScore(Int_t i){iscore=i;};
248    
249        TrkTrack*   GetTrkTrack(){return trk_track;};
250        CaloTrkVar* GetCaloTrack(){return calo_track;};
251        ToFTrkVar*  GetToFTrack(){return tof_track;};
252    
253      ClassDef(PamTrack,1);      ClassDef(PamTrack,1);
254  };  };
# Line 64  class PamTrack : public TrkTrack, public Line 260  class PamTrack : public TrkTrack, public
260   * 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.
261   */   */
262                                    
263  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 {
264    class PamLevel2 : public TObject {
265    
266   private:   private:
267    
268      TrkLevel1*   trk_l1_obj;      // ------------------------------
269          TrkLevel2*   trk_obj;      // variables to retrieve run info
270          CaloLevel2*  calo_obj;      // ------------------------------
271        Long64_t  irun;          //current run entry, relative to the TChain
272        Long64_t  irunt;         //current run entry, relative to the TTree in the TChain
273        Long64_t  runfirstentry; //first entry of run, relative to the TChain
274        Long64_t  runlastentry;  //last  entry of run, relative to the TChain
275        Long64_t  irunentry;     //position of current entry, relative to the run
276        TChain*   run_tree;
277        TTree*    run_tree_clone;
278        TChain*   sel_tree;
279        TTree*    sel_tree_clone;
280        int *irunoffset;
281    
282        // -------------------------------------------
283        // variables to store total live and dead time
284        // -------------------------------------------
285        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events  
286    
287        // ---------------------------------
288        // variables to retrieve level0 info
289        // ---------------------------------
290        TFile*      l0_file;
291        TTree*      l0_tree;
292        Int_t       iroot;      
293        Int_t       prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time
294        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
295        // --------------------------------------------
296        // variables to store Level0 info (if required)
297        // --------------------------------------------
298        EventHeader* h0_obj;
299        TrkLevel0*   trk0_obj;
300        CaloLevel0*  calo0_obj;
301    
302        // ---------------------------
303        // variables to connect the DB
304        // ---------------------------
305        TSQLServer* dbc;
306        TString host;
307        TString user;
308        TString psw;
309    
310        // ---------------------------------------------
311        // variables to retrieve and clone level1-2 info
312        // ---------------------------------------------
313        TChain* pam_tree;
314        TTree*  pam_tree_clone[NCLONES];
315    
316        // --------------------------------------------
317        // variables to store Level1-2 info
318        // --------------------------------------------
319        TrkLevel1*   trk1_obj;
320        TrkLevel2*   trk2_obj;
321        TrkHough*    trkh_obj;
322        CaloLevel1*  calo1_obj;
323        CaloLevel2*  calo2_obj;
324      ToFLevel2*   tof_obj;      ToFLevel2*   tof_obj;
325      TrigLevel2*  trig_obj;      TrigLevel2*  trig_obj;
326      S4Level2*    s4_obj;      S4Level2*    s4_obj;
327      NDLevel2*    nd_obj;      NDLevel2*    nd_obj;
328      AcLevel2*    ac_obj;      AcLevel2*    ac_obj;
329          OrbitalInfo* orb_obj;      OrbitalInfo* orb_obj;
330                GL_RUN*      run_obj;
331          GL_RUN*      run_obj;      SoftInfo*    soft_obj; // Emiliano
332        
333          TRefArray*   sorted_tracks;      // --------------------------------------------
334                // variables to handle GP
335          TChain*      Tout;      // --------------------------------------------
336                GPamela*     gp_obj;
337   public:  
338        // -----------------------------------
339        // variables to sort the track
340        // -----------------------------------
341        TClonesArray *tsorted;
342        TClonesArray *timage;
343        PamTrack    *track; //store the last required
344        TString howtosort;
345        Float_t sortthr;
346    
347        // --------------------------------------
348        // flags to load level0-1-2 trees/braches
349        // --------------------------------------
350        Bool_t CAL2; //level2 branch
351        Bool_t CAL1; //level1 branch
352        Bool_t CAL0; //yoda
353    
354        Bool_t TRK2; //level2 branch
355        Bool_t TRK1; //level1 branch
356        Bool_t TRKh; //hough branch
357        Bool_t TRK0; //yoda
358    
359        Bool_t TOF;
360        Bool_t TOF0;
361    
362        Bool_t TRG;
363        Bool_t S4;
364        Bool_t ND;
365        Bool_t AC;
366        Bool_t ORB;
367    
368        Bool_t GP; //GP simulation tree
369    
370        // -----------
371        // other flags
372        // -----------
373        Bool_t RUN;   // RunInfo tree
374        Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
375        Int_t  SELLI; // tree of selected events (there is a selection list)
376    
377          Bool_t CAL;      void Initialize();
378          Bool_t TRK;  
379          Bool_t TRG;  
380          Bool_t TOF;   public:
381          Bool_t S4;          
382          Bool_t ND;      PamLevel2(TString ddir,TString list,TString detlist);
383          Bool_t AC;      PamLevel2(TString ddir,TString list);
         Bool_t ORB;  
           
         Bool_t TRK_L1;  
   
384      PamLevel2();      PamLevel2();
385          ~PamLevel2();      ~PamLevel2(){Delete();};
386                    
387          void Clear();      void Clear();
388        void Reset();
389        void Delete();
390    
391            
392    //  =============================================
393    //  methods needed to read/write Level0-1-2 trees
394    //  =============================================
395        void    SetWhichTrees(TString);
396        void    GetWhichTrees(TFile*);
397    
398        TTree*  GetPamTree(TFile*,TString);
399        TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};    
400    
401        TList*  GetListOfLevel2Files(TString, TString);
402        Bool_t  CheckLevel2File(TString);
403        TChain* GetPamTree(TList*,TString);
404        TChain* GetPamTree(TString ddir,TString list,TString detlist){ return GetPamTree(GetListOfLevel2Files(ddir,list),detlist); };
405        TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
406        TChain* AddFriend(TString);
407    
408        TTree*  GetRunTree(TFile*);
409        TChain* GetRunTree(TList*);
410        TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
411        Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
412        Bool_t  UpdateRunInfo(TTree *,Long64_t);
413        Bool_t  UpdateRunInfo(TChain *,Long64_t);
414    
415        void    CreateCloneTrees0(TChain *fChain, TFile* file);
416        void    CreateCloneTrees(TFile* file);
417        void    CreateCloneTrees(){CreateCloneTrees(NULL);};
418        void    FillCloneTrees();
419        TTree  *GetCloneTree(TString);
420        void    WriteCloneTrees();
421    
422        void    SetBranchAddress(TTree *);
423        void    SetBranchAddress(TChain *);
424    
425    //    Int_t   GetEntry(Int_t);
426        Int_t   GetEntry(Long64_t);
427        Int_t   GetYodaEntry();
428    
429          TrkLevel1*   GetTrkLevel1()  { return trk_l1_obj;  };      void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
         TrkLevel2*   GetTrkLevel2()  { return trk_obj;  };  
         CaloLevel2*  GetCaloLevel2() { return calo_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;  };  
         GL_RUN*      GetRunInfo()    { return run_obj;  };  
           
 //  =====================================  
 //  methods needed to read Level2 trees  
 //  =====================================  
         void    SetWhichTrees(TString);  
         TList*  GetListOfLevel2Files(TString, TString);  
         Bool_t  CheckLevel2File(TString);  
         void    SetBranchAddress();  
         TTree*  GetPamTree(TFile*,TString);  
         TTree*  GetPamTree(TFile* f){return GetPamTree(f,"+ALL");};  
         TChain* GetPamTree(TList*,TString);  
         TTree*  GetRunTree(TFile*);  
         TChain* GetRunTree(TList*);  
430            
431  //  =====================================      Long64_t GetEntries(){return pam_tree->GetEntries();};
432  //  methods needed to set brach addresses      TChain* GetPamTree(){return pam_tree;};
433  //  =====================================      TChain* GetRunTree(){return run_tree;};
434          void*       GetPointerToTrk()  {return &trk_obj;};      TTree*  GetYodaTree();
435          void*       GetPointerToTrk(int s){  
436                  switch (s) {      TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things
437                          case 1:  return &trk_l1_obj;  
438                          case 2:  return &trk_obj;  //  =============
439                          default: return 0;  //  DB connection
440                  };  //  =============
441          };      void SetDBHost(TString str){host=str;};
442          void*       GetPointerToCalo() {return &calo_obj; };      void SetDBUser(TString str){user=str;};
443          void*       GetPointerToToF()  {return &tof_obj;  };      void SetDBPsw(TString str){psw=str;};
444          void*       GetPointerToTrig() {return &trig_obj; };      TSQLServer* GetSQLServer(){return dbc;};
445          void*       GetPointerToS4()   {return &s4_obj;   };      Bool_t SetDBConnection();
446          void*       GetPointerToND()   {return &nd_obj;   };  
447          void*       GetPointerToAc()   {return &ac_obj;   };  //  ========================================================
448          void*       GetPointerToOrb()  {return &orb_obj;  };  //  methods needed to retrieve subdetector general variables
449          void*       GetPointerToRun()  {return &run_obj;  };  //  ========================================================
450          void*       GetPointerTo(const char*);      void*       GetPointerTo(const char*);
451        EventHeader  *GetEventHeader()  { return h0_obj; };
452    
453    /*     TrkLevel0    *GetTrkLevel0()   { return trk0_obj; }; */
454    /*     CaloLevel0   *GetCaloLevel0()  { return calo0_obj; }; */
455        TrkLevel0    *GetTrkLevel0();
456        CaloLevel0   *GetCaloLevel0();
457    
458    
459        TrkLevel1    *GetTrkLevel1()  { return trk1_obj; };
460        TrkLevel2    *GetTrkLevel2()  { return trk2_obj; };
461        TrkHough     *GetTrkHough()   { return trkh_obj; };
462        CaloLevel1   *GetCaloLevel1() { return calo1_obj;};
463        CaloLevel2   *GetCaloLevel2() { return calo2_obj;};
464        ToFLevel2    *GetToFLevel2()  { return tof_obj;  };
465        TrigLevel2   *GetTrigLevel2() { return trig_obj; };
466        S4Level2     *GetS4Level2()   { return s4_obj;   };
467        NDLevel2     *GetNDLevel2()   { return nd_obj;   };
468        AcLevel2     *GetAcLevel2()   { return ac_obj;   };
469        OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
470        GL_RUN       *GetRunInfo()    { return run_obj;  };
471        SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
472        GPamela      *GetGPamela()    { return gp_obj; };
473    
474  //  ==============================================  //  ==============================================
475  //  methods to retrieve subdetectors stored tracks  //  methods to retrieve subdetectors stored tracks
476  //  ==============================================  //  ==============================================
477          TrkTrack*   GetTrkStoredTrack(int seqno){return TrkLevel2::GetStoredTrack(seqno); };      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
478          CaloTrkVar* GetCaloStoredTrack(int);      CaloTrkVar* GetCaloStoredTrack(int);
479          ToFTrkVar*  GetToFStoredTrack(int);      ToFTrkVar*  GetToFStoredTrack(int);
480          //      //
481          PamTrack*   GetPamTrackAlong(TrkTrack*);      PamTrack*   GetPamTrackAlong(TrkTrack*);
482  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
483  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
484                    
# Line 162  class PamLevel2 : public CaloLevel2, pub Line 486  class PamLevel2 : public CaloLevel2, pub
486  //  methods to retrieve pamela tracks  //  methods to retrieve pamela tracks
487  //  =================================  //  =================================
488      PamTrack*     GetStoredTrack(Int_t);      PamTrack*     GetStoredTrack(Int_t);
489          //      //
490          void          SortTracks(TString);          //assign value to variable sorted_tracks      void          SortTracks(TString);          //assign value to variable sorted_tracks
491          TRefArray*    GetTracks();                  //return sorted_tracks      void          SortTracks();          //assign value to variable sorted_tracks
492          PamTrack*     GetTrack(int);  //    TRefArray*    GetTracks();                  //return sorted_tracks
493        TClonesArray* GetTracks();                  //return sorted_tracks
494        PamTrack*     GetTrack(int);
495      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImage(int);
496        void          SetSortingMethod(TString how){ howtosort = how; return;};
497        TString       GetSortingMethod(){return howtosort;};
498        void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
499        Float_t       GetNucleiSortingThreshold(){ return sortthr;};
500                    
501    //  =================================
502    //  other methods
503    //  =================================
504    //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
505    
506        // -----
507        // flags methods
508        // -----
509        Bool_t IsCAL2(){return CAL2;}; //level2 branch
510        Bool_t IsCAL1(){return CAL1;}; //level1 branch
511        Bool_t IsCAL0(){return CAL0;}; //yoda
512    
513        Bool_t IsTRK2(){return TRK2;}; //level2 branch
514        Bool_t IsTRK1(){return TRK1;}; //level1 branch
515        Bool_t IsTRKh(){return TRKh;}; //hough branch
516        Bool_t IsTRK0(){return TRK0;}; //yoda
517    
518        Bool_t IsTOF(){return TOF;};
519        Bool_t IsTOF0(){return TOF0;};
520    
521        Bool_t IsTRG(){return TRG;};
522        Bool_t IsS4(){return S4;};
523        Bool_t IsND(){return ND;};
524        Bool_t IsAC(){return AC;};
525        Bool_t IsORB(){return ORB;};
526    
527        Bool_t IsRUN(){return RUN;};
528    
529        Bool_t IsGood();
530    
531        Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)
532    
533        Bool_t IsGP(){return ISGP;};
534    
535        void   SetGP(Bool_t tt){ISGP=tt;};
536    
537    
538      ClassDef(PamLevel2,1);      ClassDef(PamLevel2,1);
539  };  };
540    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.23