#ifndef CBMSTTRECOTRK_HH
#define CBMSTTRECOTRK_HH

#include "CbmTrackParam.h"
#include "TVector3.h"
#include "TVector2.h"

class CbmSttRecoTrk  : public CbmTrackParam
{
 public:
  
  CbmSttRecoTrk();
  
  CbmSttRecoTrk(Int_t trackID,Int_t eventID, TVector3 xy_par,TVector3 ptot);
  CbmSttRecoTrk(Int_t trackID,Int_t eventID, TVector3 xy_par,TVector2 z_par, Int_t fs, Double_t ptot);
  
  virtual ~CbmSttRecoTrk();
  
  /** Output to screen (not yet implemented) **/
  virtual void Print(const Option_t* opt = 0) const {;}

  /** Public method Clear
   ** Resets the flag to -1
   **/
  void Clear();

  /** Accessors **/
  Int_t GetFlag() const { return fFlag; }; 


  /** Modifiers **/
  void SetFlag (Int_t flag) { fFlag = flag; };
  // void SetFitStatus(Int_t fs);
  // Double_t GetMomentum(Int_t trackID,Int_t eventID);
  TVector3 GetxyR(){return fR0_phi0_R;}
  Double_t GetPt() {return fpt;}
  Int_t GetEvent() {return fEventID;}
  Int_t GetTrack() {return fTrackID;}
  Double_t GetMomentum() {return fptot;}		

 protected: 
  Int_t fFlag;
  

 private:
 
  Int_t fTrackID;
  Int_t fEventID;
  Int_t fitstatus;     // +1 every time the fit is successful
  TVector3 fR0_phi0_R; // (R_0, phi_0) = helix center in cilindrical coordinates; R = curvature radius
  TVector2 fdip_z0;    // fdip_z0.X() = tan(lambda), fdip_z0.Y() = z0
  Double_t fpt;        // transverse momentum Gev/c
  Double_t fptot;      // total reconstructed momentum GeV/c
  ClassDef(CbmSttRecoTrk, 1);
};

#endif