#include <TrkHitOnTrk.h>
Inheritance diagram for TrkHitOnTrk:
Public Types | |
typedef std::unary_function< TrkHitOnTrk, bool > | predicate_type |
typedef std::unary_function< TrkHitOnTrk, bool > | predicate_type |
Public Member Functions | |
virtual int | ambig () const |
virtual int | ambig () const |
virtual TrkHitOnTrk * | clone (TrkRep *parentRep, const TrkDifTraj *trkTraj=0) const =0 |
virtual TrkHitOnTrk * | clone (TrkRep *parentRep, const TrkDifTraj *trkTraj=0) const =0 |
double | fltLen () const |
double | fltLen () const |
TrkErrCode | getFitStuff (double &deltaChi) const |
TrkErrCode | getFitStuff (HepVector &derivs, double &deltaChi) const |
TrkErrCode | getFitStuff (double &deltaChi) const |
TrkErrCode | getFitStuff (HepVector &derivs, double &deltaChi) const |
const TrkRep * | getParentRep () const |
const TrkRep * | getParentRep () const |
const TrkRecoTrk * | getParentTrack () const |
const TrkRecoTrk * | getParentTrack () const |
bool | hasResidual () const |
bool | hasResidual () const |
TrkFundHit * | hit () |
const TrkFundHit * | hit () const |
TrkFundHit * | hit () |
const TrkFundHit * | hit () const |
double | hitLen () const |
double | hitLen () const |
double | hitRms () const |
double | hitRms () const |
virtual const Trajectory * | hitTraj () const =0 |
virtual const Trajectory * | hitTraj () const =0 |
bool | isActive () const |
bool | isActive () const |
bool | isUsable () const |
bool | isUsable () const |
virtual unsigned | layerNumber () const =0 |
virtual unsigned | layerNumber () const =0 |
virtual const MdcHitOnTrack * | mdcHitOnTrack () const |
virtual const MdcHitOnTrack * | mdcHitOnTrack () const |
bool | mustUse () const |
bool | mustUse () const |
bool | operator< (const TrkHitOnTrk &rhs) const |
bool | operator< (const TrkHitOnTrk &rhs) const |
bool | operator== (const TrkHitOnTrk &) const |
bool | operator== (const TrkHitOnTrk &) const |
TrkRecoTrk * | parentTrack () const |
TrkRecoTrk * | parentTrack () const |
PdtPid::PidType | particleType () const |
PdtPid::PidType | particleType () const |
virtual void | print (std::ostream &) const |
virtual void | print (std::ostream &) const |
virtual void | printAll (std::ostream &) const |
virtual void | printAll (std::ostream &) const |
bool | resid (double &resid, double &residErr, bool exclude=false) const |
double | resid (bool exclude=false) const |
bool | resid (double &resid, double &residErr, bool exclude=false) const |
double | resid (bool exclude=false) const |
double | residual () const |
double | residual () const |
void | setActivity (bool turnOn) |
void | setActivity (bool turnOn) |
virtual void | setAmbig (int newambig) |
virtual void | setAmbig (int newambig) |
void | setFltLen (double f) |
void | setFltLen (double f) |
void | setHitRms (double newRms) |
void | setHitRms (double newRms) |
void | setUsability (int usability) |
void | setUsability (int usability) |
virtual const SvtHitOnTrack * | svtHitOnTrack () const |
virtual const SvtHitOnTrack * | svtHitOnTrack () const |
virtual bool | timeAbsolute (double &time, double &error) const =0 |
virtual bool | timeAbsolute (double &time, double &error) const =0 |
virtual bool | timeResid (double &resid, double &error) const =0 |
virtual bool | timeResid (double &resid, double &error) const =0 |
TrkHitOnTrk (const TrkFundHit *, double tolerance) | |
TrkHitOnTrk (const TrkFundHit *, double tolerance) | |
const TrkDifTraj * | trkTraj () const |
const TrkDifTraj * | trkTraj () const |
double | weight () const |
double | weight () const |
virtual TrkEnums::TrkViewInfo | whatView () const =0 |
virtual TrkEnums::TrkViewInfo | whatView () const =0 |
virtual | ~TrkHitOnTrk () |
virtual | ~TrkHitOnTrk () |
Protected Member Functions | |
TrkRep * | parentRep () const |
TrkRep * | parentRep () const |
void | setHitLen (double h) |
void | setHitLen (double h) |
void | setHitResid (double newResid) |
void | setHitResid (double newResid) |
void | setUnusedHit () |
void | setUnusedHit () |
void | setUsedHit () |
void | setUsedHit () |
TrkHitOnTrk (const TrkHitOnTrk &hitToBeCopied, TrkRep *newRep, const TrkDifTraj *trkTraj=0) | |
TrkHitOnTrk (const TrkHitOnTrk &hitToBeCopied, TrkRep *newRep, const TrkDifTraj *trkTraj=0) | |
virtual TrkErrCode | updateMeasurement (const TrkDifTraj *traj, bool maintainAmbiguity)=0 |
virtual TrkErrCode | updateMeasurement (const TrkDifTraj *traj, bool maintainAmbiguity)=0 |
TrkErrCode | updatePoca (const TrkDifTraj *trkTraj, bool maintainAmbiguity) |
TrkErrCode | updatePoca (const TrkDifTraj *trkTraj, bool maintainAmbiguity) |
Protected Attributes | |
double | _hitLen |
double | _hitRms |
bool | _isActive |
int | _isUsable |
TrkRep * | _parentRep |
TrkRep * | _parentRep |
TrkPoca * | _poca |
TrkPoca * | _poca |
double | _resid |
TrkFundHit * | _theHit |
TrkFundHit * | _theHit |
double | _tolerance |
double | _trkLen |
const TrkDifTraj * | _trkTraj |
const TrkDifTraj * | _trkTraj |
Private Member Functions | |
TrkHitOnTrk & | operator= (const TrkHitOnTrk &) |
TrkHitOnTrk & | operator= (const TrkHitOnTrk &) |
TrkHitOnTrk * | setActive (bool active) |
TrkHitOnTrk * | setActive (bool active) |
TrkHitOnTrk * | setParent (TrkRep *rep) |
TrkHitOnTrk * | setParent (TrkRep *rep) |
TrkHitOnTrk (const TrkHitOnTrk &hit) | |
TrkHitOnTrk (const TrkHitOnTrk &hit) | |
Friends | |
class | TrkBase::Functors::setActive |
class | TrkBase::Functors::setParent |
class | TrkBase::Functors::updateMeasurement |
class | TrkHitOnTrkUpdater |
class | TrkHotList |
class | TrkRecoTrk |
class | TrkRep |
|
|
|
|
|
00033 : 00034 _parentRep(0), 00035 _theHit(const_cast<TrkFundHit*>(hit)), 00036 _isActive(true), 00037 _isUsable(true), 00038 //make caches invalid 00039 _hitRms(-9.e50), 00040 _trkLen(0.0), 00041 _hitLen(0.0), 00042 _trkTraj(0), 00043 _poca(0), 00044 _tolerance(tolerance) 00045 { }
|
|
00088 { 00089 delete _poca; 00090 if ( hit() != 0 && getParentRep() != 0 ) { 00091 setUnusedHit(); 00092 } 00093 }
|
|
00050 : _parentRep(newRep), 00051 _theHit(oldHit._theHit), 00052 _isActive(oldHit._isActive), 00053 _isUsable(oldHit._isUsable), 00054 _hitRms(oldHit._hitRms), 00055 _trkLen(oldHit._trkLen), 00056 _hitLen(oldHit._hitLen), 00057 _resid(9999.9), 00058 _trkTraj(0), 00059 _poca(0), 00060 _tolerance(oldHit._tolerance) 00061 { 00062 assert (0 != newRep); 00063 if( oldHit._trkTraj!=0 && trkTraj!=0 && oldHit._trkTraj == trkTraj ) { 00064 // re-use cache as traj are the same 00065 _resid=oldHit._resid; 00066 _trkTraj=trkTraj; 00067 _poca = (oldHit._poca==0 ? 0 :new TrkPoca(*oldHit._poca)); 00068 } else { 00069 double fl = oldHit.fltLen(); 00070 double dum; 00071 const TrkSimpTraj *t1= (oldHit._trkTraj==0?0:oldHit._trkTraj->localTrajectory(fl,dum)); 00072 const TrkSimpTraj *t2= (trkTraj==0?0:trkTraj->localTrajectory(fl,dum)); 00073 if( t1 != 0 && t2 != 0 && t1->parameters()->parameter() == t2->parameters()->parameter() ) { 00074 // re-use cache as traj are sufficiently equiv 00075 _resid=oldHit._resid; 00076 _trkTraj=trkTraj; 00077 _poca = (oldHit._poca==0 ? 0: new TrkPoca(*oldHit._poca)); 00078 00079 } 00080 } 00081 // Only record hots if on default TrkRep 00082 if (getParentRep()->particleType() == getParentTrack()->defaultType()) { 00083 setUsedHit(); 00084 } 00085 }
|
|
|
|
|
|
|
|
|
|
|
|
Reimplemented in MdcHitOnTrack, and MdcHitOnTrack. |
|
Reimplemented in MdcHitOnTrack, and MdcHitOnTrack. 00182 { 00183 return 0;// by default no ambiguity 00184 }
|
|
Implemented in MdcRecoHitOnTrack, TrkBmSpotOnTrk, MdcRecoHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcRecoHitOnTrack, TrkBmSpotOnTrk, MdcRecoHitOnTrack, and TrkBmSpotOnTrk. |
|
00091 {return _trkLen;}
|
|
00091 {return _trkLen;}
|
|
|
|
|
|
00272 { 00273 assert (_trkTraj == &(getParentRep()->traj())); 00274 deltaChi=_resid/hitRms(); // NOTE: use _INTERNAL_ residual 00275 return TrkErrCode(TrkErrCode::succeed); 00276 }
|
|
00249 { 00250 if (_poca==0 || _poca->status().failure()) { 00251 return TrkErrCode(TrkErrCode::fail); 00252 } 00253 // FIXME: I wish I could tell poca to NOT iterate 00254 // and ONLY compute the distance & derivatives... 00255 TrkDifPoca poca(*_trkTraj,fltLen(),*hitTraj(), hitLen(),_tolerance); 00256 if (poca.status().failure()) { 00257 return TrkErrCode(TrkErrCode::fail); 00258 } 00259 if (derivs.num_row() != 0) { 00260 poca.fetchDerivs(derivs); 00261 } else { 00262 derivs = poca.derivs(); 00263 } 00264 double sigInv = 1. / hitRms(); 00265 deltaChi = _resid * sigInv; // NOTE: use _INTERNAL_ residual 00266 derivs *= sigInv; 00267 return TrkErrCode(TrkErrCode::succeed); 00268 }
|
|
00073 {return _parentRep;}
|
|
00073 {return _parentRep;}
|
|
|
|
00153 { 00154 return getParentRep()->parentTrack(); 00155 }
|
|
00104 { return _poca != 0; }
|
|
00104 { return _poca != 0; }
|
|
00076 {return _theHit;}
|
|
00075 {return _theHit;}
|
|
00076 {return _theHit;}
|
|
00075 {return _theHit;}
|
|
00092 {return _hitLen;}
|
|
00092 {return _hitLen;}
|
|
00089 {return _hitRms;}
|
|
00089 {return _hitRms;}
|
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
|
|
00200 {return _isActive;}
|
|
|
|
00201 {return (_isUsable > 0);}
|
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Reimplemented in MdcHitOnTrack, and MdcHitOnTrack. |
|
Reimplemented in MdcHitOnTrack, and MdcHitOnTrack. 00030 {return 0;}
|
|
|
|
00202 {return (_isUsable > 1);}
|
|
|
|
|
|
|
|
|
|
|
|
00176 { 00177 return this == &rhs; 00178 }
|
|
00177 { return _parentRep;}
|
|
00177 { return _parentRep;}
|
|
|
|
00147 { 00148 return parentRep()->parentTrack(); 00149 }
|
|
|
|
00159 { 00160 return getParentRep()->particleType(); 00161 }
|
|
|
|
00132 { 00133 hit()->printAll(o); 00134 o << " hitlen " << hitLen() 00135 << " fltlen " << fltLen() 00136 << " act " << (isActive() != 0) << endl; 00137 }
|
|
|
|
00141 { 00142 print(o); 00143 }
|
|
|
|
|
|
00206 { 00207 assert(getParentRep()!=0); 00208 return getParentRep()->resid(this,resid,residErr,exclude); 00209 }
|
|
00192 { 00193 double r(-99999.9),re(-9999.9); 00194 bool s=getParentRep()->resid(this,r,re,exclude); 00195 if (!s && r<-99999.8) { 00196 #ifdef MDCPATREC_ROUTINE 00197 std::cout<<"ErrMsg(routine) " 00198 << "error calling parentRep()->residual()" << std::endl; 00199 #endif 00200 } 00201 return r; 00202 }
|
|
|
|
00213 { 00214 assert (_trkTraj == &(getParentRep()->traj())); 00215 return _resid; 00216 }
|
|
00195 { _isActive = active; return this; }
|
|
00195 { _isActive = active; return this; }
|
|
|
|
00097 { 00098 if (!isUsable() || isActive()==turnOn) return; 00099 if (getParentRep() != 0) { // needed until Rep-less HoTs go away 00100 turnOn ? parentRep()->activateHot(this) 00101 : parentRep()->deactivateHot(this); 00102 } else { 00103 _isActive = turnOn; 00104 } 00105 }
|
|
Reimplemented in MdcHitOnTrack, and MdcHitOnTrack. |
|
Reimplemented in MdcHitOnTrack, and MdcHitOnTrack. 00188 {} // by default nothing to set
|
|
00147 {_trkLen = f;}
|
|
00147 {_trkLen = f;}
|
|
00178 {_hitLen = h;}
|
|
00178 {_hitLen = h;}
|
|
00176 {_resid = newResid;}
|
|
00176 {_resid = newResid;}
|
|
00154 {_hitRms = newRms;}
|
|
00154 {_hitRms = newRms;}
|
|
00196 { _parentRep = rep; return this; }
|
|
00196 { _parentRep = rep; return this; }
|
|
|
|
00171 { 00172 if (hit() != 0) hit()->setUnusedHit(this); 00173 }
|
|
|
|
00109 { 00110 _isUsable = usability; 00111 if (isActive() && !isUsable()) { 00112 _isActive = false; 00113 if(getParentRep() != 0)parentRep()->deactivateHot(this); 00114 } 00115 if (!isActive() && mustUse()) { 00116 _isActive = true; 00117 if(getParentRep() != 0)parentRep()->activateHot(this); 00118 } 00119 }
|
|
|
|
00165 { 00166 if (hit() != 0) hit()->setUsedHit(this); 00167 }
|
|
|
|
00031 {return 0;}
|
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
00077 { return _trkTraj;}
|
|
00077 { return _trkTraj;}
|
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
|
|
00220 { 00221 _trkTraj = (trkTraj!=0?trkTraj:&(getParentRep()->traj())); 00222 if (_poca==0) { 00223 _poca = new TrkPoca(*_trkTraj,fltLen(), 00224 *hitTraj(), hitLen(),_tolerance); 00225 } else { 00226 *_poca = TrkPoca(*_trkTraj,fltLen(), 00227 *hitTraj(), hitLen(),_tolerance); 00228 } 00229 if(_poca->status().failure()) { 00230 if(isActive()){ 00231 #ifdef MDCPATREC_WARNING 00232 std::cout<<"ErrMsg(warning) " 00233 << " TrkPoca failed in TrkHitOnTrk::updatePoca" 00234 << std::endl; 00235 #endif 00236 } 00237 delete _poca; _poca=0; 00238 return TrkErrCode(TrkErrCode::fail,4); 00239 } 00240 _trkLen = _poca->flt1(); 00241 _hitLen = _poca->flt2(); 00242 double dca=_poca->doca(); 00243 if (!maintainAmb) setAmbig(dca>0?+1:-1); 00244 return TrkErrCode(TrkErrCode::succeed); 00245 }
|
|
|
|
00123 { 00124 // could be cached 00125 double rms=hitRms(); 00126 assert(rms > 0); 00127 return double(1) / ( rms * rms ); 00128 }
|
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
Implemented in MdcHitOnTrack, TrkBmSpotOnTrk, MdcHitOnTrack, and TrkBmSpotOnTrk. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|