/[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.41 by pam-fi, Wed Apr 1 09:25:36 2009 UTC revision 1.54 by pam-fi, Tue Mar 24 09:15:00 2015 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 39  Line 40 
40  #include <TChain.h>  #include <TChain.h>
41  #include <TChainElement.h>  #include <TChainElement.h>
42  #include <TFile.h>  #include <TFile.h>
43    #include <TTimeStamp.h>
44  #include <TList.h>  #include <TList.h>
45  #include <TKey.h>  #include <TKey.h>
46  #include <TSystemFile.h>  #include <TSystemFile.h>
# Line 225  class PamTrack : public TObject { Line 227  class PamTrack : public TObject {
227    
228   private:   private:
229    
230        ExtTrack*   trk_ext_track;
231      TrkTrack*   trk_track;      TrkTrack*   trk_track;
232      CaloTrkVar* calo_track;      CaloTrkVar* calo_track;
233      ToFTrkVar*  tof_track;      ToFTrkVar*  tof_track;
# Line 237  class PamTrack : public TObject { Line 240  class PamTrack : public TObject {
240    
241   public:   public:
242    
243    
244      PamTrack();      PamTrack();
245      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );      PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
246        PamTrack(ExtTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
247      PamTrack(const PamTrack&);      PamTrack(const PamTrack&);
248      ~PamTrack(){ Delete(); };      ~PamTrack(){ Delete(); };
249    
250      void Clear();      void Copy(PamTrack& ) const;
251      void Clear(Option_t *option){Clear();};  
252        void Clear(Option_t *option);
253        void Clear(){Clear("");};
254      void Delete();      void Delete();
255    
256      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 260  class PamTrack : public TObject {
260      void SetPScore(Int_t p){pscore=p;};      void SetPScore(Int_t p){pscore=p;};
261      void SetIScore(Int_t i){iscore=i;};      void SetIScore(Int_t i){iscore=i;};
262    
263        ExtTrack*   GetExtTrack(){return trk_ext_track;};
264      TrkTrack*   GetTrkTrack(){return trk_track;};      TrkTrack*   GetTrkTrack(){return trk_track;};
265      CaloTrkVar* GetCaloTrack(){return calo_track;};      CaloTrkVar* GetCaloTrack(){return calo_track;};
266      ToFTrkVar*  GetToFTrack(){return tof_track;};      ToFTrkVar*  GetToFTrack(){return tof_track;};
267      OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};      OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};
268    
269      ClassDef(PamTrack,2);      ClassDef(PamTrack,3);
270  };  };
271    
272  /**  /**
# Line 281  class PamLevel2 : public TObject { Line 289  class PamLevel2 : public TObject {
289      Long64_t  runfirstentry; //first entry of run, relative to the TChain      Long64_t  runfirstentry; //first entry of run, relative to the TChain
290      Long64_t  runlastentry;  //last  entry of run, relative to the TChain      Long64_t  runlastentry;  //last  entry of run, relative to the TChain
291      Long64_t  irunentry;     //position of current entry, relative to the run      Long64_t  irunentry;     //position of current entry, relative to the run
292        Long64_t  totrunentry;
293        Long64_t  totrunentrymax;
294        Long64_t  totrunentrymin;
295        ULong64_t    prevabstime; // abs time of previous event
296        UInt_t    prevpktnum; // packet number of previous event
297        ULong64_t    abstime; // current abs time
298        UInt_t    pktnum; // current packet number
299        Bool_t isFragment;
300        Bool_t isSync; ///< the file is fully sinchronized between L2 and L0/DB data
301            //    Bool_t hasL0EE; ///< backward compatibility flag for preselected event list
302        Long64_t il0entry; ///< entry of level0 file in selection list
303      TChain*   run_tree;      TChain*   run_tree;
304      TTree*    run_tree_clone;      TTree*    run_tree_clone;
305        TChain*   proc_tree;
306        TTree*    proc_tree_clone;
307      TChain*   sel_tree;      TChain*   sel_tree;
308      TTree*    sel_tree_clone;      TTree*    sel_tree_clone;
309      int *irunoffset;      int *irunoffset;
310      GL_TIMESYNC *gltsync; // Emiliano      GL_TIMESYNC *gltsync; // Emiliano
311        Bool_t fUpdateRunInfo; // Emiliano
312        Bool_t fUseDBinRunInfo; // Emiliano
313    
314        Bool_t fDiscarded; ///< true if any file from a list is missing
315        TString customString; ///< String passed to ProcessingInfo
316    
317      // -------------------------------------------      // -------------------------------------------
318      // variables to store total live and dead time      // variables to store total live and dead time
# Line 299  class PamLevel2 : public TObject { Line 325  class PamLevel2 : public TObject {
325      TFile*      l0_file;      TFile*      l0_file;
326      TTree*      l0_tree;      TTree*      l0_tree;
327      Int_t       iroot;      Int_t       iroot;
328      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
329        Long64_t       yprevshift; // shift needed to sync level2 data with level0 data (nested+DV_skipped packets) - run wise EM
330      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
331      // --------------------------------------------      // --------------------------------------------
332      // variables to store Level0 info (if required)      // variables to store Level0 info (if required)
# Line 330  class PamLevel2 : public TObject { Line 357  class PamLevel2 : public TObject {
357      TrkHough*    trkh_obj;      TrkHough*    trkh_obj;
358      CaloLevel1*  calo1_obj;      CaloLevel1*  calo1_obj;
359      CaloLevel2*  calo2_obj;      CaloLevel2*  calo2_obj;
360      ToFLevel2*   tof_obj;      ToFLevel2*   tof2_obj;
361      TrigLevel2*  trig_obj;      TrigLevel2*  trig_obj;
362      S4Level2*    s4_obj;      S4Level2*    s4_obj;
363      NDLevel2*    nd_obj;      NDLevel2*    nd_obj;
364      AcLevel2*    ac_obj;      AcLevel2*    ac_obj;
365      OrbitalInfo* orb_obj;      OrbitalInfo* orb2_obj;
366      GL_RUN*      run_obj;      GL_RUN*      run_obj;
367      SoftInfo*    soft_obj; // Emiliano      SoftInfo*    soft_obj; // Emiliano
368        ProcInfo*    proc_obj; //EM
369    
370    
371    
372        Int_t extAlgFlag;
373    
374        TClonesArray* trk_ext_obj;
375        TClonesArray* trk_ext_nuc_obj;
376        TClonesArray* trk_nuc_obj;
377    
378    
379        TClonesArray* calo_ext_obj;
380        TClonesArray* calo_ext_nuc_obj;
381        TClonesArray* calo_nuc_obj;
382    
383        TClonesArray* tof_ext_obj;
384        TClonesArray* tof_ext_nuc_obj;
385        TClonesArray* tof_nuc_obj;
386    
387        TClonesArray* orb_ext_obj;
388        TClonesArray* orb_ext_nuc_obj;
389        TClonesArray* orb_nuc_obj;
390    
391    
392        TrkLevel2*   trk2_nuc_obj;
393        CaloLevel2*  calo2_nuc_obj;
394        ToFLevel2*   tof2_nuc_obj;
395        OrbitalInfo* orb2_nuc_obj;
396    
397    
398    
399      // --------------------------------------------      // --------------------------------------------
400      // variables to handle GP      // variables to handle GP
# Line 349  class PamLevel2 : public TObject { Line 406  class PamLevel2 : public TObject {
406      // -----------------------------------      // -----------------------------------
407      TClonesArray *tsorted;      TClonesArray *tsorted;
408      TClonesArray *timage;      TClonesArray *timage;
409        TClonesArray *text; //array of pamtrack objects built from extended algorythm
410        TClonesArray *tsorted_nuc;
411        TClonesArray *timage_nuc;
412        TClonesArray *text_nuc; //array of pamtrack objects built from extended algorythm
413    
414      PamTrack    *track; //store the last required      PamTrack    *track; //store the last required
415      TString howtosort;      TString howtosort;
416      Float_t sortthr;      Float_t sortthr;
417    
418      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
419      Long64_t lastsorted; // the last event whose tracks has been sorted      Long64_t lastsorted; // the last event whose tracks has been sorted
420    
421    
422        Bool_t issorted_new;     // true if the current event has already been sorted with the current sort criteria
423        Long64_t lastsorted_new; // the last event whose tracks has been sorted
424    
425      // --------------------------------------      // --------------------------------------
426      // flags to load level0-1-2 trees/braches      // flags to load level0-1-2 trees/braches
427      // --------------------------------------      // --------------------------------------
# Line 378  class PamLevel2 : public TObject { Line 445  class PamLevel2 : public TObject {
445    
446      Bool_t GP; //GP simulation tree      Bool_t GP; //GP simulation tree
447    
448    
449        Bool_t EXT;//get tracks from extended tracking algorythm
450        Bool_t NUC;//get tracks from nuclei algorythm
451    
452        Bool_t PROC;//EM
453    
454        const char*  trkAlg;
455    
456      // -----------      // -----------
457      // other flags      // other flags
458      // -----------      // -----------
# Line 422  class PamLevel2 : public TObject { Line 497  class PamLevel2 : public TObject {
497      TTree*  GetRunTree(TFile*);      TTree*  GetRunTree(TFile*);
498      TChain* GetRunTree(TList*);      TChain* GetRunTree(TList*);
499      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
500      Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };      Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug
501      Bool_t  UpdateRunInfo(TTree *,Long64_t);      Bool_t  UpdateRunInfo(TTree *,Long64_t);
502      Bool_t  UpdateRunInfo(TChain *,Long64_t);      Bool_t  UpdateRunInfo(TChain *,Long64_t);
503    
504        void    SkipRunInfoUpdate(); // Emiliano
505        void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano
506        void    UseDBinRunInfo(){fUseDBinRunInfo = true;}; // [default] Emiliano
507        void    NoDBinRunInfo(){fUseDBinRunInfo = false;}; // EM
508        void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano
509    
510      void    CreateCloneTrees0(TChain *fChain, TFile* file);      void    CreateCloneTrees0(TChain *fChain, TFile* file);
511      void    CreateCloneTrees(TFile* file);      void    CreateCloneTrees(TFile* file);
512      void    CreateCloneTrees(){CreateCloneTrees(NULL);};      void    CreateCloneTrees(){CreateCloneTrees(NULL);};
# Line 441  class PamLevel2 : public TObject { Line 522  class PamLevel2 : public TObject {
522      Long64_t GetReadEntry();      Long64_t GetReadEntry();
523      Int_t   GetYodaEntry();      Int_t   GetYodaEntry();
524    
525      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
526    
527      Long64_t GetEntries(){return pam_tree->GetEntries();};      Long64_t GetEntries(){return pam_tree->GetEntries();};
528      TChain* GetPamTree(){return pam_tree;};      TChain* GetPamTree(){return pam_tree;};
529      TChain* GetRunTree(){return run_tree;};      TChain* GetRunTree(){return run_tree;};
530        TChain* GetProcessingInfoTree(){return proc_tree;};
531      TTree*  GetYodaTree();      TTree*  GetYodaTree();
532    
533      TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things      TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things
# Line 476  class PamLevel2 : public TObject { Line 558  class PamLevel2 : public TObject {
558      TrkHough     *GetTrkHough()   { return trkh_obj; };      TrkHough     *GetTrkHough()   { return trkh_obj; };
559      CaloLevel1   *GetCaloLevel1() { return calo1_obj;};      CaloLevel1   *GetCaloLevel1() { return calo1_obj;};
560      CaloLevel2   *GetCaloLevel2() { return calo2_obj;};      CaloLevel2   *GetCaloLevel2() { return calo2_obj;};
561      ToFLevel2    *GetToFLevel2()  { return tof_obj;  };      ToFLevel2    *GetToFLevel2()  { return tof2_obj;  };
562      TrigLevel2   *GetTrigLevel2() { return trig_obj; };      TrigLevel2   *GetTrigLevel2() { return trig_obj; };
563      S4Level2     *GetS4Level2()   { return s4_obj;   };      S4Level2     *GetS4Level2()   { return s4_obj;   };
564      NDLevel2     *GetNDLevel2()   { return nd_obj;   };      NDLevel2     *GetNDLevel2()   { return nd_obj;   };
565      AcLevel2     *GetAcLevel2()   { return ac_obj;   };      AcLevel2     *GetAcLevel2()   { return ac_obj;   };
566      OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };      OrbitalInfo  *GetOrbitalInfo(){ return orb2_obj;  };
567      GL_RUN       *GetRunInfo()    { return run_obj;  };      GL_RUN       *GetRunInfo()    { return run_obj;  };
568      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
569        ProcInfo     *GetProcInfo()   { return proc_obj; }; // Emiliano
570      GPamela      *GetGPamela()    { return gp_obj;   };      GPamela      *GetGPamela()    { return gp_obj;   };
571    
572  //  ==============================================  //  ==============================================
# Line 494  class PamLevel2 : public TObject { Line 577  class PamLevel2 : public TObject {
577      ToFTrkVar*  GetToFStoredTrack(int);      ToFTrkVar*  GetToFStoredTrack(int);
578      OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);      OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);
579      //      //
580      PamTrack*   GetPamTrackAlong(TrkTrack*);  //    PamTrack*   GetPamTrackAlong(TrkTrack*);
581  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
582  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
583    
584  //  =================================  //  =================================
585  //  methods to retrieve pamela tracks  //  methods to retrieve pamela tracks
586  //  =================================  //  =================================
587      PamTrack*     GetStoredTrack(Int_t);  //    PamTrack*     GetStoredTrack(Int_t);//obsolete
588      //      //
589      void          SortTracks(TString);          //assign value to variable sorted_tracks      void          SortTracks(TString);          //assign value to variable sorted_tracks
590      void          SortTracks();          //assign value to variable sorted_tracks      void          SortTracks();          //assign value to variable sorted_tracks
591  //    TRefArray*    GetTracks();                  //return sorted_tracks  //    TRefArray*    GetTracks();                  //return sorted_tracks
592      TClonesArray* GetTracks();                  //return sorted_tracks      TClonesArray* GetTracks();                  //return sorted_tracks
593      PamTrack*     GetTrack(int);      PamTrack*     GetTrackOld(int);
594      PamTrack*     GetTrackImage(int);      PamTrack*     GetTrackImageOld(int);
595      void          SetSortingMethod(TString how);      PamTrack*     GetTrack(int i){ return GetTrackOld(i); };
596        PamTrack*     GetTrackImage(int i){ return GetTrackImageOld(i); };
597       void          SetSortingMethod(TString how);
598      TString       GetSortingMethod(){return howtosort;};      TString       GetSortingMethod(){return howtosort;};
599      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
600      Float_t       GetNucleiSortingThreshold(){ return sortthr;};      Float_t       GetNucleiSortingThreshold(){ return sortthr;};
601    
602    
603    
604    
605    //  =================================
606    //  methods to handle new tracking algorythms
607    //  =================================
608    /*     Int_t GetTrackingAlgorythmFlag(){return extAlgFlag;}; */
609    /*     void  SetTrackingAlgorythm(const char * alg); */
610    /*     const char* GetTrackingAlgorythm(); */
611    
612        Int_t  GetNTracks(const char*);
613        PamTrack*     GetTrack(int,const char*);
614        PamTrack*     GetTrackImage(int, const char*);
615        Bool_t TrackHasImage(int i, const char* alg){return GetTrackImage(i,alg)!=NULL;};
616    
617    /*     Int_t  GetNTracks(){return GetNTracks(trkAlg); }; */
618    /*     PamTrack*     GetTrack(int i) { return GetTrack(i,trkAlg); }; */
619    /*     PamTrack*     GetTrackImage(int i){ return GetTrackImage(i,trkAlg); }; */
620    
621        void          SortTracksNew();          //assign value to variable sorted_tracks
622    
623    
624  //  =================================  //  =================================
625  //  other methods  //  other methods
626  //  =================================  //  =================================
# Line 546  class PamLevel2 : public TObject { Line 653  class PamLevel2 : public TObject {
653    
654      Bool_t IsRUN(){return RUN;};      Bool_t IsRUN(){return RUN;};
655    
656        Bool_t IsEXT(){return EXT;};
657        Bool_t IsNUC(){return NUC;};
658    
659      Bool_t IsGood(Bool_t strict=true);      Bool_t IsGood(Bool_t strict=true);
660    
661      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)
662      void  SetSELLI(Int_t selli){SELLI=selli;};      void  SetSELLI(Int_t selli){SELLI=selli;};
663    
664        Bool_t IsPROC(){return PROC;};//EM
665    
666        Bool_t HasDiscardedFiles(){return fDiscarded;}; ///< true if something went wrong and not all files were processed
667    
668      Bool_t IsGP(){return ISGP;};      Bool_t IsGP(){return ISGP;};
669    
670      void   SetGP(Bool_t tt){ISGP=tt;};      void   SetGP(Bool_t tt){ISGP=tt;};
671    
672      void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };      void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
673    
674        void SetCustomString(TString cus){customString=cus;};
675      ClassDef(PamLevel2,1);      TString GetCustomString(){return customString;};
676    
677        ClassDef(PamLevel2,4);
678  };  };
679    
680  //===============================================================================================  //===============================================================================================
Legend:
Removed from v.1.41  
changed lines
  Added in v.1.54
  ViewVC Help
Powered by ViewVC 1.1.23