GFAbsTrackRep.cxx

Go to the documentation of this file.

/* Copyright 2008-2010, Technische Universitaet Muenchen,
   Authors: Christian Hoeppner & Sebastian Neubert

   This file is part of GENFIT.

   GENFIT is free software: you can redistribute it and/or modify
   it under the terms of the GNU Lesser General Public License as published
   by the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.

   GENFIT is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with GENFIT.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "larreco/Genfit/GFAbsTrackRep.h"
#include "larreco/Genfit/GFException.h"
#include <iostream>

genf::GFAbsTrackRep::GFAbsTrackRep() : fDimension(5),fState(5,1), fCov(5,5), fChiSqu(0), fNdf(0), fStatusFlag(0), fInverted(false), fFirstState(5,1), fFirstCov(5,5), fLastState(5,1), fLastCov(5,5)
{
}

genf::GFAbsTrackRep::GFAbsTrackRep(int dim) : fDimension(dim), fState(dim,1), fCov(dim,dim), fChiSqu(0), fNdf(0), fStatusFlag(0), fInverted(false), fFirstState(dim,1), fFirstCov(dim,dim), fLastState(dim,1), fLastCov(dim,dim)
{
}

genf::GFAbsTrackRep::~GFAbsTrackRep() {}

double genf::GFAbsTrackRep::extrapolate(const GFDetPlane& plane){
  TMatrixT<Double_t> statePred(fDimension,1);
  TMatrixT<Double_t> covPred(fDimension,fDimension);
  double retVal = extrapolate(plane,statePred,covPred);
  setData(statePred,plane,&covPred);
  return retVal;
}

//default implentation might be overwritten, please see the doxy docu
double genf::GFAbsTrackRep::extrapolate(const GFDetPlane& plane, TMatrixT<Double_t>& statePred){
  TMatrixT<Double_t> cov(fDimension,fDimension);
  return extrapolate(plane,statePred,cov);
}

void genf::GFAbsTrackRep::Abort(std::string method){
  std::cerr << method <<  " as implemented in " << __FILE__ 
            << " was called. This means that this feature was used "
            << "in a track rep which didnt overwrite this method. "
            << std::endl << "C++ throw;" << std::endl;
  //system call abort
  throw GFException("genf::GFAbsTrackRep: " + method + "() not implemented", __LINE__, __FILE__).setFatal();
}

void genf::GFAbsTrackRep::extrapolateToPoint(const TVector3& /* point */,
                                    TVector3& /* poca */,
                                    TVector3& /* normVec */){
  Abort("extrapolateToPoca()");
}

void genf::GFAbsTrackRep::extrapolateToLine(const TVector3& /* point1 */, 
                                                                        const TVector3& /* point2 */,
                                                                        TVector3& /* poca */,
                                                                        TVector3& /* normVec */,
                                                                        TVector3& /* poca_onwire */){
  Abort("extrapolateToLine()");
}
  

void genf::GFAbsTrackRep::stepalong(double /* h */){
  Abort("stepalong()");
}

void genf::GFAbsTrackRep::getPosMomCov(const GFDetPlane& /* pl */,TVector3& /* pos */,TVector3& /* mom */,TMatrixT<Double_t>& /* cov */){
  Abort("getPosMomCov()");
}

void
genf::GFAbsTrackRep::reset(){
  std::cout<<"GFAbsTrackRep::reset"<<std::endl;
  TVector3 nullVec(0.,0.,0.);
  fRefPlane.set(nullVec,nullVec,nullVec);
  fState.Zero();
  fCov.Zero();
  fFirstState.Zero();
  fFirstCov.Zero();
  fLastState.Zero();
  fLastCov.Zero();
}

void
genf::GFAbsTrackRep::Print(std::ostream& out /* = std::cout */) const {
  out << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << std::endl;
  out <<"GFAbsTrackRep::Parameters at reference plane ";
  fRefPlane.Print(out);
  out <<"GFAbsTrackRep::State"<<std::endl;
  PrintROOTmatrix(out, fState);
  out <<"GFAbsTrackRep::Covariances"<<std::endl;
  PrintROOTmatrix(out, fCov);
  out <<"GFAbsTrackRep::chi^2"<<std::endl;
  out <<fChiSqu<<std::endl;
  out << "++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++" << std::endl;
}



//ClassImp(GFAbsTrackRep)