1 #ifndef ALICOMPARISONRES_H
2 #define ALICOMPARISONRES_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 //------------------------------------------------------------------------------
25 #include "THnSparse.h"
26 #include "AliComparisonObject.h"
28 class AliComparisonRes : public AliComparisonObject {
31 AliComparisonRes(Char_t* name, Char_t* title, Int_t analysisMode, Bool_t hptGenerator);
32 virtual ~AliComparisonRes();
38 virtual void Exec(AliMCInfo* const infoMC, AliESDRecInfo *const infoRC);
40 // Merge output objects (needed by PROOF)
41 virtual Long64_t Merge(TCollection* const list);
43 // Analyse output histograms
44 virtual void Analyse();
46 // Get analysis folder
47 virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
50 void ProcessConstrained(AliMCInfo* const infoMC, AliESDRecInfo* const infoRC);
51 void ProcessTPC(AliMCInfo* const infoMC, AliESDRecInfo* const infoRC);
52 void ProcessTPCITS(AliMCInfo* const infoMC, AliESDRecInfo* const infoRC);
54 // Create folder for analysed histograms
55 TFolder *CreateFolder(TString folder = "folderRes",TString title = "Analysed Resolution histograms");
57 // Export objects to folder
58 TFolder *ExportToFolder(TObjArray * array=0);
61 void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
62 void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
64 AliRecInfoCuts* GetAliRecInfoCuts() const {return fCutsRC;}
65 AliMCInfoCuts* GetAliMCInfoCuts() const {return fCutsMC;}
67 static TH1F* MakeResol(TH2F * his, Int_t integ, Bool_t type);
71 THnSparse *GetResolHisto() const { return fResolHisto; }
72 THnSparse *GetPullHisto() const { return fPullHisto; }
77 // 5 track parameters (details in STEER/AliExternalTrackParam.h)
80 // resolution histogram
81 THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_1pt:y:z:eta:phi:pt
84 THnSparseF *fPullHisto; //-> pull_y:pull_z:pull_phi:pull_lambda:pull_1pt:y:z:eta:phi:pt
86 // Global cuts objects
87 AliRecInfoCuts* fCutsRC; // selection cuts for reconstructed tracks
88 AliMCInfoCuts* fCutsMC; // selection cuts for MC tracks
91 TFolder *fAnalysisFolder; // folder for analysed histograms
93 AliComparisonRes(const AliComparisonRes&); // not implemented
94 AliComparisonRes& operator=(const AliComparisonRes&); // not implemented
96 ClassDef(AliComparisonRes,1);