1 #ifndef ALIPERFORMANCEEFF_H
2 #define ALIPERFORMANCEEFF_H
4 //------------------------------------------------------------------------------
5 // Class to keep information from comparison of
6 // reconstructed and MC particle tracks (TPC efficiency).
8 // Author: J.Otwinowski 04/02/2008
9 //------------------------------------------------------------------------------
29 #include "AliPerformanceObject.h"
31 class AliPerformanceEff : public AliPerformanceObject {
33 AliPerformanceEff(const Char_t* name="AliPerformanceEff",const Char_t*title="AliPerformanceEff",Int_t analysisMode=0, Bool_t hptGenerator=kFALSE);
34 virtual ~AliPerformanceEff();
40 virtual void Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent, AliESDfriend *const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);
42 // Merge output objects (needed by PROOF)
43 virtual Long64_t Merge(TCollection* const list);
45 // Analyse output histograms
46 virtual void Analyse();
48 // Get analysis folder
49 virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
51 // Create folder for analysed histograms
52 TFolder *CreateFolder(TString folder = "folderEff",TString title = "Analysed Efficiency histograms");
54 // Export objects to folder
55 TFolder *ExportToFolder(TObjArray *array=0);
58 void ProcessTPC(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent);
59 void ProcessTPCITS(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent);
60 void ProcessConstrained(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent);
61 void ProcessTPCSec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent);
63 Bool_t IsRecTPC(AliESDtrack *track);
64 Bool_t IsRecTPCITS(AliESDtrack *track);
65 Bool_t IsRecConstrained(AliESDtrack *track);
67 Bool_t IsFindable(const AliMCEvent *mcEvent, Int_t label);
68 Int_t TransformToPID(TParticle *particle);
71 void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
72 void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
75 AliRecInfoCuts* GetAliRecInfoCuts() const {return fCutsRC;}
76 AliMCInfoCuts* GetAliMCInfoCuts() const {return fCutsMC;}
78 THnSparseF* GetEffHisto() const {return fEffHisto;}
79 THnSparseF* GetEffSecHisto() const {return fEffSecHisto;}
83 static const Int_t fgkMaxClones = 3, fgkMaxFakes = 3;
86 TH1D* AddHistoEff(Int_t axis, const Char_t *name, const Char_t* vsTitle, const Int_t type, const Int_t secondary = 0);
89 THnSparseF *fEffHisto; //-> mceta:mcphi:mcpt:pid:isPrim:recStatus:findable:charge
90 THnSparseF *fEffSecHisto; //-> mceta:mcphi:mcpt:pid:isPrim:recStatus:findable:mcR:mother_phi:mother_eta:charge
92 // Global cuts objects
93 AliRecInfoCuts* fCutsRC; // selection cuts for reconstructed tracks
94 AliMCInfoCuts* fCutsMC; // selection cuts for MC tracks
97 TFolder *fAnalysisFolder; // folder for analysed histograms
99 AliPerformanceEff(const AliPerformanceEff&); // not implemented
100 AliPerformanceEff& operator=(const AliPerformanceEff&); // not implemented
102 ClassDef(AliPerformanceEff,2);