[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / ChargedHadrons / dNdPt / AlidNdPtAnalysis.h

#ifndef ALIDNDPTANALYSIS_H
#define ALIDNDPTANALYSIS_H

//------------------------------------------------------------------------------
// AlidNdPtAnalysis class used for dNdPt analysis. 
// 
// Author: J.Otwinowski 04/11/2008 
// last change: 2013-06-13 by M.Knichel
//------------------------------------------------------------------------------

class iostream;

class TFile;
class TCint;
class TProfile;
class TFolder;
class TObjArray;
class TString;
class THnSparse;

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

#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* const list);

  // Analyse output histograms 
  virtual void Analyse();

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

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

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

  // Define 0-multiplicity bin for LHC
  // background calculation
  static Bool_t IsBinZeroSPDvtx(const AliESDEvent* esdEvent);
  static Bool_t IsBinZeroTrackSPDvtx(const AliESDEvent* esdEvent);
    
  // Create folder for analysed histograms
  TFolder *CreateFolder(TString folder = "folderdNdPtAnalysis",TString title = "Analysed dNdPt histograms");

  // Set binning for Histograms (if not set default binning is used)
  void SetBinsMult(Int_t nbins, Double_t* edges) { fMultNbins = nbins; fBinsMult = CloneArray(nbins+1,edges); }
  void SetBinsPt(Int_t nbins, Double_t* edges) { fPtNbins = nbins; fBinsPt = CloneArray(nbins+1,edges); }
  void SetBinsPtCorr(Int_t nbins, Double_t* edges) { fPtCorrNbins = nbins; fBinsPtCorr = CloneArray(nbins+1,edges); }
  void SetBinsEta(Int_t nbins, Double_t* edges) { fEtaNbins = nbins; fBinsEta = CloneArray(nbins+1,edges); }
  void SetBinsZv(Int_t nbins, Double_t* edges) { fZvNbins = nbins; fBinsZv = CloneArray(nbins+1,edges); }


  // Fill histograms
  void FillHistograms(AliESDtrack *const esdTrack, AliStack *const stack, const Double_t zv, AlidNdPtHelper::TrackObject trackObj, Int_t multMB);
  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() const {return fEventMultCorrelationMatrix;}
  THnSparseF *GetTrackPtCorrelationMatrix()   const {return fTrackPtCorrelationMatrix;}

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

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

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

  //
  THnSparseF *GetRecCandleEventMatrix() const {return fRecCandleEventMatrix;}

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

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

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

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

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

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

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

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

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

  THnSparseF *GetRecTrackHist2() const {return fRecTrackHist2;}
  THnSparseF *GetEventCount() const {return fEventCount;}  

  //
  // Generic histograms to be corrected
  //
  THnSparseF *GetRecEventHist() const {return fRecEventHist;} 
  THnSparseF *GetRecTrackHist() const {return fRecTrackHist;} 
  THnSparseF* GetMCPrimTrackHist() const {return fMCPrimTrackHist;}   

private:

  // analysis folder 
  TFolder *fAnalysisFolder; // folder for analysed histograms
  Bool_t fHistogramsOn; // switch on/off filling of control histograms 

  // 
  // correlation matrices (histograms)
  //

  // event rec. track vs true track multiplicity correlation matrix 
  THnSparseF *fEventMultCorrelationMatrix; //-> multRecMult:multTrueMC:multMB

  // 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:multTrueMC (inelastic)
  THnSparseF *fGenSDEventMatrix; //-> mcZv:multTrueMC (single diffractive)
  THnSparseF *fGenDDEventMatrix; //-> mcZv:multTrueMC (single diffractive)
  THnSparseF *fGenNDEventMatrix; //-> mcZv:multTrueMC (non diffractive)
  THnSparseF *fGenNSDEventMatrix; //-> mcZv:multTrueMC (non single diffractive)

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

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


  //
  // track-event level correction 
  //

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

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

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

  //
  // track level correction 
  //

  // track rec. efficiency correction (fRecPrimTrackMatrix / fGenPrimTrackMatrix)
  THnSparseF *fGenTrackMatrix; //-> mcZv:mcPt:mcEta
  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:multRecMult
  THnSparseF *fRecMCEventHist1; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv
  THnSparseF *fRecMCEventHist2; //-> Xv-mcXv:Zv-mcZv:multMB
  THnSparseF *fRecMCEventHist3; //-> multRecMult: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 *fMCPrimTrackHist2[AlidNdPtHelper::kCutSteps]; //-> pdg: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

  // track control histograms
  THnSparseF *fRecTrackHist2;  //-> nclust:chi2:Pt:Eta:Phi

  //
  // Generic histograms to be corrected
  //
  THnSparseF *fRecEventHist;    //-> Zv:multMB  
  THnSparseF *fRecTrackHist; //-> Zv:pT:eta:multRecMult
  THnSparseF *fEventCount; //-> trig, trig + vertex, selected event
  THnSparseF *fMCPrimTrackHist; //-> Zv:mcpT:mceta:multTrueMC  

  //
  // candle events track corrections
  //
  THnSparseF *fRecCandleEventMatrix; //-> Zv:multTrueMC

  //binning for THNsparse
  Int_t	fMultNbins;
  Int_t fPtNbins;
  Int_t fPtCorrNbins;
  Int_t fEtaNbins;
  Int_t fZvNbins;
  Double_t* fBinsMult;
  Double_t* fBinsPt;
  Double_t* fBinsPtCorr;
  Double_t* fBinsEta;
  Double_t* fBinsZv;
  
  Bool_t fIsInit;

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

  ClassDef(AlidNdPtAnalysis,7);
};

#endif
Commit	Line	Data
00b81913	1	#ifndef ALIDNDPTANALYSIS_H
00b81913	2	#define ALIDNDPTANALYSIS_H
0aaa8b91	3
0aaa8b91	4	//------------------------------------------------------------------------------
00b81913	5	// AlidNdPtAnalysis class used for dNdPt analysis.
0aaa8b91	6	//
0aaa8b91	7	// Author: J.Otwinowski 04/11/2008
df1c0513	8	// last change: 2013-06-13 by M.Knichel
0aaa8b91	9	//------------------------------------------------------------------------------
0aaa8b91	10
00b81913	11	class iostream;
	12
	13	class TFile;
	14	class TCint;
0aaa8b91	15	class TProfile;
	16	class TFolder;
	17	class TObjArray;
	18	class TString;
00b81913	19	class THnSparse;
0aaa8b91	20
	21	class AliESDtrackCuts;
	22	class AliVertexerTracks;
	23	class AliESD;
	24	class AliESDfriend;
	25	class AliESDfriendTrack;
bad4ba69	26	class AlidNdPtHelper;
bad4ba69	27
0aaa8b91	28	#include "AlidNdPt.h"
0aaa8b91	29
	30	class AlidNdPtAnalysis : public AlidNdPt {
	31	public :
	32	AlidNdPtAnalysis();
	33	AlidNdPtAnalysis(Char_t* name, Char_t* title);
	34	~AlidNdPtAnalysis();
	35
	36	// Init data members
	37	virtual void Init();
	38
	39	// Process events
	40	virtual void Process(AliESDEvent const esdEvent=0, AliMCEvent const mcEvent=0);
	41
	42	// Merge output objects (needed by PROOF)
00b81913	43	virtual Long64_t Merge(TCollection* const list);
0aaa8b91	44
	45	// Analyse output histograms
	46	virtual void Analyse();
	47
	48	// Export objects to folder
c931948b	49	virtual TFolder ExportToFolder(TObjArray const array=0);
0aaa8b91	50
0aaa8b91	51	// Get analysis folder
00b81913	52	TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
0aaa8b91	53
0aaa8b91	54	// Fill control histograms
00b81913	55	void SetHistogramsOn(const Bool_t histOn=kTRUE) {fHistogramsOn = histOn;}
00b81913	56	Bool_t IsHistogramsOn() const {return fHistogramsOn;}
0aaa8b91	57
791aaf54	58	// Define 0-multiplicity bin for LHC
	59	// background calculation
	60	static Bool_t IsBinZeroSPDvtx(const AliESDEvent* esdEvent);
	61	static Bool_t IsBinZeroTrackSPDvtx(const AliESDEvent* esdEvent);
	62
0aaa8b91	63	// Create folder for analysed histograms
	64	TFolder *CreateFolder(TString folder = "folderdNdPtAnalysis",TString title = "Analysed dNdPt histograms");
	65
68f10917	66	// Set binning for Histograms (if not set default binning is used)
fc98fbb5	67	void SetBinsMult(Int_t nbins, Double_t* edges) { fMultNbins = nbins; fBinsMult = CloneArray(nbins+1,edges); }
	68	void SetBinsPt(Int_t nbins, Double_t* edges) { fPtNbins = nbins; fBinsPt = CloneArray(nbins+1,edges); }
	69	void SetBinsPtCorr(Int_t nbins, Double_t* edges) { fPtCorrNbins = nbins; fBinsPtCorr = CloneArray(nbins+1,edges); }
	70	void SetBinsEta(Int_t nbins, Double_t* edges) { fEtaNbins = nbins; fBinsEta = CloneArray(nbins+1,edges); }
	71	void SetBinsZv(Int_t nbins, Double_t* edges) { fZvNbins = nbins; fBinsZv = CloneArray(nbins+1,edges); }
68f10917	72
68f10917	73
0aaa8b91	74	// Fill histograms
fff1caf5	75	void FillHistograms(AliESDtrack const esdTrack, AliStack const stack, const Double_t zv, AlidNdPtHelper::TrackObject trackObj, Int_t multMB);
0aaa8b91	76	void FillHistograms(AliStack *const stack, Int_t label, AlidNdPtHelper::TrackObject trackObj);
	77	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);
	78
	79	// Getters
00b81913	80	THnSparseF *GetEventMultCorrelationMatrix() const {return fEventMultCorrelationMatrix;}
00b81913	81	THnSparseF *GetTrackPtCorrelationMatrix() const {return fTrackPtCorrelationMatrix;}
0aaa8b91	82
0aaa8b91	83	//
00b81913	84	THnSparseF *GetGenEventMatrix() const {return fGenEventMatrix;}
	85	THnSparseF *GetGenSDEventMatrix() const {return fGenSDEventMatrix;}
	86	THnSparseF *GetGenDDEventMatrix() const {return fGenDDEventMatrix;}
	87	THnSparseF *GetGenNDEventMatrix() const {return fGenNDEventMatrix;}
	88	THnSparseF *GetGenNSDEventMatrix() const {return fGenNSDEventMatrix;}
	89
	90	THnSparseF *GetTriggerEventMatrix() const {return fTriggerEventMatrix;}
	91	THnSparseF *GetTriggerSDEventMatrix() const {return fTriggerSDEventMatrix;}
	92	THnSparseF *GetTriggerDDEventMatrix() const {return fTriggerDDEventMatrix;}
	93	THnSparseF *GetTriggerNDEventMatrix() const {return fTriggerNDEventMatrix;}
	94	THnSparseF *GetTriggerNSDEventMatrix() const {return fTriggerNSDEventMatrix;}
	95
	96	THnSparseF *GetRecEventMatrix() const {return fRecEventMatrix;}
	97	THnSparseF *GetRecSDEventMatrix() const {return fRecSDEventMatrix;}
	98	THnSparseF *GetRecDDEventMatrix() const {return fRecDDEventMatrix;}
	99	THnSparseF *GetRecNDEventMatrix() const {return fRecNDEventMatrix;}
	100	THnSparseF *GetRecNSDEventMatrix() const {return fRecNSDEventMatrix;}
0aaa8b91	101
55468faf	102	//
	103	THnSparseF *GetRecCandleEventMatrix() const {return fRecCandleEventMatrix;}
	104
0aaa8b91	105	//
00b81913	106	THnSparseF *GetGenTrackEventMatrix() const {return fGenTrackEventMatrix;}
	107	THnSparseF *GetGenTrackSDEventMatrix() const {return fGenTrackSDEventMatrix;}
	108	THnSparseF *GetGenTrackDDEventMatrix() const {return fGenTrackDDEventMatrix;}
	109	THnSparseF *GetGenTrackNDEventMatrix() const {return fGenTrackNDEventMatrix;}
	110	THnSparseF *GetGenTrackNSDEventMatrix() const {return fGenTrackNSDEventMatrix;}
	111
	112	THnSparseF *GetTriggerTrackEventMatrix() const {return fTriggerTrackEventMatrix;}
	113	THnSparseF *GetTriggerTrackSDEventMatrix() const {return fTriggerTrackSDEventMatrix;}
	114	THnSparseF *GetTriggerTrackDDEventMatrix() const {return fTriggerTrackDDEventMatrix;}
	115	THnSparseF *GetTriggerTrackNDEventMatrix() const {return fTriggerTrackNDEventMatrix;}
	116	THnSparseF *GetTriggerTrackNSDEventMatrix() const {return fTriggerTrackNSDEventMatrix;}
	117
	118	THnSparseF *GetRecTrackEventMatrix() const {return fRecTrackEventMatrix;}
	119	THnSparseF *GetRecTrackSDEventMatrix() const {return fRecTrackSDEventMatrix;}
	120	THnSparseF *GetRecTrackDDEventMatrix() const {return fRecTrackDDEventMatrix;}
	121	THnSparseF *GetRecTrackNDEventMatrix() const {return fRecTrackNDEventMatrix;}
	122	THnSparseF *GetRecTrackNSDEventMatrix() const {return fRecTrackNSDEventMatrix;}
0aaa8b91	123
0aaa8b91	124	//
00b81913	125	THnSparseF *GetGenTrackMatrix() const {return fGenTrackMatrix;}
	126	THnSparseF *GetGenPrimTrackMatrix() const {return fGenPrimTrackMatrix;}
	127	THnSparseF *GetRecPrimTrackMatrix() const {return fRecPrimTrackMatrix;}
0aaa8b91	128
00b81913	129	THnSparseF *GetRecTrackMatrix() const {return fRecTrackMatrix;}
	130	THnSparseF *GetRecSecTrackMatrix() const {return fRecSecTrackMatrix;}
	131	THnSparseF *GetRecMultTrackMatrix() const {return fRecMultTrackMatrix;}
0aaa8b91	132
	133	//
	134	// control histograms
	135	//
00b81913	136	THnSparseF *GetMCEventHist1() const {return fMCEventHist1;}
	137	THnSparseF *GetRecEventHist1() const {return fRecEventHist1;}
	138	THnSparseF *GetRecEventHist2() const {return fRecEventHist2;}
	139	THnSparseF *GetRecMCEventHist1() const {return fRecMCEventHist1;}
	140	THnSparseF *GetRecMCEventHist2() const {return fRecMCEventHist2;}
	141	THnSparseF *GetRecMCEventHist3() const {return fRecMCEventHist3;}
0aaa8b91	142
00b81913	143	THnSparseF *GetMCTrackHist1(Int_t i) const {return fMCTrackHist1[i];}
00b81913	144	THnSparseF *GetMCPrimTrackHist1(Int_t i) const {return fMCPrimTrackHist1[i];}
7f7dd416	145	THnSparseF *GetMCPrimTrackHist2(Int_t i) const {return fMCPrimTrackHist2[i];}
00b81913	146	THnSparseF *GetMCSecTrackHist1(Int_t i) const {return fMCSecTrackHist1[i];}
0aaa8b91	147
00b81913	148	THnSparseF *GetRecTrackHist1(Int_t i) const {return fRecTrackHist1[i];}
00b81913	149	THnSparseF *GetRecTrackMultHist1(Int_t i) const {return fRecTrackMultHist1[i];}
0aaa8b91	150
00b81913	151	THnSparseF *GetRecMCTrackHist1() const {return fRecMCTrackHist1;}
00b81913	152	THnSparseF *GetMCMultRecTrackHist1() const {return fMCMultRecTrackHist1;}
0aaa8b91	153
00b81913	154	THnSparseF *GetRecTrackHist2() const {return fRecTrackHist2;}
24c88fc4	155	THnSparseF *GetEventCount() const {return fEventCount;}
f4f8cae4	156
	157	//
	158	// Generic histograms to be corrected
	159	//
	160	THnSparseF *GetRecEventHist() const {return fRecEventHist;}
	161	THnSparseF *GetRecTrackHist() const {return fRecTrackHist;}
df1c0513	162	THnSparseF* GetMCPrimTrackHist() const {return fMCPrimTrackHist;}
f4f8cae4	163
0aaa8b91	164	private:
	165
	166	// analysis folder
	167	TFolder *fAnalysisFolder; // folder for analysed histograms
00b81913	168	Bool_t fHistogramsOn; // switch on/off filling of control histograms
0aaa8b91	169
	170	//
	171	// correlation matrices (histograms)
	172	//
	173
	174	// event rec. track vs true track multiplicity correlation matrix
df1c0513	175	THnSparseF *fEventMultCorrelationMatrix; //-> multRecMult:multTrueMC:multMB
0aaa8b91	176
	177	// rec. track pt vs true track pt correlation matrix for given eta
	178	THnSparseF *fTrackPtCorrelationMatrix; //-> Pt:mcPt:mcEta
	179
	180	//
	181	// event level correction
	182	//
	183
	184	// all genertated
df1c0513	185	THnSparseF *fGenEventMatrix; //-> mcZv:multTrueMC (inelastic)
	186	THnSparseF *fGenSDEventMatrix; //-> mcZv:multTrueMC (single diffractive)
	187	THnSparseF *fGenDDEventMatrix; //-> mcZv:multTrueMC (single diffractive)
	188	THnSparseF *fGenNDEventMatrix; //-> mcZv:multTrueMC (non diffractive)
	189	THnSparseF *fGenNSDEventMatrix; //-> mcZv:multTrueMC (non single diffractive)
0aaa8b91	190
0aaa8b91	191	// trigger bias corrections (fTriggerEventMatrix / fGenEventMatrix)
df1c0513	192	THnSparseF *fTriggerEventMatrix; //-> mcZv:multTrueMC
	193	THnSparseF *fTriggerSDEventMatrix; //-> mcZv:multTrueMC
	194	THnSparseF *fTriggerDDEventMatrix; //-> mcZv:multTrueMC
	195	THnSparseF *fTriggerNDEventMatrix; //-> mcZv:multTrueMC
	196	THnSparseF *fTriggerNSDEventMatrix; //-> mcZv:multTrueMC
0aaa8b91	197
0aaa8b91	198	// event vertex rec. eff correction (fRecEventMatrix / fTriggerEventMatrix)
df1c0513	199	THnSparseF *fRecEventMatrix; //-> mcZv:multTrueMC
	200	THnSparseF *fRecSDEventMatrix; //-> mcZv:multTrueMC
	201	THnSparseF *fRecDDEventMatrix; //-> mcZv:multTrueMC
	202	THnSparseF *fRecNDEventMatrix; //-> mcZv:multTrueMC
	203	THnSparseF *fRecNSDEventMatrix; //-> mcZv:multTrueMC
0aaa8b91	204
55468faf	205
0aaa8b91	206	//
	207	// track-event level correction
	208	//
	209
df1c0513	210	THnSparseF *fGenTrackEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	211	THnSparseF *fGenTrackSDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	212	THnSparseF *fGenTrackDDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	213	THnSparseF *fGenTrackNDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	214	THnSparseF *fGenTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
0aaa8b91	215
0aaa8b91	216	// trigger bias corrections (fTriggerTrackEventMatrix / fGenTrackEventMatrix)
df1c0513	217	THnSparseF *fTriggerTrackEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	218	THnSparseF *fTriggerTrackSDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	219	THnSparseF *fTriggerTrackDDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	220	THnSparseF *fTriggerTrackNDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
	221	THnSparseF *fTriggerTrackNSDEventMatrix; //-> mcZv:mcPt:mcEta:multTrueMC
0aaa8b91	222
0aaa8b91	223	// event vertex rec. corrections (fRecTrackEventMatrix / fTriggerTrackEventMatrix)
df1c0513	224	THnSparseF *fRecTrackEventMatrix; //-> mcZv:Pt:mcEta:multTrueMC
	225	THnSparseF *fRecTrackSDEventMatrix; //-> mcZv:Pt:mcEta:multTrueMC
	226	THnSparseF *fRecTrackDDEventMatrix; //-> mcZv:Pt:mcEta:multTrueMC
	227	THnSparseF *fRecTrackNDEventMatrix; //-> mcZv:Pt:mcEta:multTrueMC
	228	THnSparseF *fRecTrackNSDEventMatrix; //-> mcZv:Pt:mcEta:multTrueMC
0aaa8b91	229
	230	//
	231	// track level correction
	232	//
	233
	234	// track rec. efficiency correction (fRecPrimTrackMatrix / fGenPrimTrackMatrix)
985bdfd5	235	THnSparseF *fGenTrackMatrix; //-> mcZv:mcPt:mcEta
0aaa8b91	236	THnSparseF *fGenPrimTrackMatrix; //-> mcZv:mcPt:mcEta
	237	THnSparseF *fRecPrimTrackMatrix; //-> mcZv:mcPt:mcEta
	238
	239	// secondary track contamination correction (fRecSecTrackMatrix / fRecTrackMatrix)
	240	THnSparseF *fRecTrackMatrix; //-> mcZv:mcPt:mcEta
	241	THnSparseF *fRecSecTrackMatrix; //-> mcZv:mcPt:mcEta
	242
	243	// multiple rec. track corrections (fRecMultTrackMatrix / fRecTrackMatrix)
	244	THnSparseF *fRecMultTrackMatrix; //-> mcZv:Pt:mcEta
	245
	246	//
	247	// ESD and MC control analysis histograms
	248	//
	249
	250	// THnSparse event histograms
	251	THnSparseF *fMCEventHist1; //-> mcXv:mcYv:mcZv
	252	THnSparseF *fRecEventHist1; //-> Xv:Yv:Zv
df1c0513	253	THnSparseF *fRecEventHist2; //-> Zv:multMB:multRecMult
0aaa8b91	254	THnSparseF *fRecMCEventHist1; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv
df1c0513	255	THnSparseF *fRecMCEventHist2; //-> Xv-mcXv:Zv-mcZv:multMB
df1c0513	256	THnSparseF *fRecMCEventHist3; //-> multRecMult:EventType (ND, DD, SD)
0aaa8b91	257
	258	// THnSparse track histograms
	259	// [0] - after charged track selection, [1] - after acceptance cuts, [2] - after esd track cuts
	260
	261	THnSparseF *fMCTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:mcPhi
	262	THnSparseF *fMCPrimTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:pid:mech:mother
7f7dd416	263	THnSparseF *fMCPrimTrackHist2[AlidNdPtHelper::kCutSteps]; //-> pdg:mech:mother
0aaa8b91	264	THnSparseF *fMCSecTrackHist1[AlidNdPtHelper::kCutSteps]; //-> mcPt:mcEta:pid:mech:mother
	265
	266	THnSparseF *fRecTrackHist1[AlidNdPtHelper::kCutSteps]; //-> Pt:Eta:Phi
bad4ba69	267	THnSparseF *fRecTrackMultHist1[AlidNdPtHelper::kCutSteps]; //-> Pt:mult
0aaa8b91	268
	269	THnSparseF *fRecMCTrackHist1; //-> mcPt:mcEta:(Pt-mcPt)/mcPt:(Eta-mcEta)
	270
	271	//multple reconstructed tracks
	272	THnSparseF *fMCMultRecTrackHist1; //-> mcPt:mcEta:pid
	273
985bdfd5	274	// track control histograms
	275	THnSparseF *fRecTrackHist2; //-> nclust:chi2:Pt:Eta:Phi
	276
f4f8cae4	277	//
	278	// Generic histograms to be corrected
	279	//
24c88fc4	280	THnSparseF *fRecEventHist; //-> Zv:multMB
df1c0513	281	THnSparseF *fRecTrackHist; //-> Zv:pT:eta:multRecMult
24c88fc4	282	THnSparseF *fEventCount; //-> trig, trig + vertex, selected event
df1c0513	283	THnSparseF *fMCPrimTrackHist; //-> Zv:mcpT:mceta:multTrueMC
f4f8cae4	284
55468faf	285	//
	286	// candle events track corrections
	287	//
df1c0513	288	THnSparseF *fRecCandleEventMatrix; //-> Zv:multTrueMC
55468faf	289
68f10917	290	//binning for THNsparse
	291	Int_t fMultNbins;
	292	Int_t fPtNbins;
	293	Int_t fPtCorrNbins;
	294	Int_t fEtaNbins;
	295	Int_t fZvNbins;
	296	Double_t* fBinsMult;
	297	Double_t* fBinsPt;
	298	Double_t* fBinsPtCorr;
	299	Double_t* fBinsEta;
	300	Double_t* fBinsZv;
	301
	302	Bool_t fIsInit;
	303
0aaa8b91	304	AlidNdPtAnalysis(const AlidNdPtAnalysis&); // not implemented
	305	AlidNdPtAnalysis& operator=(const AlidNdPtAnalysis&); // not implemented
	306
df1c0513	307	ClassDef(AlidNdPtAnalysis,7);
0aaa8b91	308	};
	309
	310	#endif