/[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.27 by mocchiut, Wed Sep 12 09:25:21 2007 UTC revision 1.35 by mocchiut, Tue Sep 30 07:45:22 2008 UTC
# Line 24  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 <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 <RunInfo.h>  #include <RunInfo.h>
31  #include <GLTables.h>  #include <GLTables.h>
32    
# Line 48  using namespace std; Line 51  using namespace std;
51    
52  using namespace std;  using namespace std;
53    
54    
55    /////////////////////////////////////////////////////////////////////////////
56    /////////////////////////////////////////////////////////////////////////////
57    /////////////////////////////////////////////////////////////////////////////
58    /////////////////////////////////////////////////////////////////////////////
59    /**
60     * \brief Class to describe GP output (to be moved to PamelaDigitizer)
61     *
62     */
63    
64    class GPamela : public TObject {
65    
66     private:
67    
68     public:
69    
70        // Declaration of leave types
71        Int_t           Irun;
72        Int_t           Ievnt;
73        Int_t           Ipa;
74        Float_t         X0;
75        Float_t         Y0;
76        Float_t         Z0;
77        Float_t         Theta;
78        Float_t         Phi;
79        Float_t         P0;
80        Int_t           Nthtof;
81        UChar_t         Ipltof[100];   //[Nthtof]
82        UChar_t         Ipaddle[100];   //[Nthtof]
83        UChar_t         Ipartof[100];   //[Nthtof]
84        Float_t         Xintof[100];   //[Nthtof]
85        Float_t         Yintof[100];   //[Nthtof]
86        Float_t         Zintof[100];   //[Nthtof]
87        Float_t         Xouttof[100];   //[Nthtof]
88        Float_t         Youttof[100];   //[Nthtof]
89        Float_t         Zouttof[100];   //[Nthtof]
90        Float_t         Ereltof[100];   //[Nthtof]
91        Float_t         Timetof[100];   //[Nthtof]
92        Float_t         Pathtof[100];   //[Nthtof]
93        Float_t         P0tof[100];   //[Nthtof]
94        Int_t           Nthcat;
95        UChar_t         Iparcat[50];   //[Nthcat]
96        UChar_t         Icat[50];   //[Nthcat]
97        Float_t         Xincat[50];   //[Nthcat]
98        Float_t         Yincat[50];   //[Nthcat]
99        Float_t         Zincat[50];   //[Nthcat]
100        Float_t         Xoutcat[50];   //[Nthcat]
101        Float_t         Youtcat[50];   //[Nthcat]
102        Float_t         Zoutcat[50];   //[Nthcat]
103        Float_t         Erelcat[50];   //[Nthcat]
104        Float_t         Timecat[50];   //[Nthcat]
105        Float_t         Pathcat[50];   //[Nthcat]
106        Float_t         P0cat[50];   //[Nthcat]
107        Int_t           Nthcas;
108        UChar_t         Iparcas[50];   //[Nthcas]
109        UChar_t         Icas[50];   //[Nthcas]
110        Float_t         Xincas[50];   //[Nthcas]
111        Float_t         Yincas[50];   //[Nthcas]
112        Float_t         Zincas[50];   //[Nthcas]
113        Float_t         Xoutcas[50];   //[Nthcas]
114        Float_t         Youtcas[50];   //[Nthcas]
115        Float_t         Zoutcas[50];   //[Nthcas]
116        Float_t         Erelcas[50];   //[Nthcas]
117        Float_t         Timecas[50];   //[Nthcas]
118        Float_t         Pathcas[50];   //[Nthcas]
119        Float_t         P0cas[50];   //[Nthcas]
120        Int_t           Nthspe;
121        UChar_t         Iparspe[100];   //[Nthspe]
122        UChar_t         Itrpb[100];   //[Nthspe]
123        UChar_t         Itrsl[100];   //[Nthspe]
124        UChar_t         Itspa[100];   //[Nthspe]
125        Float_t         Xinspe[100];   //[Nthspe]
126        Float_t         Yinspe[100];   //[Nthspe]
127        Float_t         Zinspe[100];   //[Nthspe]
128        Float_t         Xoutspe[100];   //[Nthspe]
129        Float_t         Youtspe[100];   //[Nthspe]
130        Float_t         Zoutspe[100];   //[Nthspe]
131        Float_t         Xavspe[100];   //[Nthspe]
132        Float_t         Yavspe[100];   //[Nthspe]
133        Float_t         Zavspe[100];   //[Nthspe]
134        Float_t         Erelspe[100];   //[Nthspe]
135        Float_t         Pathspe[100];   //[Nthspe]
136        Float_t         P0spe[100];   //[Nthspe]
137        UChar_t         Nxmult[100];   //[Nthspe]
138        UChar_t         Nymult[100];   //[Nthspe]
139        Int_t           Nstrpx;
140        UChar_t         Npstripx[1000];   //[Nstrpx]
141        UChar_t         Ntstripx[1000];   //[Nstrpx]
142        UShort_t        Istripx[1000];   //[Nstrpx]
143        //UInt_t        Istripx[1000];   //[Nstrpx]
144        Float_t         Qstripx[1000];   //[Nstrpx]
145        Float_t         Xstripx[1000];   //[Nstrpx]
146        Int_t           Nstrpy;
147        UChar_t         Npstripy[1000];   //[Nstrpy]
148        UChar_t         Ntstripy[1000];   //[Nstrpy]
149        UShort_t        Istripy[1000];   //[Nstrpy]
150        //UInt_t        Istripy[1000];   //[Nstrpy]
151        Float_t         Qstripy[1000];   //[Nstrpy]
152        Float_t         Ystripy[1000];   //[Nstrpy]
153        Int_t           Nthcali;
154        UChar_t         Icaplane[4224];   //[Nthcali]
155        UChar_t         Icastrip[4224];   //[Nthcali]
156        Int_t           Icamod[4224];   //[Nthcali]
157        Float_t         Enestrip[4224];   //[Nthcali]
158        Int_t           Nthcal;
159        UChar_t         Icapl[1000];   //[Nthcal]
160        UChar_t         Icasi[1000];   //[Nthcal]
161        UChar_t         Icast[1000];   //[Nthcal]
162        Float_t         Xincal[1000];   //[Nthcal]
163        Float_t         Yincal[1000];   //[Nthcal]
164        Float_t         Zincal[1000];   //[Nthcal]
165        Float_t         Erelcal[1000];   //[Nthcal]
166        Int_t           Nthnd;
167        UChar_t         Itubend[200];   //[Nthnd]
168        UChar_t         Iparnd[200];   //[Nthnd]
169        Float_t         Xinnd[200];   //[Nthnd]
170        Float_t         Yinnd[200];   //[Nthnd]
171        Float_t         Zinnd[200];   //[Nthnd]
172        Float_t         Xoutnd[200];   //[Nthnd]
173        Float_t         Youtnd[200];   //[Nthnd]
174        Float_t         Zoutnd[200];   //[Nthnd]
175        Float_t         Erelnd[200];   //[Nthnd]
176        Float_t         Timend[200];   //[Nthnd]
177        Float_t         Pathnd[200];   //[Nthnd]
178        Float_t         P0nd[200];   //[Nthnd]
179        Int_t           Nthcard;
180        UChar_t         Iparcard[100];   //[Nthcard]
181        UChar_t         Icard[100];   //[Nthcard]
182        Float_t         Xincard[100];   //[Nthcard]
183        Float_t         Yincard[100];   //[Nthcard]
184        Float_t         Zincard[100];   //[Nthcard]
185        Float_t         Xoutcard[100];   //[Nthcard]
186        Float_t         Youtcard[100];   //[Nthcard]
187        Float_t         Zoutcard[100];   //[Nthcard]
188        Float_t         Erelcard[100];   //[Nthcard]
189        Float_t         Timecard[100];   //[Nthcard]
190        Float_t         Pathcard[100];   //[Nthcard]
191        Float_t         P0card[100];   //[Nthcard]
192    
193        GPamela(){ Clear(); };
194    //    ~GPamela(){ Delete(); };
195        ~GPamela(){ Clear(); };
196    
197        void Clear();
198        void Delete();
199        void SetBranchAddress( TChain* );
200        void SetBranchAddress( TTree* t){ SetBranchAddress( (TChain*)t );};
201    
202        GPamela*    GetGPamela(){return this;};
203    
204        ClassDef(GPamela,1);
205    
206    };
207    /////////////////////////////////////////////////////////////////////////////
208    /////////////////////////////////////////////////////////////////////////////
209    /////////////////////////////////////////////////////////////////////////////
210    /////////////////////////////////////////////////////////////////////////////
211    
212    
213    
214  /**  /**
215   * \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).
216   *   *
# Line 108  class PamLevel2 : public TObject { Line 271  class PamLevel2 : public TObject {
271      // ------------------------------      // ------------------------------
272      // variables to retrieve run info      // variables to retrieve run info
273      // ------------------------------      // ------------------------------
 //    Int_t     irun;          //current run entry, relative to the TChain  
274      Long64_t  irun;          //current run entry, relative to the TChain      Long64_t  irun;          //current run entry, relative to the TChain
275  //    ULong64_t runfirstentry; //first entry of run, relative to the TChain      Long64_t  irunt;         //current run entry, relative to the TTree in the TChain
 //    ULong64_t runlastentry;  //last  entry of run, relative to the TChain  
 //    Int_t     irunentry;     //position of current entry, relative to the run  
276      Long64_t  runfirstentry; //first entry of run, relative to the TChain      Long64_t  runfirstentry; //first entry of run, relative to the TChain
277      Long64_t  runlastentry;  //last  entry of run, relative to the TChain      Long64_t  runlastentry;  //last  entry of run, relative to the TChain
278      Long64_t  irunentry;     //position of current entry, relative to the run      Long64_t  irunentry;     //position of current entry, relative to the run
# Line 120  class PamLevel2 : public TObject { Line 280  class PamLevel2 : public TObject {
280      TTree*    run_tree_clone;      TTree*    run_tree_clone;
281      TChain*   sel_tree;      TChain*   sel_tree;
282      TTree*    sel_tree_clone;      TTree*    sel_tree_clone;
283        int *irunoffset;
284    
285      // -------------------------------------------      // -------------------------------------------
286      // variables to store total live and dead time      // variables to store total live and dead time
# Line 133  class PamLevel2 : public TObject { Line 294  class PamLevel2 : public TObject {
294      TTree*      l0_tree;      TTree*      l0_tree;
295      Int_t       iroot;            Int_t       iroot;      
296      Int_t       prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time      Int_t       prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time
297        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
298      // --------------------------------------------      // --------------------------------------------
299      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
300      // --------------------------------------------      // --------------------------------------------
# Line 155  class PamLevel2 : public TObject { Line 317  class PamLevel2 : public TObject {
317      TTree*  pam_tree_clone[NCLONES];      TTree*  pam_tree_clone[NCLONES];
318    
319      // --------------------------------------------      // --------------------------------------------
320      // variables to store Level0 info (if required)      // variables to store Level1-2 info
321      // --------------------------------------------      // --------------------------------------------
322      TrkLevel1*   trk1_obj;      TrkLevel1*   trk1_obj;
323      TrkLevel2*   trk2_obj;      TrkLevel2*   trk2_obj;
# Line 171  class PamLevel2 : public TObject { Line 333  class PamLevel2 : public TObject {
333      GL_RUN*      run_obj;      GL_RUN*      run_obj;
334      SoftInfo*    soft_obj; // Emiliano      SoftInfo*    soft_obj; // Emiliano
335    
336        // --------------------------------------------
337        // variables to handle GP
338        // --------------------------------------------
339        GPamela*     gp_obj;
340    
341      // -----------------------------------      // -----------------------------------
342      // variables to sort the track      // variables to sort the track
343      // -----------------------------------      // -----------------------------------
# Line 180  class PamLevel2 : public TObject { Line 347  class PamLevel2 : public TObject {
347      TString howtosort;      TString howtosort;
348      Float_t sortthr;      Float_t sortthr;
349    
350      // -----      // --------------------------------------
351      // flags      // flags to load level0-1-2 trees/braches
352      // -----      // --------------------------------------
353      Bool_t CAL2; //level2 branch      Bool_t CAL2; //level2 branch
354      Bool_t CAL1; //level1 branch      Bool_t CAL1; //level1 branch
355      Bool_t CAL0; //yoda      Bool_t CAL0; //yoda
# Line 201  class PamLevel2 : public TObject { Line 368  class PamLevel2 : public TObject {
368      Bool_t AC;      Bool_t AC;
369      Bool_t ORB;      Bool_t ORB;
370    
371      Bool_t RUN;      Bool_t GP; //GP simulation tree
372    
373      Bool_t ISGP; //GP simulation      // -----------
374        // other flags
375        // -----------
376        Bool_t RUN;   // RunInfo tree
377        Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
378        Int_t  SELLI; // tree of selected events (there is a selection list)
379    
380        Bool_t DBG;
     Int_t  SELLI; //trees of selected events (there is a selection list)  
381    
382      void Initialize();      void Initialize();
383    
# Line 242  class PamLevel2 : public TObject { Line 413  class PamLevel2 : public TObject {
413      TTree*  GetRunTree(TFile*);      TTree*  GetRunTree(TFile*);
414      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
415      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
416        Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
417      Bool_t  UpdateRunInfo(TTree *,Long64_t);      Bool_t  UpdateRunInfo(TTree *,Long64_t);
418      Bool_t  UpdateRunInfo(TChain *,Long64_t);      Bool_t  UpdateRunInfo(TChain *,Long64_t);
     Bool_t  UpdateRunInfo(Long64_t iev);//{ return UpdateRunInfo(run_tree,iev); };  
419    
420      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
421      void    CreateCloneTrees(TFile* file);      void    CreateCloneTrees(TFile* file);
# Line 260  class PamLevel2 : public TObject { Line 431  class PamLevel2 : public TObject {
431      Int_t   GetEntry(Long64_t);      Int_t   GetEntry(Long64_t);
432      Int_t   GetYodaEntry();      Int_t   GetYodaEntry();
433    
434        void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
435        
436      Long64_t GetEntries(){return pam_tree->GetEntries();};      Long64_t GetEntries(){return pam_tree->GetEntries();};
437      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
438      TChain* GetRunTree(){return run_tree;};      TChain* GetRunTree(){return run_tree;};
# Line 301  class PamLevel2 : public TObject { Line 474  class PamLevel2 : public TObject {
474      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
475      GL_RUN       *GetRunInfo()    { return run_obj;  };      GL_RUN       *GetRunInfo()    { return run_obj;  };
476      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
477        GPamela      *GetGPamela()    { return gp_obj;   };
478    
479  //  ==============================================  //  ==============================================
480  //  methods to retrieve subdetectors stored tracks  //  methods to retrieve subdetectors stored tracks
# Line 334  class PamLevel2 : public TObject { Line 508  class PamLevel2 : public TObject {
508  //  =================================  //  =================================
509  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
510    
511    
512        Long64_t  GetRunID(){ return irun;};
513        Long64_t  GetRunEntryID(){ return irunentry; };
514    
515      // -----      // -----
516      // flags methods      // flags methods
517      // -----      // -----
# Line 365  class PamLevel2 : public TObject { Line 543  class PamLevel2 : public TObject {
543    
544      void   SetGP(Bool_t tt){ISGP=tt;};      void   SetGP(Bool_t tt){ISGP=tt;};
545    
546        void SetDebug(Bool_t set){ DBG=set; };
547    
548      ClassDef(PamLevel2,1);      ClassDef(PamLevel2,1);
549  };  };
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.23