// -------------------------------------------------------------------------
// -----                      CbmStsPoint header file                  -----
// -----                  Created 26/07/04  by V. Friese               -----
// -------------------------------------------------------------------------


/**  CbmStsPoint.h
 *@author V.Friese <v.friese@gsi.de>
 *
 * Interception of MC track with a STT detetcor. Holds in addition
 * to the base class the coordinates and momentum at the exit from
 * the active volume.
 **/


#ifndef CBMSTTPOINT_H
#define CBMSTTPOINT_H


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

using namespace std;



class CbmSttPoint : public FairMCPoint
{

 public:

  /** Default constructor **/
  CbmSttPoint();


  /** Constructor with arguments
   *@param trackID       Index of MCTrack
   *@param detID         Detector ID
   *@param pos           Coordinates at wire center of active volume [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]
   **/
  CbmSttPoint(Int_t trackID, Int_t detID, TVector3 pos, 
	      TVector3 posInLocal, TVector3 posOutLocal, 
	      TVector3 momIn, TVector3 momOut, TVector3 wireDir,
	      Double_t tof, Double_t length, Double_t eLoss, Double_t mass, TVector3 postot); // da cancellare postot


  /** Copy constructor **/
  CbmSttPoint(const CbmSttPoint& point) { *this = point; };


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

  // da cancellare
  Double_t GetXtot()  const { return fxtot; }
  Double_t GetYtot()  const { return fytot; }
  Double_t GetZtot()  const { return fztot; }
  void SetPositiontot(TVector3 postot);
  //////

  /** 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; }

  Double_t GetPxOut() const { return fPx_out; }
  Double_t GetPyOut() const { return fPy_out; }
  Double_t GetPzOut() const { return fPz_out; }

  Double_t GetXWireDirection() const { return fX_wire_dir; }
  Double_t GetYWireDirection() const { return fY_wire_dir; }
  Double_t GetZWireDirection() const { return fZ_wire_dir; }

  Double_t GetMass() const {return fMass;}

  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); }
  void WireDirection(TVector3& wire) { wire.SetXYZ(fX_wire_dir, fY_wire_dir, fZ_wire_dir); }

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

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



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

  Double32_t fPx_out, fPy_out, fPz_out;
  // wire direction
  Double32_t fX_wire_dir, fY_wire_dir, fZ_wire_dir;

  // stt1 - particle mass
  Double_t fMass;

  // da cancellare
  Double_t fxtot, fytot, fztot; 
  //////

  ClassDef(CbmSttPoint,1)

};



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

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


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

inline void CbmSttPoint::SetWireDirection(TVector3 wire) {
    fX_wire_dir = wire.Px();
    fY_wire_dir = wire.Py();
    fZ_wire_dir = wire.Pz();
}

// da cancellare
inline void CbmSttPoint::SetPositiontot(TVector3 postot){
  fxtot = postot.X();
  fytot = postot.Y();
  fztot = postot.Z();
}
////////
#endif