#include <TrkBmSpotOnTrk.h>

Inheritance diagram for TrkBmSpotOnTrk:

Public Member Functions
	TrkBmSpotOnTrk (const HepPoint3D &ip, const HepSymMatrix &size)

	~TrkBmSpotOnTrk ()

TrkBmSpotOnTrk *	clone (TrkRep , const TrkDifTraj t=0) const

virtual TrkEnums::TrkViewInfo	whatView () const

virtual unsigned	layerNumber () const

const Trajectory *	hitTraj () const

const HepPoint3D &	ip () const

virtual bool	timeResid (double &t, double &tErr) const

virtual bool	timeAbsolute (double &t, double &tErr) const

Public Member Functions inherited from TrkHitOnTrk
	TrkHitOnTrk (const TrkFundHit *, double tolerance)

virtual	~TrkHitOnTrk ()

const TrkRecoTrk *	getParentTrack () const

const TrkRep *	getParentRep () const

PdtPid::PidType	particleType () const

const TrkFundHit *	hit () const

TrkFundHit *	hit ()

const TrkDifTraj *	trkTraj () const

virtual const MdcHitOnTrack *	mdcHitOnTrack () 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

virtual int	ambig () const

virtual void	setAmbig (int newambig)

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
	TrkBmSpotOnTrk (const TrkBmSpotOnTrk &hitToBeCopied, TrkRep newRep, const TrkDifTraj trkTraj=0)

virtual TrkErrCode	updateMeasurement (const TrkDifTraj *traj, bool maintainAmbiguity)

double	GetRms ()

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)

Static Protected Member Functions
static const TrkLineTraj	FindBeamTrajectory (const HepPoint3D &point, const HepSymMatrix &error)

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 28 of file TrkBmSpotOnTrk.h.

Constructor & Destructor Documentation

◆ TrkBmSpotOnTrk() [1/2]

TrkBmSpotOnTrk::TrkBmSpotOnTrk	(	const HepPoint3D &	ip,
		const HepSymMatrix &	size )

Definition at line 28 of file TrkBmSpotOnTrk.cxx.

  : TrkHitOnTrk(0,0.5e-4), 
    _beamTraj(FindBeamTrajectory(ip,size)),
    _ip(ip), 
    _size(size)
{}

Referenced by clone().

◆ ~TrkBmSpotOnTrk()

TrkBmSpotOnTrk::~TrkBmSpotOnTrk ( )

Definition at line 44 of file TrkBmSpotOnTrk.cxx.

45{ }

◆ TrkBmSpotOnTrk() [2/2]

TrkBmSpotOnTrk::TrkBmSpotOnTrk	(	const TrkBmSpotOnTrk &	hitToBeCopied,
		TrkRep *	newRep,
		const TrkDifTraj *	trkTraj = 0 )

protected

Definition at line 36 of file TrkBmSpotOnTrk.cxx.

  : TrkHitOnTrk(hot,newRep,trkTraj), 
    _beamTraj(hot._beamTraj),
    _ip(hot.ip()),
    _size(hot._size)
{}

Member Function Documentation

◆ clone()

TrkBmSpotOnTrk * TrkBmSpotOnTrk::clone	(	TrkRep *	rep,
		const TrkDifTraj *	t = 0 ) const

virtual

Implements TrkHitOnTrk.

Definition at line 48 of file TrkBmSpotOnTrk.cxx.

{
  return  new TrkBmSpotOnTrk(*this,rep,trkTraj);
}

◆ FindBeamTrajectory()

const TrkLineTraj TrkBmSpotOnTrk::FindBeamTrajectory	(	const HepPoint3D &	point,
		const HepSymMatrix &	error )

staticprotected

Definition at line 156 of file TrkBmSpotOnTrk.cxx.

{
    int ifail;
    HepSymMatrix cover(error.inverse(ifail));
    
    if (ifail) {
#ifdef MDCPATREC_FATAL
          std::cout<<"ErrMsg(fatal) TrkLineTraj: "
          <<"Error inverting beamspot error matrix" << std::endl;
#endif
    }
    double dx = -cover.fast(3,1)/cover.fast(1,1);
    double dy = -cover.fast(3,2)/cover.fast(2,2);
    
    HepPoint3D p1 = point + Hep3Vector(-dx,-dy,-1);
    HepPoint3D p2 = point + Hep3Vector(+dx,+dy,+1);
 
    return TrkLineTraj( p1, p2 );
}

◆ GetRms()

double TrkBmSpotOnTrk::GetRms ( )

protected

Definition at line 102 of file TrkBmSpotOnTrk.cxx.

{
    //
    // Get direction
    //
    const TrkDifTraj& trkTraj = getParentRep()->traj();
    Hep3Vector trkDir = trkTraj.direction(fltLen());
    
    //
    // Get errors (assume no correlation)
    //
    double Mxx = 1.0/_size.fast(1,1);
    double Myy = 1.0/_size.fast(2,2);
    
    //
    // Normalized track directions in x/y
    //
    double vx = trkDir[0];
    double vy = trkDir[1];
    double normxy = (vx*vx + vy*vy);
    if (normxy <= 0) return 999.9;
    normxy = sqrt(normxy);
    
    vx /= normxy;
    vy /= normxy;
    
    //
    // Solve for point of least chi2
    //
    double s = vx*vy*(Mxx-Myy)/(vx*vx*Mxx + vy*vy*Myy);
    
    double dx = (-vy + s*vx);
    double dy = (+vx + s*vy);
    
    double chi2 = dx*dx*Mxx + dy*dy*Myy;
    
    return chi2 <= 0 ? 0.0 : (1.0/sqrt(chi2));
}

Referenced by updateMeasurement().

◆ hitTraj()

const Trajectory * TrkBmSpotOnTrk::hitTraj ( ) const

virtual

Implements TrkHitOnTrk.

Definition at line 79 of file TrkBmSpotOnTrk.cxx.

{
  return &_beamTraj;
}

◆ ip()

const HepPoint3D & TrkBmSpotOnTrk::ip ( ) const

Definition at line 85 of file TrkBmSpotOnTrk.cxx.

{
  return _ip;
}

◆ layerNumber()

virtual unsigned TrkBmSpotOnTrk::layerNumber ( ) const

inlinevirtual

Implements TrkHitOnTrk.

Definition at line 36 of file TrkBmSpotOnTrk.h.

36{ return 0; }

◆ timeAbsolute()

virtual bool TrkBmSpotOnTrk::timeAbsolute	(	double &	t,
		double &	tErr ) const

inlinevirtual

Implements TrkHitOnTrk.

Definition at line 43 of file TrkBmSpotOnTrk.h.

43{ return false; }

◆ timeResid()

virtual bool TrkBmSpotOnTrk::timeResid	(	double &	t,
		double &	tErr ) const

inlinevirtual

Implements TrkHitOnTrk.

Definition at line 42 of file TrkBmSpotOnTrk.h.

42{ return false; }

◆ updateMeasurement()

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

protectedvirtual

Implements TrkHitOnTrk.

Definition at line 55 of file TrkBmSpotOnTrk.cxx.

{
  TrkErrCode status=updatePoca(traj,x);
  if (status.success()) {
     assert(_poca!=0);
     setHitResid(_poca->doca());
     setHitRms(GetRms());
  } else {
#ifdef MDCPATREC_WARNING
    std::cout<<"ErrMsg(warning) TrkBmSpotOnTrk::updateMeasurement failed" << std::endl;
#endif
     setHitResid(9999.9);
     setHitRms(9999.9);
  }
  return status;
}

◆ whatView()

TrkEnums::TrkViewInfo TrkBmSpotOnTrk::whatView ( ) const

virtual

Implements TrkHitOnTrk.

Definition at line 73 of file TrkBmSpotOnTrk.cxx.

{
  return TrkEnums::xyView;
}

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

Public Member Functions

Protected Member Functions

Static Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TrkBmSpotOnTrk() [1/2]

◆ ~TrkBmSpotOnTrk()

◆ TrkBmSpotOnTrk() [2/2]

Member Function Documentation

◆ clone()

◆ FindBeamTrajectory()

◆ GetRms()

◆ hitTraj()

◆ ip()

◆ layerNumber()

◆ timeAbsolute()

◆ timeResid()

◆ updateMeasurement()

◆ whatView()