/[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.18 by mocchiut, Wed Mar 21 13:52:14 2007 UTC revision 1.31 by pam-fi, Wed Nov 14 11:14:27 2007 UTC
# Line 31  Line 31 
31  #include <TTree.h>  #include <TTree.h>
32  #include <TFriendElement.h>  #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>
# Line 38  Line 39 
39  #include <TSystemDirectory.h>  #include <TSystemDirectory.h>
40  #include <TSQLServer.h>  #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  #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   *   *
# Line 58  class PamTrack : public TObject { Line 222  class PamTrack : public TObject {
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;      bool candeleteobj;
228    
229  //    TClonesArray *Cluster;        ///< clusters (if stored)  //    TClonesArray *Cluster;        ///< clusters (if stored)
# Line 74  class PamTrack : public TObject { Line 239  class PamTrack : public TObject {
239      void Clear(Option_t *option){Clear();};      void Clear(Option_t *option){Clear();};
240      void Delete();      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;};      TrkTrack*   GetTrkTrack(){return trk_track;};
250      CaloTrkVar* GetCaloTrack(){return calo_track;};      CaloTrkVar* GetCaloTrack(){return calo_track;};
251      ToFTrkVar*  GetToFTrack(){return tof_track;};      ToFTrkVar*  GetToFTrack(){return tof_track;};
# Line 96  class PamLevel2 : public TObject { Line 268  class PamLevel2 : public TObject {
268      // ------------------------------      // ------------------------------
269      // variables to retrieve run info      // variables to retrieve run info
270      // ------------------------------      // ------------------------------
 //    Int_t     irun;          //current run entry, relative to the TChain  
271      Long64_t  irun;          //current run entry, relative to the TChain      Long64_t  irun;          //current run entry, relative to the TChain
272  //    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  
273      Long64_t  runfirstentry; //first entry of run, relative to the TChain      Long64_t  runfirstentry; //first entry of run, relative to the TChain
274      Long64_t  runlastentry;  //last  entry of run, relative to the TChain      Long64_t  runlastentry;  //last  entry of run, relative to the TChain
275      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 277  class PamLevel2 : public TObject {
277      TTree*    run_tree_clone;      TTree*    run_tree_clone;
278      TChain*   sel_tree;      TChain*   sel_tree;
279      TTree*    sel_tree_clone;      TTree*    sel_tree_clone;
280        
281        // -------------------------------------------
282        // variables to store total live and dead time
283        // -------------------------------------------
284        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events  
285    
286      // ---------------------------------      // ---------------------------------
287      // variables to retrieve level0 info      // variables to retrieve level0 info
288      // ---------------------------------      // ---------------------------------
# Line 116  class PamLevel2 : public TObject { Line 290  class PamLevel2 : public TObject {
290      TTree*      l0_tree;      TTree*      l0_tree;
291      Int_t       iroot;            Int_t       iroot;      
292      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
293        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
294      // --------------------------------------------      // --------------------------------------------
295      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
296      // --------------------------------------------      // --------------------------------------------
# Line 138  class PamLevel2 : public TObject { Line 313  class PamLevel2 : public TObject {
313      TTree*  pam_tree_clone[NCLONES];      TTree*  pam_tree_clone[NCLONES];
314    
315      // --------------------------------------------      // --------------------------------------------
316      // variables to store Level0 info (if required)      // variables to store Level1-2 info
317      // --------------------------------------------      // --------------------------------------------
318      TrkLevel1*   trk1_obj;      TrkLevel1*   trk1_obj;
319      TrkLevel2*   trk2_obj;      TrkLevel2*   trk2_obj;
# Line 154  class PamLevel2 : public TObject { Line 329  class PamLevel2 : public TObject {
329      GL_RUN*      run_obj;      GL_RUN*      run_obj;
330      SoftInfo*    soft_obj; // Emiliano      SoftInfo*    soft_obj; // Emiliano
331    
332        // --------------------------------------------
333        // variables to handle GP
334        // --------------------------------------------
335        GPamela*     gp_obj;
336    
337      // -----------------------------------      // -----------------------------------
338      // variables to sort the track      // variables to sort the track
339      // -----------------------------------      // -----------------------------------
340      TClonesArray *tsorted;      TClonesArray *tsorted;
341      TClonesArray *timage;      TClonesArray *timage;
342      PamTrack    *track; //store the last required      PamTrack    *track; //store the last required
343        TString howtosort;
344        Float_t sortthr;
345    
346        // --------------------------------------
347      // -----      // flags to load level0-1-2 trees/braches
348      // flags      // --------------------------------------
     // -----  
349      Bool_t CAL2; //level2 branch      Bool_t CAL2; //level2 branch
350      Bool_t CAL1; //level1 branch      Bool_t CAL1; //level1 branch
351      Bool_t CAL0; //yoda      Bool_t CAL0; //yoda
# Line 183  class PamLevel2 : public TObject { Line 364  class PamLevel2 : public TObject {
364      Bool_t AC;      Bool_t AC;
365      Bool_t ORB;      Bool_t ORB;
366    
367      Bool_t RUN;      Bool_t GP; //GP simulation tree
368    
369      Int_t  SELLI; //trees of selected events (there is a selection list)      // -----------
370        // other flags
371        // -----------
372        Bool_t RUN;   // RunInfo tree
373        Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
374        Int_t  SELLI; // tree of selected events (there is a selection list)
375    
376      void Initialize();      void Initialize();
377    
# Line 209  class PamLevel2 : public TObject { Line 395  class PamLevel2 : public TObject {
395      void    GetWhichTrees(TFile*);      void    GetWhichTrees(TFile*);
396    
397      TTree*  GetPamTree(TFile*,TString);      TTree*  GetPamTree(TFile*,TString);
398      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};    
399    
400      TList*  GetListOfLevel2Files(TString, TString);      TList*  GetListOfLevel2Files(TString, TString);
401      Bool_t  CheckLevel2File(TString);      Bool_t  CheckLevel2File(TString);
402      TChain* GetPamTree(TList*,TString);      TChain* GetPamTree(TList*,TString);
403      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); };
404      TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };      TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
405        TChain* AddFriend(TString);
406    
407      TTree*  GetRunTree(TFile*);      TTree*  GetRunTree(TFile*);
408      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
409      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
410        Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
411      Bool_t  UpdateRunInfo(TTree *,Long64_t);      Bool_t  UpdateRunInfo(TTree *,Long64_t);
412      Bool_t  UpdateRunInfo(TChain *,Long64_t);      Bool_t  UpdateRunInfo(TChain *,Long64_t);
     Bool_t  UpdateRunInfo(Long64_t iev){ return UpdateRunInfo(run_tree,iev); };  
413    
414      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
415      void    CreateCloneTrees(TFile* file);      void    CreateCloneTrees(TFile* file);
416        void    CreateCloneTrees(){CreateCloneTrees(NULL);};
417      void    FillCloneTrees();      void    FillCloneTrees();
418      TTree  *GetCloneTree(TString);      TTree  *GetCloneTree(TString);
419      void    WriteCloneTrees();      void    WriteCloneTrees();
# Line 236  class PamLevel2 : public TObject { Line 424  class PamLevel2 : public TObject {
424  //    Int_t   GetEntry(Int_t);  //    Int_t   GetEntry(Int_t);
425      Int_t   GetEntry(Long64_t);      Int_t   GetEntry(Long64_t);
426      Int_t   GetYodaEntry();      Int_t   GetYodaEntry();
427    
428        void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
429            
430        Long64_t GetEntries(){return pam_tree->GetEntries();};
431      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
432      TChain* GetRunTree(){return run_tree;};      TChain* GetRunTree(){return run_tree;};
433      TTree*  GetYodaTree();      TTree*  GetYodaTree();
# Line 249  class PamLevel2 : public TObject { Line 440  class PamLevel2 : public TObject {
440      void SetDBHost(TString str){host=str;};      void SetDBHost(TString str){host=str;};
441      void SetDBUser(TString str){user=str;};      void SetDBUser(TString str){user=str;};
442      void SetDBPsw(TString str){psw=str;};      void SetDBPsw(TString str){psw=str;};
443        TSQLServer* GetSQLServer(){return dbc;};
444        Bool_t SetDBConnection();
445    
446  //  ========================================================  //  ========================================================
447  //  methods needed to retrieve subdetector general variables  //  methods needed to retrieve subdetector general variables
448  //  ========================================================  //  ========================================================
449      void*       GetPointerTo(const char*);      void*       GetPointerTo(const char*);
450      EventHeader  *GetEventHeader()  { return h0_obj; };      EventHeader  *GetEventHeader()  { return h0_obj; };
451      TrkLevel0    *GetTrkLevel0()  { return trk0_obj; };  
452      CaloLevel0   *GetCaloLevel0()  { return calo0_obj; };  /*     TrkLevel0    *GetTrkLevel0()   { return trk0_obj; }; */
453    /*     CaloLevel0   *GetCaloLevel0()  { return calo0_obj; }; */
454        TrkLevel0    *GetTrkLevel0();
455        CaloLevel0   *GetCaloLevel0();
456    
457    
458      TrkLevel1    *GetTrkLevel1()  { return trk1_obj; };      TrkLevel1    *GetTrkLevel1()  { return trk1_obj; };
459      TrkLevel2    *GetTrkLevel2()  { return trk2_obj; };      TrkLevel2    *GetTrkLevel2()  { return trk2_obj; };
# Line 272  class PamLevel2 : public TObject { Line 468  class PamLevel2 : public TObject {
468      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
469      GL_RUN       *GetRunInfo()    { return run_obj;  };      GL_RUN       *GetRunInfo()    { return run_obj;  };
470      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
471        GPamela      *GetGPamela()    { return gp_obj; };
472    
473  //  ==============================================  //  ==============================================
474  //  methods to retrieve subdetectors stored tracks  //  methods to retrieve subdetectors stored tracks
# Line 290  class PamLevel2 : public TObject { Line 487  class PamLevel2 : public TObject {
487      PamTrack*     GetStoredTrack(Int_t);      PamTrack*     GetStoredTrack(Int_t);
488      //      //
489      void          SortTracks(TString);          //assign value to variable sorted_tracks      void          SortTracks(TString);          //assign value to variable sorted_tracks
490        void          SortTracks();          //assign value to variable sorted_tracks
491  //    TRefArray*    GetTracks();                  //return sorted_tracks  //    TRefArray*    GetTracks();                  //return sorted_tracks
492      TClonesArray* GetTracks();                  //return sorted_tracks      TClonesArray* GetTracks();                  //return sorted_tracks
493      PamTrack*     GetTrack(int);      PamTrack*     GetTrack(int);
494      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImage(int);
495        void          SetSortingMethod(TString how){ howtosort = how; return;};
496        TString       GetSortingMethod(){return howtosort;};
497        void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
498        Float_t       GetNucleiSortingThreshold(){ return sortthr;};
499                    
500  //  =================================  //  =================================
501  //  other methods  //  other methods
# Line 323  class PamLevel2 : public TObject { Line 525  class PamLevel2 : public TObject {
525    
526      Bool_t IsRUN(){return RUN;};      Bool_t IsRUN(){return RUN;};
527    
528        Bool_t IsGood();
529    
530      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)
531    
532        Bool_t IsGP(){return ISGP;};
533    
534        void   SetGP(Bool_t tt){ISGP=tt;};
535    
536    
537      ClassDef(PamLevel2,1);      ClassDef(PamLevel2,1);
538  };  };
539    
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.23