#ifndef PNDSTTPOINT_H
#define PNDSTTPOINT_H


#include "TObject.h"
#include "TVector3.h"
#include "FairMCPoint.h"


class PndSttPoint : public FairMCPoint 
{

 public:

  /** Default constructor **/
  PndSttPoint();


  /** Constructor with arguments
   *@param trackID       Index of MCTrack
   *@param detID         Detector ID
   *@param tubeID        Number of tube
   *@param pos           Coordinates of MC point [cm]
   *@param posInLocal    Local coordinates at entrance to active volume [cm]
   *@param posOutLocal   Local coordinates at exit of active volume [cm]
   *@param momIn         Momentum of track at entrance [GeV]
   *@param momOut        Momentum of track at exit [GeV]
   *@param tof           Time since event start [ns]
   *@param length        Track length since creation [cm]
   *@param eLoss         Energy deposit [GeV]
   **/
  PndSttPoint(Int_t trackID, Int_t detID, TVector3 pos, 
	      TVector3 posInLocal, TVector3 posOutLocal, 
	      TVector3 momIn, TVector3 momOut, 
	      Double_t tof, Double_t length, Double_t eLoss, Double_t mass);


  /** Copy constructor **/
  PndSttPoint(const PndSttPoint& point);


  /** Destructor **/
  virtual ~PndSttPoint();

  /** Accessors **/
  Double_t GetXOutLocal()  const { return fX_out_local; }
  Double_t GetYOutLocal()  const { return fY_out_local; }
  Double_t GetZOutLocal()  const { return fZ_out_local; }
  Double_t GetXInLocal()  const { return fX_in_local; }
  Double_t GetYInLocal()  const { return fY_in_local; }
  Double_t GetZInLocal()  const { return fZ_in_local; }

  // to be deleted  CHECK
  Double_t GetXtot()  const { return GetX(); }
  Double_t GetYtot()  const { return GetY(); }
  Double_t GetZtot()  const { return GetZ(); }
  
  Double_t GetPxOut() const { return fPx_out; } 
  Double_t GetPyOut() const { return fPy_out; } 
  Double_t GetPzOut() const { return fPz_out; } 

  Double_t GetMass() const {return fMass;}

  Double_t GetTrueDistance();
  
  void PositionOutLocal(TVector3& pos) { pos.SetXYZ(fX_out_local,fY_out_local,fZ_out_local); } 
  void PositionInLocal(TVector3& pos)  { pos.SetXYZ(fX_in_local,fY_in_local,fZ_in_local); }
  void MomentumOut(TVector3& mom) { mom.SetXYZ(fPx_out,fPy_out,fPz_out); } 

  /** Modifiers **/
  void SetPositionOutLocal(TVector3 pos);
  void SetPositionInLocal(TVector3 pos); 
  void SetMomentumOut(TVector3 mom);     

  // tube ID // CHECK added
  void SetTubeID(Int_t tubeid) { fTubeID = tubeid; }
  Int_t GetTubeID() { return fTubeID; }

  /** Output to screen **/
  virtual void Print(const Option_t* opt) const;



 protected:
  // exit coordinates in straw frame
  Double_t fX_out_local,  fY_out_local,  fZ_out_local;
  // entry coordinates in straw frame
  Double_t fX_in_local,  fY_in_local,  fZ_in_local;

  Double_t fPx_in; //!
  Double_t fPy_in; //!
  Double_t fPz_in; //!
  Double_t fPx_out;
  Double_t fPy_out;
  Double_t fPz_out;

  // stt1 - particle mass
  Double_t fMass;

  Int_t fTubeID; // CHECK added

  ClassDef(PndSttPoint,1)

    };



inline void PndSttPoint::SetPositionOutLocal(TVector3 pos) {
  fX_out_local = pos.X();
  fY_out_local = pos.Y();
  fZ_out_local = pos.Z();
}

inline void PndSttPoint::SetPositionInLocal(TVector3 pos) {
  fX_in_local = pos.X();
  fY_in_local = pos.Y();
  fZ_in_local = pos.Z();
}


inline void PndSttPoint::SetMomentumOut(TVector3 mom) {
  fPx_out = mom.Px();
  fPy_out = mom.Py();
  fPz_out = mom.Pz();
}

#endif