PWGPP/TPC/AliPerformanceMatch.h

   1 #ifndef ALIPERFORMANCEMATCH_H
   2 #define ALIPERFORMANCEMATCH_H
   3
   4 //------------------------------------------------------------------------------
   5 // Class keeps matching information between
   6 // central barrel detectors.
   7 //
   8 // Author: J.Otwinowski 17/10/2009
   9 // Changes by M.Knichel 22/10/2010
  10 //------------------------------------------------------------------------------
  11
  12 class TString;
  13 class TNamed;
  14 class TCanvas;
  15 class TH1F;
  16 class TH2F;
  17
  18 class AliESDVertex;
  19 class AliESDtrack;
  20 class AliMCEvent;
  21 class AliStack;
  22 class AliTrackReference;
  23 class AliESDEvent;
  24 class AliESDfriend;
  25 class AliESDfriendTrack;
  26 class AliMCEvent;
  27 class AliMCParticle;
  28 class AliMCInfoCuts;
  29 class AliRecInfoCuts;
  30 class AliExternalTrackParam;
  31
  32 #include "THnSparse.h"
  33 #include "AliPerformanceObject.h"
  34
  35 class AliPerformanceMatch : public AliPerformanceObject {
  36 public :
  37   AliPerformanceMatch(const Char_t* name="AliPerformanceMatch", const Char_t* title="AliPerformanceMatch",Int_t analysisMode=0,Bool_t hptGenerator=kFALSE);
  38   virtual ~AliPerformanceMatch();
  39
  40   // Init data members
  41   virtual void  Init();
  42
  43   // Execute analysis
  44   virtual void  Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent,AliESDfriend *const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);
  45
  46   // Merge output objects (needed by PROOF)
  47   virtual Long64_t Merge(TCollection* const list);
  48
  49   // Project Histograms store in AnalysisFolder
  50   virtual void Analyse();
  51
  52   // Analyse Projected Histograms to create Efficiency and AddToFolder
  53   //  virtual void AnalyseFinal();
  54
  55   // Get analysis folder
  56   virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
  57
  58   // Process matching
  59   void ProcessTPCITS(AliStack* const stack, AliESDtrack *const esdTrack);
  60   void ProcessTPCTRD(AliStack* const stack, AliESDtrack *const esdTrack, AliESDfriendTrack *const friendTrack);
  61   void ProcessITSTPC(Int_t trackIdx, AliESDEvent* const esdEvent, AliStack* const stack, AliESDtrack *const esdTrack);
  62   void ProcessTPCConstrain(AliStack* const stack, AliESDEvent *const esdEvent, AliESDtrack *const esdTrack); // - 01.11.2011
  63
  64   // Fill histogrrams
  65   void FillHistograms(AliESDtrack *const refParam, AliESDtrack *const param, Bool_t isRec);
  66
  67   // Create folder for analysed histograms
  68   TFolder *CreateFolder(TString folder = "folderRes",TString title = "Analysed Resolution histograms");
  69
  70   // Export objects to folder
  71   TFolder *ExportToFolder(TObjArray * array=0);
  72
  73   // Selection cuts
  74   void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
  75   void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
  76
  77   AliRecInfoCuts*  GetAliRecInfoCuts() const {return fCutsRC;}
  78   AliMCInfoCuts*   GetAliMCInfoCuts()  const {return fCutsMC;}
  79
  80   TH1F*  MakeResol(TH2F * his, Int_t integ=0, Bool_t type=kFALSE, Int_t cut=0);
  81
  82   // getters
  83   //
  84   THnSparse *GetResolHisto() const  { return fResolHisto; }
  85   THnSparse *GetPullHisto()  const  { return fPullHisto; }
  86   THnSparse *GetTrackEffHisto() const  { return fTrackingEffHisto; }
  87   THnSparse *GetTPCConstrain() const { return fTPCConstrain; }
  88
  89   TObjArray* GetHistos() const { return fFolderObj; }
  90
  91   static Bool_t GetMergeTHnSparse() { return fgMergeTHnSparse; }
  92   static void SetMergeTHnSparse(Bool_t mergeTHnSparse) {fgUseMergeTHnSparse = kTRUE; fgMergeTHnSparse = mergeTHnSparse; }
  93
  94   void SetUseHLT(Bool_t useHLT = kTRUE) {fUseHLT = useHLT;}
  95   Bool_t GetUseHLT() { return fUseHLT; }
  96
  97 private:
  98
  99   static Bool_t fgMergeTHnSparse;
 100   static Bool_t fgUseMergeTHnSparse;
 101
 102   //
 103   // Control histograms
 104   // 5 track parameters (details in STEER/AliExternalTrackParam.h)
 105   // + isRec flag to determine ITS/TRD tracking efficiency
 106   // w.r.t TPC
 107
 108   // resolution histogram
 109   THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_pt:y:z:phi:eta:pt:isRec
 110
 111   // pull histogram
 112   THnSparseF *fPullHisto;  //-> pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt:isRec
 113
 114   // tracking efficiency using ITS stand-alone tracks histogram
 115   THnSparseF *fTrackingEffHisto;  //-> has match:y:z:snp:tgl:phi:pt:ITSclusters
 116
 117   // TPC Inner constrained to global tracks - 01.11.2011
 118   THnSparseF *fTPCConstrain;  //-> pull_phi:phi,pt,eta
 119
 120
 121   TObjArray* fFolderObj; // array of analysed histograms
 122
 123   // Global cuts objects
 124   AliRecInfoCuts*  fCutsRC;      // selection cuts for reconstructed tracks
 125   AliMCInfoCuts*  fCutsMC;       // selection cuts for MC tracks
 126
 127   // analysis folder
 128   TFolder *fAnalysisFolder; // folder for analysed histograms
 129
 130   Bool_t fUseHLT; // use HLT ESD
 131
 132   AliPerformanceMatch(const AliPerformanceMatch&); // not implemented
 133   AliPerformanceMatch& operator=(const AliPerformanceMatch&); // not implemented
 134
 135   ClassDef(AliPerformanceMatch,3);
 136 };
 137
 138 #endif