[u/mrichter/AliRoot.git] / PWG0 / dNdPt / AlidNdPtAnalysis.h

#ifndef ALIDEDPTANALYSIS_H
#define ALIDEDPTANALYSIS_H

//------------------------------------------------------------------------------
// AlidNdPtAnalysis class. 
// 
// a. functionality:
// - fills analysis control histograms
// - fills generic correction matrices 
// - generates correction matrices 
//
// b. data members:
// - generic correction matrices
// - control histograms
//
// Author: J.Otwinowski 04/11/2008 
//------------------------------------------------------------------------------

class TProfile;
class TFolder;
class TObjArray;
class TString;

class AliESDtrackCuts;
class AliVertexerTracks;
class AliESD;
class AliESDfriend;
class AliESDfriendTrack;

class AlidNdPtHelper;

#include "THnSparse.h"
#include "AlidNdPt.h"

class AlidNdPtAnalysis : public AlidNdPt {
public :
  AlidNdPtAnalysis(); 
  AlidNdPtAnalysis(Char_t* name, Char_t* title);
  ~AlidNdPtAnalysis();

  // Init data members
  virtual void Init();

  // Process events
  virtual void Process(AliESDEvent *const esdEvent=0, AliMCEvent *const mcEvent=0);

  // Merge output objects (needed by PROOF) 
  virtual Long64_t Merge(TCollection* list);

  // Analyse output histograms 
  virtual void Analyse();

  // Export objects to folder
  virtual TFolder *ExportToFolder(TObjArray * array=0);

  // Get analysis folder
  TFolder* GetAnalysisFolder() {return fAnalysisFolder;}

  // Fill control histograms
  void SetHistogramsOn(Bool_t histOn=kTRUE) {fHistogramsOn = histOn;}
  Bool_t IsHistogramsOn() {return fHistogramsOn;}

  // Create folder for analysed histograms
  TFolder *CreateFolder(TString folder = "folderdNdPtAnalysis",TString title = "Analysed dNdPt histograms");

  // Fill histograms
  void FillHistograms(AliESDtrack *const esdTrack, AliStack *const stack, AlidNdPtHelper::TrackObject trackObj);
  void FillHistograms(AliStack *const stack, Int_t label, AlidNdPtHelper::TrackObject trackObj);
  void FillHistograms(TObjArray *const allChargedTracks,Int_t *const labelsAll,Int_t multAll,Int_t *const labelsAcc,Int_t multAcc,Int_t *const labelsRec,Int_t multRec);

  // Getters
  THnSparseF *GetEventMultCorrelationMatrix() {return fEventMultCorrelationMatrix;}
  THnSparseF *GetTrackPtCorrelationMatrix()   {return fTrackPtCorrelationMatrix;}

  //
  THnSparseF *GetGenEventMatrix() {return fGenEventMatrix;}
  THnSparseF *GetGenSDEventMatrix() {return fGenSDEventMatrix;}
  THnSparseF *GetGenDDEventMatrix() {return fGenDDEventMatrix;}
  THnSparseF *GetGenNDEventMatrix() {return fGenNDEventMatrix;}
  THnSparseF *GetGenNSDEventMatrix() {return fGenNSDEventMatrix;}

  THnSparseF *GetTriggerEventMatrix() {return fTriggerEventMatrix;}
  THnSparseF *GetTriggerSDEventMatrix() {return fTriggerSDEventMatrix;}
  THnSparseF *GetTriggerDDEventMatrix() {return fTriggerDDEventMatrix;}
  THnSparseF *GetTriggerNDEventMatrix() {return fTriggerNDEventMatrix;}
  THnSparseF *GetTriggerNSDEventMatrix() {return fTriggerNSDEventMatrix;}

  THnSparseF *GetRecEventMatrix() {return fRecEventMatrix;}
  THnSparseF *GetRecSDEventMatrix() {return fRecSDEventMatrix;}
  THnSparseF *GetRecDDEventMatrix() {return fRecDDEventMatrix;}
  THnSparseF *GetRecNDEventMatrix() {return fRecNDEventMatrix;}
  THnSparseF *GetRecNSDEventMatrix() {return fRecNSDEventMatrix;}

  // 
  THnSparseF *GetGenTrackEventMatrix() {return fGenTrackEventMatrix;}
  THnSparseF *GetGenTrackSDEventMatrix() {return fGenTrackSDEventMatrix;}
  THnSparseF *GetGenTrackDDEventMatrix() {return fGenTrackDDEventMatrix;}
  THnSparseF *GetGenTrackNDEventMatrix() {return fGenTrackNDEventMatrix;}
  THnSparseF *GetGenTrackNSDEventMatrix() {return fGenTrackNSDEventMatrix;}

  THnSparseF *GetTriggerTrackEventMatrix() {return fTriggerTrackEventMatrix;}
  THnSparseF *GetTriggerTrackSDEventMatrix() {return fTriggerTrackSDEventMatrix;}
  THnSparseF *GetTriggerTrackDDEventMatrix() {return fTriggerTrackDDEventMatrix;}
  THnSparseF *GetTriggerTrackNDEventMatrix() {return fTriggerTrackNDEventMatrix;}
  THnSparseF *GetTriggerTrackNSDEventMatrix() {return fTriggerTrackNSDEventMatrix;}

  THnSparseF *GetRecTrackEventMatrix() {return fRecTrackEventMatrix;}
  THnSparseF *GetRecTrackSDEventMatrix() {return fRecTrackSDEventMatrix;}
  THnSparseF *GetRecTrackDDEventMatrix() {return fRecTrackDDEventMatrix;}
  THnSparseF *GetRecTrackNDEventMatrix() {return fRecTrackNDEventMatrix;}
  THnSparseF *GetRecTrackNSDEventMatrix() {return fRecTrackNSDEventMatrix;}

  //
  THnSparseF *GetGenPrimTrackMatrix() {return fGenPrimTrackMatrix;}
  THnSparseF *GetRecPrimTrackMatrix() {return fRecPrimTrackMatrix;}

  THnSparseF *GetRecTrackMatrix() {return fRecTrackMatrix;}
  THnSparseF *GetRecSecTrackMatrix() {return fRecSecTrackMatrix;}
  THnSparseF *GetRecMultTrackMatrix() {return fRecMultTrackMatrix;}

  //
  // control histograms
  //
  THnSparseF *GetMCEventHist1() {return fMCEventHist1;}
  THnSparseF *GetRecEventHist1() {return fRecEventHist1;}
  THnSparseF *GetRecEventHist2() {return fRecEventHist2;}
  THnSparseF *GetRecMCEventHist1() {return fRecMCEventHist1;}
  THnSparseF *GetRecMCEventHist2() {return fRecMCEventHist2;}
  THnSparseF *GetRecMCEventHist3() {return fRecMCEventHist3;}

  THnSparseF *GetMCTrackHist1(Int_t i) {return fMCTrackHist1[i];}
  THnSparseF *GetMCPrimTrackHist1(Int_t i) {return fMCPrimTrackHist1[i];}
  THnSparseF *GetMCSecTrackHist1(Int_t i) {return fMCSecTrackHist1[i];}

  THnSparseF *GetRecTrackHist1(Int_t i) {return fRecTrackHist1[i];}
  THnSparseF *GetRecTrackMultHist1(Int_t i) {return fRecTrackMultHist1[i];}

  THnSparseF *GetRecMCTrackHist1() {return fRecMCTrackHist1;}
  THnSparseF *GetMCMultRecTrackHist1() {return fMCMultRecTrackHist1;}

private:

  // analysis folder 
  TFolder *fAnalysisFolder; // folder for analysed histograms

  // switch on/off filling of control histograms
  Bool_t fHistogramsOn; 

  // 
  // correlation matrices (histograms)
  //

  // event rec. track vs true track multiplicity correlation matrix 
  THnSparseF *fEventMultCorrelationMatrix; //-> mult:mult_true_tracks

  // rec. track pt vs true track pt correlation matrix for given eta
  THnSparseF *fTrackPtCorrelationMatrix; //-> Pt:mcPt:mcEta

  //
  // event level correction 
  //

  // all genertated
  THnSparseF *fGenEventMatrix; //-> mcZv:multMB (inelastic)
  THnSparseF *fGenSDEventMatrix; //-> mcZv:multMB (single diffractive)
  THnSparseF *fGenDDEventMatrix; //-> mcZv:multMB (single diffractive)
  THnSparseF *fGenNDEventMatrix; //-> mcZv:multMB (non diffractive)
  THnSparseF *fGenNSDEventMatrix; //-> mcZv:multMB (non single diffractive)

  // trigger bias corrections (fTriggerEventMatrix / fGenEventMatrix)
  THnSparseF *fTriggerEventMatrix; //-> mcZv:multMB
  THnSparseF *fTriggerSDEventMatrix; //-> mcZv:multMB
  THnSparseF *fTriggerDDEventMatrix; //-> mcZv:multMB
  THnSparseF *fTriggerNDEventMatrix; //-> mcZv:multMB
  THnSparseF *fTriggerNSDEventMatrix; //-> mcZv:multMB

  // event vertex rec. eff correction (fRecEventMatrix / fTriggerEventMatrix)
  THnSparseF *fRecEventMatrix; //-> mcZv:multMB 
  THnSparseF *fRecSDEventMatrix; //-> mcZv:multMB
  THnSparseF *fRecDDEventMatrix; //-> mcZv:multMB
  THnSparseF *fRecNDEventMatrix; //-> mcZv:multMB
  THnSparseF *fRecNSDEventMatrix; //-> mcZv:multMB

  //
  // track-event level correction 
  //

  THnSparseF *fGenTrackEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fGenTrackSDEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fGenTrackDDEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fGenTrackNDEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fGenTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta

  // trigger bias corrections (fTriggerTrackEventMatrix / fGenTrackEventMatrix)
  THnSparseF *fTriggerTrackEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fTriggerTrackSDEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fTriggerTrackDDEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fTriggerTrackNDEventMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fTriggerTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta

  // event vertex rec. corrections (fRecTrackEventMatrix / fTriggerTrackEventMatrix)
  THnSparseF *fRecTrackEventMatrix; //-> mcZv:Pt:mcEta
  THnSparseF *fRecTrackSDEventMatrix; //-> mcZv:Pt:mcEta
  THnSparseF *fRecTrackDDEventMatrix; //-> mcZv:Pt:mcEta
  THnSparseF *fRecTrackNDEventMatrix; //-> mcZv:Pt:mcEta
  THnSparseF *fRecTrackNSDEventMatrix; //-> mcZv:Pt:mcEta

  //
  // track level correction 
  //

  // track rec. efficiency correction (fRecPrimTrackMatrix / fGenPrimTrackMatrix)
  THnSparseF *fGenPrimTrackMatrix; //-> mcZv:mcPt:mcEta
  THnSparseF *fRecPrimTrackMatrix; //-> mcZv:mcPt:mcEta

  // secondary track contamination correction (fRecSecTrackMatrix / fRecTrackMatrix)
  THnSparseF *fRecTrackMatrix;    //-> mcZv:mcPt:mcEta
  THnSparseF *fRecSecTrackMatrix; //-> mcZv:mcPt:mcEta

  // multiple rec. track corrections (fRecMultTrackMatrix / fRecTrackMatrix)
  THnSparseF *fRecMultTrackMatrix; //-> mcZv:Pt:mcEta

  //
  // ESD and MC control analysis histograms
  //

  // THnSparse event histograms
  THnSparseF *fMCEventHist1;  //-> mcXv:mcYv:mcZv
  THnSparseF *fRecEventHist1; //-> Xv:Yv:Zv
  THnSparseF *fRecEventHist2; //-> Zv:multMB:mult
  THnSparseF *fRecMCEventHist1; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv
  THnSparseF *fRecMCEventHist2; //-> Xv-mcXv:Zv-mcZv:mult
  THnSparseF *fRecMCEventHist3; //-> mult:EventType (ND, DD, SD)

  // THnSparse track histograms
  // [0] - after charged track selection, [1] - after acceptance cuts, [2] - after esd track cuts

  THnSparseF *fMCTrackHist1[AlidNdPtHelper::kCutSteps];     //-> mcPt:mcEta:mcPhi
  THnSparseF *fMCPrimTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:pid:mech:mother
  THnSparseF *fMCSecTrackHist1[AlidNdPtHelper::kCutSteps];  //-> mcPt:mcEta:pid:mech:mother

  THnSparseF *fRecTrackHist1[AlidNdPtHelper::kCutSteps];     //-> Pt:Eta:Phi
  THnSparseF *fRecTrackMultHist1[AlidNdPtHelper::kCutSteps]; //-> Pt:mult

  THnSparseF *fRecMCTrackHist1; //-> mcPt:mcEta:(Pt-mcPt)/mcPt:(Eta-mcEta)

  //multple reconstructed tracks
  THnSparseF *fMCMultRecTrackHist1; //-> mcPt:mcEta:pid

  AlidNdPtAnalysis(const AlidNdPtAnalysis&); // not implemented
  AlidNdPtAnalysis& operator=(const AlidNdPtAnalysis&); // not implemented

  ClassDef(AlidNdPtAnalysis,1);
};

#endif
Commit	Line	Data
0aaa8b91	1	#ifndef ALIDEDPTANALYSIS_H
	2	#define ALIDEDPTANALYSIS_H
	3
	4	//------------------------------------------------------------------------------
	5	// AlidNdPtAnalysis class.
	6	//
	7	// a. functionality:
	8	// - fills analysis control histograms
	9	// - fills generic correction matrices
	10	// - generates correction matrices
	11	//
	12	// b. data members:
	13	// - generic correction matrices
	14	// - control histograms
	15	//
	16	// Author: J.Otwinowski 04/11/2008
	17	//------------------------------------------------------------------------------
	18
	19	class TProfile;
	20	class TFolder;
	21	class TObjArray;
	22	class TString;
	23
	24	class AliESDtrackCuts;
	25	class AliVertexerTracks;
	26	class AliESD;
	27	class AliESDfriend;
	28	class AliESDfriendTrack;
	29
bad4ba69	30	class AlidNdPtHelper;
bad4ba69	31
0aaa8b91	32	#include "THnSparse.h"
0aaa8b91	33	#include "AlidNdPt.h"
0aaa8b91	34
	35	class AlidNdPtAnalysis : public AlidNdPt {
	36	public :
	37	AlidNdPtAnalysis();
	38	AlidNdPtAnalysis(Char_t* name, Char_t* title);
	39	~AlidNdPtAnalysis();
	40
	41	// Init data members
	42	virtual void Init();
	43
	44	// Process events
	45	virtual void Process(AliESDEvent const esdEvent=0, AliMCEvent const mcEvent=0);
	46
	47	// Merge output objects (needed by PROOF)
	48	virtual Long64_t Merge(TCollection* list);
	49
	50	// Analyse output histograms
	51	virtual void Analyse();
	52
	53	// Export objects to folder
	54	virtual TFolder ExportToFolder(TObjArray array=0);
	55
	56	// Get analysis folder
	57	TFolder* GetAnalysisFolder() {return fAnalysisFolder;}
	58
	59	// Fill control histograms
	60	void SetHistogramsOn(Bool_t histOn=kTRUE) {fHistogramsOn = histOn;}
	61	Bool_t IsHistogramsOn() {return fHistogramsOn;}
	62
	63	// Create folder for analysed histograms
	64	TFolder *CreateFolder(TString folder = "folderdNdPtAnalysis",TString title = "Analysed dNdPt histograms");
	65
	66	// Fill histograms
	67	void FillHistograms(AliESDtrack const esdTrack, AliStack const stack, AlidNdPtHelper::TrackObject trackObj);
	68	void FillHistograms(AliStack *const stack, Int_t label, AlidNdPtHelper::TrackObject trackObj);
	69	void FillHistograms(TObjArray const allChargedTracks,Int_t const labelsAll,Int_t multAll,Int_t const labelsAcc,Int_t multAcc,Int_t const labelsRec,Int_t multRec);
	70
	71	// Getters
	72	THnSparseF *GetEventMultCorrelationMatrix() {return fEventMultCorrelationMatrix;}
	73	THnSparseF *GetTrackPtCorrelationMatrix() {return fTrackPtCorrelationMatrix;}
	74
	75	//
	76	THnSparseF *GetGenEventMatrix() {return fGenEventMatrix;}
	77	THnSparseF *GetGenSDEventMatrix() {return fGenSDEventMatrix;}
	78	THnSparseF *GetGenDDEventMatrix() {return fGenDDEventMatrix;}
	79	THnSparseF *GetGenNDEventMatrix() {return fGenNDEventMatrix;}
	80	THnSparseF *GetGenNSDEventMatrix() {return fGenNSDEventMatrix;}
	81
	82	THnSparseF *GetTriggerEventMatrix() {return fTriggerEventMatrix;}
	83	THnSparseF *GetTriggerSDEventMatrix() {return fTriggerSDEventMatrix;}
	84	THnSparseF *GetTriggerDDEventMatrix() {return fTriggerDDEventMatrix;}
	85	THnSparseF *GetTriggerNDEventMatrix() {return fTriggerNDEventMatrix;}
	86	THnSparseF *GetTriggerNSDEventMatrix() {return fTriggerNSDEventMatrix;}
	87
	88	THnSparseF *GetRecEventMatrix() {return fRecEventMatrix;}
	89	THnSparseF *GetRecSDEventMatrix() {return fRecSDEventMatrix;}
	90	THnSparseF *GetRecDDEventMatrix() {return fRecDDEventMatrix;}
	91	THnSparseF *GetRecNDEventMatrix() {return fRecNDEventMatrix;}
	92	THnSparseF *GetRecNSDEventMatrix() {return fRecNSDEventMatrix;}
	93
	94	//
	95	THnSparseF *GetGenTrackEventMatrix() {return fGenTrackEventMatrix;}
	96	THnSparseF *GetGenTrackSDEventMatrix() {return fGenTrackSDEventMatrix;}
	97	THnSparseF *GetGenTrackDDEventMatrix() {return fGenTrackDDEventMatrix;}
98	THnSparseF *GetGenTrackNDEventMatrix() {return fGenTrackNDEventMatrix;}
99	THnSparseF *GetGenTrackNSDEventMatrix() {return fGenTrackNSDEventMatrix;}
100
101	THnSparseF *GetTriggerTrackEventMatrix() {return fTriggerTrackEventMatrix;}
102	THnSparseF *GetTriggerTrackSDEventMatrix() {return fTriggerTrackSDEventMatrix;}
103	THnSparseF *GetTriggerTrackDDEventMatrix() {return fTriggerTrackDDEventMatrix;}
104	THnSparseF *GetTriggerTrackNDEventMatrix() {return fTriggerTrackNDEventMatrix;}
105	THnSparseF *GetTriggerTrackNSDEventMatrix() {return fTriggerTrackNSDEventMatrix;}
106
107	THnSparseF *GetRecTrackEventMatrix() {return fRecTrackEventMatrix;}
108	THnSparseF *GetRecTrackSDEventMatrix() {return fRecTrackSDEventMatrix;}
109	THnSparseF *GetRecTrackDDEventMatrix() {return fRecTrackDDEventMatrix;}
110	THnSparseF *GetRecTrackNDEventMatrix() {return fRecTrackNDEventMatrix;}
111	THnSparseF *GetRecTrackNSDEventMatrix() {return fRecTrackNSDEventMatrix;}
112
113	//
114	THnSparseF *GetGenPrimTrackMatrix() {return fGenPrimTrackMatrix;}
115	THnSparseF *GetRecPrimTrackMatrix() {return fRecPrimTrackMatrix;}
116
117	THnSparseF *GetRecTrackMatrix() {return fRecTrackMatrix;}
118	THnSparseF *GetRecSecTrackMatrix() {return fRecSecTrackMatrix;}
119	THnSparseF *GetRecMultTrackMatrix() {return fRecMultTrackMatrix;}
120
121	//
122	// control histograms
123	//
124	THnSparseF *GetMCEventHist1() {return fMCEventHist1;}
125	THnSparseF *GetRecEventHist1() {return fRecEventHist1;}
126	THnSparseF *GetRecEventHist2() {return fRecEventHist2;}
127	THnSparseF *GetRecMCEventHist1() {return fRecMCEventHist1;}
128	THnSparseF *GetRecMCEventHist2() {return fRecMCEventHist2;}
129	THnSparseF *GetRecMCEventHist3() {return fRecMCEventHist3;}
130
131	THnSparseF *GetMCTrackHist1(Int_t i) {return fMCTrackHist1[i];}
132	THnSparseF *GetMCPrimTrackHist1(Int_t i) {return fMCPrimTrackHist1[i];}
133	THnSparseF *GetMCSecTrackHist1(Int_t i) {return fMCSecTrackHist1[i];}
134
135	THnSparseF *GetRecTrackHist1(Int_t i) {return fRecTrackHist1[i];}
136	THnSparseF *GetRecTrackMultHist1(Int_t i) {return fRecTrackMultHist1[i];}
137
138	THnSparseF *GetRecMCTrackHist1() {return fRecMCTrackHist1;}
139	THnSparseF *GetMCMultRecTrackHist1() {return fMCMultRecTrackHist1;}
140
141	private:
142
143	// analysis folder
144	TFolder *fAnalysisFolder; // folder for analysed histograms
145
146	// switch on/off filling of control histograms
147	Bool_t fHistogramsOn;
148
149	//
150	// correlation matrices (histograms)
151	//
152
153	// event rec. track vs true track multiplicity correlation matrix
154	THnSparseF *fEventMultCorrelationMatrix; //-> mult:mult_true_tracks
155
156	// rec. track pt vs true track pt correlation matrix for given eta
157	THnSparseF *fTrackPtCorrelationMatrix; //-> Pt:mcPt:mcEta
158
159	//
160	// event level correction
161	//
162
163	// all genertated
bad4ba69	164	THnSparseF *fGenEventMatrix; //-> mcZv:multMB (inelastic)
	165	THnSparseF *fGenSDEventMatrix; //-> mcZv:multMB (single diffractive)
	166	THnSparseF *fGenDDEventMatrix; //-> mcZv:multMB (single diffractive)
	167	THnSparseF *fGenNDEventMatrix; //-> mcZv:multMB (non diffractive)
	168	THnSparseF *fGenNSDEventMatrix; //-> mcZv:multMB (non single diffractive)
0aaa8b91	169
0aaa8b91	170	// trigger bias corrections (fTriggerEventMatrix / fGenEventMatrix)
bad4ba69	171	THnSparseF *fTriggerEventMatrix; //-> mcZv:multMB
	172	THnSparseF *fTriggerSDEventMatrix; //-> mcZv:multMB
	173	THnSparseF *fTriggerDDEventMatrix; //-> mcZv:multMB
	174	THnSparseF *fTriggerNDEventMatrix; //-> mcZv:multMB
	175	THnSparseF *fTriggerNSDEventMatrix; //-> mcZv:multMB
0aaa8b91	176
0aaa8b91	177	// event vertex rec. eff correction (fRecEventMatrix / fTriggerEventMatrix)
bad4ba69	178	THnSparseF *fRecEventMatrix; //-> mcZv:multMB
	179	THnSparseF *fRecSDEventMatrix; //-> mcZv:multMB
	180	THnSparseF *fRecDDEventMatrix; //-> mcZv:multMB
	181	THnSparseF *fRecNDEventMatrix; //-> mcZv:multMB
	182	THnSparseF *fRecNSDEventMatrix; //-> mcZv:multMB
0aaa8b91	183
	184	//
	185	// track-event level correction
	186	//
	187
	188	THnSparseF *fGenTrackEventMatrix; //-> mcZv:mcPt:mcEta
	189	THnSparseF *fGenTrackSDEventMatrix; //-> mcZv:mcPt:mcEta
	190	THnSparseF *fGenTrackDDEventMatrix; //-> mcZv:mcPt:mcEta
	191	THnSparseF *fGenTrackNDEventMatrix; //-> mcZv:mcPt:mcEta
	192	THnSparseF *fGenTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta
	193
	194	// trigger bias corrections (fTriggerTrackEventMatrix / fGenTrackEventMatrix)
	195	THnSparseF *fTriggerTrackEventMatrix; //-> mcZv:mcPt:mcEta
	196	THnSparseF *fTriggerTrackSDEventMatrix; //-> mcZv:mcPt:mcEta
	197	THnSparseF *fTriggerTrackDDEventMatrix; //-> mcZv:mcPt:mcEta
	198	THnSparseF *fTriggerTrackNDEventMatrix; //-> mcZv:mcPt:mcEta
	199	THnSparseF *fTriggerTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta
	200
	201	// event vertex rec. corrections (fRecTrackEventMatrix / fTriggerTrackEventMatrix)
	202	THnSparseF *fRecTrackEventMatrix; //-> mcZv:Pt:mcEta
	203	THnSparseF *fRecTrackSDEventMatrix; //-> mcZv:Pt:mcEta
	204	THnSparseF *fRecTrackDDEventMatrix; //-> mcZv:Pt:mcEta
	205	THnSparseF *fRecTrackNDEventMatrix; //-> mcZv:Pt:mcEta
	206	THnSparseF *fRecTrackNSDEventMatrix; //-> mcZv:Pt:mcEta
	207
	208	//
	209	// track level correction
	210	//
	211
	212	// track rec. efficiency correction (fRecPrimTrackMatrix / fGenPrimTrackMatrix)
	213	THnSparseF *fGenPrimTrackMatrix; //-> mcZv:mcPt:mcEta
	214	THnSparseF *fRecPrimTrackMatrix; //-> mcZv:mcPt:mcEta
	215
	216	// secondary track contamination correction (fRecSecTrackMatrix / fRecTrackMatrix)
	217	THnSparseF *fRecTrackMatrix; //-> mcZv:mcPt:mcEta
	218	THnSparseF *fRecSecTrackMatrix; //-> mcZv:mcPt:mcEta
	219
	220	// multiple rec. track corrections (fRecMultTrackMatrix / fRecTrackMatrix)
	221	THnSparseF *fRecMultTrackMatrix; //-> mcZv:Pt:mcEta
	222
	223	//
	224	// ESD and MC control analysis histograms
	225	//
	226
	227	// THnSparse event histograms
	228	THnSparseF *fMCEventHist1; //-> mcXv:mcYv:mcZv
	229	THnSparseF *fRecEventHist1; //-> Xv:Yv:Zv
bad4ba69	230	THnSparseF *fRecEventHist2; //-> Zv:multMB:mult
0aaa8b91	231	THnSparseF *fRecMCEventHist1; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv
bad4ba69	232	THnSparseF *fRecMCEventHist2; //-> Xv-mcXv:Zv-mcZv:mult
bad4ba69	233	THnSparseF *fRecMCEventHist3; //-> mult:EventType (ND, DD, SD)
0aaa8b91	234
	235	// THnSparse track histograms
	236	// [0] - after charged track selection, [1] - after acceptance cuts, [2] - after esd track cuts
	237
	238	THnSparseF *fMCTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:mcPhi
	239	THnSparseF *fMCPrimTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:pid:mech:mother
	240	THnSparseF *fMCSecTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:pid:mech:mother
	241
	242	THnSparseF *fRecTrackHist1[AlidNdPtHelper::kCutSteps]; //-> Pt:Eta:Phi
bad4ba69	243	THnSparseF *fRecTrackMultHist1[AlidNdPtHelper::kCutSteps]; //-> Pt:mult
0aaa8b91	244
	245	THnSparseF *fRecMCTrackHist1; //-> mcPt:mcEta:(Pt-mcPt)/mcPt:(Eta-mcEta)
	246
	247	//multple reconstructed tracks
	248	THnSparseF *fMCMultRecTrackHist1; //-> mcPt:mcEta:pid
	249
	250	AlidNdPtAnalysis(const AlidNdPtAnalysis&); // not implemented
	251	AlidNdPtAnalysis& operator=(const AlidNdPtAnalysis&); // not implemented
	252
	253	ClassDef(AlidNdPtAnalysis,1);
	254	};
	255
	256	#endif