/[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.32 by pam-fi, Tue Nov 20 18:16:29 2007 UTC revision 1.41 by pam-fi, Wed Apr 1 09:25:36 2009 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 55  using namespace std; Line 61  using namespace std;
61  /////////////////////////////////////////////////////////////////////////////  /////////////////////////////////////////////////////////////////////////////
62  /**  /**
63   * \brief Class to describe GP output (to be moved to PamelaDigitizer)   * \brief Class to describe GP output (to be moved to PamelaDigitizer)
64   *   *
65   */   */
66    
67  class GPamela : public TObject {  class GPamela : public TObject {
# Line 211  class GPamela : public TObject { Line 217  class GPamela : public TObject {
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 {  class PamTrack : public TObject {
# Line 222  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;      Int_t pscore;
233      Int_t iscore;      Int_t iscore;
234      bool candeleteobj;      bool candeleteobj;
# Line 231  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 249  class PamTrack : public TObject { Line 256  class PamTrack : public TObject {
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  /**  /**
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 CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {  //class PamLevel2 : public CaloLevel2, public CaloLevel1, public ToFLevel2, public TrkLevel2, public TrkLevel1, public TrkHough, public TrigLevel2, public S4Level2, public NDLevel2, public AcLevel2, public OrbitalInfo {
272  class PamLevel2 : public TObject {  class PamLevel2 : public TObject {
273    
# Line 278  class PamLevel2 : public TObject { Line 286  class PamLevel2 : public TObject {
286      TChain*   sel_tree;      TChain*   sel_tree;
287      TTree*    sel_tree_clone;      TTree*    sel_tree_clone;
288      int *irunoffset;      int *irunoffset;
289        GL_TIMESYNC *gltsync; // Emiliano
290    
291      // -------------------------------------------      // -------------------------------------------
292      // variables to store total live and dead time      // variables to store total live and dead time
293      // -------------------------------------------      // -------------------------------------------
294      ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events
295    
296      // ---------------------------------      // ---------------------------------
297      // variables to retrieve level0 info      // variables to retrieve level0 info
298      // ---------------------------------      // ---------------------------------
299      TFile*      l0_file;      TFile*      l0_file;
300      TTree*      l0_tree;      TTree*      l0_tree;
301      Int_t       iroot;            Int_t       iroot;
302      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
303      Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10      Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
304      // --------------------------------------------      // --------------------------------------------
305      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
306      // --------------------------------------------      // --------------------------------------------
# Line 314  class PamLevel2 : public TObject { Line 323  class PamLevel2 : public TObject {
323      TTree*  pam_tree_clone[NCLONES];      TTree*  pam_tree_clone[NCLONES];
324    
325      // --------------------------------------------      // --------------------------------------------
326      // variables to store Level1-2 info      // variables to store Level1-2 info
327      // --------------------------------------------      // --------------------------------------------
328      TrkLevel1*   trk1_obj;      TrkLevel1*   trk1_obj;
329      TrkLevel2*   trk2_obj;      TrkLevel2*   trk2_obj;
# Line 340  class PamLevel2 : public TObject { Line 349  class PamLevel2 : public TObject {
349      // -----------------------------------      // -----------------------------------
350      TClonesArray *tsorted;      TClonesArray *tsorted;
351      TClonesArray *timage;      TClonesArray *timage;
352      PamTrack    *track; //store the last required      PamTrack    *track; //store the last required
353      TString howtosort;      TString howtosort;
354      Float_t sortthr;      Float_t sortthr;
355    
356        Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria
357        Long64_t lastsorted; // the last event whose tracks has been sorted
358      // --------------------------------------      // --------------------------------------
359      // flags to load level0-1-2 trees/braches      // flags to load level0-1-2 trees/braches
360      // --------------------------------------      // --------------------------------------
# Line 374  class PamLevel2 : public TObject { Line 385  class PamLevel2 : public TObject {
385      Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)      Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
386      Int_t  SELLI; // tree of selected events (there is a selection list)      Int_t  SELLI; // tree of selected events (there is a selection list)
387    
388        Bool_t DBG;
389    
390      void Initialize();      void Initialize();
391    
392    
393   public:   public:
394            
395        PamLevel2(TString ddir,TList *list,TString detlist);
396      PamLevel2(TString ddir,TString list,TString detlist);      PamLevel2(TString ddir,TString list,TString detlist);
397      PamLevel2(TString ddir,TString list);      PamLevel2(TString ddir,TString list);
398      PamLevel2();      PamLevel2();
399      ~PamLevel2(){Delete();};      ~PamLevel2(){Delete();};
400            
401      void Clear();      void Clear();
402      void Reset();      void Reset();
403      void Delete();      void Delete();
404    
405            
406  //  =============================================  //  =============================================
407  //  methods needed to read/write Level0-1-2 trees  //  methods needed to read/write Level0-1-2 trees
408  //  =============================================  //  =============================================
# Line 396  class PamLevel2 : public TObject { Line 410  class PamLevel2 : public TObject {
410      void    GetWhichTrees(TFile*);      void    GetWhichTrees(TFile*);
411    
412      TTree*  GetPamTree(TFile*,TString);      TTree*  GetPamTree(TFile*,TString);
413      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};          TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};
414    
415      TList*  GetListOfLevel2Files(TString, TString);      TList*  GetListOfLevel2Files(TString, TString);
416      Bool_t  CheckLevel2File(TString);      Bool_t  CheckLevel2File(TString);
# Line 424  class PamLevel2 : public TObject { Line 438  class PamLevel2 : public TObject {
438    
439  //    Int_t   GetEntry(Int_t);  //    Int_t   GetEntry(Int_t);
440      Int_t   GetEntry(Long64_t);      Int_t   GetEntry(Long64_t);
441      Int_t   GetYodaEntry();      Long64_t GetReadEntry();
442        Int_t   GetYodaEntry();
443    
444      void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable      void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
445        
446      Long64_t GetEntries(){return pam_tree->GetEntries();};      Long64_t GetEntries(){return pam_tree->GetEntries();};
447      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
448      TChain* GetRunTree(){return run_tree;};      TChain* GetRunTree(){return run_tree;};
# Line 469  class PamLevel2 : public TObject { Line 484  class PamLevel2 : public TObject {
484      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
485      GL_RUN       *GetRunInfo()    { return run_obj;  };      GL_RUN       *GetRunInfo()    { return run_obj;  };
486      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
487      GPamela      *GetGPamela()    { return gp_obj; };      GPamela      *GetGPamela()    { return gp_obj;   };
488    
489  //  ==============================================  //  ==============================================
490  //  methods to retrieve subdetectors stored tracks  //  methods to retrieve subdetectors stored tracks
# Line 477  class PamLevel2 : public TObject { Line 492  class PamLevel2 : public TObject {
492      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
493      CaloTrkVar* GetCaloStoredTrack(int);      CaloTrkVar* GetCaloStoredTrack(int);
494      ToFTrkVar*  GetToFStoredTrack(int);      ToFTrkVar*  GetToFStoredTrack(int);
495        OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);
496      //      //
497      PamTrack*   GetPamTrackAlong(TrkTrack*);      PamTrack*   GetPamTrackAlong(TrkTrack*);
498  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
499  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
500            
501  //  =================================  //  =================================
502  //  methods to retrieve pamela tracks  //  methods to retrieve pamela tracks
503  //  =================================  //  =================================
# Line 493  class PamLevel2 : public TObject { Line 509  class PamLevel2 : public TObject {
509      TClonesArray* GetTracks();                  //return sorted_tracks      TClonesArray* GetTracks();                  //return sorted_tracks
510      PamTrack*     GetTrack(int);      PamTrack*     GetTrack(int);
511      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImage(int);
512      void          SetSortingMethod(TString how){ howtosort = how; return;};      void          SetSortingMethod(TString how);
513      TString       GetSortingMethod(){return howtosort;};      TString       GetSortingMethod(){return howtosort;};
514      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
515      Float_t       GetNucleiSortingThreshold(){ return sortthr;};      Float_t       GetNucleiSortingThreshold(){ return sortthr;};
516            
517  //  =================================  //  =================================
518  //  other methods  //  other methods
519  //  =================================  //  =================================
520  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
521    
522    
523        Long64_t  GetRunID(){ return irun;};
524        Long64_t  GetRunEntryID(){ return irunentry; };
525    
526      // -----      // -----
527      // flags methods      // flags methods
528      // -----      // -----
# Line 526  class PamLevel2 : public TObject { Line 546  class PamLevel2 : public TObject {
546    
547      Bool_t IsRUN(){return RUN;};      Bool_t IsRUN(){return RUN;};
548    
549      Bool_t IsGood();      Bool_t IsGood(Bool_t strict=true);
550    
551      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)
552        void  SetSELLI(Int_t selli){SELLI=selli;};
553    
554      Bool_t IsGP(){return ISGP;};      Bool_t IsGP(){return ISGP;};
555    
556      void   SetGP(Bool_t tt){ISGP=tt;};      void   SetGP(Bool_t tt){ISGP=tt;};
557    
558        void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
559    
560    
561      ClassDef(PamLevel2,1);      ClassDef(PamLevel2,1);
562  };  };
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.23