Main Page | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

TrkCircleRep Class Reference

#include <TrkCircleRep.h>

Inheritance diagram for TrkCircleRep:

TrkSimpleRep TrkSimpleRep TrkRep TrkRep TrkRep TrkRep TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus TrkHitOnTrkUpdater TrkFit TrkFitStatus List of all members.

Public Types

typedef std::unary_function<
TrkFitStatus, bool > 		predicate_type
typedef std::unary_function<
TrkFitStatus, bool > 		predicate_type

Public Member Functions

virtual void activateHot (TrkHitOnTrk *theHot)
virtual void activateHot (TrkHitOnTrk *theHot)
template<class T>
void addHistory (std::pair< T, T > p)
template<class T>
void addHistory (T begin, T end)
virtual void addHistory (const TrkErrCode &status, const char *modulename)
template<class T>
void addHistory (std::pair< T, T > p)
template<class T>
void addHistory (T begin, T end)
virtual void addHistory (const TrkErrCode &status, const char *modulename)
virtual void addHot (TrkHitOnTrk *theHot)
virtual void addHot (TrkHitOnTrk *theHot)
virtual double arrivalTime (double fltL) const
virtual double arrivalTime (double fltL) const
history_iterator beginHistory () const
history_iterator beginHistory () const
virtual int charge () const
virtual int charge () const
double chisq () const
double chisq () const
virtual ChisqConsistency chisqConsistency () const
virtual ChisqConsistency chisqConsistency () const
virtual TrkCircleRep * clone (TrkRecoTrk *newTrack) const
virtual TrkCircleRep * clone (TrkRecoTrk *newTrack) const
virtual TrkCircleRep * cloneNewHypo (PdtPid::PidType hypo)
virtual TrkCircleRep * cloneNewHypo (PdtPid::PidType hypo)
virtual void deactivateHot (TrkHitOnTrk *theHot)
virtual void deactivateHot (TrkHitOnTrk *theHot)
virtual Hep3Vector direction (double fltL) const
virtual Hep3Vector direction (double fltL) const
virtual BesVectorErr directionErr (double fltL) const
virtual BesVectorErr directionErr (double fltL) const
virtual double endFoundRange () const
virtual double endFoundRange () const
history_iterator endHistory () const
history_iterator endHistory () const
double endValidRange () const
double endValidRange () const
virtual TrkErrCode fit ()
virtual TrkErrCode fit ()
bool fitCurrent () const
bool fitCurrent () const
const TrkErrCode & fitStatus () const
const TrkErrCode & fitStatus () const
TrkHelixFitter & fitter ()
TrkHelixFitter & fitter ()
bool fitValid () const
bool fitValid () const
virtual void getAllCovs (double fltL, HepSymMatrix &xxCov, HepSymMatrix &ppCov, HepMatrix &xpCov) const
virtual void getAllCovs (double fltL, HepSymMatrix &xxCov, HepSymMatrix &ppCov, HepMatrix &xpCov) const
virtual void getAllWeights (const HepPoint3D &pt, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const
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
virtual void getAllWeights (double fltL, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const
virtual TrkExchangePar helix (double fltLen) const
virtual TrkExchangePar helix (double fltLen) const
std::pair< history_iterator,
history_iterator > 		history () const
std::pair< history_iterator,
history_iterator > 		history () const
const std::vector< TrkHistory > & historyVector () const
const std::vector< TrkHistory > & historyVector () const
bool hitCapable () const
bool hitCapable () const
virtual const TrkHotList * hotList () const
virtual TrkHotList * hotList ()
virtual const TrkHotList * hotList () const
virtual TrkHotList * hotList ()
bool is2d () const
bool is2d () const
virtual Hep3Vector momentum (double fltL=0.) const
virtual Hep3Vector momentum (double fltL=0.) const
virtual BesVectorErr momentumErr (double fltL) const
virtual BesVectorErr momentumErr (double fltL) const
bool multScat () const
bool multScat () const
virtual const IfdKey & myKey () const
virtual const IfdKey & myKey () const
virtual int nActive () const
virtual int nActive () const
int nDof () const
int nDof () const
virtual int nMdc () const
virtual int nMdc () const
virtual int nSvt () const
virtual int nSvt () const
bool operator== (const TrkRep &)
bool operator== (const TrkRep &)
const TrkRecoTrk * parentTrack () const
TrkRecoTrk * parentTrack ()
const TrkRecoTrk * parentTrack () const
TrkRecoTrk * parentTrack ()
virtual PdtPid::PidType particleType () const
virtual PdtPid::PidType particleType () const
virtual HepPoint3D position (double fltL) const
virtual HepPoint3D position (double fltL) const
virtual BesPointErr positionErr (double fltL) const
virtual BesPointErr positionErr (double fltL) const
virtual HepMatrix posmomCov (double fltL) const
virtual HepMatrix posmomCov (double fltL) const
virtual void print (std::ostream &ostr) const
virtual void print (std::ostream &ostr) const
virtual void printAll (std::ostream &ostr) const
virtual void printAll (std::ostream &ostr) const
std::ostream & printHistory (std::ostream &os=std::cout) const
std::ostream & printHistory (std::ostream &os=std::cout) const
std::ostream & printStatus (std::ostream &os=std::cout) const
std::ostream & printStatus (std::ostream &os=std::cout) const
void printType (std::ostream &ostr) const
void printType (std::ostream &ostr) const
virtual double pt (double fltL=0.) const
virtual double pt (double fltL=0.) const
virtual void removeHot (TrkHitOnTrk *theHot)
virtual void removeHot (TrkHitOnTrk *theHot)
virtual bool resid (const TrkHitOnTrk *theHot, double &residual, double &residErr, bool exclude=false) const
virtual bool resid (const TrkHitOnTrk *theHot, double &residual, double &residErr, bool exclude=false) const
history_riterator reverseBeginHistory () const
history_riterator reverseBeginHistory () const
history_riterator reverseEndHistory () const
history_riterator reverseEndHistory () const
void set2d (bool d)
void set2d (bool d)
void setChisq (double c)
void setChisq (double c)
void setCurrent (bool c)
void setCurrent (bool c)
void setMultScat (bool m)
void setMultScat (bool m)
void setValid (bool v)
void setValid (bool v)
virtual double startFoundRange () const
virtual double startFoundRange () const
double startValidRange () const
double startValidRange () const
double trackT0 () const
double trackT0 () const
virtual const TrkDifTraj & traj () const
virtual TrkDifTraj & traj ()
virtual const TrkDifTraj & traj () const
virtual TrkDifTraj & traj ()
 TrkCircleRep (const TrkExchangePar &inPar, TrkRecoTrk *myTrack, PdtPid::PidType hypo, int nActive, int nSvt, int nDch, double chi2, double startFoundRange, double endFoundRange)
 TrkCircleRep (const TrkExchangePar &inPar, TrkRecoTrk *myTrack, PdtPid::PidType hypo, const TrkHotList *=0)
 TrkCircleRep (const TrkExchangePar &inPar, TrkRecoTrk *myTrack, PdtPid::PidType hypo, int nActive, int nSvt, int nDch, double chi2, double startFoundRange, double endFoundRange)
 TrkCircleRep (const TrkExchangePar &inPar, TrkRecoTrk *myTrack, PdtPid::PidType hypo, const TrkHotList *=0)
virtual void updateHots ()
virtual void updateHots ()
virtual bool validFlightLength (double fltL, double tolerance=0.0) const
virtual bool validFlightLength (double fltL, double tolerance=0.0) const
virtual ~TrkCircleRep ()
virtual ~TrkCircleRep ()

Protected Member Functions

TrkBase::Functors::setActive setActive (bool active) const
TrkBase::Functors::setActive setActive (bool active) const
void setActivity (TrkHitOnTrk &hot, bool active) const
void setActivity (TrkHitOnTrk &hot, bool active) const
TrkBase::Functors::setParent setParent (TrkRep *parent) const
void setParent (TrkHitOnTrk &hot, TrkRep *parent) const
TrkBase::Functors::setParent setParent (TrkRep *parent) const
void setParent (TrkHitOnTrk &hot, TrkRep *parent) const
virtual const TrkSimpTraj & simpTraj () const
virtual TrkSimpTraj & simpTraj ()
virtual const TrkSimpTraj & simpTraj () const
virtual TrkSimpTraj & simpTraj ()
 TrkCircleRep (const TrkCircleRep &right, TrkRecoTrk *trk, PdtPid::PidType)
 TrkCircleRep (const TrkCircleRep &right, TrkRecoTrk *trk, PdtPid::PidType)
const TrkContext & trkContext () const
const TrkContext & trkContext () const
TrkBase::Functors::updateMeasurement updateMeasurement (const TrkDifTraj *traj=0, bool maintainAmbiguity=false) const
TrkErrCode updateMeasurement (TrkHitOnTrk &hot, const TrkDifTraj *traj=0, bool maintainAmbiguity=false) const
TrkBase::Functors::updateMeasurement updateMeasurement (const TrkDifTraj *traj=0, bool maintainAmbiguity=false) const
TrkErrCode updateMeasurement (TrkHitOnTrk &hot, const TrkDifTraj *traj=0, bool maintainAmbiguity=false) const

Protected Attributes

std::auto_ptr< TrkHotList > _hotList
std::auto_ptr< TrkHotList > _hotList

Private Member Functions

TrkCircleRep & operator= (const TrkCircleRep &right)
TrkCircleRep & operator= (const TrkCircleRep &right)

Private Attributes

TrkCircleTraj _traj

Member Typedef Documentation

typedef std::unary_function<TrkFitStatus,bool> TrkFitStatus::predicate_type [inherited]
 

typedef std::unary_function<TrkFitStatus,bool> TrkFitStatus::predicate_type [inherited]
 

Constructor & Destructor Documentation

TrkCircleRep::TrkCircleRep (  const TrkExchangePar &  inPar,
TrkRecoTrk *  myTrack,
PdtPid::PidType  hypo,
const TrkHotList *  = 0
) 
  

00022                                                     :
00023   TrkSimpleRep(inHots, trk, hypo),
00024   _traj(inPar)
00025 {
00026   set2d(true);
00027 }

TrkCircleRep::TrkCircleRep (  const TrkExchangePar &  inPar,
TrkRecoTrk *  myTrack,
PdtPid::PidType  hypo,
int  nActive,
int  nSvt,
int  nDch,
double  chi2,
double  startFoundRange,
double  endFoundRange
) 
  

00032                                                               :
00033   TrkSimpleRep(trk, hypo, nact, nsv, ndc, chi2, stFndRng, endFndRng),
00034   _traj(inPar)
00035 {
00036   set2d(true);
00037 }

TrkCircleRep::~TrkCircleRep (   )  [virtual]
  

00048 {}

TrkCircleRep::TrkCircleRep (  const TrkCircleRep &  right,
TrkRecoTrk *  trk,
PdtPid::PidType 
)  [protected]
  

00041                                                :
00042   TrkSimpleRep(right, trk, hypo), _traj(right._traj) 
00043 {
00044     set2d(true);
00045 }

TrkCircleRep::TrkCircleRep (  const TrkExchangePar &  inPar,
TrkRecoTrk *  myTrack,
PdtPid::PidType  hypo,
const TrkHotList *  = 0
) 
 

TrkCircleRep::TrkCircleRep (  const TrkExchangePar &  inPar,
TrkRecoTrk *  myTrack,
PdtPid::PidType  hypo,
int  nActive,
int  nSvt,
int  nDch,
double  chi2,
double  startFoundRange,
double  endFoundRange
) 
 

virtual TrkCircleRep::~TrkCircleRep (   )  [virtual]
 

TrkCircleRep::TrkCircleRep (  const TrkCircleRep &  right,
TrkRecoTrk *  trk,
PdtPid::PidType 
)  [protected]
 

Member Function Documentation

virtual void TrkRep::activateHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
 

void TrkRep::activateHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
  

00155 {
00156   if(!hot->isActive()){
00157 // make sure this is my hot we're talking about
00158     if(this == hot->getParentRep()){
00159       setCurrent(false);
00160 // actually activate the hot; this is now the rep's job
00161       hot->setActive(true);
00162     }
00163   }
00164 }

template<class T>
void TrkFitStatus::addHistory (  std::pair< T, T >  p  )  [inline, inherited]
  

00059 { addHistory(p.first,p.second); }

template<class T>
void TrkFitStatus::addHistory (  T  begin,
T  end
)  [inline, inherited]
  

00057 { _history.insert(_history.end(),begin,end); }

virtual void TrkFitStatus::addHistory (  const TrkErrCode &  status,
const char *  modulename
)  [virtual, inherited]
 

template<class T>
void TrkFitStatus::addHistory (  std::pair< T, T >  p  )  [inline, inherited]
  

00059 { addHistory(p.first,p.second); }

template<class T>
void TrkFitStatus::addHistory (  T  begin,
T  end
)  [inline, inherited]
  

00057 { _history.insert(_history.end(),begin,end); }

void TrkFitStatus::addHistory (  const TrkErrCode &  status,
const char *  modulename
)  [virtual, inherited]
  

00078 {
00079 // append the module name to the error code message
00080   _history.push_back(TrkHistory(status,modulename));
00081 }

virtual void TrkRep::addHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
 

void TrkRep::addHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
  

00141 {
00142   if (newHot->isActive()) setCurrent(false);
00143   hotList()->append(newHot);
00144 }

virtual double TrkRep::arrivalTime (  double  fltL  )  const [virtual, inherited]
 

Implements TrkFit.

Reimplemented in TrkLineRep, and TrkLineRep.

double TrkRep::arrivalTime (  double  fltL  )  const [virtual, inherited]
 

Implements TrkFit.

Reimplemented in TrkLineRep, and TrkLineRep.

00193 {
00194   static double cinv = 1./Constants::c;
00195   // Initialize cache
00196   if (_betainv < 0.0) {
00197     double mass2 = Pdt::lookup(particleType())->mass();
00198     mass2 = mass2 * mass2;
00199     double ptot2 = momentum(0.).mag2();
00200     assert(ptot2 != 0.0);
00201     _betainv = sqrt( (ptot2 +  mass2)/ ptot2);
00202   }
00203   double tof = fltL * _betainv * cinv;
00204   return trackT0() + tof;
00205 }

history_iterator TrkFitStatus::beginHistory (   )  const [inline, inherited]
  

00046 { return _history.begin(); }

history_iterator TrkFitStatus::beginHistory (   )  const [inline, inherited]
  

00046 { return _history.begin(); }

virtual int TrkSimpleRep::charge (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

int TrkSimpleRep::charge (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

00118 {
00119   const BField& theField = parentTrack()->bField();
00120   return TrkMomCalculator::charge(simpTraj(), theField, 0.);
00121 }

double TrkSimpleRep::chisq (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

double TrkSimpleRep::chisq (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

00087 {
00088   // Note possible inefficiency in calling virtual simpTraj(); implementation 
00089   //    could be moved to derived classes instead.  
00090   return _chisq;
00091   //  return simpTraj().parameters()->chisquared();
00092 }

virtual ChisqConsistency TrkRep::chisqConsistency (   )  const [virtual, inherited]
 

Implements TrkFit.

ChisqConsistency TrkRep::chisqConsistency (   )  const [virtual, inherited]
 

Implements TrkFit.

00363                                {
00364   if(fitValid())
00365     return ChisqConsistency(chisq(),nDof());
00366   else
00367     return ChisqConsistency();
00368 }

virtual TrkCircleRep* TrkCircleRep::clone (  TrkRecoTrk *  newTrack  )  const [virtual]
 

Implements TrkRep.

TrkCircleRep * TrkCircleRep::clone (  TrkRecoTrk *  newTrack  )  const [virtual]
 

Implements TrkRep.

00052 {
00053   TrkCircleRep* newRep = new TrkCircleRep(*this, theTrack, this->particleType());
00054   newRep->setValid(fitValid());
00055   newRep->setCurrent(fitCurrent());
00056   return newRep;
00057         
00058 }

virtual TrkCircleRep* TrkCircleRep::cloneNewHypo (  PdtPid::PidType  hypo  )  [virtual]
 

Implements TrkRep.

TrkCircleRep * TrkCircleRep::cloneNewHypo (  PdtPid::PidType  hypo  )  [virtual]
 

Implements TrkRep.

00062 {
00063   TrkCircleRep* newRep = new TrkCircleRep(*this, parentTrack(), hypo);
00064   newRep->setValid(fitValid());
00065   newRep->setCurrent(fitCurrent());  // true for circles, anyway
00066   return newRep;
00067 }

virtual void TrkRep::deactivateHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
 

void TrkRep::deactivateHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
  

00168 {
00169   if(hot->isActive()){
00170 // make sure this is my hot we're talking about
00171     if(this == hot->getParentRep()){
00172       setCurrent(false);
00173 // actually deactivate the hot; this is now the rep's job
00174       hot->setActive(false);
00175     }
00176   }
00177 }

virtual Hep3Vector TrkRep::direction (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

Hep3Vector TrkRep::direction (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00187 {
00188   return traj().direction(fltL);
00189 }

virtual BesVectorErr TrkRep::directionErr (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

BesVectorErr TrkRep::directionErr (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00273 {
00274   static DifPoint posD;
00275   static DifVector dirD;
00276   traj().getDFInfo2(fltL, posD, dirD);
00277   BesError symErr(3);
00278   symErr.assign( dirD.errorMatrix( dirD.x.indepPar()->covariance() ));
00279   Hep3Vector dir(dirD.x.number(), dirD.y.number(), dirD.z.number());
00280   return BesVectorErr(dir, symErr);
00281 }

virtual double TrkRep::endFoundRange (   )  const [virtual, inherited]
 

Implements TrkFit.

double TrkRep::endFoundRange (   )  const [virtual, inherited]
 

Implements TrkFit.

00303 {
00304   return hotList()->endFoundRange();
00305 }

history_iterator TrkFitStatus::endHistory (   )  const [inline, inherited]
  

00047 { return _history.end(); }

history_iterator TrkFitStatus::endHistory (   )  const [inline, inherited]
  

00047 { return _history.end(); }

double TrkRep::endValidRange (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

double TrkRep::endValidRange (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

00291 {
00292   return traj().hiRange();
00293 }

virtual TrkErrCode TrkSimpleRep::fit (   )  [virtual, inherited]
 

Implements TrkRep.

TrkErrCode TrkSimpleRep::fit (   )  [virtual, inherited]
 

Implements TrkRep.

00061 {
00062   if (!hitCapable()) {
00063     setValid(false);
00064     return TrkErrCode(TrkErrCode::fail,100,
00065                       "Tried to fit without a list of hits.");
00066   }
00067   TrkErrCode error = _fitter.fit(*hotList(), simpTraj());
00068   if (error.failure()) {
00069     setValid(false);
00070   }
00071   else {
00072     setValid(true);
00073     setCurrent(true);
00074   }
00075   setChisq(_fitter.lastChisq());
00076   return error;
00077 }

bool TrkFitStatus::fitCurrent (   )  const [inline, inherited]
  

00033 {return _fitCurrent;}

bool TrkFitStatus::fitCurrent (   )  const [inline, inherited]
  

00033 {return _fitCurrent;}

const TrkErrCode& TrkFitStatus::fitStatus (   )  const [inline, inherited]
  

00043 { return _history.back().status(); }

const TrkErrCode& TrkFitStatus::fitStatus (   )  const [inline, inherited]
  

00043 { return _history.back().status(); }

TrkHelixFitter& TrkSimpleRep::fitter (  void   )  [inline, inherited]
  

00096 {return _fitter;}

TrkHelixFitter& TrkSimpleRep::fitter (   )  [inline, inherited]
  

00096 {return _fitter;}

bool TrkFitStatus::fitValid (   )  const [inline, inherited]
  

00034 {return _fitValid;}

bool TrkFitStatus::fitValid (   )  const [inline, inherited]
  

00034 {return _fitValid;}

virtual void TrkSimpleRep::getAllCovs (  double  fltL,
HepSymMatrix &  xxCov,
HepSymMatrix &  ppCov,
HepMatrix &  xpCov
)  const [virtual, inherited]
 

Implements TrkAbsFit.

void TrkSimpleRep::getAllCovs (  double  fltL,
HepSymMatrix &  xxCov,
HepSymMatrix &  ppCov,
HepMatrix &  xpCov
)  const [virtual, inherited]
 

Implements TrkAbsFit.

00148 {
00149   const BField& theField = parentTrack()->bField();
00150   TrkMomCalculator::getAllCovs(simpTraj(), theField, fltL,
00151                                xxCov,ppCov,xpCov);
00152 }

virtual void TrkSimpleRep::getAllWeights (  const HepPoint3D &  pt,
HepVector &  pos,
HepVector &  mom,
HepSymMatrix &  xxWeight,
HepSymMatrix &  ppWeight,
HepMatrix &  xpWeight
)  const [virtual, inherited]
 

virtual void TrkSimpleRep::getAllWeights (  double  fltL,
HepVector &  pos,
HepVector &  mom,
HepSymMatrix &  xxWeight,
HepSymMatrix &  ppWeight,
HepMatrix &  xpWeight
)  const [virtual, inherited]
 

Implements TrkAbsFit.

void TrkSimpleRep::getAllWeights (  const HepPoint3D &  pt,
HepVector &  pos,
HepVector &  mom,
HepSymMatrix &  xxWeight,
HepSymMatrix &  ppWeight,
HepMatrix &  xpWeight
)  const [virtual, inherited]
  

00172 {
00173   double fltL=0;
00174   TrkPoca poca(traj(),fltL, pt);
00175   fltL = poca.flt1();
00176   getAllWeights(fltL,pos,mom,xxWeight,ppWeight,xpWeight);
00177 }

void TrkSimpleRep::getAllWeights (  double  fltL,
HepVector &  pos,
HepVector &  mom,
HepSymMatrix &  xxWeight,
HepSymMatrix &  ppWeight,
HepMatrix &  xpWeight
)  const [virtual, inherited]
 

Implements TrkAbsFit.

00160 {
00161   const BField& theField = parentTrack()->bField();
00162   TrkMomCalculator::getAllWeights(simpTraj(), theField, fltL,
00163                                    pos,mom,xxWeight,ppWeight,xpWeight);
00164 }

virtual TrkExchangePar TrkCircleRep::helix (  double  fltLen  )  const [virtual]
 

Implements TrkFit.

TrkExchangePar TrkCircleRep::helix (  double  fltLen  )  const [virtual]
 

Implements TrkFit.

00071 {
00072   TrkExchangePar localPar(_traj.d0(), BesAngle(_traj.phi0()).rad(), 
00073                           _traj.omega(), 0.0, 0.0);
00074   HepSymMatrix& theCov = localPar.covariance();
00075   for (int i = 1; i < 4; i++) {
00076     for (int j = 1; j < 4; j++) {
00077       theCov(i,j) = _traj.parameters()->covariance()(i,j);
00078     }
00079   }
00080   return localPar;
00081 }

std::pair<history_iterator,history_iterator> TrkFitStatus::history (   )  const [inline, inherited]
  

00050                                                              {
00051           return std::pair<history_iterator,history_iterator>(beginHistory(),endHistory());
00052   }

std::pair<history_iterator,history_iterator> TrkFitStatus::history (   )  const [inline, inherited]
  

00050                                                              {
00051           return std::pair<history_iterator,history_iterator>(beginHistory(),endHistory());
00052   }

const std::vector<TrkHistory>& TrkFitStatus::historyVector (   )  const [inline, inherited]
  

00053 { return _history; }

const std::vector<TrkHistory>& TrkFitStatus::historyVector (   )  const [inline, inherited]
  

00053 { return _history; }

bool TrkRep::hitCapable (   )  const [inline, inherited]
  

00118 {return hotList()->hitCapable();}

bool TrkRep::hitCapable (   )  const [inline, inherited]
  

00118 {return hotList()->hitCapable();}

virtual const TrkHotList* TrkRep::hotList (   )  const [inline, virtual, inherited]
  

00110 {return _hotList.get();}

virtual TrkHotList* TrkRep::hotList (   )  [inline, virtual, inherited]
  

00109 {return _hotList.get();}

virtual const TrkHotList* TrkRep::hotList (   )  const [inline, virtual, inherited]
  

00110 {return _hotList.get();}

virtual TrkHotList* TrkRep::hotList (   )  [inline, virtual, inherited]
  

00109 {return _hotList.get();}

bool TrkFitStatus::is2d (   )  const [inline, inherited]
  

00032 {return _is2d;}

bool TrkFitStatus::is2d (   )  const [inline, inherited]
  

00032 {return _is2d;}

virtual Hep3Vector TrkSimpleRep::momentum (  double  fltL = 0.  )  const [virtual, inherited]
 

Implements TrkAbsFit.

Hep3Vector TrkSimpleRep::momentum (  double  fltL = 0.  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00104 {
00105   const BField& theField = parentTrack()->bField();
00106   return TrkMomCalculator::vecMom(simpTraj(), theField, fltL);
00107 }

virtual BesVectorErr TrkSimpleRep::momentumErr (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

BesVectorErr TrkSimpleRep::momentumErr (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00125 {
00126   const BField& theField = parentTrack()->bField();
00127   return TrkMomCalculator::errMom(simpTraj(), theField, fltL);
00128 }

bool TrkFitStatus::multScat (   )  const [inline, inherited]
  

00035 {return _multScat;}

bool TrkFitStatus::multScat (   )  const [inline, inherited]
  

00035 {return _multScat;}

virtual const IfdKey& TrkRep::myKey (   )  const [virtual, inherited]
 

Reimplemented in TrkHelixRep, and TrkHelixRep.

const IfdKey & TrkRep::myKey (   )  const [virtual, inherited]
 

Reimplemented in TrkHelixRep, and TrkHelixRep.

00315 {
00316   // This provides a default key (used to provide Rep-specific interfaces
00317   //   to TrkRecoTrk consumers).
00318   static IfdIntKey _theKey(0);
00319   return _theKey;
00320 }

virtual int TrkRep::nActive (   )  const [virtual, inherited]
 

Implements TrkFit.

int TrkRep::nActive (   )  const [virtual, inherited]
 

Implements TrkFit.

00331 {
00332   return hotList()->nActive();
00333 }

int TrkSimpleRep::nDof (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

int TrkSimpleRep::nDof (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

00081 {
00082   return nActive() - simpTraj().parameters()->nPar();
00083 }

virtual int TrkRep::nMdc (   )  const [virtual, inherited]
 

Implements TrkFit.

int TrkRep::nMdc (   )  const [virtual, inherited]
 

Implements TrkFit.

00343 {
00344   return hotList()->nMdc();
00345 }

virtual int TrkRep::nSvt (   )  const [virtual, inherited]
 

Implements TrkFit.

int TrkRep::nSvt (   )  const [virtual, inherited]
 

Implements TrkFit.

00337 {
00338   return hotList()->nSvt();
00339 }

TrkCircleRep& TrkCircleRep::operator= (  const TrkCircleRep &  right  )  [private]
 

TrkCircleRep& TrkCircleRep::operator= (  const TrkCircleRep &  right  )  [private]
 

bool TrkRep::operator== (  const TrkRep &   )  [inherited]
 

bool TrkRep::operator== (  const TrkRep &   )  [inherited]
  

00135 {
00136   return (&rhs == this);
00137 }

const TrkRecoTrk* TrkRep::parentTrack (   )  const [inline, inherited]
  

00083 {return _parentTrack;}

TrkRecoTrk* TrkRep::parentTrack (   )  [inline, inherited]
  

00082 {return _parentTrack;}

const TrkRecoTrk* TrkRep::parentTrack (   )  const [inline, inherited]
  

00083 {return _parentTrack;}

TrkRecoTrk* TrkRep::parentTrack (   )  [inline, inherited]
  

00082 {return _parentTrack;}

virtual PdtPid::PidType TrkRep::particleType (   )  const [virtual, inherited]
 

Implements TrkFit.

PdtPid::PidType TrkRep::particleType (   )  const [virtual, inherited]
 

Implements TrkFit.

00309 {
00310   return _partHypo;
00311 }

virtual HepPoint3D TrkRep::position (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

HepPoint3D TrkRep::position (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00181 {
00182   return traj().position(fltL);
00183 }

virtual BesPointErr TrkRep::positionErr (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

BesPointErr TrkRep::positionErr (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00215 {
00216   static DifPoint posD;
00217   static DifVector dirD;
00218   traj().getDFInfo2(fltL, posD, dirD);
00219   HepMatrix err = posD.errorMatrix( posD.x.indepPar()->covariance() );
00220   HepPoint3D point(posD.x.number(), posD.y.number(), posD.z.number());
00221   BesError symErr(3);
00222   symErr.assign(err);
00223 
00224   if (false) {
00225 #ifdef MDCPATREC_ROUTINE 
00226     cout<< "ErrMsg(routine) " << "Pos " 
00227             << err.num_row() << " " << err.num_col() << endl
00228             << "output:" << endl
00229       //    << err(1,1) << endl
00230       //    << err(2,1) << "  " << err(2,2) << endl
00231       //    << err(3,1) << "  " << err(3,2) << "  " << err(3,3) << endl
00232                     << "x deriv: " << endl
00233                     << posD.x.derivatives() << endl
00234                     << "y deriv: " << endl
00235                     << posD.y.derivatives() << endl
00236                     << endl;
00237 #endif
00238     //  }
00239 
00240     Hep3Vector pointDir(point.x(), point.y());
00241     double dirMag = pointDir.mag();
00242     double dist = 5.e-3;
00243     double delx = dist * point.x() / dirMag;
00244     double dely = dist * point.y() / dirMag;
00245     int ierr = 0;
00246     HepMatrix weight = err.inverse(ierr);
00247     double chisq =     weight(1,1) * delx * delx +
00248       2 * weight(2,1) * delx * dely +
00249       weight(2,2) * dely * dely;
00250 #ifdef MDCPATREC_DEBUG
00251     cout << point << endl;
00252     cout << symErr << endl;
00253     cout << "delta: " << delx << "  " << dely << endl;
00254     cout << "chisq: " << chisq << endl;
00255 #endif
00256     double phi0 = helix(fltL).phi0();
00257     delx = dist * cos(phi0);
00258     dely = dist * sin(phi0);
00259     chisq =            weight(1,1) * delx * delx +
00260                    2 * weight(2,1) * delx * dely +
00261                        weight(2,2) * dely * dely;
00262 #ifdef MDCPATREC_DEBUG
00263     cout << "delta: " << delx << "  " << dely << endl;
00264     cout << "chisq: " << chisq << endl;
00265     cout << endl << endl;
00266 #endif
00267   }
00268   return BesPointErr(point, symErr);
00269 }

virtual HepMatrix TrkSimpleRep::posmomCov (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

HepMatrix TrkSimpleRep::posmomCov (  double  fltL  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00138 {
00139   const BField& theField = parentTrack()->bField();
00140   return TrkMomCalculator::posmomCov(simpTraj(), theField, fltL);
00141 }

virtual void TrkCircleRep::print (  std::ostream &  ostr  )  const [virtual]
 

Implements TrkAbsFit.

void TrkCircleRep::print (  std::ostream &  ostr  )  const [virtual]
 

Implements TrkAbsFit.

00085 {
00086   ostr << "TrkCircleRep "
00087        << "phi0: " << BesAngle(_traj.phi0()).rad()
00088        << " om: "<<_traj.omega()
00089        << " d0: " << _traj.d0() 
00090        << " parent track:"<<parentTrack()->id()
00091        << endl;
00092        
00093 }

virtual void TrkCircleRep::printAll (  std::ostream &  ostr  )  const [virtual]
 

Implements TrkAbsFit.

void TrkCircleRep::printAll (  std::ostream &  ostr  )  const [virtual]
 

Implements TrkAbsFit.

00097 {
00098   print(ostr);
00099 }

std::ostream& TrkFitStatus::printHistory (  std::ostream &  os = std::cout  )  const [inherited]
 

ostream & TrkFitStatus::printHistory (  std::ostream &  os = std::cout  )  const [inherited]
  

00085 {
00086   //FIXME: this doesn't work because ostream isn't a std::ostream, so 
00087   //       std::ostream_iterator doesn't like this...
00088   // std::copy(_history.begin(),_history.end(),std::ostream_iterator<TrkHistory>(os,"\n"));
00089   for(history_iterator i=_history.begin();i!=_history.end();++i)
00090           os << *i << endl;
00091   return os;
00092 }

std::ostream& TrkFitStatus::printStatus (  std::ostream &  os = std::cout  )  const [inherited]
 

ostream & TrkFitStatus::printStatus (  std::ostream &  os = std::cout  )  const [inherited]
  

00070 {
00071   os << " 3-d: " << (is2d() == 0) << " ";
00072   return printHistory(os);
00073 }

void TrkFit::printType (  std::ostream &  ostr  )  const [inherited]
 

void TrkFit::printType (  std::ostream &  ostr  )  const [inherited]
  

00032 {
00033   ostr << "Particle type: " << Pdt::lookup( particleType() )->name();
00034 }

virtual double TrkSimpleRep::pt (  double  fltL = 0.  )  const [virtual, inherited]
 

Implements TrkAbsFit.

double TrkSimpleRep::pt (  double  fltL = 0.  )  const [virtual, inherited]
 

Implements TrkAbsFit.

00111 {
00112   const BField& theField = parentTrack()->bField();
00113   return TrkMomCalculator::ptMom(simpTraj(), theField, fltL);
00114 }

virtual void TrkRep::removeHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
 

void TrkRep::removeHot (  TrkHitOnTrk *  theHot  )  [virtual, inherited]
  

00148 {
00149   if(theHot->isActive()) setCurrent(false);     // fit no longer current
00150   hotList()->remove(theHot);
00151 }

virtual bool TrkRep::resid (  const TrkHitOnTrk *  theHot,
double &  residual,
double &  residErr,
bool  exclude = false
)  const [virtual, inherited]
 

Reimplemented in TrkHelixRep, and TrkHelixRep.

bool TrkRep::resid (  const TrkHitOnTrk *  theHot,
double &  residual,
double &  residErr,
bool  exclude = false
)  const [virtual, inherited]
 

Reimplemented in TrkHelixRep, and TrkHelixRep.

00351 {
00352   assert (h != 0);
00353   if (h->parentRep() != this) return false;     // HOT must belong to Rep
00354   if (!h->hasResidual()) return false;          // Residual must be available
00355   if (exclude) return false;                    // FIXME: Can't do unbiased residuals (yet!)
00356 
00357   residual=h->residual();
00358   residErr=h->hitRms();
00359   return true;
00360 }

history_riterator TrkFitStatus::reverseBeginHistory (   )  const [inline, inherited]
  

00048 { return _history.rbegin(); }

history_riterator TrkFitStatus::reverseBeginHistory (   )  const [inline, inherited]
  

00048 { return _history.rbegin(); }

history_riterator TrkFitStatus::reverseEndHistory (   )  const [inline, inherited]
  

00049 { return _history.rend(); }

history_riterator TrkFitStatus::reverseEndHistory (   )  const [inline, inherited]
  

00049 { return _history.rend(); }

void TrkFitStatus::set2d (  bool  d  )  [inline, inherited]
  

00040 { _is2d = d; };

void TrkFitStatus::set2d (  bool  d  )  [inline, inherited]
  

00040 { _is2d = d; };

TrkBase::Functors::setActive TrkHitOnTrkUpdater::setActive (  bool  active  )  const [inline, protected, inherited]
  

00056   { return TrkBase::Functors::setActive(active); }

TrkBase::Functors::setActive TrkHitOnTrkUpdater::setActive (  bool  active  )  const [inline, protected, inherited]
  

00056   { return TrkBase::Functors::setActive(active); }

void TrkHitOnTrkUpdater::setActivity (  TrkHitOnTrk &  hot,
bool  active
)  const [inline, protected, inherited]
  

00044                                                        {
00045     hot.setActive(active); }

void TrkHitOnTrkUpdater::setActivity (  TrkHitOnTrk &  hot,
bool  active
)  const [inline, protected, inherited]
  

00044                                                        {
00045     hot.setActive(active); }

void TrkSimpleRep::setChisq (  double  c  )  [inherited]
 

void TrkSimpleRep::setChisq (  double  c  )  [inherited]
  

00096 {
00097   // See inefficiency comment above.
00098   //  simpTraj().parameters()->setChisquared(c);
00099   _chisq = c;
00100 }

void TrkFitStatus::setCurrent (  bool  c  )  [inline, inherited]
  

00039 { _fitCurrent = c; }

void TrkFitStatus::setCurrent (  bool  c  )  [inline, inherited]
  

00039 { _fitCurrent = c; }

void TrkFitStatus::setMultScat (  bool  m  )  [inline, inherited]
  

00041 { _multScat = m; };

void TrkFitStatus::setMultScat (  bool  m  )  [inline, inherited]
  

00041 { _multScat = m; };

TrkBase::Functors::setParent TrkHitOnTrkUpdater::setParent (  TrkRep *  parent  )  const [inline, protected, inherited]
  

00054   { return TrkBase::Functors::setParent(parent); }

void TrkHitOnTrkUpdater::setParent (  TrkHitOnTrk &  hot,
TrkRep *  parent
)  const [inline, protected, inherited]
  

00047                                                         {
00048     hot._parentRep = parent;
00049   }

TrkBase::Functors::setParent TrkHitOnTrkUpdater::setParent (  TrkRep *  parent  )  const [inline, protected, inherited]
  

00054   { return TrkBase::Functors::setParent(parent); }

void TrkHitOnTrkUpdater::setParent (  TrkHitOnTrk &  hot,
TrkRep *  parent
)  const [inline, protected, inherited]
  

00047                                                         {
00048     hot._parentRep = parent;
00049   }

void TrkFitStatus::setValid (  bool  v  )  [inherited]
 

void TrkFitStatus::setValid (  bool  v  )  [inherited]
  

00058 {
00059   _fitValid = v;
00060   if (!v) setCurrent(false);
00061 }

virtual const TrkSimpTraj& TrkCircleRep::simpTraj (   )  const [protected, virtual]
 

Implements TrkSimpleRep.

virtual TrkSimpTraj& TrkCircleRep::simpTraj (   )  [protected, virtual]
 

Implements TrkSimpleRep.

const TrkSimpTraj & TrkCircleRep::simpTraj (   )  const [protected, virtual]
 

Implements TrkSimpleRep.

00109 {
00110   return _traj;
00111 }

TrkSimpTraj & TrkCircleRep::simpTraj (   )  [protected, virtual]
 

Implements TrkSimpleRep.

00103 {
00104   return _traj;
00105 }

virtual double TrkRep::startFoundRange (   )  const [virtual, inherited]
 

Implements TrkFit.

double TrkRep::startFoundRange (   )  const [virtual, inherited]
 

Implements TrkFit.

00297 {
00298   return hotList()->startFoundRange();
00299 }

double TrkRep::startValidRange (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

double TrkRep::startValidRange (   )  const [virtual, inherited]
 

Implements TrkAbsFit.

00285 {
00286   return traj().lowRange();
00287 }

double TrkRep::trackT0 (   )  const [inherited]
 

double TrkRep::trackT0 (   )  const [inherited]
  

00209 {
00210   return parentTrack()->trackT0();
00211 }

virtual const TrkDifTraj& TrkCircleRep::traj (   )  const [virtual]
 

Implements TrkAbsFit.

virtual TrkDifTraj& TrkCircleRep::traj (   )  [virtual]
 

const TrkDifTraj & TrkCircleRep::traj (   )  const [virtual]
 

Implements TrkAbsFit.

00120 {
00121   return _traj;
00122 }

TrkDifTraj & TrkCircleRep::traj (   )  [virtual]
  

00115 {
00116   return _traj;
00117 }

const TrkContext& TrkRep::trkContext (   )  const [protected, inherited]
 

const TrkContext& TrkRep::trkContext (   )  const [protected, inherited]
 

virtual void TrkRep::updateHots (   )  [virtual, inherited]
 

void TrkRep::updateHots (   )  [virtual, inherited]
  

00324 {
00325   setCurrent(false);
00326   hotList()->updateHots();
00327 }

TrkBase::Functors::updateMeasurement TrkHitOnTrkUpdater::updateMeasurement (  const TrkDifTraj *  traj = 0,
bool  maintainAmbiguity = false
)  const [inline, protected, inherited]
  

00052   { return TrkBase::Functors::updateMeasurement(traj,maintainAmbiguity); }

TrkErrCode TrkHitOnTrkUpdater::updateMeasurement (  TrkHitOnTrk &  hot,
const TrkDifTraj *  traj = 0,
bool  maintainAmbiguity = false
)  const [inline, protected, inherited]
  

00042     { return hot.updateMeasurement(traj,maintainAmbiguity);}

TrkBase::Functors::updateMeasurement TrkHitOnTrkUpdater::updateMeasurement (  const TrkDifTraj *  traj = 0,
bool  maintainAmbiguity = false
)  const [inline, protected, inherited]
  

00052   { return TrkBase::Functors::updateMeasurement(traj,maintainAmbiguity); }

TrkErrCode TrkHitOnTrkUpdater::updateMeasurement (  TrkHitOnTrk &  hot,
const TrkDifTraj *  traj = 0,
bool  maintainAmbiguity = false
)  const [inline, protected, inherited]
  

00042     { return hot.updateMeasurement(traj,maintainAmbiguity);}

virtual bool TrkFit::validFlightLength (  double  fltL,
double  tolerance = 0.0
)  const [virtual, inherited]
 

bool TrkFit::validFlightLength (  double  fltL,
double  tolerance = 0.0
)  const [virtual, inherited]
  

00038                                                             {
00039   return fltl+tolerance >= startValidRange() &&
00040     fltl-tolerance <= endValidRange();
00041 }

Member Data Documentation

std::auto_ptr<TrkHotList> TrkRep::_hotList [protected, inherited]
 

std::auto_ptr<TrkHotList> TrkRep::_hotList [protected, inherited]
 

TrkCircleTraj TrkCircleRep::_traj [private]
 

The documentation for this class was generated from the following files:
Generated on Wed Feb 2 19:06:07 2011 for BOSS6.5.5 by  doxygen 1.3.9.1