////////////////////////////////////////////////////////////////////////////////
//                                                                            //
// AliFemtoBasicTrackCut - the basic cut for tracks.                          //
// Cuts on particle identification, transverse momentum, rapidity, distance   //
// of closest approach to primary vertex and charge                           //
//                                                                            //
////////////////////////////////////////////////////////////////////////////////

#ifndef ALIFEMTOBASICTRACKCUT_H
#define ALIFEMTOBASICTRACKCUT_H

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

#include "AliFemtoTrackCut.h"

class AliFemtoBasicTrackCut : public AliFemtoTrackCut {

public:

  AliFemtoBasicTrackCut();
  //~mikesTrackCut();

  virtual bool Pass(const AliFemtoTrack* tr);

  virtual AliFemtoString Report();


  void SetNSigmaPion(const float& lo, const float& hi);
  void SetNSigmaKaon(const float& lo, const float& hi);
  void SetNSigmaProton(const float& lo, const float& hi);

  void SetNHits(const int& lo, const int& hi);
  void SetPt(const float& lo, const float& hi);
  void SetRapidity(const float& lo, const float& hi);
  void SetDCA(const float& lo, const float& hi);
  void SetCharge(const int& ch);


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

  int               fCharge;             // charge of the track
  float             fNSigmaPion[2];      // bounds for nsigma dEdx from pion band 
  float             fNSigmaKaon[2];      // bounds for nsigma dEdx from kaon band
  float             fNSigmaProton[2];    // bounds for nsigma dEdx from proton band
  int               fNHits[2];           // bounds for number of hits
  float             fPt[2];              // bounds for transverse momentum
  float             fRapidity[2];        // bounds for rapidity
  float             fDCA[2];             // bounds for DCA to primary vertex

  long              fNTracksPassed;      // passed tracks counter
  long              fNTracksFailed;      // falied tracks counter

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


inline void AliFemtoBasicTrackCut::SetNSigmaPion(const float& lo, const float& hi){fNSigmaPion[0]=lo; fNSigmaPion[1]=hi;}
inline void AliFemtoBasicTrackCut::SetNSigmaKaon(const float& lo, const float& hi){fNSigmaKaon[0]=lo; fNSigmaKaon[1]=hi;}
inline void AliFemtoBasicTrackCut::SetNSigmaProton(const float& lo, const float& hi){fNSigmaProton[0]=lo; fNSigmaProton[1]=hi;}

inline void AliFemtoBasicTrackCut::SetNHits(const int& lo, const int& hi){fNHits[0]=lo;fNHits[1]=hi;}
inline void AliFemtoBasicTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
inline void AliFemtoBasicTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
inline void AliFemtoBasicTrackCut::SetDCA(const float& lo,const float& hi){fDCA[0]=lo; fDCA[1]=hi;}
inline void AliFemtoBasicTrackCut::SetCharge(const int& ch){fCharge = ch;}

#endif