don't lie in the log!
[u/mrichter/AliRoot.git] / 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