#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.

114{ ; }

◆ 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.

67{ 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.

123{ 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(), HoughFinder::storeTracks(), timeResid(), and CgemMdcFitAlg::updateTracks().

◆ 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.

75 { return _ambig!=0 ? _drift[_ambig<0 ? 0:1]

76 : (_drift[1]-_drift[0])*0.5; }

Referenced by dcaToWire(), MdcTrkRecon::fillTrackList(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), timeResid(), updateMeasurement(), and CgemMdcFitAlg::updateTracks().

◆ drift() [2/2]

double MdcHitOnTrack::drift ( double dca ) const

inline

Definition at line 77 of file MdcHitOnTrack.h.

77{ 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(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ 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.

71{ 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.

115 { 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(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

◆ 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.

104{ _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(), MdcTrack::storeTrack(), HoughFinder::storeTracks(), and CgemMdcFitAlg::updateTracks().

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MdcHitOnTrack() [1/3]

◆ ~MdcHitOnTrack()

◆ MdcHitOnTrack() [2/3]

◆ MdcHitOnTrack() [3/3]

Member Function Documentation

◆ ambig()

◆ baseHit()

◆ changeBase()

◆ charge()

◆ dcaToWire()

◆ dipAngle()

◆ drift() [1/2]

◆ drift() [2/2]

◆ entranceAngle() [1/2]

◆ entranceAngle() [2/2]

◆ entranceAngleHit()

◆ fitTime()

◆ hitTraj()

◆ isBeyondEndflange()

◆ layer()

◆ layerNumber()

◆ layernumber()

◆ mdcHit()

◆ mdcHitOnTrack()

◆ rawTime()

◆ setAmbig()

◆ setT0()

◆ status()

◆ tdcIndex()

◆ timeAbsolute()

◆ timeResid()

◆ updateAmbiguity()

◆ updateMeasurement()

◆ whatView()

◆ whichView()

◆ wire()

◆ wireAmbig()