#include "AbsTrackRep.h"
#include <iostream>


AbsTrackRep::AbsTrackRep(int dim) : state(dim,1), cov(dim,dim), startState(dim,1), startCov(dim,dim) , chiSqu(0), inverted(false){
  statusFlag = 0;
}

AbsTrackRep::~AbsTrackRep() {}

void 
AbsTrackRep::init(const TVector3& startpoint, 
			  const TVector3& startmomentum, 
			  const TVector3& locBField)
{throw;}

DetPlane 
AbsTrackRep::getVirtualDetPlane(AbsRecoHit* hit) const
{throw;}


// this changes the state of the rep!!! ->
void
AbsTrackRep::extrapolate(const DetPlane& plane){
  TMatrixT<double> newstate=state;
  TMatrixT<double> newcov=cov;
  TMatrixT<double> jac;
  extrapolate(plane,newstate);//,newcov,jac);
  state=newstate;
  cov=newcov;
  setReferencePlane(plane);
}

void
AbsTrackRep::reset(){
  std::cout<<"AbsTrackRep::reset"<<std::endl;
  state.Zero();
  cov.Zero();
  startState.Zero();
  startCov.Zero();
}

void
AbsTrackRep::Print() const {
  std::cout << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << std::endl;
  std::cout<<"AbsTrackRep::Parameters at reference plane ";
  _refPlane.Print();
  std::cout<<"AbsTrackRep::State"<<std::endl;
  state.Print();
  std::cout<<"AbsTrackRep::Covariances"<<std::endl;
  cov.Print();
  std::cout<<"AbsTrackRep::chi^2"<<std::endl;
  std::cout<<chiSqu<<std::endl;
  std::cout << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << std::endl;
}

ClassImp(AbsTrackRep)