//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      a GEANE (sd-system)  track representation
//      (q/p, v',w',v,w) 
//      (v,w) refers to DetPlane system
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Sebastian Neubert    TUM            (original author)
//      Andrea Fontana       INFN
//
//
//-----------------------------------------------------------

#ifndef GeaneTRACKREP_HH
#define GeaneTRACKREP_HH

// Base Class Headers ----------------
#include "AbsTrackRep.h"
#include "CbmTrackParP.h"

// Collaborating Class Headers -------
#include <ostream> // remove if you do not need streaming op
#include "TVectorT.h"


// Collaborating Class Declarations --
class CbmGeanePro;


class GeaneTrackRep : public AbsTrackRep {
public:

  // Constructors/Destructors ---------
  GeaneTrackRep();
  GeaneTrackRep(CbmGeanePro* geane,const TVector3& pos,
			     const TVector3& mom,
			     const TVector3& poserr,
			     const TVector3& momerr,
			     double q,
			     DetPlane& plane);

  //GeaneTrackRep(const GeaneTrackRep& rep); // copy ctor

  virtual ~GeaneTrackRep();


  virtual AbsTrackRep* clone(){std::cout<<"cloning GeaneTrackRep"<<std::endl;return new GeaneTrackRep(*this);}

  virtual AbsTrackRep* prototype(){return new GeaneTrackRep();}

  // Operators
  friend std::ostream& operator<< (std::ostream& s, const GeaneTrackRep& me);

  // Accessors -----------------------
  virtual int getDim() const {return 5;}

  // Modifiers
  void setReferencePlane(const DetPlane& pl) {_refPlane=pl;}

  // Operations ----------------------
  virtual void predict(const DetPlane&,TMatrixT<double>& statePred,
		       TMatrixT<double>& covPred,TMatrixT<double>& jacobian);
  virtual void extrapolate(const DetPlane&, TMatrixT<double>& statePred);
  //virtual void extrapolate(const DetPlane&, 
  //			   const TMatrixT<double>& stateFrom, 
  //			   TMatrixT<double>& stateResult);

  virtual void extrapolate(const DetPlane&, 
			   TMatrixT<double>& statePred,
			   TMatrixT<double>& covPred,
			   TMatrixT<double>& jacobian);
  virtual void stepalong(double h){}

 virtual DetPlane getVirtualDetPlane(AbsRecoHit* hit) const {return DetPlane();}



  virtual TVector3 getPos(const DetPlane&) ;
  virtual TVector3 getMom(const DetPlane&) ;
  virtual double getCharge()const {return state[0][0] > 0 ? 1.: -1.;}
  
private:

  // Private Data Members ------------
  CbmGeanePro* _geane;
  // Private Methods -----------------
  
  // calculate jacobian of extrapolation
  //void Jacobian(const DetPlane& pl,
  //	const TMatrixT<double>& statePred,
  //	TMatrixT<double>& jacResult);

 public:
  ClassDef(GeaneTrackRep,1)

};


#endif

//--------------------------------------------------------------
// $Log$
//--------------------------------------------------------------