/[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.44 by mocchiut, Mon Dec 12 16:11:41 2011 UTC revision 1.51 by mocchiut, Fri Oct 17 14:03:45 2014 UTC
# Line 32  Line 32 
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>
35    #include <ExtTrack.h>
36    
37    
38  #include <TTree.h>  #include <TTree.h>
# Line 225  class PamTrack : public TObject { Line 226  class PamTrack : public TObject {
226    
227   private:   private:
228    
229        ExtTrack*   trk_ext_track;
230      TrkTrack*   trk_track;      TrkTrack*   trk_track;
231      CaloTrkVar* calo_track;      CaloTrkVar* calo_track;
232      ToFTrkVar*  tof_track;      ToFTrkVar*  tof_track;
# Line 239  class PamTrack : public TObject { Line 241  class PamTrack : public TObject {
241    
242      PamTrack();      PamTrack();
243      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
244        PamTrack(ExtTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
245      PamTrack(const PamTrack&);      PamTrack(const PamTrack&);
246      ~PamTrack(){ Delete(); };      ~PamTrack(){ Delete(); };
247    
248      void Clear();      void Clear(Option_t *option);
249      void Clear(Option_t *option){Clear();};      void Clear(){Clear("");};
250      void Delete();      void Delete();
251    
252      Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};      Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
# Line 253  class PamTrack : public TObject { Line 256  class PamTrack : public TObject {
256      void SetPScore(Int_t p){pscore=p;};      void SetPScore(Int_t p){pscore=p;};
257      void SetIScore(Int_t i){iscore=i;};      void SetIScore(Int_t i){iscore=i;};
258    
259        ExtTrack*   GetExtTrack(){return trk_ext_track;};
260      TrkTrack*   GetTrkTrack(){return trk_track;};      TrkTrack*   GetTrkTrack(){return trk_track;};
261      CaloTrkVar* GetCaloTrack(){return calo_track;};      CaloTrkVar* GetCaloTrack(){return calo_track;};
262      ToFTrkVar*  GetToFTrack(){return tof_track;};      ToFTrkVar*  GetToFTrack(){return tof_track;};
263      OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};      OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};
264    
265      ClassDef(PamTrack,2);      ClassDef(PamTrack,3);
266  };  };
267    
268  /**  /**
# Line 289  class PamLevel2 : public TObject { Line 293  class PamLevel2 : public TObject {
293      ULong64_t    abstime; // current abs time      ULong64_t    abstime; // current abs time
294      UInt_t    pktnum; // current packet number      UInt_t    pktnum; // current packet number
295      Bool_t isFragment;      Bool_t isFragment;
296        Bool_t isSync; ///< the file is fully sinchronized between L2 and L0/DB data
297            //    Bool_t hasL0EE; ///< backward compatibility flag for preselected event list
298        Long64_t il0entry; ///< entry of level0 file in selection list
299      TChain*   run_tree;      TChain*   run_tree;
300      TTree*    run_tree_clone;      TTree*    run_tree_clone;
301      TChain*   sel_tree;      TChain*   sel_tree;
# Line 296  class PamLevel2 : public TObject { Line 303  class PamLevel2 : public TObject {
303      int *irunoffset;      int *irunoffset;
304      GL_TIMESYNC *gltsync; // Emiliano      GL_TIMESYNC *gltsync; // Emiliano
305      Bool_t fUpdateRunInfo; // Emiliano      Bool_t fUpdateRunInfo; // Emiliano
306        Bool_t fUseDBinRunInfo; // Emiliano
307    
308        Bool_t fDiscarded; ///< true if any file from a list is missing
309    
310      // -------------------------------------------      // -------------------------------------------
311      // variables to store total live and dead time      // variables to store total live and dead time
# Line 308  class PamLevel2 : public TObject { Line 318  class PamLevel2 : public TObject {
318      TFile*      l0_file;      TFile*      l0_file;
319      TTree*      l0_tree;      TTree*      l0_tree;
320      Int_t       iroot;      Int_t       iroot;
321      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
322        Long64_t       yprevshift; // shift needed to sync level2 data with level0 data (nested+DV_skipped packets) - run wise EM
323      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
324      // --------------------------------------------      // --------------------------------------------
325      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
# Line 339  class PamLevel2 : public TObject { Line 350  class PamLevel2 : public TObject {
350      TrkHough*    trkh_obj;      TrkHough*    trkh_obj;
351      CaloLevel1*  calo1_obj;      CaloLevel1*  calo1_obj;
352      CaloLevel2*  calo2_obj;      CaloLevel2*  calo2_obj;
353      ToFLevel2*   tof_obj;      ToFLevel2*   tof2_obj;
354      TrigLevel2*  trig_obj;      TrigLevel2*  trig_obj;
355      S4Level2*    s4_obj;      S4Level2*    s4_obj;
356      NDLevel2*    nd_obj;      NDLevel2*    nd_obj;
357      AcLevel2*    ac_obj;      AcLevel2*    ac_obj;
358      OrbitalInfo* orb_obj;      OrbitalInfo* orb2_obj;
359      GL_RUN*      run_obj;      GL_RUN*      run_obj;
360      SoftInfo*    soft_obj; // Emiliano      SoftInfo*    soft_obj; // Emiliano
361    
362    
363    
364        Int_t extAlgFlag;
365    
366        TClonesArray* trk_ext_obj;
367        TClonesArray* trk_ext_nuc_obj;
368        TClonesArray* trk_nuc_obj;
369    
370    
371        TClonesArray* calo_ext_obj;
372        TClonesArray* calo_ext_nuc_obj;
373        TClonesArray* calo_nuc_obj;
374    
375        TClonesArray* tof_ext_obj;
376        TClonesArray* tof_ext_nuc_obj;
377        TClonesArray* tof_nuc_obj;
378    
379        TClonesArray* orb_ext_obj;
380        TClonesArray* orb_ext_nuc_obj;
381        TClonesArray* orb_nuc_obj;
382    
383    
384        TrkLevel2*   trk2_nuc_obj;
385        CaloLevel2*  calo2_nuc_obj;
386        ToFLevel2*   tof2_nuc_obj;
387        OrbitalInfo* orb2_nuc_obj;
388    
389    
390    
391      // --------------------------------------------      // --------------------------------------------
392      // variables to handle GP      // variables to handle GP
393      // --------------------------------------------      // --------------------------------------------
# Line 358  class PamLevel2 : public TObject { Line 398  class PamLevel2 : public TObject {
398      // -----------------------------------      // -----------------------------------
399      TClonesArray *tsorted;      TClonesArray *tsorted;
400      TClonesArray *timage;      TClonesArray *timage;
401        TClonesArray *text; //array of pamtrack objects built from extended algorythm
402        TClonesArray *tsorted_nuc;
403        TClonesArray *timage_nuc;
404        TClonesArray *text_nuc; //array of pamtrack objects built from extended algorythm
405    
406      PamTrack    *track; //store the last required      PamTrack    *track; //store the last required
407      TString howtosort;      TString howtosort;
408      Float_t sortthr;      Float_t sortthr;
409    
410      Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria      Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria
411      Long64_t lastsorted; // the last event whose tracks has been sorted      Long64_t lastsorted; // the last event whose tracks has been sorted
412    
413    
414        Bool_t issorted_new;     // true if the current event has already been sorted with the current sort criteria
415        Long64_t lastsorted_new; // the last event whose tracks has been sorted
416    
417      // --------------------------------------      // --------------------------------------
418      // flags to load level0-1-2 trees/braches      // flags to load level0-1-2 trees/braches
419      // --------------------------------------      // --------------------------------------
# Line 387  class PamLevel2 : public TObject { Line 437  class PamLevel2 : public TObject {
437    
438      Bool_t GP; //GP simulation tree      Bool_t GP; //GP simulation tree
439    
440    
441        Bool_t EXT;//get tracks from extended tracking algorythm
442        Bool_t NUC;//get tracks from nuclei algorythm
443    
444        const char*  trkAlg;
445    
446      // -----------      // -----------
447      // other flags      // other flags
448      // -----------      // -----------
# Line 432  class PamLevel2 : public TObject { Line 488  class PamLevel2 : public TObject {
488      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
489      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
490      Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug      Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug
     Bool_t  UpdateRunInfo_bckw_compat(Long64_t);///< backward compatibility method  
491      Bool_t  UpdateRunInfo(TTree *,Long64_t);      Bool_t  UpdateRunInfo(TTree *,Long64_t);
492      Bool_t  UpdateRunInfo(TChain *,Long64_t);      Bool_t  UpdateRunInfo(TChain *,Long64_t);
493    
494      void    SkipRunInfoUpdate(); // Emiliano      void    SkipRunInfoUpdate(); // Emiliano
495      void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano      void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano
496        void    UseDBinRunInfo(){fUseDBinRunInfo = true;}; // [default] Emiliano
497        void    NoDBinRunInfo(){fUseDBinRunInfo = false;}; // EM
498      void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano      void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano
499    
500      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
# Line 455  class PamLevel2 : public TObject { Line 512  class PamLevel2 : public TObject {
512      Long64_t GetReadEntry();      Long64_t GetReadEntry();
513      Int_t   GetYodaEntry();      Int_t   GetYodaEntry();
514    
515      void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable      void SetMaxShift(Int_t sh); // EMILIANO method to change the maxshift variable
516    
517      Long64_t GetEntries(){return pam_tree->GetEntries();};      Long64_t GetEntries(){return pam_tree->GetEntries();};
518      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
# Line 490  class PamLevel2 : public TObject { Line 547  class PamLevel2 : public TObject {
547      TrkHough     *GetTrkHough()   { return trkh_obj; };      TrkHough     *GetTrkHough()   { return trkh_obj; };
548      CaloLevel1   *GetCaloLevel1() { return calo1_obj;};      CaloLevel1   *GetCaloLevel1() { return calo1_obj;};
549      CaloLevel2   *GetCaloLevel2() { return calo2_obj;};      CaloLevel2   *GetCaloLevel2() { return calo2_obj;};
550      ToFLevel2    *GetToFLevel2()  { return tof_obj;  };      ToFLevel2    *GetToFLevel2()  { return tof2_obj;  };
551      TrigLevel2   *GetTrigLevel2() { return trig_obj; };      TrigLevel2   *GetTrigLevel2() { return trig_obj; };
552      S4Level2     *GetS4Level2()   { return s4_obj;   };      S4Level2     *GetS4Level2()   { return s4_obj;   };
553      NDLevel2     *GetNDLevel2()   { return nd_obj;   };      NDLevel2     *GetNDLevel2()   { return nd_obj;   };
554      AcLevel2     *GetAcLevel2()   { return ac_obj;   };      AcLevel2     *GetAcLevel2()   { return ac_obj;   };
555      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };      OrbitalInfo  *GetOrbitalInfo(){ return orb2_obj;  };
556      GL_RUN       *GetRunInfo()    { return run_obj;  };      GL_RUN       *GetRunInfo()    { return run_obj;  };
557      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
558      GPamela      *GetGPamela()    { return gp_obj;   };      GPamela      *GetGPamela()    { return gp_obj;   };
# Line 521  class PamLevel2 : public TObject { Line 578  class PamLevel2 : public TObject {
578      void          SortTracks();          //assign value to variable sorted_tracks      void          SortTracks();          //assign value to variable sorted_tracks
579  //    TRefArray*    GetTracks();                  //return sorted_tracks  //    TRefArray*    GetTracks();                  //return sorted_tracks
580      TClonesArray* GetTracks();                  //return sorted_tracks      TClonesArray* GetTracks();                  //return sorted_tracks
581      PamTrack*     GetTrack(int);      PamTrack*     GetTrackOld(int);
582      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImageOld(int);
583      void          SetSortingMethod(TString how);      PamTrack*     GetTrack(int i){ return GetTrackOld(i); };
584        PamTrack*     GetTrackImage(int i){ return GetTrackImageOld(i); };
585       void          SetSortingMethod(TString how);
586      TString       GetSortingMethod(){return howtosort;};      TString       GetSortingMethod(){return howtosort;};
587      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
588      Float_t       GetNucleiSortingThreshold(){ return sortthr;};      Float_t       GetNucleiSortingThreshold(){ return sortthr;};
589    
590    
591    
592    
593    //  =================================
594    //  methods to handle new tracking algorythms
595    //  =================================
596    /*     Int_t GetTrackingAlgorythmFlag(){return extAlgFlag;}; */
597    /*     void  SetTrackingAlgorythm(const char * alg); */
598    /*     const char* GetTrackingAlgorythm(); */
599    
600        Int_t  GetNTracks(const char*);
601        PamTrack*     GetTrack(int,const char*);
602        PamTrack*     GetTrackImage(int, const char*);
603    
604    
605    /*     Int_t  GetNTracks(){return GetNTracks(trkAlg); }; */
606    /*     PamTrack*     GetTrack(int i) { return GetTrack(i,trkAlg); }; */
607    /*     PamTrack*     GetTrackImage(int i){ return GetTrackImage(i,trkAlg); }; */
608    
609        void          SortTracksNew();          //assign value to variable sorted_tracks
610    
611    
612  //  =================================  //  =================================
613  //  other methods  //  other methods
614  //  =================================  //  =================================
# Line 560  class PamLevel2 : public TObject { Line 641  class PamLevel2 : public TObject {
641    
642      Bool_t IsRUN(){return RUN;};      Bool_t IsRUN(){return RUN;};
643    
644        Bool_t IsEXT(){return EXT;};
645        Bool_t IsNUC(){return NUC;};
646    
647      Bool_t IsGood(Bool_t strict=true);      Bool_t IsGood(Bool_t strict=true);
648    
649      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)
650      void  SetSELLI(Int_t selli){SELLI=selli;};      void  SetSELLI(Int_t selli){SELLI=selli;};
651    
652        Bool_t HasDiscardedFiles(){return fDiscarded;}; ///< true if something went wrong and not all files were processed
653    
654      Bool_t IsGP(){return ISGP;};      Bool_t IsGP(){return ISGP;};
655    
656      void   SetGP(Bool_t tt){ISGP=tt;};      void   SetGP(Bool_t tt){ISGP=tt;};
657    
658      void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };      void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
659    
660    
661      ClassDef(PamLevel2,2);      ClassDef(PamLevel2,3);
662  };  };
663    
664  //===============================================================================================  //===============================================================================================
Legend:
Removed from v.1.44  
changed lines
  Added in v.1.51
  ViewVC Help
Powered by ViewVC 1.1.23