Reconstruction/MdcPatRec/TrkBase/include/TrkBase/TrkCompTrk.h

Go to the documentation of this file.

#ifndef TRKCOMPTRK_HH
#define TRKCOMPTRK_HH
 
#include "GaudiKernel/DataObject.h"
#include "TrkBase/TrkAbsFit.h"
#include <memory>
 
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Vector.h"
#include <iosfwd>
 
using CLHEP::Hep3Vector;
using CLHEP::HepMatrix;
using CLHEP::HepSymMatrix;
class IfdKey;
class TrkExchangePar;
class TrkVolume;
class TrkDifTraj;
class TrkSimpTraj;
namespace MdcPatRec {
  class BField;
}
 
class TrkCompTrk : public TrkAbsFit, public DataObject {
public:
  //  friend class TrkCompMaker;
 
  //****************
  // Global track quantities:
  //****************
  int                      charge() const;
  double                   chisq() const;
  int                      nDof() const;
  const TrkDifTraj&        traj() const;
  const MdcPatRec::BField& bField() const { return *_bf; }
 
  //****************
  // Information about track at a given position (flight length)
  //****************
 
  HepPoint3D   position( double fltL ) const;
  Hep3Vector   direction( double fltL ) const;
  Hep3Vector   momentum( double fltL = 0. ) const;
  double       pt( double fltL = 0. ) const;
  BesPointErr  positionErr( double fltL ) const;
  BesVectorErr directionErr( double fltL ) const;
  BesVectorErr momentumErr( double fltL ) const;
 
  // maybe this isn't needed?
  //  TrkExchangePar     helix(double fltL)        const;
 
  //****************
  // Valid flight range
  //****************
  //  validRange = trajectory's valid range = flight length range
  //  over which the track can be used.
  double startValidRange() const;
  double endValidRange() 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;
 
  //****************
  // Printing
  //****************
  virtual void print( std::ostream& ) const;
  virtual void printAll( std::ostream& ) const;
 
  //****************
  // Constructors and such
  //****************
  // Constructor from parameters
  TrkCompTrk( const BesPointErr& pos, const BesVectorErr& mom, const HepMatrix& xpCov,
              int charge, double chisq, int nDoF, const MdcPatRec::BField* bf );
  // Copy constructor (leaves original unchanged):
  TrkCompTrk( const TrkCompTrk& right );
  // Destructor
  virtual ~TrkCompTrk();
  const TrkCompTrk& operator=( const TrkCompTrk& right );
 
private:
  std::auto_ptr<TrkSimpTraj> _traj;
  const MdcPatRec::BField*   _bf;
  double                     _chisq;
  int                        _charge;
  int                        _nDof;
};
 
#endif