[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / AliFemtoUser / AliFemtoKKTrackCut.h

///////////////////////////////////////////////////////////////////////////
//                                                                       //
// AliFemtoESDTrackCut: A basic track cut that used information from     //
// ALICE ESD to accept or reject the track.                              //  
// Enables the selection on charge, transverse momentum, rapidity,       //
// pid probabilities, number of ITS and TPC clusters                     //
// Author: Marek Chojnacki (WUT), mchojnacki@knf.pw.edu.pl               //
//                                                                       //
///////////////////////////////////////////////////////////////////////////
#include "AliESDtrackCuts.h"

#ifndef ALIFEMTOKKTRACKCUT_H
#define ALIFEMTOKKTRACKCUT_H

//#ifndef StMaker_H
//#include "StMaker.h"
//#endif

#include "AliESDtrackCuts.h" //for enum with ITS layers
#include "AliFemtoTrackCut.h"


class AliFemtoKKTrackCut : public AliFemtoTrackCut 
{
  public:

  enum PIDMethodType {knSigma=0, kContour=1};
  typedef enum PIDMethodType ReadPIDMethodType; 

  AliFemtoKKTrackCut();
  virtual ~AliFemtoKKTrackCut();

  virtual bool Pass(const AliFemtoTrack* aTrack);

  virtual AliFemtoString Report();
  virtual TList *ListSettings();
  virtual AliFemtoParticleType Type(){return hbtTrack;}

  void SetPt(const float& lo, const float& hi);
  void SetRapidity(const float& lo, const float& hi);
  void SetEta(const float& lo, const float& hi);
  void SetCharge(const int& ch);
  void SetPidProbElectron(const float& lo, const float& hi);
  void SetPidProbPion(const float& lo, const float& hi);
  void SetPidProbKaon(const float& lo, const float& hi);
  void SetPidProbProton(const float& lo, const float& hi);
  void SetPidProbMuon(const float& lo, const float& hi);
  void SetLabel(const bool& flag);
  void SetStatus(const long& w);
  void SetminTPCclsF(const short& s);
  void SetminTPCncls(const short& s);
  void SetminITScls(const int& s);
  void SetRemoveKinks(const bool& flag);
  void SetRemoveITSFake(const bool& flag);
  void SetMaxITSChiNdof(const float& maxchi);
  void SetMaxTPCChiNdof(const float& maxchi);
  void SetMaxSigmaToVertex(const float& maxsig);
  void SetMaxImpactXY(const float& maximpxy);
  void SetMaxImpactZ(const float& maximpz);
  void SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow);
  void SetMostProbablePion();
  void SetMostProbableKaon();
  void SetMostProbableProton();
  void SetNoMostProbable(); 
  void SetPIDMethod(ReadPIDMethodType newMethod);
  void SetClusterRequirementITS(AliESDtrackCuts::Detector det, AliESDtrackCuts::ITSClusterRequirement req = AliESDtrackCuts::kOff);

  void SetMomRangeTOFpidIs(const float& minp, const float& maxp);
  void SetMomRangeTPCpidIs(const float& minp, const float& maxp);
  void SetMomRangeITSpidIs(const float& minp, const float& maxp);

  //ml
 void SetNsigmaTPCle250(Double_t);
 void SetNsigmaTPC250_400(Double_t);
 void SetNsigmaTPC400_450(Double_t);
 void SetNsigmaTPC450_500(Double_t);
 void SetNsigmaTPCge500(Double_t);

 void SetNsigmaTOF500_800(Double_t);
 void SetNsigmaTOF800_1000(Double_t);
 void SetNsigmaTOFge1000(Double_t);


 private:   // here are the quantities I want to cut on...

  int               fCharge;             // particle charge
  float             fPt[2];              // bounds for transverse momentum
  float             fRapidity[2];        // bounds for rapidity
  float             fEta[2];             // bounds for pseudorapidity
  float             fPidProbElectron[2]; // bounds for electron probability
  float             fPidProbPion[2];     // bounds for pion probability
  float             fPidProbKaon[2];     // bounds for kaon probability
  float             fPidProbProton[2];   // bounds for proton probability
  float             fPidProbMuon[2];     // bounds for muon probability 

  AliESDtrackCuts::ITSClusterRequirement fCutClusterRequirementITS[3];  // detailed ITS cluster requirements for (SPD, SDD, SSD) - from AliESDtrackcuts!
  bool              fLabel;              // if true label<0 will not pass throught 
  long              fStatus;             // staus flag
  ReadPIDMethodType fPIDMethod;          // which PID mehod to use. 0 - nsgima, 1 - contour 

  //ml
 Double_t fNsigmaTPCle250;
 Double_t fNsigmaTPC250_400;
 Double_t fNsigmaTPC400_450;
 Double_t fNsigmaTPC450_500;
 Double_t fNsigmaTPCge500;

 Double_t fNsigmaTOF500_800;
 Double_t fNsigmaTOF800_1000;
 Double_t fNsigmaTOFge1000;


  short             fminTPCclsF;         // min number of findable clusters in the TPC
  short             fminTPCncls;         // min number of clusters in the TPC
  int               fminITScls;          // min number of clusters assigned in the ITS 
  float             fMaxITSchiNdof;      // maximum allowed chi2/ndof for ITS clusters
  float             fMaxTPCchiNdof;      // maximum allowed chi2/ndof for TPC clusters
  float             fMaxSigmaToVertex;   // maximum allowed sigma to primary vertex
  long              fNTracksPassed;      // passed tracks count
  long              fNTracksFailed;      // failed tracks count
  bool              fRemoveKinks;        // if true particles with any kink label will not pass
  bool              fRemoveITSFake;      // if true particles with ITS fake flag will not pass
  int               fMostProbable;       // this particle type is required to be most probable

  float             fMaxImpactXY;        // Max XY impact parameter
  float             fMaxImpactZ;         // Max Z impact parameter

  float             fMaxImpactXYPtOff;   // Max XY DCA Pt dependent offset
  float             fMaxImpactXYPtNrm;   // Max XY DCA Pt dependent normalization
  float             fMaxImpactXYPtPow;   // Max XY DCA Pt dependent power

  float             fMinPforTOFpid;  // momentum from which TOF PID is requested
  float             fMaxPforTOFpid;  // momentum till which TOF PID is requested
  float             fMinPforTPCpid;  // momentum from which TPC PID is requested
  float             fMaxPforTPCpid;  // momentum till which TPC PID is requested
  float             fMinPforITSpid;  // momentum from which ITS PID is requested
  float             fMaxPforITSpid;  // momentum till which ITS PID is requested

  float PidFractionElectron(float mom) const;
  float PidFractionPion(float mom) const;
  float PidFractionKaon(float mom) const;
  float PidFractionProton(float mom) const;

  bool IsPionTPCdEdx(float mom, float dEdx);
  bool IsKaonTPCdEdx(float mom, float dEdx);
  bool IsProtonTPCdEdx(float mom, float dEdx);

  bool IsPionTOFTime(float mom, float ttof);
  bool IsKaonTOFTime(float mom, float ttof);
  bool IsProtonTOFTime(float mom, float ttof);

  bool IsKaonTPCdEdxNSigma(float mom, float nsigma);
  bool IsKaonTOFNSigma(float mom, float nsigma);
  bool IsKaonNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
  bool IsPionNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
  bool IsProtonNSigma(float mom, float nsigmaTPC, float nsigmaTOF);

  Bool_t CheckITSClusterRequirement(AliESDtrackCuts::ITSClusterRequirement req, Bool_t clusterL1, Bool_t clusterL2); //the same as in AliESDtrackCuts


#ifdef __ROOT__ 
  ClassDef(AliFemtoKKTrackCut, 1)
#endif
    };


inline void AliFemtoKKTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
inline void AliFemtoKKTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
inline void AliFemtoKKTrackCut::SetEta(const float& lo,const float& hi){fEta[0]=lo; fEta[1]=hi;}
inline void AliFemtoKKTrackCut::SetCharge(const int& ch){fCharge = ch;}
inline void AliFemtoKKTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
inline void AliFemtoKKTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
inline void AliFemtoKKTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
inline void AliFemtoKKTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
inline void AliFemtoKKTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
inline void AliFemtoKKTrackCut::SetLabel(const bool& flag){fLabel=flag;}
inline void AliFemtoKKTrackCut::SetStatus(const long& status){fStatus=status;}
inline void AliFemtoKKTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
inline void AliFemtoKKTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
inline void AliFemtoKKTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
inline void AliFemtoKKTrackCut::SetMostProbablePion() { fMostProbable = 2; }
inline void AliFemtoKKTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
inline void AliFemtoKKTrackCut::SetMostProbableProton() { fMostProbable = 4; }
inline void AliFemtoKKTrackCut::SetNoMostProbable() { fMostProbable = 0; }
inline void AliFemtoKKTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
inline void AliFemtoKKTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
inline void AliFemtoKKTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
inline void AliFemtoKKTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; }
inline void AliFemtoKKTrackCut::SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow) { fMaxImpactXYPtOff = maxoff; fMaxImpactXYPtNrm = maxnrm; fMaxImpactXYPtPow = maxpow; }
inline void AliFemtoKKTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; }


#endif
Commit	Line	Data
9bba4e99	1	///////////////////////////////////////////////////////////////////////////
	2	// //
	3	// AliFemtoESDTrackCut: A basic track cut that used information from //
	4	// ALICE ESD to accept or reject the track. //
	5	// Enables the selection on charge, transverse momentum, rapidity, //
	6	// pid probabilities, number of ITS and TPC clusters //
	7	// Author: Marek Chojnacki (WUT), mchojnacki@knf.pw.edu.pl //
	8	// //
	9	///////////////////////////////////////////////////////////////////////////
	10	#include "AliESDtrackCuts.h"
	11
	12	#ifndef ALIFEMTOKKTRACKCUT_H
	13	#define ALIFEMTOKKTRACKCUT_H
	14
	15	//#ifndef StMaker_H
	16	//#include "StMaker.h"
	17	//#endif
	18
	19	#include "AliESDtrackCuts.h" //for enum with ITS layers
	20	#include "AliFemtoTrackCut.h"
	21
	22
	23	class AliFemtoKKTrackCut : public AliFemtoTrackCut
	24	{
	25	public:
	26
	27	enum PIDMethodType {knSigma=0, kContour=1};
	28	typedef enum PIDMethodType ReadPIDMethodType;
	29
	30	AliFemtoKKTrackCut();
	31	virtual ~AliFemtoKKTrackCut();
	32
	33	virtual bool Pass(const AliFemtoTrack* aTrack);
	34
	35	virtual AliFemtoString Report();
	36	virtual TList *ListSettings();
	37	virtual AliFemtoParticleType Type(){return hbtTrack;}
	38
	39	void SetPt(const float& lo, const float& hi);
	40	void SetRapidity(const float& lo, const float& hi);
	41	void SetEta(const float& lo, const float& hi);
	42	void SetCharge(const int& ch);
	43	void SetPidProbElectron(const float& lo, const float& hi);
	44	void SetPidProbPion(const float& lo, const float& hi);
	45	void SetPidProbKaon(const float& lo, const float& hi);
	46	void SetPidProbProton(const float& lo, const float& hi);
	47	void SetPidProbMuon(const float& lo, const float& hi);
	48	void SetLabel(const bool& flag);
	49	void SetStatus(const long& w);
	50	void SetminTPCclsF(const short& s);
	51	void SetminTPCncls(const short& s);
	52	void SetminITScls(const int& s);
	53	void SetRemoveKinks(const bool& flag);
	54	void SetRemoveITSFake(const bool& flag);
	55	void SetMaxITSChiNdof(const float& maxchi);
	56	void SetMaxTPCChiNdof(const float& maxchi);
	57	void SetMaxSigmaToVertex(const float& maxsig);
	58	void SetMaxImpactXY(const float& maximpxy);
	59	void SetMaxImpactZ(const float& maximpz);
	60	void SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow);
	61	void SetMostProbablePion();
	62	void SetMostProbableKaon();
	63	void SetMostProbableProton();
	64	void SetNoMostProbable();
65	void SetPIDMethod(ReadPIDMethodType newMethod);
66	void SetClusterRequirementITS(AliESDtrackCuts::Detector det, AliESDtrackCuts::ITSClusterRequirement req = AliESDtrackCuts::kOff);
67
68	void SetMomRangeTOFpidIs(const float& minp, const float& maxp);
69	void SetMomRangeTPCpidIs(const float& minp, const float& maxp);
70	void SetMomRangeITSpidIs(const float& minp, const float& maxp);
71
8a345524	72	//ml
	73	void SetNsigmaTPCle250(Double_t);
	74	void SetNsigmaTPC250_400(Double_t);
	75	void SetNsigmaTPC400_450(Double_t);
	76	void SetNsigmaTPC450_500(Double_t);
	77	void SetNsigmaTPCge500(Double_t);
	78
	79	void SetNsigmaTOF500_800(Double_t);
	80	void SetNsigmaTOF800_1000(Double_t);
	81	void SetNsigmaTOFge1000(Double_t);
	82
	83
	84
9bba4e99	85	private: // here are the quantities I want to cut on...
	86
	87	int fCharge; // particle charge
	88	float fPt[2]; // bounds for transverse momentum
	89	float fRapidity[2]; // bounds for rapidity
	90	float fEta[2]; // bounds for pseudorapidity
	91	float fPidProbElectron[2]; // bounds for electron probability
	92	float fPidProbPion[2]; // bounds for pion probability
	93	float fPidProbKaon[2]; // bounds for kaon probability
	94	float fPidProbProton[2]; // bounds for proton probability
	95	float fPidProbMuon[2]; // bounds for muon probability
	96
	97	AliESDtrackCuts::ITSClusterRequirement fCutClusterRequirementITS[3]; // detailed ITS cluster requirements for (SPD, SDD, SSD) - from AliESDtrackcuts!
	98	bool fLabel; // if true label<0 will not pass throught
	99	long fStatus; // staus flag
	100	ReadPIDMethodType fPIDMethod; // which PID mehod to use. 0 - nsgima, 1 - contour
	101
8a345524	102	//ml
	103	Double_t fNsigmaTPCle250;
	104	Double_t fNsigmaTPC250_400;
	105	Double_t fNsigmaTPC400_450;
	106	Double_t fNsigmaTPC450_500;
	107	Double_t fNsigmaTPCge500;
	108
	109	Double_t fNsigmaTOF500_800;
	110	Double_t fNsigmaTOF800_1000;
	111	Double_t fNsigmaTOFge1000;
	112
	113
	114
9bba4e99	115	short fminTPCclsF; // min number of findable clusters in the TPC
	116	short fminTPCncls; // min number of clusters in the TPC
	117	int fminITScls; // min number of clusters assigned in the ITS
	118	float fMaxITSchiNdof; // maximum allowed chi2/ndof for ITS clusters
	119	float fMaxTPCchiNdof; // maximum allowed chi2/ndof for TPC clusters
	120	float fMaxSigmaToVertex; // maximum allowed sigma to primary vertex
	121	long fNTracksPassed; // passed tracks count
	122	long fNTracksFailed; // failed tracks count
	123	bool fRemoveKinks; // if true particles with any kink label will not pass
	124	bool fRemoveITSFake; // if true particles with ITS fake flag will not pass
	125	int fMostProbable; // this particle type is required to be most probable
	126
	127	float fMaxImpactXY; // Max XY impact parameter
	128	float fMaxImpactZ; // Max Z impact parameter
	129
	130	float fMaxImpactXYPtOff; // Max XY DCA Pt dependent offset
	131	float fMaxImpactXYPtNrm; // Max XY DCA Pt dependent normalization
	132	float fMaxImpactXYPtPow; // Max XY DCA Pt dependent power
	133
	134	float fMinPforTOFpid; // momentum from which TOF PID is requested
	135	float fMaxPforTOFpid; // momentum till which TOF PID is requested
	136	float fMinPforTPCpid; // momentum from which TPC PID is requested
	137	float fMaxPforTPCpid; // momentum till which TPC PID is requested
	138	float fMinPforITSpid; // momentum from which ITS PID is requested
	139	float fMaxPforITSpid; // momentum till which ITS PID is requested
	140
	141	float PidFractionElectron(float mom) const;
	142	float PidFractionPion(float mom) const;
	143	float PidFractionKaon(float mom) const;
	144	float PidFractionProton(float mom) const;
	145
	146	bool IsPionTPCdEdx(float mom, float dEdx);
	147	bool IsKaonTPCdEdx(float mom, float dEdx);
	148	bool IsProtonTPCdEdx(float mom, float dEdx);
	149
	150	bool IsPionTOFTime(float mom, float ttof);
	151	bool IsKaonTOFTime(float mom, float ttof);
	152	bool IsProtonTOFTime(float mom, float ttof);
	153
	154	bool IsKaonTPCdEdxNSigma(float mom, float nsigma);
	155	bool IsKaonTOFNSigma(float mom, float nsigma);
	156	bool IsKaonNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
	157	bool IsPionNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
	158	bool IsProtonNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
	159
	160	Bool_t CheckITSClusterRequirement(AliESDtrackCuts::ITSClusterRequirement req, Bool_t clusterL1, Bool_t clusterL2); //the same as in AliESDtrackCuts
	161
	162
	163	#ifdef __ROOT__
	164	ClassDef(AliFemtoKKTrackCut, 1)
	165	#endif
	166	};
	167
	168
	169	inline void AliFemtoKKTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
	170	inline void AliFemtoKKTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
	171	inline void AliFemtoKKTrackCut::SetEta(const float& lo,const float& hi){fEta[0]=lo; fEta[1]=hi;}
	172	inline void AliFemtoKKTrackCut::SetCharge(const int& ch){fCharge = ch;}
	173	inline void AliFemtoKKTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
	174	inline void AliFemtoKKTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
	175	inline void AliFemtoKKTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
	176	inline void AliFemtoKKTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
	177	inline void AliFemtoKKTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
	178	inline void AliFemtoKKTrackCut::SetLabel(const bool& flag){fLabel=flag;}
179	inline void AliFemtoKKTrackCut::SetStatus(const long& status){fStatus=status;}
180	inline void AliFemtoKKTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
181	inline void AliFemtoKKTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
182	inline void AliFemtoKKTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
183	inline void AliFemtoKKTrackCut::SetMostProbablePion() { fMostProbable = 2; }
184	inline void AliFemtoKKTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
185	inline void AliFemtoKKTrackCut::SetMostProbableProton() { fMostProbable = 4; }
186	inline void AliFemtoKKTrackCut::SetNoMostProbable() { fMostProbable = 0; }
187	inline void AliFemtoKKTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
188	inline void AliFemtoKKTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
189	inline void AliFemtoKKTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
190	inline void AliFemtoKKTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; }
191	inline void AliFemtoKKTrackCut::SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow) { fMaxImpactXYPtOff = maxoff; fMaxImpactXYPtNrm = maxnrm; fMaxImpactXYPtPow = maxpow; }
192	inline void AliFemtoKKTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; }
193
194
8a345524	195
8a345524	196
9bba4e99	197	#endif
9bba4e99	198