SlTrackRep.h

Go to the documentation of this file.

#ifndef SLTRACKREP_HH
#define SLTRACKREP_HH
#include "larreco/Genfit/GFAbsTrackRep.h"

namespace genf {

class SlTrackRep : public GFAbsTrackRep {
public:

  // Constructors/Destructors ---------
  SlTrackRep();
  SlTrackRep(const TMatrixT<double>&, const TMatrixT<double>&);
  SlTrackRep(const TMatrixT<double>&, const TMatrixT<double>&,const double);
  SlTrackRep(const GFDetPlane&,const TMatrixT<double>&, const TMatrixT<double>&);
  SlTrackRep(const TVector3& pos, const TVector3& dir);


  virtual ~SlTrackRep();


  virtual GFAbsTrackRep* clone() const {return new SlTrackRep(*this);}
  virtual GFAbsTrackRep* prototype()const{return new SlTrackRep();}

  void setReferencePlane(const GFDetPlane& pl) {fRefPlane=pl;}

  // Operations ----------------------

  virtual double extrapolate(const GFDetPlane&,
                           TMatrixT<double>& statePred,
                           TMatrixT<double>& covPred);
  virtual double extrapolate(const GFDetPlane&,
                           TMatrixT<double>& statePred);


  void extrapolateToPoint(const TVector3& pos,
                         TVector3& poca,
                         TVector3& dirInPoca);

  void extrapolateToLine(const TVector3& point1,
                         const TVector3& point2,
                         TVector3& poca,
                         TVector3& dirInPoca,
                         TVector3& poca_onwire);


  virtual TVector3 getPos(const GFDetPlane&) ;
  virtual TVector3 getMom(const GFDetPlane&);

  virtual void getPosMom(const GFDetPlane&,TVector3& pos,TVector3& mom) ;
  virtual double getCharge()const {return 0;}

  void switchDirection(){_backw=-_backw;}

private:

  // Private Data Members ------------
  int _backw; // (-1,0,1) -> (backward prop,decide myself,forward)

  // public:
  //ClassDef(SlTrackRep,1);

};
} // namespace genf

#endif