/[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.40 by mocchiut, Fri Dec 12 14:11:47 2008 UTC revision 1.47 by mocchiut, Wed Mar 12 16:37:26 2014 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  #include <InclinationInfo.h> // EMILIANO, now this library is included in libDarthVader.so
28                           // and the header installed in $PAM_INC. If you get the                           // and the header installed in $PAM_INC. If you get the
29                           // error "Cannot find InclinationInfo.h" update your DV installation                           // error "Cannot find InclinationInfo.h" update your DV installation
30  #include <OrientationInfo.h> // EMILIANO, now this library is included in libDarthVader.so  #include <OrientationInfo.h> // EMILIANO, now this library is included in libDarthVader.so
31                           // and the header installed in $PAM_INC. If you get the                           // and the header installed in $PAM_INC. If you get the
32                           // error "Cannot find InclinationInfo.h" update your DV installation                           // error "Cannot find InclinationInfo.h" update your DV installation
33  #include <RunInfo.h>  #include <RunInfo.h>
34  #include <GLTables.h>  #include <GLTables.h>
# Line 61  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 217  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 264  class PamTrack : public TObject { Line 264  class PamTrack : public TObject {
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 281  class PamLevel2 : public TObject { Line 281  class PamLevel2 : public TObject {
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
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        Bool_t isSync; ///< the file is fully sinchronized between L2 and L0/DB data
293            //    Bool_t hasL0EE; ///< backward compatibility flag for preselected event list
294        Long64_t il0entry; ///< entry of level0 file in selection list
295      TChain*   run_tree;      TChain*   run_tree;
296      TTree*    run_tree_clone;      TTree*    run_tree_clone;
297      TChain*   sel_tree;      TChain*   sel_tree;
298      TTree*    sel_tree_clone;      TTree*    sel_tree_clone;
299      int *irunoffset;      int *irunoffset;
300      GL_TIMESYNC *gltsync; // Emiliano      GL_TIMESYNC *gltsync; // Emiliano
301        Bool_t fUpdateRunInfo; // Emiliano
302        Bool_t fUseDBinRunInfo; // Emiliano
303    
304      // -------------------------------------------      // -------------------------------------------
305      // variables to store total live and dead time      // variables to store total live and dead time
306      // -------------------------------------------      // -------------------------------------------
307      ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events        ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events
308    
309      // ---------------------------------      // ---------------------------------
310      // variables to retrieve level0 info      // variables to retrieve level0 info
311      // ---------------------------------      // ---------------------------------
312      TFile*      l0_file;      TFile*      l0_file;
313      TTree*      l0_tree;      TTree*      l0_tree;
314      Int_t       iroot;            Int_t       iroot;
315      Int_t       prevshift; // to remember if any shift is needed to synchronize l0 and l2 data instead of calculating any time      Long64_t       prevshift; // shift needed to sync level2 data with the DB (nested+DV_skipped packets) - global position in level2 file EM
316      Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10      Long64_t       yprevshift; // shift needed to sync level2 data with level0 data (nested+DV_skipped packets) - run wise EM
317        Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
318      // --------------------------------------------      // --------------------------------------------
319      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
320      // --------------------------------------------      // --------------------------------------------
# Line 323  class PamLevel2 : public TObject { Line 337  class PamLevel2 : public TObject {
337      TTree*  pam_tree_clone[NCLONES];      TTree*  pam_tree_clone[NCLONES];
338    
339      // --------------------------------------------      // --------------------------------------------
340      // variables to store Level1-2 info      // variables to store Level1-2 info
341      // --------------------------------------------      // --------------------------------------------
342      TrkLevel1*   trk1_obj;      TrkLevel1*   trk1_obj;
343      TrkLevel2*   trk2_obj;      TrkLevel2*   trk2_obj;
# Line 349  class PamLevel2 : public TObject { Line 363  class PamLevel2 : public TObject {
363      // -----------------------------------      // -----------------------------------
364      TClonesArray *tsorted;      TClonesArray *tsorted;
365      TClonesArray *timage;      TClonesArray *timage;
366      PamTrack    *track; //store the last required      PamTrack    *track; //store the last required
367      TString howtosort;      TString howtosort;
368      Float_t sortthr;      Float_t sortthr;
369    
370        Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria
371        Long64_t lastsorted; // the last event whose tracks has been sorted
372      // --------------------------------------      // --------------------------------------
373      // flags to load level0-1-2 trees/braches      // flags to load level0-1-2 trees/braches
374      // --------------------------------------      // --------------------------------------
# Line 389  class PamLevel2 : public TObject { Line 405  class PamLevel2 : public TObject {
405    
406    
407   public:   public:
408            
409      PamLevel2(TString ddir,TList *list,TString detlist);      PamLevel2(TString ddir,TList *list,TString detlist);
410      PamLevel2(TString ddir,TString list,TString detlist);      PamLevel2(TString ddir,TString list,TString detlist);
411      PamLevel2(TString ddir,TString list);      PamLevel2(TString ddir,TString list);
412      PamLevel2();      PamLevel2();
413      ~PamLevel2(){Delete();};      ~PamLevel2(){Delete();};
414            
415      void Clear();      void Clear();
416      void Reset();      void Reset();
417      void Delete();      void Delete();
418    
419            
420  //  =============================================  //  =============================================
421  //  methods needed to read/write Level0-1-2 trees  //  methods needed to read/write Level0-1-2 trees
422  //  =============================================  //  =============================================
# Line 408  class PamLevel2 : public TObject { Line 424  class PamLevel2 : public TObject {
424      void    GetWhichTrees(TFile*);      void    GetWhichTrees(TFile*);
425    
426      TTree*  GetPamTree(TFile*,TString);      TTree*  GetPamTree(TFile*,TString);
427      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};          TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};
428    
429      TList*  GetListOfLevel2Files(TString, TString);      TList*  GetListOfLevel2Files(TString, TString);
430      Bool_t  CheckLevel2File(TString);      Bool_t  CheckLevel2File(TString);
# Line 420  class PamLevel2 : public TObject { Line 436  class PamLevel2 : public TObject {
436      TTree*  GetRunTree(TFile*);      TTree*  GetRunTree(TFile*);
437      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
438      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
439      Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };      Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug
440      Bool_t  UpdateRunInfo(TTree *,Long64_t);      Bool_t  UpdateRunInfo(TTree *,Long64_t);
441      Bool_t  UpdateRunInfo(TChain *,Long64_t);      Bool_t  UpdateRunInfo(TChain *,Long64_t);
442    
443        void    SkipRunInfoUpdate(); // Emiliano
444        void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano
445        void    UseDBinRunInfo(){fUseDBinRunInfo = true;}; // [default] Emiliano
446        void    NoDBinRunInfo(){fUseDBinRunInfo = false;}; // EM
447        void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano
448    
449      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
450      void    CreateCloneTrees(TFile* file);      void    CreateCloneTrees(TFile* file);
451      void    CreateCloneTrees(){CreateCloneTrees(NULL);};      void    CreateCloneTrees(){CreateCloneTrees(NULL);};
# Line 436  class PamLevel2 : public TObject { Line 458  class PamLevel2 : public TObject {
458    
459  //    Int_t   GetEntry(Int_t);  //    Int_t   GetEntry(Int_t);
460      Int_t   GetEntry(Long64_t);      Int_t   GetEntry(Long64_t);
461      Int_t   GetYodaEntry();      Long64_t GetReadEntry();
462        Int_t   GetYodaEntry();
463    
464        void SetMaxShift(Int_t sh); // EMILIANO method to change the maxshift variable
465    
     void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable  
       
466      Long64_t GetEntries(){return pam_tree->GetEntries();};      Long64_t GetEntries(){return pam_tree->GetEntries();};
467      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
468      TChain* GetRunTree(){return run_tree;};      TChain* GetRunTree(){return run_tree;};
# Line 494  class PamLevel2 : public TObject { Line 517  class PamLevel2 : public TObject {
517      PamTrack*   GetPamTrackAlong(TrkTrack*);      PamTrack*   GetPamTrackAlong(TrkTrack*);
518  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
519  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
520            
521  //  =================================  //  =================================
522  //  methods to retrieve pamela tracks  //  methods to retrieve pamela tracks
523  //  =================================  //  =================================
# Line 506  class PamLevel2 : public TObject { Line 529  class PamLevel2 : public TObject {
529      TClonesArray* GetTracks();                  //return sorted_tracks      TClonesArray* GetTracks();                  //return sorted_tracks
530      PamTrack*     GetTrack(int);      PamTrack*     GetTrack(int);
531      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImage(int);
532      void          SetSortingMethod(TString how){ howtosort = how; return;};      void          SetSortingMethod(TString how);
533      TString       GetSortingMethod(){return howtosort;};      TString       GetSortingMethod(){return howtosort;};
534      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
535      Float_t       GetNucleiSortingThreshold(){ return sortthr;};      Float_t       GetNucleiSortingThreshold(){ return sortthr;};
536            
537  //  =================================  //  =================================
538  //  other methods  //  other methods
539  //  =================================  //  =================================
# Line 546  class PamLevel2 : public TObject { Line 569  class PamLevel2 : public TObject {
569      Bool_t IsGood(Bool_t strict=true);      Bool_t IsGood(Bool_t strict=true);
570    
571      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)
572      void  SetSELLI(Int_t selli){SELLI=selli;};      void  SetSELLI(Int_t selli){SELLI=selli;};
573    
574      Bool_t IsGP(){return ISGP;};      Bool_t IsGP(){return ISGP;};
575    
# Line 554  class PamLevel2 : public TObject { Line 577  class PamLevel2 : public TObject {
577    
578      void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };      void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
579    
580      ClassDef(PamLevel2,1);  
581        ClassDef(PamLevel2,2);
582  };  };
583    
584  //===============================================================================================  //===============================================================================================
Legend:
Removed from v.1.40  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.23