#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 ()

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

virtual const SvtHitOnTrack *	svtHitOnTrack () const

bool	isActive () const

bool	isUsable () const

bool	mustUse () const

virtual TrkEnums::TrkViewInfo	whatView () const =0

virtual unsigned	layerNumber () const =0

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

virtual const Trajectory *	hitTraj () const =0

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

virtual bool	timeResid (double &resid, double &error) const =0

virtual bool	timeAbsolute (double &time, double &error) const =0

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 ()

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

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)
{}

◆ ~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()