[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemtoUser / AliFemtoESDTrackCut.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               //
//                                                                       //
///////////////////////////////////////////////////////////////////////////

#ifndef ALIFEMTOESDTRACKCUT_H
#define ALIFEMTOESDTRACKCUT_H

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

#include "AliFemtoTrackCut.h"

class AliFemtoESDTrackCut : public AliFemtoTrackCut 
{

 public:
  AliFemtoESDTrackCut();
  virtual ~AliFemtoESDTrackCut();

  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 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 SetMostProbablePion();
  void SetMostProbableKaon();
  void SetMostProbableProton();
  void SetNoMostProbable(); 

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

 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 
  bool              fLabel;              // if true label<0 will not pass throught 
  long              fStatus;             // staus flag

  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
  int               fMostProbable;       // this particle type is required to be most probable

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

  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);

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


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

#endif
Commit	Line	Data
d92ed900	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	///////////////////////////////////////////////////////////////////////////
67427ff7	10
d92ed900	11	#ifndef ALIFEMTOESDTRACKCUT_H
d92ed900	12	#define ALIFEMTOESDTRACKCUT_H
67427ff7	13
	14	//#ifndef StMaker_H
	15	//#include "StMaker.h"
	16	//#endif
	17
65423af9	18	#include "AliFemtoTrackCut.h"
67427ff7	19
	20	class AliFemtoESDTrackCut : public AliFemtoTrackCut
	21	{
	22
d92ed900	23	public:
d92ed900	24	AliFemtoESDTrackCut();
cc5faabc	25	virtual ~AliFemtoESDTrackCut();
67427ff7	26
d92ed900	27	virtual bool Pass(const AliFemtoTrack* aTrack);
67427ff7	28
d92ed900	29	virtual AliFemtoString Report();
3a74204a	30	virtual TList *ListSettings();
0b3bd1ac	31	virtual AliFemtoParticleType Type(){return hbtTrack;}
67427ff7	32
d92ed900	33	void SetPt(const float& lo, const float& hi);
d92ed900	34	void SetRapidity(const float& lo, const float& hi);
c5c1b5e8	35	void SetEta(const float& lo, const float& hi);
cc5faabc	36	void SetCharge(const int& ch);
d92ed900	37	void SetPidProbElectron(const float& lo, const float& hi);
	38	void SetPidProbPion(const float& lo, const float& hi);
	39	void SetPidProbKaon(const float& lo, const float& hi);
	40	void SetPidProbProton(const float& lo, const float& hi);
	41	void SetPidProbMuon(const float& lo, const float& hi);
	42	void SetLabel(const bool& flag);
	43	void SetStatus(const long& w);
	44	void SetminTPCclsF(const short& s);
0b3bd1ac	45	void SetminTPCncls(const short& s);
d92ed900	46	void SetminITScls(const int& s);
cc5faabc	47	void SetRemoveKinks(const bool& flag);
0b3bd1ac	48	void SetMaxITSChiNdof(const float& maxchi);
	49	void SetMaxTPCChiNdof(const float& maxchi);
	50	void SetMaxSigmaToVertex(const float& maxsig);
ab74d4d3	51	void SetMaxImpactXY(const float& maximpxy);
ab74d4d3	52	void SetMaxImpactZ(const float& maximpz);
cc5faabc	53	void SetMostProbablePion();
	54	void SetMostProbableKaon();
	55	void SetMostProbableProton();
	56	void SetNoMostProbable();
	57
e3f643f0	58	void SetMomRangeTOFpidIs(const float& minp, const float& maxp);
	59	void SetMomRangeTPCpidIs(const float& minp, const float& maxp);
	60	void SetMomRangeITSpidIs(const float& minp, const float& maxp);
	61
d92ed900	62	private: // here are the quantities I want to cut on...
67427ff7	63
d92ed900	64	int fCharge; // particle charge
	65	float fPt[2]; // bounds for transverse momentum
	66	float fRapidity[2]; // bounds for rapidity
c5c1b5e8	67	float fEta[2]; // bounds for pseudorapidity
d92ed900	68	float fPidProbElectron[2]; // bounds for electron probability
	69	float fPidProbPion[2]; // bounds for pion probability
	70	float fPidProbKaon[2]; // bounds for kaon probability
	71	float fPidProbProton[2]; // bounds for proton probability
	72	float fPidProbMuon[2]; // bounds for muon probability
	73	bool fLabel; // if true label<0 will not pass throught
	74	long fStatus; // staus flag
db2d5b4b	75
d92ed900	76	short fminTPCclsF; // min number of findable clusters in the TPC
0b3bd1ac	77	short fminTPCncls; // min number of clusters in the TPC
d92ed900	78	int fminITScls; // min number of clusters assigned in the ITS
0b3bd1ac	79	float fMaxITSchiNdof; // maximum allowed chi2/ndof for ITS clusters
	80	float fMaxTPCchiNdof; // maximum allowed chi2/ndof for TPC clusters
	81	float fMaxSigmaToVertex; // maximum allowed sigma to primary vertex
d92ed900	82	long fNTracksPassed; // passed tracks count
d92ed900	83	long fNTracksFailed; // failed tracks count
cc5faabc	84	bool fRemoveKinks; // if true particles with any kink label will not pass
	85	int fMostProbable; // this particle type is required to be most probable
	86
ab74d4d3	87	float fMaxImpactXY; // Max XY impact parameter
	88	float fMaxImpactZ; // Max Z impact parameter
	89
e3f643f0	90	float fMinPforTOFpid; // momentum from which TOF PID is requested
	91	float fMaxPforTOFpid; // momentum till which TOF PID is requested
	92	float fMinPforTPCpid; // momentum from which TPC PID is requested
	93	float fMaxPforTPCpid; // momentum till which TPC PID is requested
	94	float fMinPforITSpid; // momentum from which ITS PID is requested
	95	float fMaxPforITSpid; // momentum till which ITS PID is requested
	96
cc5faabc	97	float PidFractionElectron(float mom) const;
	98	float PidFractionPion(float mom) const;
	99	float PidFractionKaon(float mom) const;
	100	float PidFractionProton(float mom) const;
67427ff7	101
c5c1b5e8	102	bool IsPionTPCdEdx(float mom, float dEdx);
	103	bool IsKaonTPCdEdx(float mom, float dEdx);
	104
67427ff7	105	#ifdef __ROOT__
	106	ClassDef(AliFemtoESDTrackCut, 1)
	107	#endif
d92ed900	108	};
67427ff7	109
	110
	111	inline void AliFemtoESDTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
	112	inline void AliFemtoESDTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
c5c1b5e8	113	inline void AliFemtoESDTrackCut::SetEta(const float& lo,const float& hi){fEta[0]=lo; fEta[1]=hi;}
67427ff7	114	inline void AliFemtoESDTrackCut::SetCharge(const int& ch){fCharge = ch;}
	115	inline void AliFemtoESDTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
	116	inline void AliFemtoESDTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
	117	inline void AliFemtoESDTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
	118	inline void AliFemtoESDTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
	119	inline void AliFemtoESDTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
	120	inline void AliFemtoESDTrackCut::SetLabel(const bool& flag){fLabel=flag;}
	121	inline void AliFemtoESDTrackCut::SetStatus(const long& status){fStatus=status;}
	122	inline void AliFemtoESDTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
0b3bd1ac	123	inline void AliFemtoESDTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
67427ff7	124	inline void AliFemtoESDTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
cc5faabc	125	inline void AliFemtoESDTrackCut::SetMostProbablePion() { fMostProbable = 2; }
	126	inline void AliFemtoESDTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
	127	inline void AliFemtoESDTrackCut::SetMostProbableProton() { fMostProbable = 4; }
	128	inline void AliFemtoESDTrackCut::SetNoMostProbable() { fMostProbable = 0; }
0b3bd1ac	129	inline void AliFemtoESDTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
	130	inline void AliFemtoESDTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
	131	inline void AliFemtoESDTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
ab74d4d3	132	inline void AliFemtoESDTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; }
ab74d4d3	133	inline void AliFemtoESDTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; }
67427ff7	134
	135	#endif
	136