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

///////////////////////////////////////////////////////////////////////////
//                                                                       //
// AliFemtoQATrackCut: 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 AliFemtoQATrackCut_H
#define AliFemtoQATrackCut_H

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

#include "AliFemtoTrackCut.h"

class AliFemtoQATrackCut : public AliFemtoTrackCut 
{

 public:
  AliFemtoQATrackCut();
  virtual ~AliFemtoQATrackCut();

  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 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 SetminTPCchiNdof(const float& s);
  void SetRemoveKinks(const bool& flag);
  void SetMaxTPCncls(const short& s);
  void SetMaxITSChiNdof(const float& maxchi);
  void SetMaxTPCChiNdof(const float& maxchi);
  void SetMaxSigmaToVertex(const float& maxsig);
  void SetMostProbablePion();
  void SetMostProbableKaon();
  void SetMostProbableProton();
  void SetNoMostProbable(); 
  
  void SetTPCnclsExclusionSwitch(const bool& excluSwitch);
  void SetTPCnclsExclusion(const short& lo, const short& hi);
  void SetTPCchiNdofExclusionSwitch(const bool& excluSwitch);
  void SetTPCchiNdofExclusion(const float& lo, const float& hi);

 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             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             fminTPCchiNdof;      // min allowed chi2/ndof for TPC clusters
  short             fMaxTPCncls;         // maximum allowed clusters in the TPC
  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
  
  bool    fTPCnclsExclusionSwitch;       // turn on/off TPCncls exclusion zone (true=on)
  short   fTPCnclsExclusion[2];          // lower and upper limit of TPCncls QA exclusion zone
  bool    fTPCchiNdofExclusionSwitch;        // turn on/off TPCchi exclusion zone (true=on)
  float   fTPCchiNdofExclusion[2];           // lower and upper limit of TPCchi QA exclusion zone

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

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


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

inline void AliFemtoQATrackCut::SetTPCnclsExclusionSwitch(const bool& excluSwitch) { fTPCnclsExclusionSwitch = excluSwitch; }
inline void AliFemtoQATrackCut::SetTPCnclsExclusion(const short& lo, const short& hi) {fTPCnclsExclusion[0] = lo; fTPCnclsExclusion[1] = hi;}
inline void AliFemtoQATrackCut::SetTPCchiNdofExclusionSwitch(const bool& excluSwitch) { fTPCchiNdofExclusionSwitch = excluSwitch; }
inline void AliFemtoQATrackCut::SetTPCchiNdofExclusion(const float& lo, const float& hi) {fTPCchiNdofExclusion[0] = lo; fTPCchiNdofExclusion[1] = hi;}

#endif
Commit	Line	Data
76ce4b5b	1	///////////////////////////////////////////////////////////////////////////
	2	// //
	3	// AliFemtoQATrackCut: 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
	11	#ifndef AliFemtoQATrackCut_H
	12	#define AliFemtoQATrackCut_H
	13
	14	//#ifndef StMaker_H
	15	//#include "StMaker.h"
	16	//#endif
	17
	18	#include "AliFemtoTrackCut.h"
	19
	20	class AliFemtoQATrackCut : public AliFemtoTrackCut
	21	{
	22
	23	public:
	24	AliFemtoQATrackCut();
	25	virtual ~AliFemtoQATrackCut();
	26
	27	virtual bool Pass(const AliFemtoTrack* aTrack);
	28
	29	virtual AliFemtoString Report();
	30	virtual TList *ListSettings();
	31	virtual AliFemtoParticleType Type(){return hbtTrack;}
	32
	33	void SetPt(const float& lo, const float& hi);
	34	void SetRapidity(const float& lo, const float& hi);
	35	void SetCharge(const int& ch);
	36	void SetPidProbElectron(const float& lo, const float& hi);
	37	void SetPidProbPion(const float& lo, const float& hi);
	38	void SetPidProbKaon(const float& lo, const float& hi);
	39	void SetPidProbProton(const float& lo, const float& hi);
	40	void SetPidProbMuon(const float& lo, const float& hi);
	41	void SetLabel(const bool& flag);
	42	void SetStatus(const long& w);
	43	void SetminTPCclsF(const short& s);
	44	void SetminTPCncls(const short& s);
	45	void SetminITScls(const int& s);
	46	void SetminTPCchiNdof(const float& s);
	47	void SetRemoveKinks(const bool& flag);
	48	void SetMaxTPCncls(const short& s);
	49	void SetMaxITSChiNdof(const float& maxchi);
	50	void SetMaxTPCChiNdof(const float& maxchi);
	51	void SetMaxSigmaToVertex(const float& maxsig);
	52	void SetMostProbablePion();
	53	void SetMostProbableKaon();
	54	void SetMostProbableProton();
	55	void SetNoMostProbable();
	56
	57	void SetTPCnclsExclusionSwitch(const bool& excluSwitch);
	58	void SetTPCnclsExclusion(const short& lo, const short& hi);
	59	void SetTPCchiNdofExclusionSwitch(const bool& excluSwitch);
	60	void SetTPCchiNdofExclusion(const float& lo, const float& hi);
	61
	62	private: // here are the quantities I want to cut on...
	63
	64	int fCharge; // particle charge
65	float fPt[2]; // bounds for transverse momentum
66	float fRapidity[2]; // bounds for rapidity
67	float fPidProbElectron[2]; // bounds for electron probability
68	float fPidProbPion[2]; // bounds for pion probability
69	float fPidProbKaon[2]; // bounds for kaon probability
70	float fPidProbProton[2]; // bounds for proton probability
71	float fPidProbMuon[2]; // bounds for muon probability
72	bool fLabel; // if true label<0 will not pass throught
73	long fStatus; // staus flag
74
75	short fminTPCclsF; // min number of findable clusters in the TPC
76	short fminTPCncls; // min number of clusters in the TPC
77	int fminITScls; // min number of clusters assigned in the ITS
78	float fminTPCchiNdof; // min allowed chi2/ndof for TPC clusters
79	short fMaxTPCncls; // maximum allowed clusters in the TPC
80	float fMaxITSchiNdof; // maximum allowed chi2/ndof for ITS clusters
81	float fMaxTPCchiNdof; // maximum allowed chi2/ndof for TPC clusters
82	float fMaxSigmaToVertex; // maximum allowed sigma to primary vertex
83	long fNTracksPassed; // passed tracks count
84	long fNTracksFailed; // failed tracks count
85	bool fRemoveKinks; // if true particles with any kink label will not pass
86	int fMostProbable; // this particle type is required to be most probable
87
88	bool fTPCnclsExclusionSwitch; // turn on/off TPCncls exclusion zone (true=on)
89	short fTPCnclsExclusion[2]; // lower and upper limit of TPCncls QA exclusion zone
90	bool fTPCchiNdofExclusionSwitch; // turn on/off TPCchi exclusion zone (true=on)
91	float fTPCchiNdofExclusion[2]; // lower and upper limit of TPCchi QA exclusion zone
92
93	float PidFractionElectron(float mom) const;
94	float PidFractionPion(float mom) const;
95	float PidFractionKaon(float mom) const;
96	float PidFractionProton(float mom) const;
97
98	#ifdef __ROOT__
99	ClassDef(AliFemtoQATrackCut, 1)
100	#endif
101	};
102
103
104	inline void AliFemtoQATrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
105	inline void AliFemtoQATrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
106	inline void AliFemtoQATrackCut::SetCharge(const int& ch){fCharge = ch;}
107	inline void AliFemtoQATrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
108	inline void AliFemtoQATrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
109	inline void AliFemtoQATrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
110	inline void AliFemtoQATrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
111	inline void AliFemtoQATrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
112	inline void AliFemtoQATrackCut::SetLabel(const bool& flag){fLabel=flag;}
113	inline void AliFemtoQATrackCut::SetStatus(const long& status){fStatus=status;}
114	inline void AliFemtoQATrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
115	inline void AliFemtoQATrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
116	inline void AliFemtoQATrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
117	inline void AliFemtoQATrackCut::SetminTPCchiNdof(const float& s){fminTPCchiNdof = s;}
118	inline void AliFemtoQATrackCut::SetMostProbablePion() { fMostProbable = 2; }
119	inline void AliFemtoQATrackCut::SetMostProbableKaon() { fMostProbable = 3; }
120	inline void AliFemtoQATrackCut::SetMostProbableProton() { fMostProbable = 4; }
121	inline void AliFemtoQATrackCut::SetNoMostProbable() { fMostProbable = 0; }
122	inline void AliFemtoQATrackCut::SetMaxTPCncls(const short& s){fMaxTPCncls=s;}
123	inline void AliFemtoQATrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
124	inline void AliFemtoQATrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
125	inline void AliFemtoQATrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
126
127	inline void AliFemtoQATrackCut::SetTPCnclsExclusionSwitch(const bool& excluSwitch) { fTPCnclsExclusionSwitch = excluSwitch; }
128	inline void AliFemtoQATrackCut::SetTPCnclsExclusion(const short& lo, const short& hi) {fTPCnclsExclusion[0] = lo; fTPCnclsExclusion[1] = hi;}
129	inline void AliFemtoQATrackCut::SetTPCchiNdofExclusionSwitch(const bool& excluSwitch) { fTPCchiNdofExclusionSwitch = excluSwitch; }
130	inline void AliFemtoQATrackCut::SetTPCchiNdofExclusion(const float& lo, const float& hi) {fTPCchiNdofExclusion[0] = lo; fTPCchiNdofExclusion[1] = hi;}
131
132	#endif
133