1 #ifndef ALIPERFORMANCERES_H
2 #define ALIPERFORMANCERES_H
4 //------------------------------------------------------------------------------
5 // Class to keep information from comparison of
6 // reconstructed and MC particle tracks (TPC resolution).
8 // Author: J.Otwinowski 04/02/2008
9 //------------------------------------------------------------------------------
21 class AliTrackReference;
24 class AliESDfriendTrack;
30 #include "THnSparse.h"
31 #include "AliPerformanceObject.h"
33 class AliPerformanceRes : public AliPerformanceObject {
36 AliPerformanceRes(Char_t* name, Char_t* title, Int_t analysisMode, Bool_t hptGenerator);
37 virtual ~AliPerformanceRes();
43 virtual void Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent, AliESDfriend *const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);
45 // Merge output objects (needed by PROOF)
46 virtual Long64_t Merge(TCollection* const list);
48 // Analyse output histograms
49 virtual void Analyse();
51 // Get analysis folder
52 virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
55 void ProcessConstrained(AliStack* const stack, AliESDtrack *const esdTrack);
56 void ProcessTPC(AliStack* const stack, AliESDtrack *const esdTrack);
57 void ProcessTPCITS(AliStack* const stack, AliESDtrack *const esdTrack);
58 void ProcessInnerTPC(AliMCEvent *const mcEvent, AliESDtrack *const esdTrack);
59 void ProcessOuterTPC(AliMCEvent *const mcEvent, AliESDtrack *const esdTrack, AliESDfriendTrack *const friendTrack);
61 AliTrackReference *GetFirstTPCTrackRef(AliMCParticle *mcParticle);
62 AliTrackReference *GetLastTPCTrackRef(AliMCParticle *mcParticle);
64 // Create folder for analysed histograms
65 TFolder *CreateFolder(TString folder = "folderRes",TString title = "Analysed Resolution histograms");
67 // Export objects to folder
68 TFolder *ExportToFolder(TObjArray * array=0);
71 void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
72 void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
74 AliRecInfoCuts* GetAliRecInfoCuts() const {return fCutsRC;}
75 AliMCInfoCuts* GetAliMCInfoCuts() const {return fCutsMC;}
77 TH1F* MakeResol(TH2F * his, Int_t integ=0, Bool_t type=kFALSE, Int_t cut=0);
81 THnSparse *GetResolHisto() const { return fResolHisto; }
82 THnSparse *GetPullHisto() const { return fPullHisto; }
86 // 5 track parameters (details in STEER/AliExternalTrackParam.h)
89 // resolution histogram
90 THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_pt:y:z:phi:eta:pt
93 //THnSparseF *fPullHisto; //-> pull_y:pull_z:pull_phi:pull_lambda:pull_1pt:y:z:eta:phi:pt
94 THnSparseF *fPullHisto; //-> pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt
96 // Global cuts objects
97 AliRecInfoCuts* fCutsRC; // selection cuts for reconstructed tracks
98 AliMCInfoCuts* fCutsMC; // selection cuts for MC tracks
101 TFolder *fAnalysisFolder; // folder for analysed histograms
103 AliPerformanceRes(const AliPerformanceRes&); // not implemented
104 AliPerformanceRes& operator=(const AliPerformanceRes&); // not implemented
106 ClassDef(AliPerformanceRes,1);