#include "AliESDtrackCuts.h"
#ifndef ALIFEMTOKKTRACKCUT_H
#define ALIFEMTOKKTRACKCUT_H
#include "AliESDtrackCuts.h"
#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);
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:
int fCharge;
float fPt[2];
float fRapidity[2];
float fEta[2];
float fPidProbElectron[2];
float fPidProbPion[2];
float fPidProbKaon[2];
float fPidProbProton[2];
float fPidProbMuon[2];
AliESDtrackCuts::ITSClusterRequirement fCutClusterRequirementITS[3];
bool fLabel;
long fStatus;
ReadPIDMethodType fPIDMethod;
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;
short fminTPCncls;
int fminITScls;
float fMaxITSchiNdof;
float fMaxTPCchiNdof;
float fMaxSigmaToVertex;
long fNTracksPassed;
long fNTracksFailed;
bool fRemoveKinks;
bool fRemoveITSFake;
int fMostProbable;
float fMaxImpactXY;
float fMaxImpactZ;
float fMaxImpactXYPtOff;
float fMaxImpactXYPtNrm;
float fMaxImpactXYPtPow;
float fMinPforTOFpid;
float fMaxPforTOFpid;
float fMinPforTPCpid;
float fMaxPforTPCpid;
float fMinPforITSpid;
float fMaxPforITSpid;
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);
#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