/[PAMELA software]/PamelaLevel2/inc/PamLevel2.h

Diff of /PamelaLevel2/inc/PamLevel2.h

Parent Directory | Revision Log | View Patch Patch

-revision 1.32 by pam-fi,
Tue Nov 20 18:16:29 2007 UTC
+revision 1.54 by pam-fi,
Tue Mar 24 09:15:00 2015 UTC
 Line 24
  #include <NDLevel2.h>
  #include <AcLevel2.h>
  #include <OrbitalInfo.h>
+ #include <InclinationInfo.h> // EMILIANO, now this library is included in libDarthVader.so
+                          // and the header installed in $PAM_INC. If you get the
+                          // error "Cannot find InclinationInfo.h" update your DV installation
+ #include <OrientationInfo.h> // EMILIANO, now this library is included in libDarthVader.so
+                          // and the header installed in $PAM_INC. If you get the
+                          // error "Cannot find InclinationInfo.h" update your DV installation
  #include <RunInfo.h>
  #include <GLTables.h>
+ #include <ExtTrack.h>
  #include <TTree.h>
-Line 33
+Line 40
  #include <TChain.h>
  #include <TChainElement.h>
  #include <TFile.h>
+ #include <TTimeStamp.h>
  #include <TList.h>
  #include <TKey.h>
  #include <TSystemFile.h>
-Line 55 
 using namespace std;
+Line 63 
 using namespace std;
  /////////////////////////////////////////////////////////////////////////////
  /**
   * \brief Class to describe GP output (to be moved to PamelaDigitizer)
   *
   */
  class GPamela : public TObject {
-Line 211 
 class GPamela : public TObject {
+Line 219 
 class GPamela : public TObject {
  /**
   * \brief Class to describe particle tracks in the apparatus, including tracker and track-related variables from other detectors (calorimeter and TOF).
   *
   * This class inherites from TrkTrack, CaloTrkVar ecc...
   *
   */
  //class PamTrack : public TrkTrack, public CaloTrkVar, public ToFTrkVar {
  class PamTrack : public TObject {
   private:
+     ExtTrack*   trk_ext_track;
      TrkTrack*   trk_track;
      CaloTrkVar* calo_track;
      ToFTrkVar*  tof_track;
+     OrbitalInfoTrkVar*  orb_track;
      Int_t pscore;
      Int_t iscore;
      bool candeleteobj;
-Line 230 
 class PamTrack : public TObject {
+Line 240 
 class PamTrack : public TObject {
   public:
      PamTrack();
-     PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar* );
+     PamTrack(TrkTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
+     PamTrack(ExtTrack* , CaloTrkVar* , ToFTrkVar*, OrbitalInfoTrkVar* );
      PamTrack(const PamTrack&);
      ~PamTrack(){ Delete(); };
-     void Clear();
+     void Copy(PamTrack& ) const;
-     void Clear(Option_t *option){Clear();};
+     void Clear(Option_t *option);
+     void Clear(){Clear("");};
      void Delete();
      Bool_t IsSolved(){if ((pscore-iscore)>0) return true; return false;};
-Line 246 
 class PamTrack : public TObject {
+Line 260 
 class PamTrack : public TObject {
      void SetPScore(Int_t p){pscore=p;};
      void SetIScore(Int_t i){iscore=i;};
+     ExtTrack*   GetExtTrack(){return trk_ext_track;};
      TrkTrack*   GetTrkTrack(){return trk_track;};
      CaloTrkVar* GetCaloTrack(){return calo_track;};
      ToFTrkVar*  GetToFTrack(){return tof_track;};
+     OrbitalInfoTrkVar*  GetOrbitalInfoTrack(){return orb_track;};
-     ClassDef(PamTrack,1);
+     ClassDef(PamTrack,3);
  };
  /**
   * \brief Class to describe PAMELA LEVEL2 data.
   *
   * This class inherites from TrkLevel2, CaloLevel2 ecc... and it overrides the TrkLevel2:GetTrack(int it) and TrkLevel2:GetTrackImage(int it) methods.
   * This allows to solve the tracker ambiguity (that is, to choose between track images) by using also information from other detectors.
   */
  //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 TObject {
-Line 273 
 class PamLevel2 : public TObject {
+Line 289 
 class PamLevel2 : public TObject {
      Long64_t  runfirstentry; //first entry of run, relative to the TChain
      Long64_t  runlastentry;  //last  entry of run, relative to the TChain
      Long64_t  irunentry;     //position of current entry, relative to the run
+     Long64_t  totrunentry;
+     Long64_t  totrunentrymax;
+     Long64_t  totrunentrymin;
+     ULong64_t    prevabstime; // abs time of previous event
+     UInt_t    prevpktnum; // packet number of previous event
+     ULong64_t    abstime; // current abs time
+     UInt_t    pktnum; // current packet number
+     Bool_t isFragment;
+     Bool_t isSync; ///< the file is fully sinchronized between L2 and L0/DB data
+         //    Bool_t hasL0EE; ///< backward compatibility flag for preselected event list
+     Long64_t il0entry; ///< entry of level0 file in selection list
      TChain*   run_tree;
      TTree*    run_tree_clone;
+     TChain*   proc_tree;
+     TTree*    proc_tree_clone;
      TChain*   sel_tree;
      TTree*    sel_tree_clone;
      int *irunoffset;
+     GL_TIMESYNC *gltsync; // Emiliano
+     Bool_t fUpdateRunInfo; // Emiliano
+     Bool_t fUseDBinRunInfo; // Emiliano
+     Bool_t fDiscarded; ///< true if any file from a list is missing
+     TString customString; ///< String passed to ProcessingInfo
      // -------------------------------------------
      // variables to store total live and dead time
      // -------------------------------------------
      ULong64_t  totdltime[3];// 0=live 1=dead 2=n.events
      // ---------------------------------
      // variables to retrieve level0 info
      // ---------------------------------
      TFile*      l0_file;
      TTree*      l0_tree;
      Int_t       iroot;
-     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
-     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
+     Int_t       maxshift; // EMILIANO maximum discrepancy between L2 and L0 data (the bigger the better concering synch but the processing speed drops) default = 10
      // --------------------------------------------
      // variables to store Level0 info (if required)
      // --------------------------------------------
-Line 314 
 class PamLevel2 : public TObject {
+Line 350 
 class PamLevel2 : public TObject {
      TTree*  pam_tree_clone[NCLONES];
      // --------------------------------------------
      // variables to store Level1-2 info
      // --------------------------------------------
      TrkLevel1*   trk1_obj;
      TrkLevel2*   trk2_obj;
      TrkHough*    trkh_obj;
      CaloLevel1*  calo1_obj;
      CaloLevel2*  calo2_obj;
-     ToFLevel2*   tof_obj;
+     ToFLevel2*   tof2_obj;
      TrigLevel2*  trig_obj;
      S4Level2*    s4_obj;
      NDLevel2*    nd_obj;
      AcLevel2*    ac_obj;
-     OrbitalInfo* orb_obj;
+     OrbitalInfo* orb2_obj;
      GL_RUN*      run_obj;
      SoftInfo*    soft_obj; // Emiliano
+     ProcInfo*    proc_obj; //EM
+     Int_t extAlgFlag;
+     TClonesArray* trk_ext_obj;
+     TClonesArray* trk_ext_nuc_obj;
+     TClonesArray* trk_nuc_obj;
+     TClonesArray* calo_ext_obj;
+     TClonesArray* calo_ext_nuc_obj;
+     TClonesArray* calo_nuc_obj;
+     TClonesArray* tof_ext_obj;
+     TClonesArray* tof_ext_nuc_obj;
+     TClonesArray* tof_nuc_obj;
+     TClonesArray* orb_ext_obj;
+     TClonesArray* orb_ext_nuc_obj;
+     TClonesArray* orb_nuc_obj;
+     TrkLevel2*   trk2_nuc_obj;
+     CaloLevel2*  calo2_nuc_obj;
+     ToFLevel2*   tof2_nuc_obj;
+     OrbitalInfo* orb2_nuc_obj;
      // --------------------------------------------
      // variables to handle GP
-Line 340 
 class PamLevel2 : public TObject {
+Line 406 
 class PamLevel2 : public TObject {
      // -----------------------------------
      TClonesArray *tsorted;
      TClonesArray *timage;
-     PamTrack    *track; //store the last required
+     TClonesArray *text; //array of pamtrack objects built from extended algorythm
+     TClonesArray *tsorted_nuc;
+     TClonesArray *timage_nuc;
+     TClonesArray *text_nuc; //array of pamtrack objects built from extended algorythm
+     PamTrack    *track; //store the last required
      TString howtosort;
      Float_t sortthr;
+     Bool_t issorted;     // true if the current event has already been sorted with the current sort criteria
+     Long64_t lastsorted; // the last event whose tracks has been sorted
+     Bool_t issorted_new;     // true if the current event has already been sorted with the current sort criteria
+     Long64_t lastsorted_new; // the last event whose tracks has been sorted
      // --------------------------------------
      // flags to load level0-1-2 trees/braches
      // --------------------------------------
-Line 367 
 class PamLevel2 : public TObject {
+Line 445 
 class PamLevel2 : public TObject {
      Bool_t GP; //GP simulation tree
+     Bool_t EXT;//get tracks from extended tracking algorythm
+     Bool_t NUC;//get tracks from nuclei algorythm
+     Bool_t PROC;//EM
+     const char*  trkAlg;
      // -----------
      // other flags
      // -----------
-Line 374 
 class PamLevel2 : public TObject {
+Line 460 
 class PamLevel2 : public TObject {
      Bool_t ISGP;  // GP simulation (...sort of patch, to skip OBT check, since OrbitalInfo is sometime missing)
      Int_t  SELLI; // tree of selected events (there is a selection list)
+     Bool_t DBG;
      void Initialize();
   public:
+     PamLevel2(TString ddir,TList *list,TString detlist);
      PamLevel2(TString ddir,TString list,TString detlist);
      PamLevel2(TString ddir,TString list);
      PamLevel2();
      ~PamLevel2(){Delete();};
      void Clear();
      void Reset();
      void Delete();
  //  =============================================
  //  methods needed to read/write Level0-1-2 trees
  //  =============================================
-Line 396 
 class PamLevel2 : public TObject {
+Line 485 
 class PamLevel2 : public TObject {
      void    GetWhichTrees(TFile*);
      TTree*  GetPamTree(TFile*,TString);
      TTree*  GetPamTree(TFile* f){return GetPamTree(f,"");};
      TList*  GetListOfLevel2Files(TString, TString);
      Bool_t  CheckLevel2File(TString);
-Line 408 
 class PamLevel2 : public TObject {
+Line 497 
 class PamLevel2 : public TObject {
      TTree*  GetRunTree(TFile*);
      TChain* GetRunTree(TList*);
      TChain* GetRunTree(TString ddir,TString list){ return GetRunTree(GetListOfLevel2Files(ddir,list)); };
-     Bool_t  UpdateRunInfo(Long64_t);//{ return UpdateRunInfo(run_tree,iev); };
+     Bool_t  UpdateRunInfo(Long64_t);///< new implementation to solve Mori's bug
      Bool_t  UpdateRunInfo(TTree *,Long64_t);
      Bool_t  UpdateRunInfo(TChain *,Long64_t);
+     void    SkipRunInfoUpdate(); // Emiliano
+     void    DoRunInfoUpdate(){fUpdateRunInfo = true;}; // [default] Emiliano
+     void    UseDBinRunInfo(){fUseDBinRunInfo = true;}; // [default] Emiliano
+     void    NoDBinRunInfo(){fUseDBinRunInfo = false;}; // EM
+     void    NoDBconnections(){SkipRunInfoUpdate();}; // Emiliano
      void    CreateCloneTrees0(TChain *fChain, TFile* file);
      void    CreateCloneTrees(TFile* file);
      void    CreateCloneTrees(){CreateCloneTrees(NULL);};
-Line 424 
 class PamLevel2 : public TObject {
+Line 519 
 class PamLevel2 : public TObject {
  //    Int_t   GetEntry(Int_t);
      Int_t   GetEntry(Long64_t);
-     Int_t   GetYodaEntry();
+     Long64_t GetReadEntry();
+     Int_t   GetYodaEntry();
+     void SetMaxShift(Int_t sh); // EMILIANO method to change the maxshift variable
-     void SetMaxShift(Int_t sh){ maxshift = sh; return;}; // EMILIANO method to change the maxshift variable
      Long64_t GetEntries(){return pam_tree->GetEntries();};
      TChain* GetPamTree(){return pam_tree;};
      TChain* GetRunTree(){return run_tree;};
+     TChain* GetProcessingInfoTree(){return proc_tree;};
      TTree*  GetYodaTree();
      TFile*  GetYodaFile(){return l0_file;}; // can be useful to know the name of level0 file we are using or similar things
-Line 461 
 class PamLevel2 : public TObject {
+Line 558 
 class PamLevel2 : public TObject {
      TrkHough     *GetTrkHough()   { return trkh_obj; };
      CaloLevel1   *GetCaloLevel1() { return calo1_obj;};
      CaloLevel2   *GetCaloLevel2() { return calo2_obj;};
-     ToFLevel2    *GetToFLevel2()  { return tof_obj;  };
+     ToFLevel2    *GetToFLevel2()  { return tof2_obj;  };
      TrigLevel2   *GetTrigLevel2() { return trig_obj; };
      S4Level2     *GetS4Level2()   { return s4_obj;   };
      NDLevel2     *GetNDLevel2()   { return nd_obj;   };
      AcLevel2     *GetAcLevel2()   { return ac_obj;   };
-     OrbitalInfo  *GetOrbitalInfo(){ return orb_obj;  };
+     OrbitalInfo  *GetOrbitalInfo(){ return orb2_obj;  };
      GL_RUN       *GetRunInfo()    { return run_obj;  };
      SoftInfo     *GetSoftInfo()   { return soft_obj; }; // Emiliano
-     GPamela      *GetGPamela()    { return gp_obj; };
+     ProcInfo     *GetProcInfo()   { return proc_obj; }; // Emiliano
+     GPamela      *GetGPamela()    { return gp_obj;   };
  //  ==============================================
  //  methods to retrieve subdetectors stored tracks
-Line 477 
 class PamLevel2 : public TObject {
+Line 575 
 class PamLevel2 : public TObject {
      TrkTrack*   GetTrkStoredTrack(int seqno){return trk2_obj->TrkLevel2::GetStoredTrack(seqno); };
      CaloTrkVar* GetCaloStoredTrack(int);
      ToFTrkVar*  GetToFStoredTrack(int);
+     OrbitalInfoTrkVar*  GetOrbitalInfoStoredTrack(int);
      //
-     PamTrack*   GetPamTrackAlong(TrkTrack*);
+ //    PamTrack*   GetPamTrackAlong(TrkTrack*);
  //      PamTrack*   GetPamTrackAlong(CaloTrkVar*);      //se serve
  //      PamTrack*   GetPamTrackAlong(ToFTrkVar*);      //se serve
  //  =================================
  //  methods to retrieve pamela tracks
  //  =================================
-     PamTrack*     GetStoredTrack(Int_t);
+ //    PamTrack*     GetStoredTrack(Int_t);//obsolete
      //
      void          SortTracks(TString);          //assign value to variable sorted_tracks
      void          SortTracks();          //assign value to variable sorted_tracks
  //    TRefArray*    GetTracks();                  //return sorted_tracks
      TClonesArray* GetTracks();                  //return sorted_tracks
-     PamTrack*     GetTrack(int);
+     PamTrack*     GetTrackOld(int);
-     PamTrack*     GetTrackImage(int);
+     PamTrack*     GetTrackImageOld(int);
-     void          SetSortingMethod(TString how){ howtosort = how; return;};
+     PamTrack*     GetTrack(int i){ return GetTrackOld(i); };
+     PamTrack*     GetTrackImage(int i){ return GetTrackImageOld(i); };
+    void          SetSortingMethod(TString how);
      TString       GetSortingMethod(){return howtosort;};
      void          SetNucleiSortingThreshold(Float_t nucleithr){ sortthr = nucleithr; return;};
      Float_t       GetNucleiSortingThreshold(){ return sortthr;};
+ //  =================================
+ //  methods to handle new tracking algorythms
+ //  =================================
+ /*     Int_t GetTrackingAlgorythmFlag(){return extAlgFlag;}; */
+ /*     void  SetTrackingAlgorythm(const char * alg); */
+ /*     const char* GetTrackingAlgorythm(); */
+     Int_t  GetNTracks(const char*);
+     PamTrack*     GetTrack(int,const char*);
+     PamTrack*     GetTrackImage(int, const char*);
+     Bool_t TrackHasImage(int i, const char* alg){return GetTrackImage(i,alg)!=NULL;};
+ /*     Int_t  GetNTracks(){return GetNTracks(trkAlg); }; */
+ /*     PamTrack*     GetTrack(int i) { return GetTrack(i,trkAlg); }; */
+ /*     PamTrack*     GetTrackImage(int i){ return GetTrackImage(i,trkAlg); }; */
+     void          SortTracksNew();          //assign value to variable sorted_tracks
  //  =================================
  //  other methods
  //  =================================
  //    Int_t ntrk(){ return trk2_obj->TrkLevel2::ntrk(); };
+     Long64_t  GetRunID(){ return irun;};
+     Long64_t  GetRunEntryID(){ return irunentry; };
      // -----
      // flags methods
      // -----
-Line 526 
 class PamLevel2 : public TObject {
+Line 653 
 class PamLevel2 : public TObject {
      Bool_t IsRUN(){return RUN;};
-     Bool_t IsGood();
+     Bool_t IsEXT(){return EXT;};
+     Bool_t IsNUC(){return NUC;};
+     Bool_t IsGood(Bool_t strict=true);
      Int_t  IsSELLI(){return SELLI;}; //trees of selected events (there is a selection list)
+     void  SetSELLI(Int_t selli){SELLI=selli;};
+     Bool_t IsPROC(){return PROC;};//EM
+     Bool_t HasDiscardedFiles(){return fDiscarded;}; ///< true if something went wrong and not all files were processed
      Bool_t IsGP(){return ISGP;};
      void   SetGP(Bool_t tt){ISGP=tt;};
+     void SetDebug(Bool_t set){ DBG=set; TrkParams::SetVerboseMode(); };
-     ClassDef(PamLevel2,1);
+     void SetCustomString(TString cus){customString=cus;};
+     TString GetCustomString(){return customString;};
+     ClassDef(PamLevel2,4);
  };
  //===============================================================================================

 Legend:



Removed from v.1.32
 


changed lines


 
Added in v.1.54
 Legend:



Removed from v.1.32
 


changed lines


 
Added in v.1.54
-Removed from v.1.32
+Added in v.1.54

	ViewVC Help
Powered by ViewVC 1.1.23