MdcHitOnTrack Class Reference

#include <MdcHitOnTrack.h>

Inheritance diagram for MdcHitOnTrack:

Public Member Functions
	MdcHitOnTrack (const TrkFundHit &fundHit, const MdcHit &baseHit, int ambig, double fittime)
virtual	~MdcHitOnTrack ()
double	entranceAngleHit () const
double	entranceAngle () const
double	entranceAngle (const HepPoint3D pos, const Hep3Vector dir) const
double	dipAngle () const
int	ambig () const
int	wireAmbig () const
double	fitTime () const
double	drift () const
double	drift (double dca) const
double	dcaToWire () const
double	rawTime () const
virtual const Trajectory *	hitTraj () const
virtual const MdcHitOnTrack *	mdcHitOnTrack () const
virtual bool	timeResid (double &t, double &tErr) const
virtual bool	timeAbsolute (double &t, double &tErr) const
virtual const MdcHit *	mdcHit () const
virtual unsigned	tdcIndex () const =0
virtual unsigned	status () const =0
int	wire () const
const MdcLayer *	layer () const
int	layernumber () const
unsigned	layerNumber () const
int	whichView () const
double	charge () const
TrkEnums::TrkViewInfo	whatView () const
void	setAmbig (int a)
void	setT0 (double t0)
Public Member Functions inherited from TrkHitOnTrk
	TrkHitOnTrk (const TrkFundHit *, double tolerance)
virtual	~TrkHitOnTrk ()
virtual TrkHitOnTrk *	clone (TrkRep *parentRep, const TrkDifTraj *trkTraj=0) const =0
const TrkRecoTrk *	getParentTrack () const
const TrkRep *	getParentRep () const
PdtPid::PidType	particleType () const
const TrkFundHit *	hit () const
TrkFundHit *	hit ()
const TrkDifTraj *	trkTraj () const
virtual const SvtHitOnTrack *	svtHitOnTrack () const
bool	isActive () const
bool	isUsable () const
bool	mustUse () const
double	hitRms () const
double	weight () const
double	fltLen () const
double	hitLen () const
bool	operator== (const TrkHitOnTrk &) const
bool	operator< (const TrkHitOnTrk &rhs) const
bool	hasResidual () const
TrkErrCode	getFitStuff (HepVector &derivs, double &deltaChi) const
TrkErrCode	getFitStuff (double &deltaChi) const
double	resid (bool exclude=false) const
bool	resid (double &resid, double &residErr, bool exclude=false) const
double	residual () const
TrkRecoTrk *	parentTrack () const
void	setActivity (bool turnOn)
void	setUsability (int usability)
void	setFltLen (double f)
void	setHitRms (double newRms)
virtual void	print (std::ostream &) const
virtual void	printAll (std::ostream &) const

Protected Member Functions
	MdcHitOnTrack (const TrkFundHit *baseHit, int ambig, double fittime, int layer, int wire)
	MdcHitOnTrack (const MdcHitOnTrack &hitToBeCopied, TrkRep *newRep, const TrkDifTraj *trkTraj, const MdcHit *hb=0)
bool	isBeyondEndflange () const
bool	updateAmbiguity (double dca)
virtual TrkErrCode	updateMeasurement (const TrkDifTraj *traj, bool maintainAmbiguity)
const MdcHit *	baseHit () const
void	changeBase (MdcHit *newBase)
Protected Member Functions inherited from TrkHitOnTrk
	TrkHitOnTrk (const TrkHitOnTrk &hitToBeCopied, TrkRep *newRep, const TrkDifTraj *trkTraj=0)
void	setHitResid (double newResid)
TrkRep *	parentRep () const
void	setHitLen (double h)
void	setUsedHit ()
void	setUnusedHit ()
TrkErrCode	updatePoca (const TrkDifTraj *trkTraj, bool maintainAmbiguity)

Additional Inherited Members
Public Types inherited from TrkHitOnTrk
typedef std::unary_function< TrkHitOnTrk, bool >	predicate_type
Protected Attributes inherited from TrkHitOnTrk
TrkRep *	_parentRep
TrkFundHit *	_theHit
bool	_isActive
int	_isUsable
double	_hitRms
double	_trkLen
double	_hitLen
double	_resid
const TrkDifTraj *	_trkTraj
TrkPoca *	_poca
double	_tolerance

Detailed Description

Definition at line 40 of file MdcHitOnTrack.h.

Constructor & Destructor Documentation

MdcHitOnTrack() [1/3]

MdcHitOnTrack::MdcHitOnTrack	(	const TrkFundHit &	fundHit,
		const MdcHit &	baseHit,
		int	ambig,
		double	fittime )

Definition at line 47 of file MdcHitOnTrack.cxx.

: TrkHitOnTrk(&fundHit,10.e-4),
  _ambig(ambig),
  _hitTraj(baseHit.hitTraj()),  
  _fitTime(baseHit.rawTime()-t0*1e-9),
  _dHit(&baseHit)
{
  // need to flag somehow that that we haven't computed things yet...
  //   now, we know that _drift[0] is for ambig <0, and _drift[1] is ambig >0
  //   and _drift is signed according to ambig. Thus _drift[0] should be -tive
  //   and _drift[1] should be +tive. To indicate this are not yet initialized,
  //   put 'impossible' values here...
  _drift[0] = +9999;
  _drift[1] = -9999;
  setHitResid(-21212121.0);
  setHitRms( 0.02 );
  setHitLen(0.5 * baseHit.layer()->zLength());
  setFltLen(0.);
  _startLen = hitTraj()->lowRange() - 5.;
  _endLen   = hitTraj()->hiRange()  + 5.;
}

~MdcHitOnTrack()

MdcHitOnTrack::~MdcHitOnTrack ( )

virtual

Definition at line 113 of file MdcHitOnTrack.cxx.

{ ; }

MdcHitOnTrack() [2/3]

MdcHitOnTrack::MdcHitOnTrack ( const TrkFundHit * baseHit, int ambig, double fittime, int layer, int wire )

protected

MdcHitOnTrack() [3/3]

MdcHitOnTrack::MdcHitOnTrack ( const MdcHitOnTrack & hitToBeCopied, TrkRep * newRep, const TrkDifTraj * trkTraj, const MdcHit * hb = 0 )

protected

Definition at line 99 of file MdcHitOnTrack.cxx.

: TrkHitOnTrk(hot,newRep,trkTraj)
{
  _ambig = hot._ambig;
  _hitTraj  = hot._hitTraj;
  _fitTime  = hot._fitTime;
  _drift[0] = hot._drift[0];
  _drift[1] = hot._drift[1];
  _startLen = hot._startLen;
  _endLen   = hot._endLen;
  _dHit = (hb==0?hot._dHit:hb);
}

Member Function Documentation

ambig()

int MdcHitOnTrack::ambig ( ) const

inlinevirtual

Reimplemented from TrkHitOnTrk.

Definition at line 67 of file MdcHitOnTrack.h.

{ return _ambig; }  // wire wrt track direction

Referenced by MdcTrackList::pickHits(), updateAmbiguity(), updateMeasurement(), and wireAmbig().

baseHit()

const MdcHit * MdcHitOnTrack::baseHit ( ) const

inlineprotected

Definition at line 123 of file MdcHitOnTrack.h.

{ return _dHit; }

Referenced by MdcHitOnTrack().

changeBase()

void MdcHitOnTrack::changeBase ( MdcHit * newBase )

protected

Definition at line 399 of file MdcHitOnTrack.cxx.

{
  _dHit = newBase;
}

charge()

double MdcHitOnTrack::charge

(

)

const

Definition at line 378 of file MdcHitOnTrack.cxx.

{
  return _dHit->charge();
}

dcaToWire()

double MdcHitOnTrack::dcaToWire

(

)

const

Definition at line 123 of file MdcHitOnTrack.cxx.

{
  double dca = -9999.;
  if ( getParentRep() == 0 ) {
    //     cout << "no parent rep" << endl;
    return dca;
  }
  // WARNING: cannot use the internal _poca, as it lags one iteration
  //          behind the fit... therfore use _EXTERNAL_ residual
  if (isActive())  {  // FIXME: currently can only use 'resid()' if isActive..
    dca = resid()+drift();
  } else {
    TrkPoca poca(getParentRep()->traj(), fltLen(), *hitTraj(), hitLen(),
    _tolerance);
    if (poca.status().success()) dca = poca.doca();
  }
  return dca;
}

Referenced by MdcTrkRecon::fillTrackList(), MdcTrack::storeTrack(), and timeResid().

dipAngle()

double MdcHitOnTrack::dipAngle

(

)

const

Definition at line 227 of file MdcHitOnTrack.cxx.

{
  return getParentRep()==0?0:Constants::pi/2-getParentRep()->traj().direction(fltLen()).theta();
}

drift() [1/2]

double MdcHitOnTrack::drift ( ) const

inline

Definition at line 75 of file MdcHitOnTrack.h.

                             { return _ambig!=0 ? _drift[_ambig<0 ? 0:1] 
      : (_drift[1]-_drift[0])*0.5; }

Referenced by dcaToWire(), MdcTrkRecon::fillTrackList(), MdcTrack::storeTrack(), timeResid(), and updateMeasurement().

drift() [2/2]

double MdcHitOnTrack::drift ( double dca ) const

inline

Definition at line 77 of file MdcHitOnTrack.h.

{ return _drift[dca<0?0:1]; }

entranceAngle() [1/2]

double MdcHitOnTrack::entranceAngle

(

)

const

Definition at line 172 of file MdcHitOnTrack.cxx.

{
  static Hep3Vector dir;
  static HepPoint3D pos;
  if (getParentRep() == 0) return 0.;
  getParentRep()->traj().getInfo(fltLen(), pos, dir);
 
  return entranceAngle(pos, dir);
}

Referenced by entranceAngle(), MdcTrkRecon::fillTrackList(), and MdcTrack::storeTrack().

entranceAngle() [2/2]

double MdcHitOnTrack::entranceAngle	(	const HepPoint3D	pos,
		const Hep3Vector	dir ) const

Definition at line 183 of file MdcHitOnTrack.cxx.

{
  double angle = EntranceAngle(dir.phi() - _dHit->phi(pos.z()));
  //std::cout<< "eAngle("<<layernumber()<<","<<wire()<<") dir.phi()  "<<dir.phi()*180./3.14<<" hit phi "<<_dHit->phi(pos.z())*180./3.14<<" eAngle "<<angle*180./3.14<<" degree "<<angle<<std::endl;
 
  //std::cout<< __FILE__ << "   " << __LINE__ << "   ("<<layernumber()<<","<<wire()<<") "<<
  //" phiPoca "<<dir.phi()*180/3.14<< " phiWire "<<_dHit->phi(pos.z())*180/3.14<<" z "<<pos.z()<<
  //" dPhiz "<<_dHit->wire()->dPhizDC(pos.z())*180/3.14<< " eAngle "<<angle*180/3.14<< " angle "<<(dir.phi() - _dHit->phi(pos.z()))*180./3.14<<std::endl;
  /*
     if(g_tupleHit && fabs(angle)>0.0001){
     int layer = layernumber();
     int wireId = wire();
     g_hitLayer = layer;
     g_hitWire = wireId;
 
     int lrCalib=2;
     if (ambig()==1) lrCalib = 0;
     else if (ambig()==-1) lrCalib = 1;
     g_hitAmbig = lrCalib;
     g_hitAmbigMc = haveDigiMcLr[layer][wireId];
     g_hitEAngle = angle*180./3.14;
     g_hitEAngleMc = haveDigiEAngle[layer][wireId]*180./3.14;
     g_hitZ = pos.z();
     g_hitZMc = haveDigiZ[layer][wireId];
     g_hitDrift = _drift[ambig()];
     g_hitDriftMc = haveDigiDrift[layer][wireId];
     g_hitTkIdMc = haveDigiTk[layer][wireId];
     g_hitPhiPoca = dir.phi()*180./3.41;
     g_hitPhiHit = _dHit->phi(pos.z())*180./3.41;
     g_hitPhiHit0 = _dHit->phi()*180./3.41;
     g_hitPhiHitDel = _dHit->wire()->dPhiz()*180./3.41;
     g_tupleHit->write();
     }
     */
  return angle;
}

entranceAngleHit()

double MdcHitOnTrack::entranceAngleHit

(

)

const

Definition at line 161 of file MdcHitOnTrack.cxx.

{
  static Hep3Vector dir;
  static HepPoint3D pos;
  if (getParentRep() == 0) return 0.;
  getParentRep()->traj().getInfo(fltLen(), pos, dir);
 
  return BesAngle(dir.phi() - pos.phi());
}

fitTime()

double MdcHitOnTrack::fitTime ( ) const

inline

Definition at line 71 of file MdcHitOnTrack.h.

{ return _fitTime; }

hitTraj()

const Trajectory * MdcHitOnTrack::hitTraj ( ) const

virtual

Implements TrkHitOnTrk.

Definition at line 384 of file MdcHitOnTrack.cxx.

{
  return _hitTraj;
}

Referenced by dcaToWire(), and MdcHitOnTrack().

isBeyondEndflange()

bool MdcHitOnTrack::isBeyondEndflange ( ) const

inlineprotected

Definition at line 114 of file MdcHitOnTrack.h.

    { return (hitLen() < _startLen || hitLen() > _endLen); }

Referenced by updateMeasurement().

layer()

const MdcLayer * MdcHitOnTrack::layer

(

)

const

Definition at line 354 of file MdcHitOnTrack.cxx.

{
  return _dHit->layer();
}

Referenced by MdcTrkRecon::fillTrackList(), MdcTrackList::pickHits(), and MdcTrack::storeTrack().

layerNumber()

unsigned MdcHitOnTrack::layerNumber ( ) const

virtual

Implements TrkHitOnTrk.

Definition at line 221 of file MdcHitOnTrack.cxx.

{
  return layernumber();
}

layernumber()

int MdcHitOnTrack::layernumber

(

)

const

Definition at line 348 of file MdcHitOnTrack.cxx.

{
  return _dHit->layernumber();
}

Referenced by layerNumber().

mdcHit()

const MdcHit * MdcHitOnTrack::mdcHit ( ) const

virtual

Reimplemented in MdcRecoHitOnTrack.

Definition at line 390 of file MdcHitOnTrack.cxx.

{
  return 0;
}

Referenced by MdcTrackList::dropMultiHotInLayer(), and wireAmbig().

mdcHitOnTrack()

const MdcHitOnTrack * MdcHitOnTrack::mdcHitOnTrack ( ) const

virtual

Reimplemented from TrkHitOnTrk.

Definition at line 155 of file MdcHitOnTrack.cxx.

{
  return this;
}

rawTime()

double MdcHitOnTrack::rawTime

(

)

const

Definition at line 372 of file MdcHitOnTrack.cxx.

{
  return _dHit->rawTime();
}

Referenced by MdcTrkRecon::fillTrackList().

setAmbig()

void MdcHitOnTrack::setAmbig ( int a )

inlinevirtual

Reimplemented from TrkHitOnTrk.

Definition at line 104 of file MdcHitOnTrack.h.

{ _ambig = a<0?-1:a>0?1:0; }

Referenced by updateAmbiguity().

setT0()

void MdcHitOnTrack::setT0

(

double

t0

)

Definition at line 117 of file MdcHitOnTrack.cxx.

{ 
  _fitTime= _dHit->rawTime()-t0*1e-9; 
}

status()

virtual unsigned MdcHitOnTrack::status ( ) const

pure virtual

Implemented in MdcRecoHitOnTrack.

Referenced by updateMeasurement().

tdcIndex()

virtual unsigned MdcHitOnTrack::tdcIndex ( ) const

pure virtual

Implemented in MdcRecoHitOnTrack.

timeAbsolute()

bool MdcHitOnTrack::timeAbsolute ( double & t, double & tErr ) const

virtual

Implements TrkHitOnTrk.

Definition at line 330 of file MdcHitOnTrack.cxx.

{
  double tresid(-1.0);
  if(timeResid(tresid,tErr)){
    // add back the track time
    t = tresid + getParentRep()->parentTrack()->trackT0();
    return true;
  } else
    return false;
}

timeResid()

bool MdcHitOnTrack::timeResid ( double & t, double & tErr ) const

virtual

Implements TrkHitOnTrk.

Definition at line 320 of file MdcHitOnTrack.cxx.

{
  double v = driftVelocity();
  if (v <= 0) return false;
  t = (fabs(drift())-fabs(dcaToWire()))/v;
  tErr= hitRms()/v;
  return true;
}

Referenced by timeAbsolute().

updateAmbiguity()

bool MdcHitOnTrack::updateAmbiguity ( double dca )

protected

Definition at line 143 of file MdcHitOnTrack.cxx.

{
  if (dca < 0 && ambig() >= 0) {
    setAmbig(-1); return isActive();
  } else if (dca > 0 && ambig() <= 0) {
    setAmbig(1); return isActive();
  } else {
    return false;
  }
}

Referenced by updateMeasurement().

updateMeasurement()

TrkErrCode MdcHitOnTrack::updateMeasurement ( const TrkDifTraj * traj, bool maintainAmbiguity )

protectedvirtual

Implements TrkHitOnTrk.

Definition at line 233 of file MdcHitOnTrack.cxx.

{
  TrkErrCode status=updatePoca(traj,maintainAmb);
  if (status.failure()) {
#ifdef MDCPATREC_DEBUG
    std::cout<<" ErrMsg(warning) " << "MdcHitOnTrack::updateMeasurement failed " << status << std::endl;
#endif
    return status;
  }
  assert (_poca!=0);
  double dca=_poca->doca();
  bool forceIteration = (maintainAmb&&ambig()!=0)?false:updateAmbiguity(dca);
  //std::cout<< __FILE__ << "   " << __LINE__ << " maintainAmb "<< maintainAmb
  //<< " maintain&& ambig "<<(maintainAmb&&ambig()!=0)
  //<< " forceIteration  "<<forceIteration<<std::endl;
  assert(ambig()!=0);
  // Check for hits beyond end plates.  !!Turn off hit if it is == temp. hack
  if (isBeyondEndflange()) setUsability(false);
  if (forceIteration || !driftCurrent() ) {
    updateCorrections(); // force recomputation of drift for current ambig(), setting of hitRms
    forceIteration=true;
  }
  setHitResid(dca-drift());
  return !forceIteration?status:
    TrkErrCode(TrkErrCode::succeed, 11, "Ambiguity flipped");
}

whatView()

TrkEnums::TrkViewInfo MdcHitOnTrack::whatView ( ) const

virtual

Implements TrkHitOnTrk.

Definition at line 342 of file MdcHitOnTrack.cxx.

{
  return _dHit->whatView();
}

whichView()

int MdcHitOnTrack::whichView

(

)

const

Definition at line 366 of file MdcHitOnTrack.cxx.

{
  return _dHit->whichView();
}

wire()

int MdcHitOnTrack::wire

(

)

const

Definition at line 360 of file MdcHitOnTrack.cxx.

{
  return _dHit->wirenumber();
}

wireAmbig()

int MdcHitOnTrack::wireAmbig

(

)

const

Definition at line 415 of file MdcHitOnTrack.cxx.

                                   {
  // hit wrt the wire location
 
  //return fabs(entranceAngle())<Constants::pi/2?ambig():-ambig();
  const TrkRep* tkRep = getParentRep();
  static Hep3Vector dir;
  static HepPoint3D pos;
  if (getParentRep() == 0) return 0.;
  getParentRep()->traj().getInfo(fltLen(), pos, dir);
 
  double wireAmb = ambig();
  if (mdcHit()->isCosmicFit()){
    HepPoint3D poca = tkRep->position(0.);
    if ( pos.y() > poca.y()){
      wireAmb = -1.*_ambig;//yzhang 2012-07-17 
      //std::cout<<"MdcHitOnTrack CosmicFit  up ambig *-1"<<std::endl;
    }
  } 
 
  return wireAmb;
}

Referenced by MdcTrkRecon::fillTrackList(), and MdcTrack::storeTrack().

---

The documentation for this class was generated from the following files:

• MdcHitOnTrack.h
• MdcHitOnTrack.cxx