Reconstruction/MdcPatRec/TrkFitter/include/TrkFitter/TrkSimpleRep.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: TrkSimpleRep.h,v 1.2 2006/03/28 01:03:35 zhangy Exp $
//
// Description:
//      Base class for simple TrkReps: helix, line, circle.
//
// Environment:
//      Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//      Steve Schaffner
//
//------------------------------------------------------------------------
 
#ifndef TRKSIMPLEREP_HH
#define TRKSIMPLEREP_HH
 
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "TrkBase/TrkRep.h"
#include "TrkFitter/TrkHelixFitter.h"
 
class TrkSimpTraj;
using CLHEP::Hep3Vector;
using CLHEP::HepMatrix;
using CLHEP::HepSymMatrix;
using CLHEP::HepVector;
 
class TrkSimpleRep : public TrkRep {
public:
  //******************************************
  // Constructors and such
  //******************************************
  // 2 ctors, one for reps with hits, one for reps without hits
  TrkSimpleRep( const TrkHotList*, TrkRecoTrk*, PdtPid::PidType );
  TrkSimpleRep( TrkRecoTrk*, PdtPid::PidType, int nActive, int nSvt, int nDch, double chi2,
                double startFoundRange, double endFoundRange );
  // Copy ctor:
  TrkSimpleRep( const TrkSimpleRep& right, TrkRecoTrk* trk, PdtPid::PidType );
  virtual ~TrkSimpleRep();
 
  //******************************************
  // Global quantities:
  //******************************************
  int         nDof() const;
  double      chisq() const;
  void        setChisq( double c );
  virtual int charge() const;
 
  //********************************************
  // Information about track at a given position
  //********************************************
  virtual Hep3Vector   momentum( double fltL = 0. ) const;
  virtual double       pt( double fltL = 0. ) const;
  virtual BesVectorErr momentumErr( double fltL ) const;
 
  // Interface to vertexing algorithms (M.Bondioli 7/17/98)
  // covariance matrices of the track at fixed flight length
  virtual HepMatrix posmomCov( double fltL ) const;
  virtual void      getAllCovs( double fltL, HepSymMatrix& xxCov, HepSymMatrix& ppCov,
                                HepMatrix& xpCov ) const;
  // accessors to 2nd derivative of chi2 wrt x and p.
  // x0 and p0 are filled with the pos and 3mom around which expansion
  // takes place, whilst Weights are filled with the 2nd deriv of chi2.
  // so that:
  //  dx=x-x0   dp=p-p0
  //
  //  chi2(x,p)=0.5 dx^t*Wxx*dx + dx^t*Wxp*dp + 0.5 dp^t*Wpp*dp
  // where:
  //  pos and mom are 3-dim Vectors,
  //  xxWeight ppWeight and xpWeight are 3 by 3 matrices
  //
  virtual void getAllWeights( double fltL, HepVector& pos, HepVector& mom,
                              HepSymMatrix& xxWeight, HepSymMatrix& ppWeight,
                              HepMatrix& xpWeight ) const;
  virtual void getAllWeights( const HepPoint3D& pt, HepVector& pos, HepVector& mom,
                              HepSymMatrix& xxWeight, HepSymMatrix& ppWeight,
                              HepMatrix& xpWeight ) const;
 
  //******************************************
  // Fitting
  //******************************************
  virtual TrkErrCode fit();
  TrkHelixFitter&    fitter() { return _fitter; }
 
protected:
  TrkSimpleRep&              operator=( const TrkSimpleRep& right );
  virtual TrkSimpTraj&       simpTraj()       = 0;
  virtual const TrkSimpTraj& simpTraj() const = 0;
 
private:
  TrkHelixFitter _fitter;
  double         _chisq;
};
#endif