/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/MdcPatRec/TrkFitter/TrkFitter-00-01-11/src/TrkSimpleRep.cxx

Go to the documentation of this file.

// File and Version Information:
//      $Id: TrkSimpleRep.cxx,v 1.1.1.1 2005/04/21 06:26:56 maqm Exp $
//
// Description:
//      implementation for TrkSimpleHelix
//
// Author(s): Steve Schaffner
//------------------------------------------------------------------------
//#include "BaBar/BaBar.h"
#include "TrkFitter/TrkSimpleRep.h"
#include "TrkBase/TrkHitOnTrk.h"
#include "TrkBase/TrkSimpTraj.h"
#include "TrkBase/TrkParams.h"
#include "TrkBase/TrkPoca.h"
#include "TrkBase/TrkMomCalculator.h"
#include "TrkBase/TrkRecoTrk.h"
#include "BField/BField.h"//yzhang debug
class BField;
#include "TrkBase/TrkErrCode.h"
#include "MdcRecoUtil/BesVectorErr.h"
#include "MdcRecoUtil/BesPointErr.h"
#include "CLHEP/Vector/ThreeVector.h"

TrkSimpleRep::TrkSimpleRep(const TrkHotList* inHots,
                           TrkRecoTrk* trk, PdtPid::PidType hypo) :
  TrkRep(inHots, trk, hypo)
{
  setChisq(-1.);
}

TrkSimpleRep::TrkSimpleRep(TrkRecoTrk* trk, PdtPid::PidType hypo, int nact,
                           int nsv, int ndc, double chi2,
                           double stFndRng, double endFndRng) :
  TrkRep(trk, hypo, nact, nsv, ndc, stFndRng, endFndRng)
{
  setChisq(chi2);
}

// copy ctor
TrkSimpleRep::TrkSimpleRep(const TrkSimpleRep& right,
                           TrkRecoTrk* trk, PdtPid::PidType hypo)  :
  TrkRep(right, trk, hypo), _fitter(right._fitter), _chisq(right._chisq)
{
  _hotList.reset( right.hotList()->clone(this) );
}

TrkSimpleRep::~TrkSimpleRep()
{
}

TrkSimpleRep&
TrkSimpleRep::operator=(const TrkSimpleRep &right)
{
  _chisq = right._chisq;
  TrkRep::operator=(right);
  return *this;
}

TrkErrCode
TrkSimpleRep::fit()
{
  if (!hitCapable()) {
    setValid(false);
    return TrkErrCode(TrkErrCode::fail,100,
                      "Tried to fit without a list of hits.");
  }
  TrkErrCode error = _fitter.fit(*hotList(), simpTraj());
  if (error.failure()) {
    setValid(false);
  }
  else {
    setValid(true);
    setCurrent(true);
  }
  setChisq(_fitter.lastChisq());
  return error;
}

int
TrkSimpleRep::nDof() const
{
  return nActive() - simpTraj().parameters()->nPar();
}

double
TrkSimpleRep::chisq() const
{
  // Note possible inefficiency in calling virtual simpTraj(); implementation 
  //    could be moved to derived classes instead.  
  return _chisq;
  //  return simpTraj().parameters()->chisquared();
}

void
TrkSimpleRep::setChisq(double c)
{
  // See inefficiency comment above.
  //  simpTraj().parameters()->setChisquared(c);
  _chisq = c;
}

Hep3Vector
TrkSimpleRep::momentum(double fltL) const
{
  const BField& theField = parentTrack()->bField();
  return TrkMomCalculator::vecMom(simpTraj(), theField, fltL);
}

double
TrkSimpleRep::pt(double fltL) const
{
  const BField& theField = parentTrack()->bField();
  return TrkMomCalculator::ptMom(simpTraj(), theField, fltL);
}

int
TrkSimpleRep::charge() const
{
  const BField& theField = parentTrack()->bField();
  return TrkMomCalculator::charge(simpTraj(), theField, 0.);
}

BesVectorErr
TrkSimpleRep::momentumErr(double fltL) const
{
  const BField& theField = parentTrack()->bField();
  return TrkMomCalculator::errMom(simpTraj(), theField, fltL);
}

/*
 * added on 7/17/98 to accomplish vertexing interface (M.Bondioli)
 * implementation temporarily needed to assure no compile crashes
 * elsewhere.
 */

HepMatrix
TrkSimpleRep::posmomCov(double fltL) const
{
  const BField& theField = parentTrack()->bField();
  return TrkMomCalculator::posmomCov(simpTraj(), theField, fltL);
}

void
TrkSimpleRep::getAllCovs(double fltL,
                        HepSymMatrix& xxCov,
                        HepSymMatrix& ppCov,
                        HepMatrix&    xpCov) const
{
  const BField& theField = parentTrack()->bField();
  TrkMomCalculator::getAllCovs(simpTraj(), theField, fltL,
                               xxCov,ppCov,xpCov);
}

void TrkSimpleRep::getAllWeights(double fltL,
                                 HepVector& pos,
                                 HepVector& mom,
                                 HepSymMatrix& xxWeight,
                                 HepSymMatrix& ppWeight,
                                 HepMatrix&    xpWeight) const
{
  const BField& theField = parentTrack()->bField();
  TrkMomCalculator::getAllWeights(simpTraj(), theField, fltL,
                                   pos,mom,xxWeight,ppWeight,xpWeight);
}

void TrkSimpleRep::getAllWeights(const HepPoint3D& pt,
                                 HepVector& pos,
                                 HepVector& mom,
                                 HepSymMatrix& xxWeight,
                                 HepSymMatrix& ppWeight,
                                 HepMatrix&    xpWeight)const
{
  double fltL=0;
  TrkPoca poca(traj(),fltL, pt);
  fltL = poca.flt1();
  getAllWeights(fltL,pos,mom,xxWeight,ppWeight,xpWeight);
}

---