/[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.24 by pam-fi, Fri May 25 13:38:50 2007 UTC revision 1.39 by mocchiut, Wed Dec 10 11:56:31 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 <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    
# Line 31  Line 37 
37  #include <TTree.h>  #include <TTree.h>
38  #include <TFriendElement.h>  #include <TFriendElement.h>
39  #include <TChain.h>  #include <TChain.h>
40    #include <TChainElement.h>
41  #include <TFile.h>  #include <TFile.h>
42  #include <TList.h>  #include <TList.h>
43  #include <TKey.h>  #include <TKey.h>
# Line 38  Line 45 
45  #include <TSystemDirectory.h>  #include <TSystemDirectory.h>
46  #include <TSQLServer.h>  #include <TSQLServer.h>
47    
48    #include <stdlib.h>
49  #include <iostream>  #include <iostream>
50    #include <iomanip>
51    using namespace std;
52    
53  #define NCLONES 10  #define NCLONES 10
54    
55  using namespace std;  using namespace std;
56    
57    
58    /////////////////////////////////////////////////////////////////////////////
59    /////////////////////////////////////////////////////////////////////////////
60    /////////////////////////////////////////////////////////////////////////////
61    /////////////////////////////////////////////////////////////////////////////
62    /**
63     * \brief Class to describe GP output (to be moved to PamelaDigitizer)
64     *
65     */
66    
67    class GPamela : public TObject {
68    
69     private:
70    
71     public:
72    
73        // Declaration of leave types
74        Int_t           Irun;
75        Int_t           Ievnt;
76        Int_t           Ipa;
77        Float_t         X0;
78        Float_t         Y0;
79        Float_t         Z0;
80        Float_t         Theta;
81        Float_t         Phi;
82        Float_t         P0;
83        Int_t           Nthtof;
84        UChar_t         Ipltof[100];   //[Nthtof]
85        UChar_t         Ipaddle[100];   //[Nthtof]
86        UChar_t         Ipartof[100];   //[Nthtof]
87        Float_t         Xintof[100];   //[Nthtof]
88        Float_t         Yintof[100];   //[Nthtof]
89        Float_t         Zintof[100];   //[Nthtof]
90        Float_t         Xouttof[100];   //[Nthtof]
91        Float_t         Youttof[100];   //[Nthtof]
92        Float_t         Zouttof[100];   //[Nthtof]
93        Float_t         Ereltof[100];   //[Nthtof]
94        Float_t         Timetof[100];   //[Nthtof]
95        Float_t         Pathtof[100];   //[Nthtof]
96        Float_t         P0tof[100];   //[Nthtof]
97        Int_t           Nthcat;
98        UChar_t         Iparcat[50];   //[Nthcat]
99        UChar_t         Icat[50];   //[Nthcat]
100        Float_t         Xincat[50];   //[Nthcat]
101        Float_t         Yincat[50];   //[Nthcat]
102        Float_t         Zincat[50];   //[Nthcat]
103        Float_t         Xoutcat[50];   //[Nthcat]
104        Float_t         Youtcat[50];   //[Nthcat]
105        Float_t         Zoutcat[50];   //[Nthcat]
106        Float_t         Erelcat[50];   //[Nthcat]
107        Float_t         Timecat[50];   //[Nthcat]
108        Float_t         Pathcat[50];   //[Nthcat]
109        Float_t         P0cat[50];   //[Nthcat]
110        Int_t           Nthcas;
111        UChar_t         Iparcas[50];   //[Nthcas]
112        UChar_t         Icas[50];   //[Nthcas]
113        Float_t         Xincas[50];   //[Nthcas]
114        Float_t         Yincas[50];   //[Nthcas]
115        Float_t         Zincas[50];   //[Nthcas]
116        Float_t         Xoutcas[50];   //[Nthcas]
117        Float_t         Youtcas[50];   //[Nthcas]
118        Float_t         Zoutcas[50];   //[Nthcas]
119        Float_t         Erelcas[50];   //[Nthcas]
120        Float_t         Timecas[50];   //[Nthcas]
121        Float_t         Pathcas[50];   //[Nthcas]
122        Float_t         P0cas[50];   //[Nthcas]
123        Int_t           Nthspe;
124        UChar_t         Iparspe[100];   //[Nthspe]
125        UChar_t         Itrpb[100];   //[Nthspe]
126        UChar_t         Itrsl[100];   //[Nthspe]
127        UChar_t         Itspa[100];   //[Nthspe]
128        Float_t         Xinspe[100];   //[Nthspe]
129        Float_t         Yinspe[100];   //[Nthspe]
130        Float_t         Zinspe[100];   //[Nthspe]
131        Float_t         Xoutspe[100];   //[Nthspe]
132        Float_t         Youtspe[100];   //[Nthspe]
133        Float_t         Zoutspe[100];   //[Nthspe]
134        Float_t         Xavspe[100];   //[Nthspe]
135        Float_t         Yavspe[100];   //[Nthspe]
136        Float_t         Zavspe[100];   //[Nthspe]
137        Float_t         Erelspe[100];   //[Nthspe]
138        Float_t         Pathspe[100];   //[Nthspe]
139        Float_t         P0spe[100];   //[Nthspe]
140        UChar_t         Nxmult[100];   //[Nthspe]
141        UChar_t         Nymult[100];   //[Nthspe]
142        Int_t           Nstrpx;
143        UChar_t         Npstripx[1000];   //[Nstrpx]
144        UChar_t         Ntstripx[1000];   //[Nstrpx]
145        UShort_t        Istripx[1000];   //[Nstrpx]
146        //UInt_t        Istripx[1000];   //[Nstrpx]
147        Float_t         Qstripx[1000];   //[Nstrpx]
148        Float_t         Xstripx[1000];   //[Nstrpx]
149        Int_t           Nstrpy;
150        UChar_t         Npstripy[1000];   //[Nstrpy]
151        UChar_t         Ntstripy[1000];   //[Nstrpy]
152        UShort_t        Istripy[1000];   //[Nstrpy]
153        //UInt_t        Istripy[1000];   //[Nstrpy]
154        Float_t         Qstripy[1000];   //[Nstrpy]
155        Float_t         Ystripy[1000];   //[Nstrpy]
156        Int_t           Nthcali;
157        UChar_t         Icaplane[4224];   //[Nthcali]
158        UChar_t         Icastrip[4224];   //[Nthcali]
159        Int_t           Icamod[4224];   //[Nthcali]
160        Float_t         Enestrip[4224];   //[Nthcali]
161        Int_t           Nthcal;
162        UChar_t         Icapl[1000];   //[Nthcal]
163        UChar_t         Icasi[1000];   //[Nthcal]
164        UChar_t         Icast[1000];   //[Nthcal]
165        Float_t         Xincal[1000];   //[Nthcal]
166        Float_t         Yincal[1000];   //[Nthcal]
167        Float_t         Zincal[1000];   //[Nthcal]
168        Float_t         Erelcal[1000];   //[Nthcal]
169        Int_t           Nthnd;
170        UChar_t         Itubend[200];   //[Nthnd]
171        UChar_t         Iparnd[200];   //[Nthnd]
172        Float_t         Xinnd[200];   //[Nthnd]
173        Float_t         Yinnd[200];   //[Nthnd]
174        Float_t         Zinnd[200];   //[Nthnd]
175        Float_t         Xoutnd[200];   //[Nthnd]
176        Float_t         Youtnd[200];   //[Nthnd]
177        Float_t         Zoutnd[200];   //[Nthnd]
178        Float_t         Erelnd[200];   //[Nthnd]
179        Float_t         Timend[200];   //[Nthnd]
180        Float_t         Pathnd[200];   //[Nthnd]
181        Float_t         P0nd[200];   //[Nthnd]
182        Int_t           Nthcard;
183        UChar_t         Iparcard[100];   //[Nthcard]
184        UChar_t         Icard[100];   //[Nthcard]
185        Float_t         Xincard[100];   //[Nthcard]
186        Float_t         Yincard[100];   //[Nthcard]
187        Float_t         Zincard[100];   //[Nthcard]
188        Float_t         Xoutcard[100];   //[Nthcard]
189        Float_t         Youtcard[100];   //[Nthcard]
190        Float_t         Zoutcard[100];   //[Nthcard]
191        Float_t         Erelcard[100];   //[Nthcard]
192        Float_t         Timecard[100];   //[Nthcard]
193        Float_t         Pathcard[100];   //[Nthcard]
194        Float_t         P0card[100];   //[Nthcard]
195    
196        GPamela(){ Clear(); };
197    //    ~GPamela(){ Delete(); };
198        ~GPamela(){ Clear(); };
199    
200        void Clear();
201        void Delete();
202        void SetBranchAddress( TChain* );
203        void SetBranchAddress( TTree* t){ SetBranchAddress( (TChain*)t );};
204    
205        GPamela*    GetGPamela(){return this;};
206    
207        ClassDef(GPamela,1);
208    
209    };
210    /////////////////////////////////////////////////////////////////////////////
211    /////////////////////////////////////////////////////////////////////////////
212    /////////////////////////////////////////////////////////////////////////////
213    /////////////////////////////////////////////////////////////////////////////
214    
215    
216    
217  /**  /**
218   * \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).
219   *   *
# Line 58  class PamTrack : public TObject { Line 228  class PamTrack : public TObject {
228      TrkTrack*   trk_track;      TrkTrack*   trk_track;
229      CaloTrkVar* calo_track;      CaloTrkVar* calo_track;
230      ToFTrkVar*  tof_track;      ToFTrkVar*  tof_track;
231        OrbitalInfoTrkVar*  orb_track;
232        Int_t pscore;
233        Int_t iscore;
234      bool candeleteobj;      bool candeleteobj;
235    
236  //    TClonesArray *Cluster;        ///< clusters (if stored)  //    TClonesArray *Cluster;        ///< clusters (if stored)
# Line 66  class PamTrack : public TObject { Line 238  class PamTrack : public TObject {
238   public:   public:
239    
240      PamTrack();      PamTrack();
241      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
242      PamTrack(const PamTrack&);      PamTrack(const PamTrack&);
243      ~PamTrack(){ Delete(); };      ~PamTrack(){ Delete(); };
244    
# Line 74  class PamTrack : public TObject { Line 246  class PamTrack : public TObject {
246      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();};
247      void Delete();      void Delete();
248    
249        Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
250        Int_t GetScore(){return (pscore-iscore);};
251        Int_t GetPScore(){return (pscore);};
252        Int_t GetIScore(){return (iscore);};
253        void SetPScore(Int_t p){pscore=p;};
254        void SetIScore(Int_t i){iscore=i;};
255    
256      TrkTrack*   GetTrkTrack(){return trk_track;};      TrkTrack*   GetTrkTrack(){return trk_track;};
257      CaloTrkVar* GetCaloTrack(){return calo_track;};      CaloTrkVar* GetCaloTrack(){return calo_track;};
258      ToFTrkVar*  GetToFTrack(){return tof_track;};      ToFTrkVar*  GetToFTrack(){return tof_track;};
259        OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};
260    
261      ClassDef(PamTrack,1);      ClassDef(PamTrack,2);
262  };  };
263    
264  /**  /**
# Line 96  class PamLevel2 : public TObject { Line 276  class PamLevel2 : public TObject {
276      // ------------------------------      // ------------------------------
277      // variables to retrieve run info      // variables to retrieve run info
278      // ------------------------------      // ------------------------------
 //    Int_t     irun;          //current run entry, relative to the TChain  
279      Long64_t  irun;          //current run entry, relative to the TChain      Long64_t  irun;          //current run entry, relative to the TChain
280  //    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  
281      Long64_t  runfirstentry; //first entry of run, relative to the TChain      Long64_t  runfirstentry; //first entry of run, relative to the TChain
282      Long64_t  runlastentry;  //last  entry of run, relative to the TChain      Long64_t  runlastentry;  //last  entry of run, relative to the TChain
283      Long64_t  irunentry;     //position of current entry, relative to the run      Long64_t  irunentry;     //position of current entry, relative to the run
# Line 108  class PamLevel2 : public TObject { Line 285  class PamLevel2 : public TObject {
285      TTree*    run_tree_clone;      TTree*    run_tree_clone;
286      TChain*   sel_tree;      TChain*   sel_tree;
287      TTree*    sel_tree_clone;      TTree*    sel_tree_clone;
288            int *irunoffset;
289    
290        // -------------------------------------------
291        // variables to store total live and dead time
292        // -------------------------------------------
293        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events  
294    
295      // ---------------------------------      // ---------------------------------
296      // variables to retrieve level0 info      // variables to retrieve level0 info
297      // ---------------------------------      // ---------------------------------
# Line 116  class PamLevel2 : public TObject { Line 299  class PamLevel2 : public TObject {
299      TTree*      l0_tree;      TTree*      l0_tree;
300      Int_t       iroot;            Int_t       iroot;      
301      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
302        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
303      // --------------------------------------------      // --------------------------------------------
304      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
305      // --------------------------------------------      // --------------------------------------------
# Line 138  class PamLevel2 : public TObject { Line 322  class PamLevel2 : public TObject {
322      TTree*  pam_tree_clone[NCLONES];      TTree*  pam_tree_clone[NCLONES];
323    
324      // --------------------------------------------      // --------------------------------------------
325      // variables to store Level0 info (if required)      // variables to store Level1-2 info
326      // --------------------------------------------      // --------------------------------------------
327      TrkLevel1*   trk1_obj;      TrkLevel1*   trk1_obj;
328      TrkLevel2*   trk2_obj;      TrkLevel2*   trk2_obj;
# Line 154  class PamLevel2 : public TObject { Line 338  class PamLevel2 : public TObject {
338      GL_RUN*      run_obj;      GL_RUN*      run_obj;
339      SoftInfo*    soft_obj; // Emiliano      SoftInfo*    soft_obj; // Emiliano
340    
341        // --------------------------------------------
342        // variables to handle GP
343        // --------------------------------------------
344        GPamela*     gp_obj;
345    
346      // -----------------------------------      // -----------------------------------
347      // variables to sort the track      // variables to sort the track
348      // -----------------------------------      // -----------------------------------
# Line 163  class PamLevel2 : public TObject { Line 352  class PamLevel2 : public TObject {
352      TString howtosort;      TString howtosort;
353      Float_t sortthr;      Float_t sortthr;
354    
355      // -----      // --------------------------------------
356      // flags      // flags to load level0-1-2 trees/braches
357      // -----      // --------------------------------------
358      Bool_t CAL2; //level2 branch      Bool_t CAL2; //level2 branch
359      Bool_t CAL1; //level1 branch      Bool_t CAL1; //level1 branch
360      Bool_t CAL0; //yoda      Bool_t CAL0; //yoda
# Line 184  class PamLevel2 : public TObject { Line 373  class PamLevel2 : public TObject {
373      Bool_t AC;      Bool_t AC;
374      Bool_t ORB;      Bool_t ORB;
375    
376      Bool_t RUN;      Bool_t GP; //GP simulation tree
377    
378        // -----------
379        // other flags
380        // -----------
381        Bool_t RUN;   // RunInfo tree
382        Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
383        Int_t  SELLI; // tree of selected events (there is a selection list)
384    
385      Int_t  SELLI; //trees of selected events (there is a selection list)      Bool_t DBG;
386    
387      void Initialize();      void Initialize();
388    
389    
390   public:   public:
391                    
392        PamLevel2(TString ddir,TList *list,TString detlist);
393      PamLevel2(TString ddir,TString list,TString detlist);      PamLevel2(TString ddir,TString list,TString detlist);
394      PamLevel2(TString ddir,TString list);      PamLevel2(TString ddir,TString list);
395      PamLevel2();      PamLevel2();
# Line 210  class PamLevel2 : public TObject { Line 407  class PamLevel2 : public TObject {
407      void    GetWhichTrees(TFile*);      void    GetWhichTrees(TFile*);
408    
409      TTree*  GetPamTree(TFile*,TString);      TTree*  GetPamTree(TFile*,TString);
410      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};    
411    
412      TList*  GetListOfLevel2Files(TString, TString);      TList*  GetListOfLevel2Files(TString, TString);
413      Bool_t  CheckLevel2File(TString);      Bool_t  CheckLevel2File(TString);
414      TChain* GetPamTree(TList*,TString);      TChain* GetPamTree(TList*,TString);
415      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); };
416      TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };      TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
417        TChain* AddFriend(TString);
418    
419      TTree*  GetRunTree(TFile*);      TTree*  GetRunTree(TFile*);
420      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
421      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
422        Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
423      Bool_t  UpdateRunInfo(TTree *,Long64_t);      Bool_t  UpdateRunInfo(TTree *,Long64_t);
424      Bool_t  UpdateRunInfo(TChain *,Long64_t);      Bool_t  UpdateRunInfo(TChain *,Long64_t);
     Bool_t  UpdateRunInfo(Long64_t iev);//{ return UpdateRunInfo(run_tree,iev); };  
425    
426      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
427      void    CreateCloneTrees(TFile* file);      void    CreateCloneTrees(TFile* file);
# Line 239  class PamLevel2 : public TObject { Line 437  class PamLevel2 : public TObject {
437      Int_t   GetEntry(Long64_t);      Int_t   GetEntry(Long64_t);
438      Int_t   GetYodaEntry();      Int_t   GetYodaEntry();
439    
440        void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
441        
442      Long64_t GetEntries(){return pam_tree->GetEntries();};      Long64_t GetEntries(){return pam_tree->GetEntries();};
443      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
444      TChain* GetRunTree(){return run_tree;};      TChain* GetRunTree(){return run_tree;};
# Line 280  class PamLevel2 : public TObject { Line 480  class PamLevel2 : public TObject {
480      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
481      GL_RUN       *GetRunInfo()    { return run_obj;  };      GL_RUN       *GetRunInfo()    { return run_obj;  };
482      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
483        GPamela      *GetGPamela()    { return gp_obj;   };
484    
485  //  ==============================================  //  ==============================================
486  //  methods to retrieve subdetectors stored tracks  //  methods to retrieve subdetectors stored tracks
# Line 287  class PamLevel2 : public TObject { Line 488  class PamLevel2 : public TObject {
488      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
489      CaloTrkVar* GetCaloStoredTrack(int);      CaloTrkVar* GetCaloStoredTrack(int);
490      ToFTrkVar*  GetToFStoredTrack(int);      ToFTrkVar*  GetToFStoredTrack(int);
491        OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);
492      //      //
493      PamTrack*   GetPamTrackAlong(TrkTrack*);      PamTrack*   GetPamTrackAlong(TrkTrack*);
494  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
# Line 313  class PamLevel2 : public TObject { Line 515  class PamLevel2 : public TObject {
515  //  =================================  //  =================================
516  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
517    
518    
519        Long64_t  GetRunID(){ return irun;};
520        Long64_t  GetRunEntryID(){ return irunentry; };
521    
522      // -----      // -----
523      // flags methods      // flags methods
524      // -----      // -----
# Line 336  class PamLevel2 : public TObject { Line 542  class PamLevel2 : public TObject {
542    
543      Bool_t IsRUN(){return RUN;};      Bool_t IsRUN(){return RUN;};
544    
545      Bool_t IsGood();      Bool_t IsGood(Bool_t strict=true);
546    
547      Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)      Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)
548    
549        Bool_t IsGP(){return ISGP;};
550    
551        void   SetGP(Bool_t tt){ISGP=tt;};
552    
553        void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
554    
555      ClassDef(PamLevel2,1);      ClassDef(PamLevel2,1);
556  };  };
557    
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.39
  ViewVC Help
Powered by ViewVC 1.1.23