/[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.1 by pam-fi, Fri Jun 16 16:43:55 2006 UTC revision 1.44 by mocchiut, Mon Dec 12 16:11:41 2011 UTC
# Line 7  Line 7 
7  #ifndef pamlevel2_h  #ifndef pamlevel2_h
8  #define pamlevel2_h  #define pamlevel2_h
9    
10  #include <TClonesArray.h>  #include <TRefArray.h>
11    
12    #include <TrkStruct.h>
13    #include <TrkParams.h>
14    #include <TrkLevel0.h>
15    #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>
23  #include <S4Level2.h>  #include <S4Level2.h>
24  #include <NDLevel2.h>  #include <NDLevel2.h>
25  #include <AcLevel2.h>  #include <AcLevel2.h>
26    #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>
34    #include <GLTables.h>
35    
36    
37  #include <TTree.h>  #include <TTree.h>
38    #include <TFriendElement.h>
39    #include <TChain.h>
40    #include <TChainElement.h>
41  #include <TFile.h>  #include <TFile.h>
42    #include <TList.h>
43    #include <TKey.h>
44    #include <TSystemFile.h>
45    #include <TSystemDirectory.h>
46    #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
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   *   *
220   * This class inherites from TrkTrack, CaloTrkVar ecc...   * This class inherites from TrkTrack, CaloTrkVar ecc...
221   *   *
222   */   */
223  class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {  //class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {
224    class PamTrack : public TObject {
225    
226   private:   private:
227    
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;
235    
236    //    TClonesArray *Cluster;        ///< clusters (if stored)
237    
238   public:   public:
239    
240      PamTrack();      PamTrack();
241        PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
242        PamTrack(const PamTrack&);
243        ~PamTrack(){ Delete(); };
244    
245        void Clear();
246        void Clear(Option_t *option){Clear();};
247        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;};
257        CaloTrkVar* GetCaloTrack(){return calo_track;};
258        ToFTrkVar*  GetToFTrack(){return tof_track;};
259        OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};
260    
261      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );      ClassDef(PamTrack,2);
   
     ClassDef(PamTrack,1);  
262  };  };
263    
264  /**  /**
265   * \brief Class to describe PAMELA LEVEL2 data.   * \brief Class to describe PAMELA LEVEL2 data.
266   *   *
267   * This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.   * This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.
268   * 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.
269   */   */
270                    
271  class PamLevel2 : public CaloLevel2, public TrkLevel2, public ToFLevel2, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2 {  //class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {
272    class PamLevel2 : public TObject {
273    
274   private:   private:
275    
276      TrkLevel2*  trk_obj;      // ------------------------------
277      CaloLevel2* calo_obj;      // variables to retrieve run info
278      ToFLevel2*  tof_obj;      // ------------------------------
279      TrigLevel2* trig_obj;      Long64_t  irun;          //current run entry, relative to the TChain
280      S4Level2*   s4_obj;      Long64_t  irunt;         //current run entry, relative to the TTree in the TChain
281      NDLevel2*   nd_obj;      Long64_t  runfirstentry; //first entry of run, relative to the TChain
282      AcLevel2*   ac_obj;      Long64_t  runlastentry;  //last  entry of run, relative to the TChain
283            Long64_t  irunentry;     //position of current entry, relative to the run
284        Long64_t  totrunentry;
285        Long64_t  totrunentrymax;
286        Long64_t  totrunentrymin;
287        ULong64_t    prevabstime; // abs time of previous event
288        UInt_t    prevpktnum; // packet number of previous event
289        ULong64_t    abstime; // current abs time
290        UInt_t    pktnum; // current packet number
291        Bool_t isFragment;
292        TChain*   run_tree;
293        TTree*    run_tree_clone;
294        TChain*   sel_tree;
295        TTree*    sel_tree_clone;
296        int *irunoffset;
297        GL_TIMESYNC *gltsync; // Emiliano
298        Bool_t fUpdateRunInfo; // Emiliano
299    
300        // -------------------------------------------
301        // variables to store total live and dead time
302        // -------------------------------------------
303        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events
304    
305        // ---------------------------------
306        // variables to retrieve level0 info
307        // ---------------------------------
308        TFile*      l0_file;
309        TTree*      l0_tree;
310        Int_t       iroot;
311        Int_t       prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time
312        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
313        // --------------------------------------------
314        // variables to store Level0 info (if required)
315        // --------------------------------------------
316        EventHeader* h0_obj;
317        TrkLevel0*   trk0_obj;
318        CaloLevel0*  calo0_obj;
319    
320        // ---------------------------
321        // variables to connect the DB
322        // ---------------------------
323        TSQLServer* dbc;
324        TString host;
325        TString user;
326        TString psw;
327    
328        // ---------------------------------------------
329        // variables to retrieve and clone level1-2 info
330        // ---------------------------------------------
331        TChain* pam_tree;
332        TTree*  pam_tree_clone[NCLONES];
333    
334        // --------------------------------------------
335        // variables to store Level1-2 info
336        // --------------------------------------------
337        TrkLevel1*   trk1_obj;
338        TrkLevel2*   trk2_obj;
339        TrkHough*    trkh_obj;
340        CaloLevel1*  calo1_obj;
341        CaloLevel2*  calo2_obj;
342        ToFLevel2*   tof_obj;
343        TrigLevel2*  trig_obj;
344        S4Level2*    s4_obj;
345        NDLevel2*    nd_obj;
346        AcLevel2*    ac_obj;
347        OrbitalInfo* orb_obj;
348        GL_RUN*      run_obj;
349        SoftInfo*    soft_obj; // Emiliano
350    
351        // --------------------------------------------
352        // variables to handle GP
353        // --------------------------------------------
354        GPamela*     gp_obj;
355    
356        // -----------------------------------
357        // variables to sort the track
358        // -----------------------------------
359        TClonesArray *tsorted;
360        TClonesArray *timage;
361        PamTrack    *track; //store the last required
362        TString howtosort;
363        Float_t sortthr;
364    
365        Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria
366        Long64_t lastsorted; // the last event whose tracks has been sorted
367        // --------------------------------------
368        // flags to load level0-1-2 trees/braches
369        // --------------------------------------
370        Bool_t CAL2; //level2 branch
371        Bool_t CAL1; //level1 branch
372        Bool_t CAL0; //yoda
373    
374        Bool_t TRK2; //level2 branch
375        Bool_t TRK1; //level1 branch
376        Bool_t TRKh; //hough branch
377        Bool_t TRK0; //yoda
378    
379        Bool_t TOF;
380        Bool_t TOF0;
381    
382        Bool_t TRG;
383        Bool_t S4;
384        Bool_t ND;
385        Bool_t AC;
386        Bool_t ORB;
387    
388        Bool_t GP; //GP simulation tree
389    
390        // -----------
391        // other flags
392        // -----------
393        Bool_t RUN;   // RunInfo tree
394        Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
395        Int_t  SELLI; // tree of selected events (there is a selection list)
396    
397        Bool_t DBG;
398    
399        void Initialize();
400    
401    
402   public:   public:
403    
404        PamLevel2(TString ddir,TList *list,TString detlist);
405        PamLevel2(TString ddir,TString list,TString detlist);
406        PamLevel2(TString ddir,TString list);
407      PamLevel2();      PamLevel2();
408        ~PamLevel2(){Delete();};
409    
410        void Clear();
411        void Reset();
412        void Delete();
413    
414    
415    //  =============================================
416    //  methods needed to read/write Level0-1-2 trees
417    //  =============================================
418        void    SetWhichTrees(TString);
419        void    GetWhichTrees(TFile*);
420    
421        TTree*  GetPamTree(TFile*,TString);
422        TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};
423    
424        TList*  GetListOfLevel2Files(TString, TString);
425        Bool_t  CheckLevel2File(TString);
426        TChain* GetPamTree(TList*,TString);
427        TChain* GetPamTree(TString ddir,TString list,TString detlist){ return GetPamTree(GetListOfLevel2Files(ddir,list),detlist); };
428        TChain* GetPamTree(TString ddir,TString list){ return GetPamTree(GetListOfLevel2Files(ddir,list),""); };
429        TChain* AddFriend(TString);
430    
431        TTree*  GetRunTree(TFile*);
432        TChain* GetRunTree(TList*);
433        TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
434        Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug
435        Bool_t  UpdateRunInfo_bckw_compat(Long64_t);///< backward compatibility method
436        Bool_t  UpdateRunInfo(TTree *,Long64_t);
437        Bool_t  UpdateRunInfo(TChain *,Long64_t);
438    
439        void    SkipRunInfoUpdate(); // Emiliano
440        void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano
441        void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano
442    
443        void    CreateCloneTrees0(TChain *fChain, TFile* file);
444        void    CreateCloneTrees(TFile* file);
445        void    CreateCloneTrees(){CreateCloneTrees(NULL);};
446        void    FillCloneTrees();
447        TTree  *GetCloneTree(TString);
448        void    WriteCloneTrees();
449    
450        void    SetBranchAddress(TTree *);
451        void    SetBranchAddress(TChain *);
452    
453    //    Int_t   GetEntry(Int_t);
454        Int_t   GetEntry(Long64_t);
455        Long64_t GetReadEntry();
456        Int_t   GetYodaEntry();
457    
458        void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
459    
460        Long64_t GetEntries(){return pam_tree->GetEntries();};
461        TChain* GetPamTree(){return pam_tree;};
462        TChain* GetRunTree(){return run_tree;};
463        TTree*  GetYodaTree();
464    
465        TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things
466    
467    //  =============
468    //  DB connection
469    //  =============
470        void SetDBHost(TString str){host=str;};
471        void SetDBUser(TString str){user=str;};
472        void SetDBPsw(TString str){psw=str;};
473        TSQLServer* GetSQLServer(){return dbc;};
474        Bool_t SetDBConnection();
475    
476    //  ========================================================
477    //  methods needed to retrieve subdetector general variables
478    //  ========================================================
479        void*       GetPointerTo(const char*);
480        EventHeader  *GetEventHeader()  { return h0_obj; };
481    
482    /*     TrkLevel0    *GetTrkLevel0()   { return trk0_obj; }; */
483    /*     CaloLevel0   *GetCaloLevel0()  { return calo0_obj; }; */
484        TrkLevel0    *GetTrkLevel0();
485        CaloLevel0   *GetCaloLevel0();
486    
487    
488        TrkLevel1    *GetTrkLevel1()  { return trk1_obj; };
489        TrkLevel2    *GetTrkLevel2()  { return trk2_obj; };
490        TrkHough     *GetTrkHough()   { return trkh_obj; };
491        CaloLevel1   *GetCaloLevel1() { return calo1_obj;};
492        CaloLevel2   *GetCaloLevel2() { return calo2_obj;};
493        ToFLevel2    *GetToFLevel2()  { return tof_obj;  };
494        TrigLevel2   *GetTrigLevel2() { return trig_obj; };
495        S4Level2     *GetS4Level2()   { return s4_obj;   };
496        NDLevel2     *GetNDLevel2()   { return nd_obj;   };
497        AcLevel2     *GetAcLevel2()   { return ac_obj;   };
498        OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
499        GL_RUN       *GetRunInfo()    { return run_obj;  };
500        SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
501        GPamela      *GetGPamela()    { return gp_obj;   };
502    
503    //  ==============================================
504    //  methods to retrieve subdetectors stored tracks
505    //  ==============================================
506        TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
507        CaloTrkVar* GetCaloStoredTrack(int);
508        ToFTrkVar*  GetToFStoredTrack(int);
509        OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);
510        //
511        PamTrack*   GetPamTrackAlong(TrkTrack*);
512    //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
513    //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
514    
     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;};  
       
 //  =================================  
 //  methods needed to ... well I am not shure they are needed  
 //  =================================  
     void*       GetPointerToTrk()  {return &trk_obj;};  
     void*       GetPointerToCalo() {return &calo_obj;};  
     void*       GetPointerToToF()  {return &tof_obj;};  
     void*       GetPointerToTrig() {return &trig_obj;};  
     void*       GetPointerToS4() {return &s4_obj;};  
     void*       GetPointerToND() {return &nd_obj;};  
     void*       GetPointerToAc() {return &ac_obj;};  
   
     TTree* LoadPamTrees(TFile*);  
       
       
515  //  =================================  //  =================================
516  //  methods to retrieve pamela tracks  //  methods to retrieve pamela tracks
517  //  =================================  //  =================================
518      PamTrack*     GetStoredTrack(Int_t itrk);      PamTrack*     GetStoredTrack(Int_t);
519      TClonesArray* GetTracks();      //
520        void          SortTracks(TString);          //assign value to variable sorted_tracks
521        void          SortTracks();          //assign value to variable sorted_tracks
522    //    TRefArray*    GetTracks();                  //return sorted_tracks
523        TClonesArray* GetTracks();                  //return sorted_tracks
524      PamTrack*     GetTrack(int);      PamTrack*     GetTrack(int);
525      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImage(int);
526      ClassDef(PamLevel2,1);      void          SetSortingMethod(TString how);
527        TString       GetSortingMethod(){return howtosort;};
528        void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
529        Float_t       GetNucleiSortingThreshold(){ return sortthr;};
530    
531    //  =================================
532    //  other methods
533    //  =================================
534    //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
535    
536    
537        Long64_t  GetRunID(){ return irun;};
538        Long64_t  GetRunEntryID(){ return irunentry; };
539    
540        // -----
541        // flags methods
542        // -----
543        Bool_t IsCAL2(){return CAL2;}; //level2 branch
544        Bool_t IsCAL1(){return CAL1;}; //level1 branch
545        Bool_t IsCAL0(){return CAL0;}; //yoda
546    
547        Bool_t IsTRK2(){return TRK2;}; //level2 branch
548        Bool_t IsTRK1(){return TRK1;}; //level1 branch
549        Bool_t IsTRKh(){return TRKh;}; //hough branch
550        Bool_t IsTRK0(){return TRK0;}; //yoda
551    
552        Bool_t IsTOF(){return TOF;};
553        Bool_t IsTOF0(){return TOF0;};
554    
555        Bool_t IsTRG(){return TRG;};
556        Bool_t IsS4(){return S4;};
557        Bool_t IsND(){return ND;};
558        Bool_t IsAC(){return AC;};
559        Bool_t IsORB(){return ORB;};
560    
561        Bool_t IsRUN(){return RUN;};
562    
563        Bool_t IsGood(Bool_t strict=true);
564    
565        Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)
566        void  SetSELLI(Int_t selli){SELLI=selli;};
567    
568        Bool_t IsGP(){return ISGP;};
569    
570        void   SetGP(Bool_t tt){ISGP=tt;};
571    
572        void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
573    
574    
575        ClassDef(PamLevel2,2);
576  };  };
577    
578  //===============================================================================================  //===============================================================================================
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.44
  ViewVC Help
Powered by ViewVC 1.1.23