1 #ifndef ALIPERFORMANCEMC_H
2 #define ALIPERFORMANCEMC_H
4 //------------------------------------------------------------------------------
5 // Class to keep information for MC tracks
6 // to check propagation algorithm, B-field and material budget.
8 // Author: J.Otwinowski 09/06/2009
9 //------------------------------------------------------------------------------
22 class AliTrackReference;
25 class AliESDfriendTrack;
30 class AliExternalTrackParam;
32 #include "THnSparse.h"
33 #include "AliPerformanceObject.h"
35 class AliPerformanceMC : public AliPerformanceObject {
37 AliPerformanceMC(const Char_t* name="AliPerformanceMC", const Char_t* title="AliPerformanceMC",Int_t analysisMode=0,Bool_t hptGenerator=kFALSE);
38 virtual ~AliPerformanceMC();
44 virtual void Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent, AliESDfriend *const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);
46 // Merge output objects (needed by PROOF)
47 virtual Long64_t Merge(TCollection* const list);
49 // Analyse output histograms
50 virtual void Analyse();
52 // Get analysis folder
53 virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
56 void ProcessTPC(AliTrackReference* const refIn, TParticle* const particle);
57 void ProcessInnerTPC(AliTrackReference* const refIn, AliTrackReference* const refOut, TParticle* const particle);
58 void ProcessOuterTPCExt(TParticle *const part, TClonesArray * const trefs);
59 AliExternalTrackParam * MakeTrack(const AliTrackReference* ref, TParticle *const part);
61 // Create folder for analysed histograms
62 TFolder *CreateFolder(TString folder = "folderMC", TString title = "Analysed MC histograms");
64 // Export objects to folder
65 TFolder *ExportToFolder(TObjArray * array=0);
68 void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
69 void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
71 AliRecInfoCuts* GetAliRecInfoCuts() const {return fCutsRC;}
72 AliMCInfoCuts* GetAliMCInfoCuts() const {return fCutsMC;}
74 TH1F* MakeResol(TH2F * his, Int_t integ=0, Bool_t type=kFALSE, Int_t cut=0);
78 THnSparse *GetResolHisto() const { return fResolHisto; }
79 THnSparse *GetPullHisto() const { return fPullHisto; }
84 // 5 track parameters (details in STEER/AliExternalTrackParam.h)
87 // resolution histogram
88 THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_pt:y:z:phi:eta:pt
91 THnSparseF *fPullHisto; //-> pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt
93 // Global cuts objects
94 AliRecInfoCuts* fCutsRC; // selection cuts for reconstructed tracks
95 AliMCInfoCuts* fCutsMC; // selection cuts for MC tracks
98 TFolder *fAnalysisFolder; // folder for analysed histograms
100 AliPerformanceMC(const AliPerformanceMC&); // not implemented
101 AliPerformanceMC& operator=(const AliPerformanceMC&); // not implemented
103 ClassDef(AliPerformanceMC,1);