[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / Kinks / AliAnalysisPionKinksMCESD.h

#ifndef AliAnalysisPionKinksMCESD_h
#define AliAnalysisPionKinksMCESD_h


/*  See cxx source for full Copyright notice */

//-----------------------------------------------------------------
//                 AliAnalysisPionKinksMCESD class
//         This task is an example of an analysis task
//                  for kink topology Study
//          Authors: Eftychios Cheiladakis, under the supervision
//           of Martha's Spyropoulou-Stassinaki   
//          Physics Department of Athens University
//------------------------------------------------------------------
class AliESDEvent;
class TF1;
class TH1;
class TH2;
class TH3;
class TParticle;
class TParticle;
class AliESDtrackCuts;
class AliPhysicsSelection;
class AliPIDResponse;

#include "AliAnalysisTaskSE.h"

class AliAnalysisPionKinksMCESD : public AliAnalysisTaskSE {
	public:
		AliAnalysisPionKinksMCESD(const char *name = "AliAnalysisPionKinksMCESD");
		virtual ~AliAnalysisPionKinksMCESD() {}
  
		virtual void   UserCreateOutputObjects();
		virtual void   UserExec(Option_t *option);
		virtual void   Terminate(Option_t *);

  		const AliESDVertex *GetEventVertex(AliESDEvent* esd);
		Double_t Energy(AliESDtrack* track) const;
		Double_t GetMCRapidity(TParticle* particle) const;
		Double_t fuRapidity(AliESDtrack* track) const;
		Double_t MCPQt(AliMCEvent* mcEvent, Int_t iMC, TParticle* MCdaughter) const;
		Double_t fuMCKinkAngle(AliMCEvent* mcEvent, Int_t iMC, TParticle* MCdaughter, Bool_t degrees) const;
		Bool_t IsGoodTrack(AliESDtrack* ESDTrack) const;
		Bool_t IsPrimaryTrack(AliESDtrack* ESDTrack) const;
		Bool_t IsRealPionKink(AliESDkink* kink, AliStack* MCstack) const;

		// Set limits & cuts
		void SetMulCut(Int_t low, Int_t up){fLowMulcut=low; fUpMulcut=up;}	  
		void SetPtCut(Double_t PtCut){cLowPt=PtCut;}
		void SetRapidityLimits(Double_t RapidityLim){cRapidityLim=RapidityLim;}
		void SetRadiusRange(Double_t LowR, Double_t UpR){cLowR=LowR; cUpR=UpR;}
		void SetZRange(Double_t LowZ, Double_t UpZ){cLowZ=LowZ; cUpZ=UpZ;}		
		void SetKinkAngleCut(Double_t LowAngle){cLowKinkAngle=LowAngle;}
		void SetQtRange(Double_t LowQt, Double_t UpQt){cLowQt=LowQt; cUpQt=UpQt;}
		void SetInvMassRange(Double_t LowInvMass, Double_t UpInvMass){cLowInvMass=LowInvMass; cUpInvMass=UpInvMass;}
		void SetSigmaCut(Double_t SigmaCut){cSigmaCut=SigmaCut;}
		void SetPdgCodes(Int_t kaon, Int_t pion, Int_t muon, Int_t electron){cPdgKaon=kaon; cPdgPion=pion; cPdgMuon=muon; cPdgElectron=electron;}
		void SetMasses(Double_t KaonMass, Double_t PionMass, Double_t MuonMass, Double_t ElectronMass){cKaonMass=KaonMass; cPionMass=PionMass; cMuonMass=MuonMass; cElectronMass=ElectronMass;}

		// Set histos limits
		void SetMultHistos(Int_t BinsMult, Int_t LowMult, Int_t UpMult){nBinsMult=BinsMult; hLowPt=LowMult; hUpMult=UpMult;}
		void SetPtHistos(Int_t BinsPt, Double_t LowPt, Double_t UpPt){nBinsPt=BinsPt; hLowPt=LowPt; hUpPt=UpPt;}   
		void SetEtaHistos(Int_t BinsEta, Double_t LowEta, Double_t UpEta){nBinsEta=BinsEta; hLowEta=LowEta; hUpEta=UpEta;}
		void SetQtHistos(Int_t BinsQt, Double_t LowQt, Double_t UpQt){nBinsQt=BinsQt; hLowQt=LowQt; hUpQt=UpQt;}
		void SetPdgHistos(Int_t BinsPdg, Double_t LowPdg, Double_t UpPdg){nBinsPdg=BinsPdg; hLowPdg=LowPdg; hUpPdg=UpPdg;}
		void SetPdg2Histos(Int_t BinsPdg2, Double_t LowPdg2, Double_t UpPdg2){nBinsPdg2=BinsPdg2; hLowPdg2=LowPdg2; hUpPdg2=UpPdg2;}
		void SetUIDHistos(Int_t BinsUID, Double_t LowUID, Double_t UpUID){nBinsUID=BinsUID; hLowUID=LowUID; hUpUID=UpUID;}
		void SetRHistos(Int_t BinsR, Double_t LowR, Double_t UpR){nBinsR=BinsR; hLowR=LowR; hUpR=UpR;}
		void SetZHistos(Int_t BinsZ, Double_t LowZ, Double_t UpZ){nBinsZ=BinsZ; hLowZ=LowZ; hUpZ=UpZ;}
		void SetXYHistos(Int_t BinsXY, Double_t LowXY, Double_t UpXY){nBinsXY=BinsXY; hLowXY=LowXY; hUpXY=UpXY;}
		void SetAngleHistos(Int_t BinsAngle, Double_t LowAngle, Double_t UpAngle){nBinsAngle=BinsAngle; hLowAngle=LowAngle; hUpAngle=UpAngle;}
		void SetZVHistos(Int_t BinsZV, Double_t LowZV, Double_t UpZV){nBinsZV=BinsZV; hLowZV=LowZV; hUpZV=UpZV;}
		void SetXYVHistos(Int_t BinsXYV, Double_t LowXYV, Double_t UpXYV){nBinsXYV=BinsXYV; hLowXYV=LowXYV; hUpXYV=UpXYV;}
		void SetInvMassHistos(Int_t BinsInvMass, Double_t LowInvMass, Double_t UpInvMass){nBinsInvMass=BinsInvMass; hLowInvMass=LowInvMass; hUpInvMass=UpInvMass;}
		void SetdEdxHistos(Int_t BinsdEdx, Double_t LowdEdx, Double_t UpdEdx){nBinsdEdx=BinsdEdx; hLowdEdx=LowdEdx; hUpdEdx=UpdEdx;}
	private:
		TF1* fMaxKinkAngKmu; 
		TF1* fMaxKinkAngPimu;

		TH1F* hMCMult;
		TH1F* hMCMultPrim;
		TH1F* hMCPtAll;
		TH1F* hMCEtaAll;
		TH1F* hMCPtPrim; 
		TH1F* hMCEtaPrim; 
		TH1F* hMCPt; 
		TH1F* hMCEta;  
		TH1F* hMCPdg; 
		TH1F* hMCMultPiPlus; 
		TH1F* hMCPtPiPlus; 
		TH1F* hMCEtaPiPlus; 
		TH1F* hMCRapidityPiPlus;
		TH1F* hMCNDaughtersPlus;
		TH1F* hMCRadPiDauPlus; 
		TH1F* hMCKinkPosZPlus;
		TH1F* hMCUIDPiDauPlus; 
		TH1F* hMCPdgPiNonDecayedPlus; 
		TH1F* hMCPdgPiDauPlus;
		TH1F* hMCQtPlus; 
		TH1F* hMCKinkAnglePlus;
		TH2F* hMCPKinkAngPlus;
		TH2F* hMCPdgCodemdPlus; 
		TH2F* hMCPtmdPlus; 
		TH1F* hMCPtPimuonPlus;
		TH1F* hMCEtaPimuonPlus;
		TH1F* hMCRapidityPimuonPlus;
		TH1F* hMCQtPimuonPlus; 
		TH1F* hMCPKinkAngPimuonPlus;
		TH1F* hMCPtPiotherPlus;
		TH1F* hMCEtaPiotherPlus;
		TH1F* hMCRapidityPiotherPlus;
		TH1F* hMCQtPiotherPlus; 
		TH1F* hMCPKinkAngPiotherPlus;
		TH1F* hMCMultPiMinus; 
		TH1F* hMCPtPiMinus; 
		TH1F* hMCEtaPiMinus; 
		TH1F* hMCRapidityPiMinus;
		TH1F* hMCNDaughtersMinus;
		TH1F* hMCRadPiDauMinus; 
		TH1F* hMCKinkPosZMinus;
		TH1F* hMCUIDPiDauMinus; 
		TH1F* hMCPdgPiNonDecayedMinus; 
		TH1F* hMCPdgPiDauMinus;
		TH1F* hMCQtMinus; 
		TH1F* hMCKinkAngleMinus;
		TH2F* hMCPKinkAngMinus;
		TH2F* hMCPdgCodemdMinus; 
		TH2F* hMCPtmdMinus; 
		TH1F* hMCPtPimuonMinus;
		TH1F* hMCEtaPimuonMinus;
		TH1F* hMCRapidityPimuonMinus;
		TH1F* hMCQtPimuonMinus; 
		TH1F* hMCPKinkAngPimuonMinus;
		TH1F* hMCPtPiotherMinus;
		TH1F* hMCEtaPiotherMinus;
		TH1F* hMCRapidityPiotherMinus;
		TH1F* hMCQtPiotherMinus; 
		TH1F* hMCPKinkAngPiotherMinus;

		TH1F* hMult;
		TH1F* hAcceptedMult;
		TH1F* hMultPS;
		TH3F* hvtx;
		TH2F* hvtxy;
		TH2F* hvtyz;
		TH2F* hvtxz;
		TH1F* hMultPSV;
		TH1F* hPtAll;
		TH1F* hEtaAll;
		TH3F* hTrackPos;
		TH2F* hTrackPosxy;
		TH2F* hTrackPosyz;
		TH2F* hTrackPosxz;
		//TH1F* hTPCchi2clusters;
		//TH1F* hdcaToVertexXY;
		//TH1F* hdcaToVertexZ;
		TH1F* hMultPrim;
		TH1F* hPtPrim;
		TH1F* hEtaPrim;
		TH3F* hPrimTrackPos;
		TH2F* hPrimTrackPosxy;
		TH2F* hPrimTrackPosyz;
		TH2F* hPrimTrackPosxz;
		TH1F* hPt;
		TH1F* hEta;
		//TH1F* hRapidity;
		TH1F* hPtKink;
		TH1F* hEtaKink;
		TH1F* hRapidityKink;
		TH2F* hPmP;
		TH1F* hKinkPosRTPCclusters1;
		TH2F* hKinkPosRTPCclusters2;
		TH1F* hQt;
		TH1F* hKinkAngle;
		TH1F* hDCAkink;
		TH2F* hPmKinkAng;
		TH2F* hKinkPosXY;
		TH2F* hKinkPosZY;
		TH2F* hKinkPosZR;
		TH1F* hKinkPosR;
		TH1F* hKinkPosZ;
		TH2F* hKinkPosZMCKinkPosZ;
		TH2F* hPdgCodemd;
		TH2F* hPmd;
		TH1F* hMinvPimu;
		TH1F* hUIDKinkDau;
		TH2F* hdEdx;
		TH1F* hPtKinkFake;
		TH1F* hEtaKinkFake;
		TH1F* hRapidityKinkFake;
		TH2F* hPmPFake;
		TH1F* hKinkPosRTPCclusters1Fake;
		TH2F* hKinkPosRTPCclusters2Fake;
		TH1F* hQtFake;
		TH1F* hKinkAngleFake;
		TH1F* hDCAkinkFake;
		TH2F* hPmKinkAngFake;
		TH2F* hKinkPosXYFake;
		TH2F* hKinkPosZYFake;
		TH2F* hKinkPosZRFake;
		TH1F* hKinkPosRFake;
		TH1F* hKinkPosZFake;
		TH2F* hKinkPosZMCKinkPosZFake;
		TH2F* hPdgCodemdFake;
		TH2F* hPmdFake;
		TH1F* hMinvPimuFake;
		TH1F* hUIDKinkDauFake;
		TH2F* hdEdxFake;
		TH1F* hPtPosRSelected;
		TH2F* hPdgCodemdZRejected;
		TH1F* hPtZSelected;
		TH2F* hPdgCodemdAngRejected;
		TH1F* hPtAngSelected;
		TH2F* hPdgCodemdPmRejected;
		TH1F* hPtPmSelected;
		TH2F* hPdgCodemdQtLowRejected;
		TH1F* hPtGoodKink; 
		TH1F* hEtaGoodKink; 
		TH1F* hRapidityGoodKink; 
		TH1F* hQtGoodKink; 
		TH2F* hPmGoodKinkAng;  
		TH2F* hPdgCodemdGoodKink; 
		TH2F* hPmdGoodKink;
		TH1F* hUIDGoodKinkDau;
		TH2F* hdEdxGoodKink;
		TH1F* hUIDPiDauPlus;
		TH1F* hMultPiPlus;
		TH1F* hPtPiPlus;
		TH1F* hEtaPiPlus;
		TH1F* hRapidityPiPlus;
		TH1F* hQtPiPlus;
		TH1F* hKinkAnglePiPlus;
		TH2F* hPmKinkAngPiPlus;
		TH2F* hKinkPosXYPiPlus;
		TH2F* hKinkPosZRPiPlus;
		TH1F* hKinkPosRPiPlus;
		TH1F* hDCAkinkPiPlus;
		TH2F* hPdgCodemdPiPlus;
		TH2F* hPmdPiPlus;
		TH2F* hdEdxPiPlus;
		TH1F* hQtPimuPlus;
		TH2F* hPmKinkAngPimuPlus;
		TH1F* hQtPiotherPlus;
		TH2F* hPmKinkAngPiotherPlus; 
		TH2F* hPdgCodemdPiotherPlus;
		TH1F* hUIDPiDauMinus;
		TH1F* hMultPiMinus;
		TH1F* hPtPiMinus;
		TH1F* hEtaPiMinus;
		TH1F* hRapidityPiMinus;
		TH1F* hQtPiMinus;
		TH1F* hKinkAnglePiMinus;
		TH2F* hPmKinkAngPiMinus;
		TH2F* hKinkPosXYPiMinus;
		TH2F* hKinkPosZRPiMinus;
		TH1F* hKinkPosRPiMinus;
		TH1F* hDCAkinkPiMinus;
		TH2F* hPdgCodemdPiMinus;
		TH2F* hPmdPiMinus;
		TH2F* hdEdxPiMinus;
		TH1F* hQtPimuMinus;
		TH2F* hPmKinkAngPimuMinus;
		TH1F* hQtPiotherMinus;
		TH2F* hPmKinkAngPiotherMinus; 
		TH2F* hPdgCodemdPiotherMinus;
		TH2F* hPdgCodemdQtRejected;
		TH1F* hPtQtSelected;
		TH2F* hPdgCodemdMaxAngRejected;
		TH1F* hPtMaxAngSelected;
		TH2F* hPdgCodemdRTPCclustersRejected;
		TH1F* hPtRTPCclustersSelected;
		TH2F* hRTPCclustersRTPCclustersSelected;
		TH2F* hPdgCodemdMinvRejected;
		TH1F* hPtSelected; 
		TH1F* hEtaSelected; 
		TH1F* hRapiditySelected; 
		TH1F* hQtSelected; 
		TH1F* hKinkAngleSelected;
		TH1F* hDCAkinkSelected;
		TH2F* hPmKinkAngSelected; 
		TH2F* hKinkPosXYSelected; 
		TH2F* hKinkPosZRSelected; 
		TH1F* hKinkPosRSelected; 
		TH2F* hPdgCodemdSelected; 
		TH2F* hPmdSelected;
		TH1F* hMinvPimuSelected;  
		TH1F* hUIDKinkDauSelected; 
		TH2F* hdEdxSelected; 
		TH1F* hPtSelectedFake; 
		TH1F* hEtaSelectedFake; 
		TH1F* hRapiditySelectedFake; 
		TH1F* hQtSelectedFake; 
		TH1F* hKinkAngleSelectedFake;
		TH1F* hDCAkinkSelectedFake;
		TH2F* hPmKinkAngSelectedFake;
		TH2F* hKinkPosXYSelectedFake;
		TH2F* hKinkPosZRSelectedFake;
		TH1F* hKinkPosRSelectedFake;
		TH2F* hPmdSelectedFake;
		TH1F* hMinvPimuSelectedFake; 
		TH2F* hdEdxSelectedFake;
		TH2F* hPdgCodemddEdxRejected;
		TH1F* hPtPiSelected; 
		TH1F* hEtaPiSelected; 
		TH1F* hRapidityPiSelected; 
		TH1F* hQtPiSelected; 
		TH1F* hKinkAnglePiSelected;
		TH1F* hDCAkinkPiSelected;
		TH2F* hPmKinkAngPiSelected; 
		TH1F* hKinkPosRTPCclusters1PiSelected;
		TH2F* hKinkPosRTPCclusters2PiSelected;
		TH2F* hKinkPosXYPiSelected; 
		TH2F* hKinkPosZRPiSelected; 
		TH1F* hKinkPosRPiSelected; 
		TH1F* hKinkPosZPiSelected; 
		TH2F* hPmPPiSelected; 
		TH2F* hPdgCodemdPiSelected; 
		TH2F* hPmdPiSelected;
		TH1F* hMinvPimuPiSelected;  
		TH1F* hUIDKinkDauPiSelected; 
		TH2F* hdEdxPiSelected;

		TH1F* hPtPiSelectedPlus; //plus
		TH1F* hEtaPiSelectedPlus; 
		TH1F* hRapidityPiSelectedPlus; 
		TH1F* hQtPiSelectedPlus; 
		TH1F* hKinkAnglePiSelectedPlus;
		TH1F* hDCAkinkPiSelectedPlus;
		TH2F* hPmKinkAngPiSelectedPlus; 
		TH2F* hKinkPosXYPiSelectedPlus; 
		TH2F* hKinkPosZRPiSelectedPlus;  
		TH2F* hPdgCodemdPiSelectedPlus; 
		TH2F* hPmdPiSelectedPlus;
		TH1F* hMinvPimuPiSelectedPlus;  
		TH1F* hUIDPiDaumuSelectedPlus; 
		TH2F* hdEdxPiSelectedPlus;
		TH1F* hPtPrimPiKinksPlus; 
		TH1F* hEtaPrimPiKinksPlus; 
		TH1F* hRapidityPrimPiKinksPlus;
		TH1F* hPtSecondPiKinksPlus; 
		TH1F* hEtaSecondPiKinksPlus; 
		TH1F* hRapiditySecondPiKinksPlus;
		TH1F* hPtNonPiKinksPlus; 
		TH1F* hEtaNonPiKinksPlus; 
		TH1F* hRapidityNonPiKinksPlus; 
		TH2F* hPdgCodemdNonPiKinksPlus; 
		TH1F* hPtPiSelectedMinus;//minus
		TH1F* hEtaPiSelectedMinus; 
		TH1F* hRapidityPiSelectedMinus; 
		TH1F* hQtPiSelectedMinus; 
		TH1F* hKinkAnglePiSelectedMinus;
		TH1F* hDCAkinkPiSelectedMinus;
		TH2F* hPmKinkAngPiSelectedMinus; 
		TH2F* hKinkPosXYPiSelectedMinus; 
		TH2F* hKinkPosZRPiSelectedMinus;  
		TH2F* hPdgCodemdPiSelectedMinus; 
		TH2F* hPmdPiSelectedMinus;
		TH1F* hMinvPimuPiSelectedMinus;  
		TH1F* hUIDPiDaumuSelectedMinus; 
		TH2F* hdEdxPiSelectedMinus;
		TH1F* hPtPrimPiKinksMinus; 
		TH1F* hEtaPrimPiKinksMinus; 
		TH1F* hRapidityPrimPiKinksMinus;
		TH1F* hPtSecondPiKinksMinus; 
		TH1F* hEtaSecondPiKinksMinus; 
		TH1F* hRapiditySecondPiKinksMinus;
		TH1F* hPtNonPiKinksMinus; 
		TH1F* hEtaNonPiKinksMinus; 
		TH1F* hRapidityNonPiKinksMinus; 
		TH2F* hPdgCodemdNonPiKinksMinus;// reconstruction histograms
		
		TList* fListOfHistos;

		// Limits and cuts
		Int_t fLowMulcut;
		Int_t fUpMulcut;
		
		Double_t cLowPt;
		Double_t cRapidityLim;
		Double_t cLowR, cUpR;
		Double_t cLowZ, cUpZ;
		Double_t cLowKinkAngle;
		Double_t cLowQt, cUpQt;
		Double_t cLowInvMass, cUpInvMass;
		Double_t cSigmaCut;
		Int_t cPdgKaon, cPdgPion, cPdgMuon, cPdgElectron;
		Double_t cKaonMass, cPionMass, cMuonMass, cElectronMass;

		// Histos limits
		Int_t nBinsMult, hLowMult, hUpMult;
		Int_t nBinsPt;
		Double_t hLowPt, hUpPt;
		Int_t nBinsEta;
		Double_t hLowEta, hUpEta;
		Int_t nBinsQt;
		Double_t hLowQt, hUpQt;
		Int_t nBinsPdg;
		Double_t hLowPdg, hUpPdg;
		Int_t nBinsPdg2;
		Double_t hLowPdg2, hUpPdg2;
		Int_t nBinsUID;
		Double_t hLowUID, hUpUID;
		Int_t nBinsR;
		Double_t hLowR, hUpR;
		Int_t nBinsZ;
		Double_t hLowZ, hUpZ;
		Int_t nBinsXY;
		Double_t hLowXY, hUpXY;
		Int_t nBinsAngle;
		Double_t hLowAngle, hUpAngle;
		Int_t nBinsZV;
		Double_t hLowZV, hUpZV;
		Int_t nBinsXYV;
		Double_t hLowXYV, hUpXYV;
		Int_t nBinsInvMass;
		Double_t hLowInvMass, hUpInvMass;
		Int_t nBinsdEdx;
		Double_t hLowdEdx, hUpdEdx;

		AliPIDResponse *fPIDResponse;     //! PID response object
		AliESDtrackCuts* fMaxDCAtoVtxCut;
		AliESDtrackCuts* fTrackCuts;

		
		AliAnalysisPionKinksMCESD(const AliAnalysisPionKinksMCESD&); 
		AliAnalysisPionKinksMCESD& operator=(const AliAnalysisPionKinksMCESD&);
  
	ClassDef(AliAnalysisPionKinksMCESD, 1); 
};
#endif
Commit	Line	Data
bd4be5f8	1	#ifndef AliAnalysisPionKinksMCESD_h
	2	#define AliAnalysisPionKinksMCESD_h
	3
	4
30687e36	5	/* See cxx source for full Copyright notice */
	6
	7	//-----------------------------------------------------------------
	8	// AliAnalysisPionKinksMCESD class
	9	// This task is an example of an analysis task
	10	// for kink topology Study
	11	// Authors: Eftychios Cheiladakis, under the supervision
	12	// of Martha's Spyropoulou-Stassinaki
	13	// Physics Department of Athens University
	14	//------------------------------------------------------------------
bd4be5f8	15	class AliESDEvent;
	16	class TF1;
	17	class TH1;
	18	class TH2;
	19	class TH3;
	20	class TParticle;
	21	class TParticle;
	22	class AliESDtrackCuts;
	23	class AliPhysicsSelection;
	24	class AliPIDResponse;
	25
	26	#include "AliAnalysisTaskSE.h"
	27
	28	class AliAnalysisPionKinksMCESD : public AliAnalysisTaskSE {
	29	public:
	30	AliAnalysisPionKinksMCESD(const char *name = "AliAnalysisPionKinksMCESD");
	31	virtual ~AliAnalysisPionKinksMCESD() {}
	32
	33	virtual void UserCreateOutputObjects();
	34	virtual void UserExec(Option_t *option);
	35	virtual void Terminate(Option_t *);
	36
	37	const AliESDVertex GetEventVertex(AliESDEvent esd);
	38	Double_t Energy(AliESDtrack* track) const;
	39	Double_t GetMCRapidity(TParticle* particle) const;
	40	Double_t fuRapidity(AliESDtrack* track) const;
	41	Double_t MCPQt(AliMCEvent* mcEvent, Int_t iMC, TParticle* MCdaughter) const;
	42	Double_t fuMCKinkAngle(AliMCEvent* mcEvent, Int_t iMC, TParticle* MCdaughter, Bool_t degrees) const;
	43	Bool_t IsGoodTrack(AliESDtrack* ESDTrack) const;
	44	Bool_t IsPrimaryTrack(AliESDtrack* ESDTrack) const;
	45	Bool_t IsRealPionKink(AliESDkink* kink, AliStack* MCstack) const;
	46
	47	// Set limits & cuts
	48	void SetMulCut(Int_t low, Int_t up){fLowMulcut=low; fUpMulcut=up;}
	49	void SetPtCut(Double_t PtCut){cLowPt=PtCut;}
	50	void SetRapidityLimits(Double_t RapidityLim){cRapidityLim=RapidityLim;}
	51	void SetRadiusRange(Double_t LowR, Double_t UpR){cLowR=LowR; cUpR=UpR;}
	52	void SetZRange(Double_t LowZ, Double_t UpZ){cLowZ=LowZ; cUpZ=UpZ;}
	53	void SetKinkAngleCut(Double_t LowAngle){cLowKinkAngle=LowAngle;}
	54	void SetQtRange(Double_t LowQt, Double_t UpQt){cLowQt=LowQt; cUpQt=UpQt;}
	55	void SetInvMassRange(Double_t LowInvMass, Double_t UpInvMass){cLowInvMass=LowInvMass; cUpInvMass=UpInvMass;}
	56	void SetSigmaCut(Double_t SigmaCut){cSigmaCut=SigmaCut;}
	57	void SetPdgCodes(Int_t kaon, Int_t pion, Int_t muon, Int_t electron){cPdgKaon=kaon; cPdgPion=pion; cPdgMuon=muon; cPdgElectron=electron;}
	58	void SetMasses(Double_t KaonMass, Double_t PionMass, Double_t MuonMass, Double_t ElectronMass){cKaonMass=KaonMass; cPionMass=PionMass; cMuonMass=MuonMass; cElectronMass=ElectronMass;}
	59
	60	// Set histos limits
	61	void SetMultHistos(Int_t BinsMult, Int_t LowMult, Int_t UpMult){nBinsMult=BinsMult; hLowPt=LowMult; hUpMult=UpMult;}
	62	void SetPtHistos(Int_t BinsPt, Double_t LowPt, Double_t UpPt){nBinsPt=BinsPt; hLowPt=LowPt; hUpPt=UpPt;}
	63	void SetEtaHistos(Int_t BinsEta, Double_t LowEta, Double_t UpEta){nBinsEta=BinsEta; hLowEta=LowEta; hUpEta=UpEta;}
	64	void SetQtHistos(Int_t BinsQt, Double_t LowQt, Double_t UpQt){nBinsQt=BinsQt; hLowQt=LowQt; hUpQt=UpQt;}
	65	void SetPdgHistos(Int_t BinsPdg, Double_t LowPdg, Double_t UpPdg){nBinsPdg=BinsPdg; hLowPdg=LowPdg; hUpPdg=UpPdg;}
	66	void SetPdg2Histos(Int_t BinsPdg2, Double_t LowPdg2, Double_t UpPdg2){nBinsPdg2=BinsPdg2; hLowPdg2=LowPdg2; hUpPdg2=UpPdg2;}
	67	void SetUIDHistos(Int_t BinsUID, Double_t LowUID, Double_t UpUID){nBinsUID=BinsUID; hLowUID=LowUID; hUpUID=UpUID;}
	68	void SetRHistos(Int_t BinsR, Double_t LowR, Double_t UpR){nBinsR=BinsR; hLowR=LowR; hUpR=UpR;}
	69	void SetZHistos(Int_t BinsZ, Double_t LowZ, Double_t UpZ){nBinsZ=BinsZ; hLowZ=LowZ; hUpZ=UpZ;}
	70	void SetXYHistos(Int_t BinsXY, Double_t LowXY, Double_t UpXY){nBinsXY=BinsXY; hLowXY=LowXY; hUpXY=UpXY;}
	71	void SetAngleHistos(Int_t BinsAngle, Double_t LowAngle, Double_t UpAngle){nBinsAngle=BinsAngle; hLowAngle=LowAngle; hUpAngle=UpAngle;}
	72	void SetZVHistos(Int_t BinsZV, Double_t LowZV, Double_t UpZV){nBinsZV=BinsZV; hLowZV=LowZV; hUpZV=UpZV;}
	73	void SetXYVHistos(Int_t BinsXYV, Double_t LowXYV, Double_t UpXYV){nBinsXYV=BinsXYV; hLowXYV=LowXYV; hUpXYV=UpXYV;}
	74	void SetInvMassHistos(Int_t BinsInvMass, Double_t LowInvMass, Double_t UpInvMass){nBinsInvMass=BinsInvMass; hLowInvMass=LowInvMass; hUpInvMass=UpInvMass;}
	75	void SetdEdxHistos(Int_t BinsdEdx, Double_t LowdEdx, Double_t UpdEdx){nBinsdEdx=BinsdEdx; hLowdEdx=LowdEdx; hUpdEdx=UpdEdx;}
	76	private:
	77	TF1* fMaxKinkAngKmu;
	78	TF1* fMaxKinkAngPimu;
79
80	TH1F* hMCMult;
81	TH1F* hMCMultPrim;
82	TH1F* hMCPtAll;
83	TH1F* hMCEtaAll;
84	TH1F* hMCPtPrim;
85	TH1F* hMCEtaPrim;
86	TH1F* hMCPt;
87	TH1F* hMCEta;
88	TH1F* hMCPdg;
89	TH1F* hMCMultPiPlus;
90	TH1F* hMCPtPiPlus;
91	TH1F* hMCEtaPiPlus;
92	TH1F* hMCRapidityPiPlus;
93	TH1F* hMCNDaughtersPlus;
94	TH1F* hMCRadPiDauPlus;
95	TH1F* hMCKinkPosZPlus;
96	TH1F* hMCUIDPiDauPlus;
97	TH1F* hMCPdgPiNonDecayedPlus;
98	TH1F* hMCPdgPiDauPlus;
99	TH1F* hMCQtPlus;
100	TH1F* hMCKinkAnglePlus;
101	TH2F* hMCPKinkAngPlus;
102	TH2F* hMCPdgCodemdPlus;
103	TH2F* hMCPtmdPlus;
104	TH1F* hMCPtPimuonPlus;
105	TH1F* hMCEtaPimuonPlus;
106	TH1F* hMCRapidityPimuonPlus;
107	TH1F* hMCQtPimuonPlus;
108	TH1F* hMCPKinkAngPimuonPlus;
109	TH1F* hMCPtPiotherPlus;
110	TH1F* hMCEtaPiotherPlus;
111	TH1F* hMCRapidityPiotherPlus;
112	TH1F* hMCQtPiotherPlus;
113	TH1F* hMCPKinkAngPiotherPlus;
114	TH1F* hMCMultPiMinus;
115	TH1F* hMCPtPiMinus;
116	TH1F* hMCEtaPiMinus;
117	TH1F* hMCRapidityPiMinus;
118	TH1F* hMCNDaughtersMinus;
119	TH1F* hMCRadPiDauMinus;
120	TH1F* hMCKinkPosZMinus;
121	TH1F* hMCUIDPiDauMinus;
122	TH1F* hMCPdgPiNonDecayedMinus;
123	TH1F* hMCPdgPiDauMinus;
124	TH1F* hMCQtMinus;
125	TH1F* hMCKinkAngleMinus;
126	TH2F* hMCPKinkAngMinus;
127	TH2F* hMCPdgCodemdMinus;
128	TH2F* hMCPtmdMinus;
129	TH1F* hMCPtPimuonMinus;
130	TH1F* hMCEtaPimuonMinus;
131	TH1F* hMCRapidityPimuonMinus;
132	TH1F* hMCQtPimuonMinus;
133	TH1F* hMCPKinkAngPimuonMinus;
134	TH1F* hMCPtPiotherMinus;
135	TH1F* hMCEtaPiotherMinus;
136	TH1F* hMCRapidityPiotherMinus;
137	TH1F* hMCQtPiotherMinus;
138	TH1F* hMCPKinkAngPiotherMinus;
139
140	TH1F* hMult;
141	TH1F* hAcceptedMult;
142	TH1F* hMultPS;
143	TH3F* hvtx;
144	TH2F* hvtxy;
145	TH2F* hvtyz;
146	TH2F* hvtxz;
147	TH1F* hMultPSV;
148	TH1F* hPtAll;
149	TH1F* hEtaAll;
150	TH3F* hTrackPos;
151	TH2F* hTrackPosxy;
152	TH2F* hTrackPosyz;
153	TH2F* hTrackPosxz;
154	//TH1F* hTPCchi2clusters;
155	//TH1F* hdcaToVertexXY;
156	//TH1F* hdcaToVertexZ;
157	TH1F* hMultPrim;
158	TH1F* hPtPrim;
159	TH1F* hEtaPrim;
160	TH3F* hPrimTrackPos;
161	TH2F* hPrimTrackPosxy;
162	TH2F* hPrimTrackPosyz;
163	TH2F* hPrimTrackPosxz;
164	TH1F* hPt;
165	TH1F* hEta;
166	//TH1F* hRapidity;
167	TH1F* hPtKink;
168	TH1F* hEtaKink;
169	TH1F* hRapidityKink;
170	TH2F* hPmP;
171	TH1F* hKinkPosRTPCclusters1;
172	TH2F* hKinkPosRTPCclusters2;
173	TH1F* hQt;
174	TH1F* hKinkAngle;
175	TH1F* hDCAkink;
176	TH2F* hPmKinkAng;
177	TH2F* hKinkPosXY;
178	TH2F* hKinkPosZY;
179	TH2F* hKinkPosZR;
180	TH1F* hKinkPosR;
181	TH1F* hKinkPosZ;
182	TH2F* hKinkPosZMCKinkPosZ;
183	TH2F* hPdgCodemd;
184	TH2F* hPmd;
185	TH1F* hMinvPimu;
186	TH1F* hUIDKinkDau;
187	TH2F* hdEdx;
188	TH1F* hPtKinkFake;
189	TH1F* hEtaKinkFake;
190	TH1F* hRapidityKinkFake;
191	TH2F* hPmPFake;
192	TH1F* hKinkPosRTPCclusters1Fake;
193	TH2F* hKinkPosRTPCclusters2Fake;
194	TH1F* hQtFake;
195	TH1F* hKinkAngleFake;
196	TH1F* hDCAkinkFake;
197	TH2F* hPmKinkAngFake;
198	TH2F* hKinkPosXYFake;
199	TH2F* hKinkPosZYFake;
200	TH2F* hKinkPosZRFake;
201	TH1F* hKinkPosRFake;
202	TH1F* hKinkPosZFake;
203	TH2F* hKinkPosZMCKinkPosZFake;
204	TH2F* hPdgCodemdFake;
205	TH2F* hPmdFake;
206	TH1F* hMinvPimuFake;
207	TH1F* hUIDKinkDauFake;
208	TH2F* hdEdxFake;
209	TH1F* hPtPosRSelected;
210	TH2F* hPdgCodemdZRejected;
211	TH1F* hPtZSelected;
212	TH2F* hPdgCodemdAngRejected;
213	TH1F* hPtAngSelected;
214	TH2F* hPdgCodemdPmRejected;
215	TH1F* hPtPmSelected;
216	TH2F* hPdgCodemdQtLowRejected;
217	TH1F* hPtGoodKink;
218	TH1F* hEtaGoodKink;
219	TH1F* hRapidityGoodKink;
220	TH1F* hQtGoodKink;
221	TH2F* hPmGoodKinkAng;
222	TH2F* hPdgCodemdGoodKink;
223	TH2F* hPmdGoodKink;
224	TH1F* hUIDGoodKinkDau;
225	TH2F* hdEdxGoodKink;
226	TH1F* hUIDPiDauPlus;
227	TH1F* hMultPiPlus;
228	TH1F* hPtPiPlus;
229	TH1F* hEtaPiPlus;
230	TH1F* hRapidityPiPlus;
231	TH1F* hQtPiPlus;
232	TH1F* hKinkAnglePiPlus;
233	TH2F* hPmKinkAngPiPlus;
234	TH2F* hKinkPosXYPiPlus;
235	TH2F* hKinkPosZRPiPlus;
236	TH1F* hKinkPosRPiPlus;
237	TH1F* hDCAkinkPiPlus;
238	TH2F* hPdgCodemdPiPlus;
239	TH2F* hPmdPiPlus;
240	TH2F* hdEdxPiPlus;
241	TH1F* hQtPimuPlus;
242	TH2F* hPmKinkAngPimuPlus;
243	TH1F* hQtPiotherPlus;
244	TH2F* hPmKinkAngPiotherPlus;
245	TH2F* hPdgCodemdPiotherPlus;
246	TH1F* hUIDPiDauMinus;
247	TH1F* hMultPiMinus;
248	TH1F* hPtPiMinus;
249	TH1F* hEtaPiMinus;
250	TH1F* hRapidityPiMinus;
251	TH1F* hQtPiMinus;
252	TH1F* hKinkAnglePiMinus;
253	TH2F* hPmKinkAngPiMinus;
254	TH2F* hKinkPosXYPiMinus;
255	TH2F* hKinkPosZRPiMinus;
256	TH1F* hKinkPosRPiMinus;
257	TH1F* hDCAkinkPiMinus;
258	TH2F* hPdgCodemdPiMinus;
259	TH2F* hPmdPiMinus;
260	TH2F* hdEdxPiMinus;
261	TH1F* hQtPimuMinus;
262	TH2F* hPmKinkAngPimuMinus;
263	TH1F* hQtPiotherMinus;
264	TH2F* hPmKinkAngPiotherMinus;
265	TH2F* hPdgCodemdPiotherMinus;
266	TH2F* hPdgCodemdQtRejected;
267	TH1F* hPtQtSelected;
268	TH2F* hPdgCodemdMaxAngRejected;
269	TH1F* hPtMaxAngSelected;
270	TH2F* hPdgCodemdRTPCclustersRejected;
271	TH1F* hPtRTPCclustersSelected;
272	TH2F* hRTPCclustersRTPCclustersSelected;
273	TH2F* hPdgCodemdMinvRejected;
274	TH1F* hPtSelected;
275	TH1F* hEtaSelected;
276	TH1F* hRapiditySelected;
277	TH1F* hQtSelected;
278	TH1F* hKinkAngleSelected;
279	TH1F* hDCAkinkSelected;
280	TH2F* hPmKinkAngSelected;
281	TH2F* hKinkPosXYSelected;
282	TH2F* hKinkPosZRSelected;
283	TH1F* hKinkPosRSelected;
284	TH2F* hPdgCodemdSelected;
285	TH2F* hPmdSelected;
286	TH1F* hMinvPimuSelected;
287	TH1F* hUIDKinkDauSelected;
288	TH2F* hdEdxSelected;
289	TH1F* hPtSelectedFake;
290	TH1F* hEtaSelectedFake;
291	TH1F* hRapiditySelectedFake;
292	TH1F* hQtSelectedFake;
293	TH1F* hKinkAngleSelectedFake;
294	TH1F* hDCAkinkSelectedFake;
295	TH2F* hPmKinkAngSelectedFake;
296	TH2F* hKinkPosXYSelectedFake;
297	TH2F* hKinkPosZRSelectedFake;
298	TH1F* hKinkPosRSelectedFake;
299	TH2F* hPmdSelectedFake;
300	TH1F* hMinvPimuSelectedFake;
301	TH2F* hdEdxSelectedFake;
302	TH2F* hPdgCodemddEdxRejected;
303	TH1F* hPtPiSelected;
304	TH1F* hEtaPiSelected;
305	TH1F* hRapidityPiSelected;
306	TH1F* hQtPiSelected;
307	TH1F* hKinkAnglePiSelected;
308	TH1F* hDCAkinkPiSelected;
309	TH2F* hPmKinkAngPiSelected;
310	TH1F* hKinkPosRTPCclusters1PiSelected;
311	TH2F* hKinkPosRTPCclusters2PiSelected;
312	TH2F* hKinkPosXYPiSelected;
313	TH2F* hKinkPosZRPiSelected;
314	TH1F* hKinkPosRPiSelected;
315	TH1F* hKinkPosZPiSelected;
316	TH2F* hPmPPiSelected;
317	TH2F* hPdgCodemdPiSelected;
318	TH2F* hPmdPiSelected;
319	TH1F* hMinvPimuPiSelected;
320	TH1F* hUIDKinkDauPiSelected;
321	TH2F* hdEdxPiSelected;
322
323	TH1F* hPtPiSelectedPlus; //plus
324	TH1F* hEtaPiSelectedPlus;
325	TH1F* hRapidityPiSelectedPlus;
326	TH1F* hQtPiSelectedPlus;
327	TH1F* hKinkAnglePiSelectedPlus;
328	TH1F* hDCAkinkPiSelectedPlus;
329	TH2F* hPmKinkAngPiSelectedPlus;
330	TH2F* hKinkPosXYPiSelectedPlus;
331	TH2F* hKinkPosZRPiSelectedPlus;
332	TH2F* hPdgCodemdPiSelectedPlus;
333	TH2F* hPmdPiSelectedPlus;
334	TH1F* hMinvPimuPiSelectedPlus;
335	TH1F* hUIDPiDaumuSelectedPlus;
336	TH2F* hdEdxPiSelectedPlus;
337	TH1F* hPtPrimPiKinksPlus;
338	TH1F* hEtaPrimPiKinksPlus;
339	TH1F* hRapidityPrimPiKinksPlus;
340	TH1F* hPtSecondPiKinksPlus;
341	TH1F* hEtaSecondPiKinksPlus;
342	TH1F* hRapiditySecondPiKinksPlus;
343	TH1F* hPtNonPiKinksPlus;
344	TH1F* hEtaNonPiKinksPlus;
345	TH1F* hRapidityNonPiKinksPlus;
346	TH2F* hPdgCodemdNonPiKinksPlus;
347	TH1F* hPtPiSelectedMinus;//minus
348	TH1F* hEtaPiSelectedMinus;
349	TH1F* hRapidityPiSelectedMinus;
350	TH1F* hQtPiSelectedMinus;
351	TH1F* hKinkAnglePiSelectedMinus;
352	TH1F* hDCAkinkPiSelectedMinus;
353	TH2F* hPmKinkAngPiSelectedMinus;
354	TH2F* hKinkPosXYPiSelectedMinus;
355	TH2F* hKinkPosZRPiSelectedMinus;
356	TH2F* hPdgCodemdPiSelectedMinus;
357	TH2F* hPmdPiSelectedMinus;
358	TH1F* hMinvPimuPiSelectedMinus;
359	TH1F* hUIDPiDaumuSelectedMinus;
360	TH2F* hdEdxPiSelectedMinus;
361	TH1F* hPtPrimPiKinksMinus;
362	TH1F* hEtaPrimPiKinksMinus;
363	TH1F* hRapidityPrimPiKinksMinus;
364	TH1F* hPtSecondPiKinksMinus;
365	TH1F* hEtaSecondPiKinksMinus;
366	TH1F* hRapiditySecondPiKinksMinus;
367	TH1F* hPtNonPiKinksMinus;
368	TH1F* hEtaNonPiKinksMinus;
369	TH1F* hRapidityNonPiKinksMinus;
370	TH2F* hPdgCodemdNonPiKinksMinus;// reconstruction histograms
371
372	TList* fListOfHistos;
373
374	// Limits and cuts
375	Int_t fLowMulcut;
376	Int_t fUpMulcut;
377
378	Double_t cLowPt;
379	Double_t cRapidityLim;
380	Double_t cLowR, cUpR;
381	Double_t cLowZ, cUpZ;
382	Double_t cLowKinkAngle;
383	Double_t cLowQt, cUpQt;
384	Double_t cLowInvMass, cUpInvMass;
385	Double_t cSigmaCut;
386	Int_t cPdgKaon, cPdgPion, cPdgMuon, cPdgElectron;
387	Double_t cKaonMass, cPionMass, cMuonMass, cElectronMass;
388
389	// Histos limits
390	Int_t nBinsMult, hLowMult, hUpMult;
391	Int_t nBinsPt;
392	Double_t hLowPt, hUpPt;
393	Int_t nBinsEta;
394	Double_t hLowEta, hUpEta;
395	Int_t nBinsQt;
396	Double_t hLowQt, hUpQt;
397	Int_t nBinsPdg;
398	Double_t hLowPdg, hUpPdg;
399	Int_t nBinsPdg2;
400	Double_t hLowPdg2, hUpPdg2;
401	Int_t nBinsUID;
402	Double_t hLowUID, hUpUID;
403	Int_t nBinsR;
404	Double_t hLowR, hUpR;
405	Int_t nBinsZ;
406	Double_t hLowZ, hUpZ;
407	Int_t nBinsXY;
408	Double_t hLowXY, hUpXY;
409	Int_t nBinsAngle;
410	Double_t hLowAngle, hUpAngle;
411	Int_t nBinsZV;
412	Double_t hLowZV, hUpZV;
413	Int_t nBinsXYV;
414	Double_t hLowXYV, hUpXYV;
415	Int_t nBinsInvMass;
416	Double_t hLowInvMass, hUpInvMass;
417	Int_t nBinsdEdx;
418	Double_t hLowdEdx, hUpdEdx;
419
420	AliPIDResponse *fPIDResponse; //! PID response object
421	AliESDtrackCuts* fMaxDCAtoVtxCut;
422	AliESDtrackCuts* fTrackCuts;
423
424
425	AliAnalysisPionKinksMCESD(const AliAnalysisPionKinksMCESD&);
426	AliAnalysisPionKinksMCESD& operator=(const AliAnalysisPionKinksMCESD&);
427
428	ClassDef(AliAnalysisPionKinksMCESD, 1);
429	};
430	#endif
431
432