#include <AlignTrack.h>

Inheritance diagram for Trk::AlignTrack:

Collaboration diagram for Trk::AlignTrack:

Public Types
enum	AlignTrackType { Unknown = 0, Original = 1, NormalRefitted = 2, BeamspotConstrained = 3, VertexConstrained = 4, VertexMember = 5, NTrackTypes }

Public Member Functions
	AlignTrack (const Track &associatedTrack)
	constructor takes the track containing the original TrackStateOnSurface objects More...

	AlignTrack (const Track &associatedTrack, AlignTSOSCollection *alignTSOSCollection)
	constructor takes collection of alignTSOS and the track containing the original TrackStateOnSurface objects More...

	AlignTrack (const AlignTrack &alignTrack)
	copy constructor More...

virtual	~AlignTrack ()
	destructor More...

AlignTrack &	operator= (const AlignTrack &alignTrack)
	assignment operator More...

void	setAlignTSOSCollection (AlignTSOSCollection *atsosColl)
	sets collection of AlignTSOS More...

const AlignTSOSCollection *	alignTSOSCollection () const
	returns collection of alignTSOS More...

AlignTSOSCollection::const_iterator	firstAtsos () const
	retrieve iterator pointer to first element in collection More...

AlignTSOSCollection::const_iterator	lastAtsos () const
	returns iterator pointer to last element in collection More...

AlignTSOSCollection::iterator	firstAtsos ()
	retrieve non-const iterator pointer to first element in collection More...

AlignTSOSCollection::iterator	lastAtsos ()
	returns non-const iterator pointer to last element in collection More...

const Track *	originalTrack () const
	retrieve pointer to original track More...

void	setOriginalTrack (const Track *track)
	set pointer to original track More...

AlignTrackType	type () const
	get and set the refit type More...

void	setType (AlignTrackType type)

void	doFindPerigee () const
	For the AlignTrack, if there is pseudo-measurement in the TSOS collection, the perigee will not alway be in the first position. More...

bool	isSLTrack () const
	method to determine whether a straight line track or not More...

void	dump (MsgStream &msg)
	dump align track information More...

const Amg::SymMatrixX *	localErrorMatrix () const
	local error matrix, calculated by AlignTrack by calling atsos->hitDistanceVar() More...

const Amg::SymMatrixX *	localErrorMatrixInv () const
	inverse local error matrix, calculated by AlignTrack by calling atsos->hitDistanceVar() More...

int	nAlignTSOSMeas () const
	number of alignTSOS (including scatterers if included on AlignTrack More...

const std::vector< AlignModuleDerivatives > *	derivatives () const
	The Amg::VectorX is a vector of first-derivatives of the alignTSOS on the alignTrack w.r.t. More...

void	setDerivatives (std::vector< AlignModuleDerivatives > *vec)

const std::vector< AlignModuleDerivatives > *	derivativeErr () const
	The Amg::VectorX is a vector of errors in first-derivatives of the alignTSOS on the alignTrack w.r.t. More...

void	setDerivativeErr (std::vector< AlignModuleDerivatives > *vec)

const std::vector< std::pair< AlignModule , std::vector< double > > >	actualSecondDerivatives () const
	The Amg::VectorX is a vector of first-derivatives of the alignTSOS on the alignTrack w.r.t. More...

void	setActualSecondDerivatives (std::vector< std::pair< AlignModule , std::vector< double > > > vec)

const Amg::VectorX *	residualVector () const
	Vector of residuals of the alignTSOS on the alignTrack, to be set by AlignTrackDresser. More...

void	setResidualVector (Amg::VectorX *vec)

const Amg::SymMatrixX *	weightMatrix () const
	Weight matrix is W from Si alignment (see Eqn. More...

void	setWeightMatrix (Amg::SymMatrixX *mat)

const Amg::SymMatrixX *	weightMatrixFirstDeriv () const
	First deriv weight matrix can be either W from Si alignment (see Eqn. More...

void	setWeightMatrixFirstDeriv (Amg::SymMatrixX *mat)

double	chi2 () const
	locally calculated chi2 More...

void	setChi2 (double chi2)

double	chi2dof () const
	locally calculated chi2 More...

void	setChi2dof (double chi2dof)

const Amg::MatrixX *	derivativeMatrix () const
	set and get derivative matrix More...

void	setDerivativeMatrix (const Amg::MatrixX *matrix)

const Amg::SymMatrixX *	fullCovarianceMatrix () const
	set and get full covariance matrix More...

void	setFullCovarianceMatrix (const Amg::SymMatrixX *matrix)

double	trackAlignParamQuality (int i) const
	set and get array (one entry per alignment parameter) of fit chi2 of chi2 vs. More...

void	setTrackAlignParamQuality (int i, double val)

const Trk::Track *	trackWithoutScattering () const
	returns track with ScatteringAngle pointers all set to zero (used for refit by iPat) More...

const AlignVertex *	getVtx () const
	set and get pointer to the associated vertex More...

AlignVertex *	getVtx ()

void	setVtx (AlignVertex *const vtx)

bool	refitD0 () const
	get refit flags More...

bool	refitZ0 () const

bool	refitPhi () const

bool	refitTheta () const

bool	refitQovP () const

void	setRefitD0 (bool flag)
	set refit flags More...

void	setRefitZ0 (bool flag)

void	setRefitPhi (bool flag)

void	setRefitTheta (bool flag)

void	setRefitQovP (bool flag)

bool	isValid () const
	returns true if the track has non-nullptr fitQuality and non-empty Trk::TrackStates More...

const FitQuality *	fitQuality () const
	return a pointer to the fit quality const-overload More...

FitQuality *	fitQuality ()
	return a pointer to the fit quality non-const overload More...

void	setFitQuality (std::unique_ptr< FitQuality > quality)
	set FitQuality. More...

const Trk::TrackStates *	trackStateOnSurfaces () const
	return a pointer to a const DataVector of const TrackStateOnSurfaces. More...

Trk::TrackStates *	trackStateOnSurfaces ()
	return a pointer to a DataVector of const TrackStateOnSurfaces. More...

void	setTrackStateOnSurfaces (std::unique_ptr< Trk::TrackStates > input)
	Set the TrackStateOnSurfaces. More...

const TrackInfo &	info () const
	Returns a const ref to info of a const tracks. More...

TrackInfo &	info ()
	returns a ref to the info. More...

void	setInfo (const TrackInfo &input)
	set the info. More...

const Trk::TrackSummary *	trackSummary () const
	Returns a pointer to the const Trk::TrackSummary owned by this const track (could be nullptr) More...

Trk::TrackSummary *	trackSummary ()
	Returns a pointer to the Trk::TrackSummary owned by this track (could be nullptr) More...

void	setTrackSummary (std::unique_ptr< Trk::TrackSummary > input)
	Set the track summary. More...

const Perigee *	perigeeParameters () const
	return Perigee. More...

const DataVector< const TrackParameters > *	trackParameters () const
	Return a pointer to a vector of TrackParameters. More...

const DataVector< const MeasurementBase > *	measurementsOnTrack () const
	return a pointer to a vector of MeasurementBase (NOT including any that come from outliers). More...

const DataVector< const MeasurementBase > *	outliersOnTrack () const
	return a pointer to a vector of MeasurementBase, which represent outliers (i.e. More...

void	resetCaches ()
	reset all caches More...

Static Public Member Functions
static bool	isSLTrack (const Track *track)
	static method to determine whether a track is straight line or not More...

static void	dumpTrackInfo (const Track &track, MsgStream &msg)
	dump track information More...

static void	dumpLessTrackInfo (const Track &track, MsgStream &msg)
	dump less track information More...

static std::size_t	numberOfInstantiations ()

Static Public Attributes
static std::atomic_size_t	s_numberOfInstantiations

Protected Types
typedef Trk::TrackStates::const_iterator	TSoS_iterator

Protected Member Functions
void	findPerigee () const
	Find perigee in the vector of track parameters. More...

void	copyHelper (const Track &rhs)
	Helper method to factor common part of copy ctor and copy assignment. More...

Protected Attributes
std::unique_ptr< TrackStates >	m_trackStateVector = nullptr
	TrackStateOnSurface. More...

CxxUtils::CachedValue< DataVector< const TrackParameters > >	m_cachedParameterVector {}
	A vector of TrackParameters: these can be any of the classes that derive from Trk::TrackParameters, for example, Perigee, MeasuredPerigee, AtaCylinder etc. More...

CxxUtils::CachedValue< DataVector< const MeasurementBase > >	m_cachedMeasurementVector {}
	A vector of MeasurementBase: these objects represent the "hits" on the track (but not outliers - see m_cachedOutlierVector) More...

CxxUtils::CachedValue< DataVector< const MeasurementBase > >	m_cachedOutlierVector {}
	These objects represent the "outliers" on the track. More...

CxxUtils::CachedValue< const Perigee * >	m_perigeeParameters {}
	A pointer to the Track's Perigee parameters. More...

std::unique_ptr< FitQuality >	m_fitQuality {nullptr}
	A pointer to the Track's FitQuality. More...

std::unique_ptr< Trk::TrackSummary >	m_trackSummary {nullptr}
	Datamember to cache the TrackSummary. More...

Trk::TrackInfo	m_trackInfo {}
	This is a class which stores the identity of where the track was created, fitted, which properties the reconstruction had. More...

Private Member Functions
	AlignTrack ()

void	setAlignParameters ()
	sets derivatives of alignment parameters More...

void	findPerigeeImpl () const
	find PerigeeImpl. More...

Private Attributes
const Track *	m_originalTrack
	point to original track More...

AlignTrackType	m_type
	refit type of this alignTrack More...

int	m_nAlignTSOSMeas
	variables More...

AlignTSOSCollection *	m_alignTSOSCollection {}
	collection of alignTSOS that make up the alignTrack More...

Amg::SymMatrixX *	m_localErrorMat
	matrix containing local error matrices along diagonal More...

Amg::SymMatrixX *	m_localErrorMatInv
	inverse matrix containing local error matrices along diagonal More...

const Amg::MatrixX *	m_derivativeMatrix
	matrix containing derivative matrix from track fitter More...

const Amg::SymMatrixX *	m_fullCovarianceMatrix
	matrix containing full covariance matrix from track fitter More...

std::vector< AlignModuleDerivatives > *	m_derivatives

std::vector< AlignModuleDerivatives > *	m_derivativeErr

std::vector< std::pair< AlignModule , std::vector< double > > >	m_actualSecondDerivatives

Amg::VectorX *	m_residuals

Amg::MatrixX *	m_weights

Amg::MatrixX *	m_weightsFirstDeriv

double	m_chi2

double	m_chi2dof

AlignVertex *	m_vtx
	pointer to the associated vertex More...

bool	m_refitD0

bool	m_refitZ0

bool	m_refitPhi

bool	m_refitTheta

bool	m_refitQovP

double *	m_trackAlignParamQuality

CxxUtils::CachedUniquePtr< Trk::Track >	m_trackWithoutScattering
	describes the quality of contribution to track to alignment parameter (used by shifting derivatives) More...

Detailed Description

Definition at line 41 of file AlignTrack.h.

Member Typedef Documentation

◆ TSoS_iterator

typedef Trk::TrackStates::const_iterator Trk::Track::TSoS_iterator

protectedinherited

Definition at line 228 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

Member Enumeration Documentation

◆ AlignTrackType

enum Trk::AlignTrack::AlignTrackType

Enumerator
Unknown	default type
Original	not refitted, just copy constructed from original Track
NormalRefitted	normally refitted, without adding any pseudo-measurement
BeamspotConstrained	refitted with beamspot constraint
VertexConstrained	refitted with vertex constraint
VertexMember	Refitted with 3DoF from the VTX.
NTrackTypes	total number of track types

Definition at line 45 of file AlignTrack.h.

                         {
         Unknown = 0,               
         Original = 1,              
         NormalRefitted = 2,        
         BeamspotConstrained = 3,   
         VertexConstrained = 4,     
         VertexMember = 5,          
  
         NTrackTypes                
     };

Constructor & Destructor Documentation

◆ AlignTrack() [1/4]

AlignTrack::AlignTrack ( const Track & associatedTrack )

constructor takes the track containing the original TrackStateOnSurface objects

Definition at line 17 of file AlignTrack.cxx.

     : Track(associatedTrack)
     , m_originalTrack(nullptr)
     , m_type(Unknown)
     , m_nAlignTSOSMeas(0)
     , m_alignTSOSCollection(nullptr)
     , m_localErrorMat(nullptr)
     , m_localErrorMatInv(nullptr)
     , m_derivativeMatrix(nullptr)
     , m_fullCovarianceMatrix(nullptr)
     , m_derivatives(nullptr)
     , m_derivativeErr(nullptr)
     , m_actualSecondDerivatives(nullptr)
     , m_residuals(nullptr)
     , m_weights(nullptr)
     , m_weightsFirstDeriv(nullptr)
     , m_chi2(0.)
     , m_chi2dof(0.)
     , m_trackAlignParamQuality(new double[6])
     , m_trackWithoutScattering()
   {
     for (int i=0;i<6;i++) m_trackAlignParamQuality[i]=0.;
     m_vtx=nullptr;
     m_refitD0=true;
     m_refitZ0=true;
     m_refitPhi=true;
     m_refitTheta=true;
     m_refitQovP=true;
   }

◆ AlignTrack() [2/4]

AlignTrack::AlignTrack	(	const Track &	associatedTrack,
		AlignTSOSCollection *	alignTSOSCollection
	)

constructor takes collection of alignTSOS and the track containing the original TrackStateOnSurface objects

Definition at line 48 of file AlignTrack.cxx.

     : Track(associatedTrack)
     , m_originalTrack(nullptr)
     , m_type(Unknown)
     , m_nAlignTSOSMeas(0)
     , m_localErrorMat(nullptr)
     , m_localErrorMatInv(nullptr)
     , m_derivativeMatrix(nullptr)
     , m_fullCovarianceMatrix(nullptr)
     , m_derivatives(nullptr)
     , m_derivativeErr(nullptr)
     , m_actualSecondDerivatives(nullptr)
     , m_residuals(nullptr)
     , m_weights(nullptr)
     , m_weightsFirstDeriv(nullptr)
     , m_chi2(0.)
     , m_chi2dof(0.)
     , m_trackAlignParamQuality(new double[6])
     , m_trackWithoutScattering()
   {
     for (int i=0;i<6;i++) m_trackAlignParamQuality[i]=0.;
     setAlignTSOSCollection(alignTSOSCollection);
     m_vtx=nullptr;
     m_refitD0=true;
     m_refitZ0=true;
     m_refitPhi=true;
     m_refitTheta=true;
     m_refitQovP=true;
   }

◆ AlignTrack() [3/4]

AlignTrack::AlignTrack ( const AlignTrack & alignTrack )

copy constructor

Definition at line 80 of file AlignTrack.cxx.

     : Track(atrack)
     , m_originalTrack(atrack.m_originalTrack)
     , m_type(atrack.m_type)
     , m_nAlignTSOSMeas(atrack.m_nAlignTSOSMeas)
     , m_localErrorMat(atrack.m_localErrorMat ?
         new Amg::SymMatrixX(*(atrack.m_localErrorMat)) : nullptr)
     , m_localErrorMatInv(atrack.m_localErrorMatInv ?
         new Amg::SymMatrixX(*(atrack.m_localErrorMatInv)) : nullptr)
     , m_derivativeMatrix(atrack.m_derivativeMatrix ?
         new Amg::MatrixX(*(atrack.m_derivativeMatrix)) : nullptr)
     , m_fullCovarianceMatrix(atrack.m_fullCovarianceMatrix ?
         new Amg::SymMatrixX(*(atrack.m_fullCovarianceMatrix)) : nullptr)
     , m_derivatives(atrack.m_derivatives ?
         new std::vector<AlignModuleDerivatives>(*(atrack.m_derivatives)) : nullptr)
     , m_derivativeErr(atrack.m_derivativeErr ?
         new std::vector<AlignModuleDerivatives>(*(atrack.m_derivativeErr)) : nullptr)
     , m_actualSecondDerivatives(atrack.m_actualSecondDerivatives ?
         new std::vector<std::pair<AlignModule*,std::vector<double> > > (*(atrack.m_actualSecondDerivatives)) : nullptr)
     , m_residuals(atrack.m_residuals ?
         new Amg::VectorX(*(atrack.m_residuals)) : nullptr)
     , m_weights(atrack.m_weights ?
         new Amg::SymMatrixX(*(atrack.m_weights)) : nullptr)
     , m_weightsFirstDeriv(atrack.m_weightsFirstDeriv ?
         new Amg::SymMatrixX(*(atrack.m_weightsFirstDeriv)) : nullptr)
     , m_chi2(atrack.m_chi2)
     , m_chi2dof(atrack.m_chi2dof)
     , m_trackAlignParamQuality(new double[6])
     , m_trackWithoutScattering()
   {
     if (atrack.m_trackWithoutScattering) {
       m_trackWithoutScattering.set(std::make_unique<Trk::Track>(*(atrack.m_trackWithoutScattering)));
     }
  
     for (int i=0;i<6;i++)
       m_trackAlignParamQuality[i] = atrack.m_trackAlignParamQuality[i];
  
     // copy AlignTSOSCollection
     if(atrack.m_alignTSOSCollection!=nullptr) {
       AlignTSOSCollection * aTSOScoll = new AlignTSOSCollection;
       aTSOScoll->reserve(atrack.m_alignTSOSCollection->size());
       AlignTSOSIt itAtsos_end = atrack.m_alignTSOSCollection->end();
       for(AlignTSOSIt itAtsos = atrack.m_alignTSOSCollection->begin(); itAtsos!=itAtsos_end; ++itAtsos) {
         assert(*itAtsos!=0); // check that is defined.
         AlignTSOS * atsos = new AlignTSOS(**itAtsos);
         aTSOScoll->push_back(atsos);
       }
       m_alignTSOSCollection=aTSOScoll;
     }
     else
       m_alignTSOSCollection=nullptr;
  
     findPerigee();
     m_vtx=nullptr;
     m_refitD0=true;
     m_refitZ0=true;
     m_refitPhi=true;
     m_refitTheta=true;
     m_refitQovP=true;
   }

◆ ~AlignTrack()

AlignTrack::~AlignTrack ( )

virtual

destructor

Definition at line 220 of file AlignTrack.cxx.

   {
     delete m_alignTSOSCollection;     m_alignTSOSCollection=nullptr;
     delete m_residuals;               m_residuals=nullptr;
     delete m_weights;                 m_weights=nullptr;
     delete m_weightsFirstDeriv;       m_weightsFirstDeriv=nullptr;
  
     delete m_derivatives;             m_derivatives=nullptr;
     delete m_derivativeErr;           m_derivativeErr=nullptr;
     delete m_actualSecondDerivatives; m_actualSecondDerivatives=nullptr;
  
     delete m_localErrorMat;           m_localErrorMat=nullptr;
     delete m_localErrorMatInv;        m_localErrorMatInv=nullptr;
  
     delete m_fullCovarianceMatrix;    m_fullCovarianceMatrix=nullptr;
     delete m_derivativeMatrix;        m_derivativeMatrix=nullptr;
  
     delete [] m_trackAlignParamQuality;
   }

◆ AlignTrack() [4/4]

Trk::AlignTrack::AlignTrack ( )

private

Member Function Documentation

◆ actualSecondDerivatives()

const std::vector<std::pair<AlignModule*,std::vector<double> > >* Trk::AlignTrack::actualSecondDerivatives ( ) const

inline

The Amg::VectorX is a vector of first-derivatives of the alignTSOS on the alignTrack w.r.t.

a particular alignment parameter. There is one Amg::VectorX in the vector for each alignment parameter of each align module. Set by DerivCalcTool.

Definition at line 146 of file AlignTrack.h.

146 { return m_actualSecondDerivatives; }

◆ alignTSOSCollection()

const AlignTSOSCollection * AlignTrack::alignTSOSCollection ( ) const

inline

returns collection of alignTSOS

Definition at line 267 of file AlignTrack.h.

267 { return m_alignTSOSCollection; }

◆ chi2()

double Trk::AlignTrack::chi2 ( ) const

inline

locally calculated chi2

Definition at line 165 of file AlignTrack.h.

165 { return m_chi2; }

◆ chi2dof()

double Trk::AlignTrack::chi2dof ( ) const

inline

locally calculated chi2

Definition at line 169 of file AlignTrack.h.

169 { return m_chi2dof; }

◆ copyHelper()

void Track::copyHelper ( const Track & rhs )

protectedinherited

Helper method to factor common part of copy ctor and copy assignment.

Definition at line 56 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

                                              {
   // set the author to be that of the Track being copied.
   m_trackInfo = rhs.m_trackInfo;
  
   // create & copy other variables if available
   if (rhs.fitQuality() != nullptr) {
     m_fitQuality = std::make_unique<Trk::FitQuality>(*(rhs.m_fitQuality));
   }
   // create & copy other variables
   if (rhs.trackSummary() != nullptr) {
     m_trackSummary = std::make_unique<Trk::TrackSummary>(*(rhs.m_trackSummary));
   }
   // Create the TrackStateVector and the perigeeParameters
  
   if (rhs.m_trackStateVector != nullptr) {
     m_trackStateVector = std::make_unique<TrackStates>();
     m_trackStateVector->reserve(rhs.m_trackStateVector->size());
  
     TSoS_iterator itTSoSEnd = rhs.m_trackStateVector->cend();
     for (TSoS_iterator itTSoS = rhs.m_trackStateVector->cbegin();
          itTSoS != itTSoSEnd; ++itTSoS) {
       assert(*itTSoS != nullptr);  // check that is defined.
       // clone and store
       TrackStateOnSurface* tsos = (**itTSoS).clone();
       m_trackStateVector->push_back(tsos);
       // Check if this a perigee so we can already cache it
       if (tsos != nullptr && tsos->type(TrackStateOnSurface::Perigee)) {
         const Trk::Perigee* perigee = nullptr;
         const Trk::TrackParameters* tp = tsos->trackParameters();
         if (tp && tp->type() == Trk::AtaSurface &&
             tp->surfaceType() == Trk::SurfaceType::Perigee) {
           perigee = static_cast<const Trk::Perigee*>(tp);
         }
         if (perigee != nullptr) {
           m_perigeeParameters.store(perigee);  // Now they will be valid
         }
       }
     }
   }
 }

◆ derivativeErr()

const std::vector<AlignModuleDerivatives>* Trk::AlignTrack::derivativeErr ( ) const

inline

The Amg::VectorX is a vector of errors in first-derivatives of the alignTSOS on the alignTrack w.r.t.

a particular alignment parameter. There is one Amg::VectorX in the vector for each alignment parameter of each align module. Set by DerivCalcTool.

Definition at line 140 of file AlignTrack.h.

140 { return m_derivativeErr; }

◆ derivativeMatrix()

const Amg::MatrixX* Trk::AlignTrack::derivativeMatrix ( ) const

inline

set and get derivative matrix

Definition at line 173 of file AlignTrack.h.

173 { return m_derivativeMatrix; }

◆ derivatives()

const std::vector<AlignModuleDerivatives>* Trk::AlignTrack::derivatives ( ) const

inline

The Amg::VectorX is a vector of first-derivatives of the alignTSOS on the alignTrack w.r.t.

a particular alignment parameter. There is one Amg::VectorX in the vector for each alignment parameter of each align module. Set by DerivCalcTool.

Definition at line 133 of file AlignTrack.h.

133 { return m_derivatives; }

◆ doFindPerigee()

void Trk::AlignTrack::doFindPerigee ( ) const

inline

For the AlignTrack, if there is pseudo-measurement in the TSOS collection, the perigee will not alway be in the first position.

We should explicitly do the perigee finding. Fortunately, base Class already has a protected function which loop over the TSOS colleciton to do perigee search and assigin the perigee parameters at the perigee. Just call it

Definition at line 104 of file AlignTrack.h.

104 { findPerigee();}

◆ dump()

void AlignTrack::dump ( MsgStream & msg )

dump align track information

Definition at line 316 of file AlignTrack.cxx.

   {
     msg<<"dumping AlignTrack"<<endmsg;
     double chi2=0.;
     int imeas(1);
     if (m_alignTSOSCollection) {
       int natsos = m_alignTSOSCollection->size();
       for (int iatsos=0;iatsos<natsos;iatsos++) {
         const AlignTSOS* atsos=(*m_alignTSOSCollection)[iatsos];
         std::vector<Residual>::const_iterator itRes     = atsos->firstResidual();
         std::vector<Residual>::const_iterator itRes_end = atsos->lastResidual();
         for ( ; itRes != itRes_end; ++itRes,++imeas) {
           double resNorm = itRes->residualNorm();
           chi2 += resNorm*resNorm;
           //msg<<"resNorm="<<resNorm
           //   <<", errorMat("<<imeas<<")="<<(*m_localErrorMat)(imeas,imeas)
           //   <<endmsg;
         }
         msg<<"iatsos "<<iatsos<<", chi2="<<chi2<<*atsos<<endmsg;
       }
     }
     msg<<"total chi2: "<<m_chi2<<endmsg;
   }

◆ dumpLessTrackInfo()

void AlignTrack::dumpLessTrackInfo	(	const Track &	track,
		MsgStream &	msg
	)

static

dump less track information

Definition at line 276 of file AlignTrack.cxx.

   {
     int ntsos(0);
     for (const TrackStateOnSurface* tsos : *track.trackStateOnSurfaces()) {
       msg<<"ntsos "<<ntsos<<":"<<", type "<<tsos->dumpType();
       //msg << " perigee center of this TSOS: "<< tsos->trackParameters()->associatedSurface()->center() << endmsg;
  
       if (tsos->type(TrackStateOnSurface::Perigee)) 
         msg << ", Perigee"<<endmsg;
  
       else if (tsos->type(TrackStateOnSurface::Outlier))
         msg << ", Outlier"<<endmsg;
  
       else if ( !tsos->type(TrackStateOnSurface::Scatterer) &&
                 !tsos->type(TrackStateOnSurface::InertMaterial)) {
           msg << "," << ( tsos->fitQualityOnSurface() )<<endmsg;
       }
  
       else {
  
         const MaterialEffectsBase*     meb = tsos->materialEffectsOnTrack();
         const MaterialEffectsOnTrack* meot = dynamic_cast<const MaterialEffectsOnTrack*>(meb);
  
         if (meot) {
           msg<<", meot";
           if (meot->scatteringAngles()) 
             msg<<", have angles"<<endmsg;
           else
             msg<<", no angles"<<endmsg;
         }
         else if (tsos->type(TrackStateOnSurface::InertMaterial))
           msg<<", InertMaterial"<<endmsg;
         else
           msg << ", hmm... no material effects on track!"<<endmsg;
       }
       ++ntsos;
     }
   }

◆ dumpTrackInfo()

void AlignTrack::dumpTrackInfo	(	const Track &	track,
		MsgStream &	msg
	)

static

dump track information

Definition at line 269 of file AlignTrack.cxx.

   {
     for (const TrackStateOnSurface* tsos : *track.trackStateOnSurfaces())
       msg<<*tsos;
   }

◆ findPerigee()

void Track::findPerigee ( ) const

protectedinherited

Find perigee in the vector of track parameters.

It can be used to lazy-init the m_perigeeParameters

Definition at line 123 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

                                  {
   if (!m_perigeeParameters.isValid()) {
     findPerigeeImpl();
   }
 }

◆ findPerigeeImpl()

void Track::findPerigeeImpl ( ) const

privateinherited

find PerigeeImpl.

Assumes that Perigee parameters are currently inValid.

Definition at line 129 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

                                      {
   // loop through all passed parameters and, if there is a at Perigee in there,
   // assign it to Perigee parameters. There should never be more
   // than one perigee type.
   // Note that there can be other objects, like VertexOnTrack measurements, with
   // params at a Perigee surface, thus the  TSoS check.
  
   const Trk::Perigee* tmpPerigeeParameters = nullptr;
   if (!m_trackStateVector) {
     return;
   }
   Trk::TrackStates::const_iterator it =
       m_trackStateVector->cbegin();
   Trk::TrackStates::const_iterator itEnd =
       m_trackStateVector->cend();
   for (; it != itEnd; ++it) {
     if ((*it)->type(TrackStateOnSurface::Perigee)) {
       const Trk::TrackParameters* tp = (*it)->trackParameters();
       if (tp && tp->type() == Trk::AtaSurface &&
           tp->surfaceType() == Trk::SurfaceType::Perigee) {
         tmpPerigeeParameters = static_cast<const Trk::Perigee*>(tp);
       }
  
       if (tmpPerigeeParameters != nullptr) {
         break;  // found perigee so stop loop.
       }
     }
   }
   // set to value and valid
   if (tmpPerigeeParameters) {
     m_perigeeParameters.set(tmpPerigeeParameters);
   }
 }

◆ firstAtsos() [1/2]

AlignTSOSCollection::iterator AlignTrack::firstAtsos ( )

inline

retrieve non-const iterator pointer to first element in collection

Definition at line 282 of file AlignTrack.h.

282 { return m_alignTSOSCollection->begin(); }

◆ firstAtsos() [2/2]

AlignTSOSCollection::const_iterator AlignTrack::firstAtsos ( ) const

inline

retrieve iterator pointer to first element in collection

Definition at line 279 of file AlignTrack.h.

279 { return m_alignTSOSCollection->cbegin(); }

◆ fitQuality() [1/2]

FitQuality* Trk::Track::fitQuality ( )

inherited

return a pointer to the fit quality non-const overload

◆ fitQuality() [2/2]

const FitQuality* Trk::Track::fitQuality ( ) const

inherited

return a pointer to the fit quality const-overload

◆ fullCovarianceMatrix()

const Amg::SymMatrixX* Trk::AlignTrack::fullCovarianceMatrix ( ) const

inline

set and get full covariance matrix

Definition at line 177 of file AlignTrack.h.

177 { return m_fullCovarianceMatrix; }

◆ getVtx() [1/2]

AlignVertex* Trk::AlignTrack::getVtx ( )

inline

Definition at line 189 of file AlignTrack.h.

189 { return m_vtx; }

◆ getVtx() [2/2]

const AlignVertex* Trk::AlignTrack::getVtx ( ) const

inline

set and get pointer to the associated vertex

Definition at line 188 of file AlignTrack.h.

188 { return m_vtx; }

◆ info() [1/2]

TrackInfo& Trk::Track::info ( )

inherited

returns a ref to the info.

non-const overload

◆ info() [2/2]

const TrackInfo& Trk::Track::info ( ) const

inherited

Returns a const ref to info of a const tracks.

◆ isSLTrack() [1/2]

bool AlignTrack::isSLTrack ( ) const

method to determine whether a straight line track or not

Definition at line 263 of file AlignTrack.cxx.

   {
     return isSLTrack(this);
   }

◆ isSLTrack() [2/2]

bool AlignTrack::isSLTrack ( const Track * track )

static

static method to determine whether a track is straight line or not

Definition at line 241 of file AlignTrack.cxx.

   {
     const Perigee* pp = track->perigeeParameters();
     if( pp ){
       
       
       const AmgSymMatrix(5)* covMatrix = pp->covariance();
       if( covMatrix ){
         // sum covariance terms of momentum, use it to determine whether fit was SL fit
         double momCov = 0.;
         for( int i=0;i<5;++i )
           momCov += fabs( (*covMatrix)(4,i));
         for( int i=0;i<4;++i )
           momCov += fabs( (*covMatrix)(i,4));
         if( momCov < 1e-10 )
           return true;
       }
     }
     return false;
   }

◆ isValid()

bool Trk::Track::isValid ( ) const

inherited

returns true if the track has non-nullptr fitQuality and non-empty Trk::TrackStates

◆ lastAtsos() [1/2]

AlignTSOSCollection::iterator AlignTrack::lastAtsos ( )

inline

returns non-const iterator pointer to last element in collection

Definition at line 283 of file AlignTrack.h.

283 { return m_alignTSOSCollection->end(); }

◆ lastAtsos() [2/2]

AlignTSOSCollection::const_iterator AlignTrack::lastAtsos ( ) const

inline

returns iterator pointer to last element in collection

Definition at line 280 of file AlignTrack.h.

280 { return m_alignTSOSCollection->cend(); }

◆ localErrorMatrix()

const Amg::SymMatrixX* Trk::AlignTrack::localErrorMatrix ( ) const

inline

local error matrix, calculated by AlignTrack by calling atsos->hitDistanceVar()

Definition at line 122 of file AlignTrack.h.

122 { return m_localErrorMat; }

◆ localErrorMatrixInv()

const Amg::SymMatrixX* Trk::AlignTrack::localErrorMatrixInv ( ) const

inline

inverse local error matrix, calculated by AlignTrack by calling atsos->hitDistanceVar()

Definition at line 125 of file AlignTrack.h.

125 { return m_localErrorMatInv; }

◆ measurementsOnTrack()

const DataVector< const Trk::MeasurementBase > * Track::measurementsOnTrack ( ) const

inherited

return a pointer to a vector of MeasurementBase (NOT including any that come from outliers).

This DataVector is lazily created by this method and cached.

Returns: Pointer to a DV of MeasurementBase. The MeasurementBases are not owned by the DV (it is a view)

Definition at line 178 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

           {
   if (!m_trackStateVector) {
     return nullptr;
   }
  
   // We only need to do work if not valid.
   if (!m_cachedMeasurementVector.isValid()) {
     // create new DataVector which DOES NOT OWN ELEMENTS .
     DataVector<const Trk::MeasurementBase> tmpMeasurementVector(
         SG::VIEW_ELEMENTS);
     // for measurements on track it is very likely that #(meas) ~ #(TSOS)->
     // reserve(#(TSOS))
     tmpMeasurementVector.reserve(m_trackStateVector->size());
  
     TSoS_iterator itTSoSEnd = m_trackStateVector->cend();
     for (TSoS_iterator itTSoS = m_trackStateVector->cbegin(); itTSoS != itTSoSEnd;
          ++itTSoS) {
       if (!(*itTSoS)->type(TrackStateOnSurface::Outlier)) {
         const Trk::MeasurementBase* rot = (*itTSoS)->measurementOnTrack();
         // does it have a measurement ?
         if (rot != nullptr) {
           tmpMeasurementVector.push_back(rot);
         }
       }
     }
     m_cachedMeasurementVector.set(std::move(tmpMeasurementVector));
   }
  
   return m_cachedMeasurementVector.ptr();
 }

◆ nAlignTSOSMeas()

int Trk::AlignTrack::nAlignTSOSMeas ( ) const

inline

number of alignTSOS (including scatterers if included on AlignTrack

Definition at line 128 of file AlignTrack.h.

128 { return m_nAlignTSOSMeas; }

◆ numberOfInstantiations()

static std::size_t Trk::ObjectCounter< Trk::Track >::numberOfInstantiations

inlinestaticinherited

Definition at line 25 of file TrkObjectCounter.h.

   {
 #ifndef NDEBUG
     return s_numberOfInstantiations.load();
 #endif
     return 0;
   }

◆ operator=()

AlignTrack & AlignTrack::operator= ( const AlignTrack & alignTrack )

assignment operator

Definition at line 142 of file AlignTrack.cxx.

   {
     if (this!=&atrack) {
  
       // assign Track content
       Track::operator=(atrack);
  
       // now fill the extra stuff
       m_originalTrack  = atrack.m_originalTrack;
       m_type           = atrack.m_type;
       m_nAlignTSOSMeas = atrack.m_nAlignTSOSMeas;
       m_chi2           = atrack.m_chi2;
       m_chi2dof        = atrack.m_chi2dof;
       m_vtx            = atrack.m_vtx;
       m_refitD0        = atrack.m_refitD0;
       m_refitZ0        = atrack.m_refitZ0;
       m_refitPhi       = atrack.m_refitPhi;
       m_refitTheta     = atrack.m_refitTheta;
       m_refitQovP      = atrack.m_refitQovP;
  
       for (int i=0;i<6;i++)
         m_trackAlignParamQuality[i] = atrack.m_trackAlignParamQuality[i];
  
       // copy AlignTSOSCollection
       delete m_alignTSOSCollection;
       if(atrack.m_alignTSOSCollection!=nullptr) {
         AlignTSOSCollection * aTSOScoll = new AlignTSOSCollection;
         aTSOScoll->reserve(atrack.m_alignTSOSCollection->size());
         AlignTSOSIt itAtsos_end = atrack.m_alignTSOSCollection->end();
         for(AlignTSOSIt itAtsos = atrack.m_alignTSOSCollection->begin(); itAtsos!=itAtsos_end; ++itAtsos) {
           assert(*itAtsos!=0); // check that is defined.
           AlignTSOS * atsos = new AlignTSOS(**itAtsos);
           aTSOScoll->push_back(atsos);
         }
         m_alignTSOSCollection=aTSOScoll;
       }
       else
         m_alignTSOSCollection=nullptr;
  
       findPerigee();
  
       // fill matrices, vectors, etc.
       delete m_derivatives;
       m_derivatives = atrack.m_derivatives ?
           new std::vector<AlignModuleDerivatives>(*(atrack.m_derivatives)) : nullptr;
       delete m_derivativeErr;
       m_derivativeErr = atrack.m_derivativeErr ?
           new std::vector<AlignModuleDerivatives>(*(atrack.m_derivativeErr)) : nullptr;
       delete m_actualSecondDerivatives;
       m_actualSecondDerivatives = atrack.m_actualSecondDerivatives ?
           new std::vector<std::pair<AlignModule*,std::vector<double> > > (*(atrack.m_actualSecondDerivatives)) : nullptr;
       delete m_localErrorMat;
       m_localErrorMat = atrack.m_localErrorMat ?
           new Amg::SymMatrixX(*(atrack.m_localErrorMat)) : nullptr;
       delete m_localErrorMatInv;
       m_localErrorMatInv = atrack.m_localErrorMatInv ?
           new Amg::SymMatrixX(*(atrack.m_localErrorMatInv)) : nullptr;
       delete m_derivativeMatrix;
       m_derivativeMatrix = atrack.m_derivativeMatrix ?
           new Amg::MatrixX(*(atrack.m_derivativeMatrix)) : nullptr;
       delete m_fullCovarianceMatrix;
       m_fullCovarianceMatrix = atrack.m_fullCovarianceMatrix ?
           new Amg::SymMatrixX(*(atrack.m_fullCovarianceMatrix)) : nullptr;
       delete m_residuals;
       m_residuals = atrack.m_residuals ?
           new Amg::VectorX(*(atrack.m_residuals)) : nullptr;
       delete m_weights;
       m_weights = atrack.m_weights ?
           new Amg::SymMatrixX(*(atrack.m_weights)) : nullptr;
       delete m_weightsFirstDeriv;
       m_weightsFirstDeriv = atrack.m_weightsFirstDeriv ?
           new Amg::SymMatrixX(*(atrack.m_weightsFirstDeriv)) : nullptr;
     }
  
     return *this;
   }

◆ originalTrack()

const Track* Trk::AlignTrack::originalTrack ( ) const

inline

retrieve pointer to original track

Definition at line 90 of file AlignTrack.h.

90 { return m_originalTrack; }

◆ outliersOnTrack()

const DataVector< const Trk::MeasurementBase > * Track::outliersOnTrack ( ) const

inherited

return a pointer to a vector of MeasurementBase, which represent outliers (i.e.

measurements not used in the track fit). This DataVector is created lazily by this method and then cached.

Returns: Pointer to a DV of MeasurementBase, representing outliers. The MeasurementBases are not owned by the DV (it is a view)

Definition at line 210 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

           {
   if (!m_trackStateVector) {
     return nullptr;
   }
   // We only need to do work if not valid
   if (!m_cachedOutlierVector.isValid()) {
     // create new DataVector which DOES NOT OWN ELEMENTS .
     DataVector<const Trk::MeasurementBase> tmpOutlierVector(SG::VIEW_ELEMENTS);
     TSoS_iterator itTSoSEnd = m_trackStateVector->cend();
     for (TSoS_iterator itTSoS = m_trackStateVector->cbegin(); itTSoS != itTSoSEnd;
          ++itTSoS) {
       if ((*itTSoS)->type(TrackStateOnSurface::Outlier)) {
         const Trk::MeasurementBase* rot = (*itTSoS)->measurementOnTrack();
         assert(rot != nullptr);
         tmpOutlierVector.push_back(rot);
       }
     }
     m_cachedOutlierVector.set(std::move(tmpOutlierVector));
   }
   return m_cachedOutlierVector.ptr();
 }

◆ perigeeParameters()

const Trk::Perigee * Track::perigeeParameters ( ) const

inherited

return Perigee.

Can be nullptr if no perigee parameters were assigned to the Track.

This method performs lazy initialization and caches the result.

PLEASE NOTE! if there is more than one Perigee in trackStateOnSurfaces (which there shouldn't be!!), only the first one will be returned by Trk::Track::perigeeParameters Although the Perigee is just a type of TrackParameter, it has a dedicated method because of the specific physics interest

Definition at line 163 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

                                                     {
   if (!m_perigeeParameters.isValid()) {
     // findPerigee performs the setting of the parameters
     // i.e does the CachedValue set
     findPerigeeImpl();
   }
  
   // Return payload if valid
   if (m_perigeeParameters.isValid()) {
     return *(m_perigeeParameters.ptr());
   }
  
   return nullptr;
 }

◆ refitD0()

bool Trk::AlignTrack::refitD0 ( ) const

inline

get refit flags

Definition at line 193 of file AlignTrack.h.

193 { return m_refitD0; }

◆ refitPhi()

bool Trk::AlignTrack::refitPhi ( ) const

inline

Definition at line 195 of file AlignTrack.h.

195 { return m_refitPhi; }

◆ refitQovP()

bool Trk::AlignTrack::refitQovP ( ) const

inline

Definition at line 197 of file AlignTrack.h.

197 { return m_refitQovP; }

◆ refitTheta()

bool Trk::AlignTrack::refitTheta ( ) const

inline

Definition at line 196 of file AlignTrack.h.

196 { return m_refitTheta; }

◆ refitZ0()

bool Trk::AlignTrack::refitZ0 ( ) const

inline

Definition at line 194 of file AlignTrack.h.

194 { return m_refitZ0; }

◆ resetCaches()

void Trk::Track::resetCaches ( )

inherited

reset all caches

◆ residualVector()

const Amg::VectorX* Trk::AlignTrack::residualVector ( ) const

inline

Vector of residuals of the alignTSOS on the alignTrack, to be set by AlignTrackDresser.

Owned by AlignTrack.

Definition at line 150 of file AlignTrack.h.

150 { return m_residuals; }

◆ setActualSecondDerivatives()

void Trk::AlignTrack::setActualSecondDerivatives ( std::vector< std::pair< AlignModule *, std::vector< double > > > * vec )

inline

Definition at line 147 of file AlignTrack.h.

147 { m_actualSecondDerivatives=vec; }

◆ setAlignParameters()

void Trk::AlignTrack::setAlignParameters ( )

private

sets derivatives of alignment parameters

◆ setAlignTSOSCollection()

void AlignTrack::setAlignTSOSCollection ( AlignTSOSCollection * atsosColl )

sets collection of AlignTSOS

Definition at line 341 of file AlignTrack.cxx.

   {
     m_alignTSOSCollection=atsosColl;
  
     // calculate m_nAlignTSOSMeas
     int natsos = m_alignTSOSCollection->size();
     for (int iatsos=0;iatsos<natsos;iatsos++)
       m_nAlignTSOSMeas += (*m_alignTSOSCollection)[iatsos]->nResDim();
  
     // calculate V and Vinv    
     m_localErrorMat    = new Amg::SymMatrixX(m_nAlignTSOSMeas,m_nAlignTSOSMeas); 
     m_localErrorMat->setZero();
     m_localErrorMatInv = new Amg::SymMatrixX(m_nAlignTSOSMeas,m_nAlignTSOSMeas);
     m_localErrorMatInv->setZero();
  
     int matrixindex(0);
     AlignTSOSCollection::const_iterator itAtsos     = firstAtsos();
     AlignTSOSCollection::const_iterator itAtsos_end = lastAtsos();
     for ( ; itAtsos != itAtsos_end; ++itAtsos) {
       const AlignTSOS * atsos = *itAtsos;
       std::vector<Residual>::const_iterator itRes     = atsos->firstResidual();
       std::vector<Residual>::const_iterator itRes_end = atsos->lastResidual();
       for ( ; itRes != itRes_end; ++itRes) {
         (*m_localErrorMat)(matrixindex,matrixindex) = itRes->errSq();
         (*m_localErrorMatInv)(matrixindex,matrixindex) = 1./itRes->errSq();
         matrixindex++;
       }
     }
   }

◆ setChi2()

void Trk::AlignTrack::setChi2 ( double chi2 )

inline

Definition at line 166 of file AlignTrack.h.

166 { m_chi2=chi2; }

◆ setChi2dof()

void Trk::AlignTrack::setChi2dof ( double chi2dof )

inline

Definition at line 170 of file AlignTrack.h.

170 { m_chi2dof=chi2dof; }

◆ setDerivativeErr()

void Trk::AlignTrack::setDerivativeErr ( std::vector< AlignModuleDerivatives > * vec )

inline

Definition at line 141 of file AlignTrack.h.

141 { m_derivativeErr=vec; }

◆ setDerivativeMatrix()

void AlignTrack::setDerivativeMatrix ( const Amg::MatrixX * matrix )

inline

Definition at line 269 of file AlignTrack.h.

                                                                        { 
     if (!m_derivativeMatrix) delete m_derivativeMatrix; 
     m_derivativeMatrix=new Amg::MatrixX(*matrix); 
   }

◆ setDerivatives()

void Trk::AlignTrack::setDerivatives ( std::vector< AlignModuleDerivatives > * vec )

inline

Definition at line 134 of file AlignTrack.h.

134 { m_derivatives=vec; }

◆ setFitQuality()

void Trk::Track::setFitQuality ( std::unique_ptr< FitQuality > quality )

inherited

set FitQuality.

◆ setFullCovarianceMatrix()

void AlignTrack::setFullCovarianceMatrix ( const Amg::SymMatrixX * matrix )

inline

Definition at line 274 of file AlignTrack.h.

                                                                               { 
     if (!m_fullCovarianceMatrix) delete m_fullCovarianceMatrix; 
     m_fullCovarianceMatrix=new Amg::SymMatrixX(*matrix); 
   }

◆ setInfo()

void Trk::Track::setInfo ( const TrackInfo & input )

inherited

set the info.

◆ setOriginalTrack()

void Trk::AlignTrack::setOriginalTrack ( const Track * track )

inline

set pointer to original track

Definition at line 93 of file AlignTrack.h.

93 { m_originalTrack = track; }

◆ setRefitD0()

void Trk::AlignTrack::setRefitD0 ( bool flag )

inline

set refit flags

Definition at line 200 of file AlignTrack.h.

200 { m_refitD0 = flag; }

◆ setRefitPhi()

void Trk::AlignTrack::setRefitPhi ( bool flag )

inline

Definition at line 202 of file AlignTrack.h.

202 { m_refitPhi = flag; }

◆ setRefitQovP()

void Trk::AlignTrack::setRefitQovP ( bool flag )

inline

Definition at line 204 of file AlignTrack.h.

204 { m_refitQovP = flag; }

◆ setRefitTheta()

void Trk::AlignTrack::setRefitTheta ( bool flag )

inline

Definition at line 203 of file AlignTrack.h.

203 { m_refitTheta = flag; }

◆ setRefitZ0()

void Trk::AlignTrack::setRefitZ0 ( bool flag )

inline

Definition at line 201 of file AlignTrack.h.

201 { m_refitZ0 = flag; }

◆ setResidualVector()

void Trk::AlignTrack::setResidualVector ( Amg::VectorX * vec )

inline

Definition at line 151 of file AlignTrack.h.

151 { m_residuals=vec;}

◆ setTrackAlignParamQuality()

void Trk::AlignTrack::setTrackAlignParamQuality	(	int	i,
		double	val
	)

inline

Definition at line 182 of file AlignTrack.h.

182 { m_trackAlignParamQuality[i]=val; }

◆ setTrackStateOnSurfaces()

void Trk::Track::setTrackStateOnSurfaces ( std::unique_ptr< Trk::TrackStates > input )

inherited

Set the TrackStateOnSurfaces.

◆ setTrackSummary()

void Trk::Track::setTrackSummary ( std::unique_ptr< Trk::TrackSummary > input )

inherited

Set the track summary.

◆ setType()

void Trk::AlignTrack::setType ( AlignTrackType type )

inline

Definition at line 97 of file AlignTrack.h.

97 {m_type = type;}

◆ setVtx()

void Trk::AlignTrack::setVtx ( AlignVertex *const vtx )

inline

Definition at line 190 of file AlignTrack.h.

190 { m_vtx = vtx; }

◆ setWeightMatrix()

void Trk::AlignTrack::setWeightMatrix ( Amg::SymMatrixX * mat )

inline

Definition at line 157 of file AlignTrack.h.

157 { m_weights=mat; }

◆ setWeightMatrixFirstDeriv()

void Trk::AlignTrack::setWeightMatrixFirstDeriv ( Amg::SymMatrixX * mat )

inline

Definition at line 162 of file AlignTrack.h.

162 { m_weightsFirstDeriv=mat; }

◆ trackAlignParamQuality()

double Trk::AlignTrack::trackAlignParamQuality ( int i ) const

inline

set and get array (one entry per alignment parameter) of fit chi2 of chi2 vs.

alignment params (used for numerical derivatives)

Definition at line 181 of file AlignTrack.h.

181 { return m_trackAlignParamQuality[i]; }

◆ trackParameters()

const DataVector< const Trk::TrackParameters > * Track::trackParameters ( ) const

inherited

Return a pointer to a vector of TrackParameters.

It is created Lazily by this method and then cached.

Returns: Pointer to a DV of TrackParameters, or 0. The TrackParameters are not owned by the DV (it is a view)

Definition at line 97 of file Tracking/TrkEvent/TrkTrack/src/Track.cxx.

           {
  
   if (!m_trackStateVector) {
     return nullptr;
   }
   // Do work only if it is not valid.
   if (!m_cachedParameterVector.isValid()) {
     // create cached parameter vector (which DOES NOT OWN ELEMENTS)
     DataVector<const Trk::TrackParameters> tmp_ParameterVector(
         SG::VIEW_ELEMENTS);
     tmp_ParameterVector.reserve(m_trackStateVector->size());
     TSoS_iterator itTSoSEnd = m_trackStateVector->cend();
     for (TSoS_iterator itTSoS = m_trackStateVector->cbegin(); itTSoS != itTSoSEnd;
          ++itTSoS) {
       const TrackParameters* trackParameters = (*itTSoS)->trackParameters();
       // check to make sure that the TrackParameters exists first
       if (trackParameters != nullptr) {
         tmp_ParameterVector.push_back(trackParameters);
       }
     }
     m_cachedParameterVector.set(std::move(tmp_ParameterVector));
   }
   return m_cachedParameterVector.ptr();
 }

◆ trackStateOnSurfaces() [1/2]

Trk::TrackStates* Trk::Track::trackStateOnSurfaces ( )

inherited

return a pointer to a DataVector of const TrackStateOnSurfaces.

non-const overload

◆ trackStateOnSurfaces() [2/2]

const Trk::TrackStates* Trk::Track::trackStateOnSurfaces ( ) const

inherited

return a pointer to a const DataVector of const TrackStateOnSurfaces.

const overload

◆ trackSummary() [1/2]

Trk::TrackSummary* Trk::Track::trackSummary ( )

inherited

Returns a pointer to the Trk::TrackSummary owned by this track (could be nullptr)

◆ trackSummary() [2/2]

const Trk::TrackSummary* Trk::Track::trackSummary ( ) const

inherited

Returns a pointer to the const Trk::TrackSummary owned by this const track (could be nullptr)

◆ trackWithoutScattering()

const Trk::Track * AlignTrack::trackWithoutScattering ( ) const

returns track with ScatteringAngle pointers all set to zero (used for refit by iPat)

Definition at line 423 of file AlignTrack.cxx.

   {
     if (not m_trackWithoutScattering) {
  
       const Trk::TrackStates* states = this->trackStateOnSurfaces();
       if (!states) return nullptr;
       
       // loop over TSOSs
       Trk::TrackStates::const_iterator tsit = states->begin();
       Trk::TrackStates::const_iterator tsit_end = states->end();
       
       // This is the list of new TSOS 
       auto newTrackStateOnSurfaces = std::make_unique<Trk::TrackStates>();
       newTrackStateOnSurfaces->reserve( states->size() );
       
       for (; tsit!=tsit_end ; ++tsit) {
         auto newMeas  = (*tsit)->measurementOnTrack() ? (*tsit)->measurementOnTrack()->uniqueClone() : nullptr;
         auto newPars  = (*tsit)->trackParameters() ? (*tsit)->trackParameters()->uniqueClone() : nullptr;
         auto newFitQoS= (*tsit)->fitQualityOnSurface();
         auto meb      = (*tsit)->materialEffectsOnTrack() ? (*tsit)->materialEffectsOnTrack()->uniqueClone() : nullptr;
   
         if (meb) {
           //meot is just used as observer, not owner, so can safely duplicate the pointer
           const auto *meot=dynamic_cast<const Trk::MaterialEffectsOnTrack*>(meb.get());
           if (meot) {
             double tinX0=meot->thicknessInX0();
             std::unique_ptr<Trk::EnergyLoss> eLoss =
               meot->energyLoss()
                 ? std::unique_ptr<Trk::EnergyLoss>(meot->energyLoss()->clone())
                 : nullptr;
             const Trk::Surface& surf = meot->associatedSurface();
             std::bitset<MaterialEffectsBase::NumberOfMaterialEffectsTypes> typeMaterial;
             if (eLoss) typeMaterial.set(MaterialEffectsBase::EnergyLossEffects);
             Trk::MaterialEffectsOnTrack* newmeot=
                 new Trk::MaterialEffectsOnTrack(tinX0,std::nullopt,std::move(eLoss),surf,typeMaterial);
             meb.reset(newmeot);
           }
         }
   
         std::bitset<TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes> typePattern;
         for (int i=0;i<(int)Trk::TrackStateOnSurface::NumberOfTrackStateOnSurfaceTypes;i++) {
           if ((*tsit)->type(Trk::TrackStateOnSurface::TrackStateOnSurfaceType(i)))
             typePattern.set(i);
         }
         const Trk::TrackStateOnSurface* newTsos =
           new Trk::TrackStateOnSurface(newFitQoS,
                                        std::move(newMeas),
                                        std::move(newPars),
                                        std::move(meb),
                                        typePattern);
         newTrackStateOnSurfaces->push_back(newTsos);
       }
       
       m_trackWithoutScattering.set(std::make_unique<Trk::Track>( this->info(), std::move(newTrackStateOnSurfaces), 
                                                                  this->fitQuality() ? 
                                                                  this->fitQuality()->uniqueClone() : nullptr ));
     }
     return m_trackWithoutScattering.get();
   }

◆ type()

AlignTrackType Trk::AlignTrack::type ( ) const

inline

get and set the refit type

Definition at line 96 of file AlignTrack.h.

96 {return m_type;}

◆ weightMatrix()

const Amg::SymMatrixX* Trk::AlignTrack::weightMatrix ( ) const

inline

Weight matrix is W from Si alignment (see Eqn.

10 of INDET-PUB-2005-002).
It is related to the residual covariance matrix from TRT alignment by W= V^{-1}RV (R is defined in Eqn. 18 of ATL-INDET-PUB-2007-009).

Definition at line 156 of file AlignTrack.h.

156 { return m_weights; }

◆ weightMatrixFirstDeriv()

const Amg::SymMatrixX* Trk::AlignTrack::weightMatrixFirstDeriv ( ) const

inline

First deriv weight matrix can be either W from Si alignment (see Eqn.

10 of INDET-PUB-2005-002) or the inverse local error matrix

Definition at line 161 of file AlignTrack.h.

161 { return m_weightsFirstDeriv; }

Member Data Documentation

◆ m_actualSecondDerivatives

std::vector<std::pair<AlignModule *,std::vector<double> > >* Trk::AlignTrack::m_actualSecondDerivatives

private

Definition at line 239 of file AlignTrack.h.

◆ m_alignTSOSCollection

AlignTSOSCollection* Trk::AlignTrack::m_alignTSOSCollection {}

private

collection of alignTSOS that make up the alignTrack

Definition at line 223 of file AlignTrack.h.

◆ m_cachedMeasurementVector

CxxUtils::CachedValue<DataVector<const MeasurementBase> > Trk::Track::m_cachedMeasurementVector {}

protectedinherited

A vector of MeasurementBase: these objects represent the "hits" on the track (but not outliers - see m_cachedOutlierVector)

It is created in the return method by looping over all Trk::TrackStateOnSurface adding their pointers to the payload of m_cachedMeasurementVector

Definition at line 271 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_cachedOutlierVector

CxxUtils::CachedValue<DataVector<const MeasurementBase> > Trk::Track::m_cachedOutlierVector {}

protectedinherited

These objects represent the "outliers" on the track.

It is created in the return method by looping over all Trk::TrackStateOnSurface adding their pointers to the payload of m_cachedRioVector

Definition at line 280 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_cachedParameterVector

CxxUtils::CachedValue<DataVector<const TrackParameters> > Trk::Track::m_cachedParameterVector {}

protectedinherited

A vector of TrackParameters: these can be any of the classes that derive from Trk::TrackParameters, for example, Perigee, MeasuredPerigee, AtaCylinder etc.

It is created in the return method by looping over all Trk::TrackStateOnSurface adding their pointers to the payload of m_cachedParameterVector

Definition at line 260 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_chi2

double Trk::AlignTrack::m_chi2

private

Definition at line 245 of file AlignTrack.h.

◆ m_chi2dof

double Trk::AlignTrack::m_chi2dof

private

Definition at line 246 of file AlignTrack.h.

◆ m_derivativeErr

std::vector<AlignModuleDerivatives >* Trk::AlignTrack::m_derivativeErr

private

Definition at line 238 of file AlignTrack.h.

◆ m_derivativeMatrix

const Amg::MatrixX* Trk::AlignTrack::m_derivativeMatrix

private

matrix containing derivative matrix from track fitter

Definition at line 232 of file AlignTrack.h.

◆ m_derivatives

std::vector<AlignModuleDerivatives >* Trk::AlignTrack::m_derivatives

private

Definition at line 237 of file AlignTrack.h.

◆ m_fitQuality

std::unique_ptr<FitQuality> Trk::Track::m_fitQuality {nullptr}

protectedinherited

A pointer to the Track's FitQuality.

Definition at line 293 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_fullCovarianceMatrix

const Amg::SymMatrixX* Trk::AlignTrack::m_fullCovarianceMatrix

private

matrix containing full covariance matrix from track fitter

Definition at line 235 of file AlignTrack.h.

◆ m_localErrorMat

Amg::SymMatrixX* Trk::AlignTrack::m_localErrorMat

private

matrix containing local error matrices along diagonal

Definition at line 226 of file AlignTrack.h.

◆ m_localErrorMatInv

Amg::SymMatrixX* Trk::AlignTrack::m_localErrorMatInv

private

inverse matrix containing local error matrices along diagonal

Definition at line 229 of file AlignTrack.h.

◆ m_nAlignTSOSMeas

int Trk::AlignTrack::m_nAlignTSOSMeas

private

variables

Definition at line 220 of file AlignTrack.h.

◆ m_originalTrack

const Track* Trk::AlignTrack::m_originalTrack

private

point to original track

Definition at line 214 of file AlignTrack.h.

◆ m_perigeeParameters

CxxUtils::CachedValue<const Perigee*> Trk::Track::m_perigeeParameters {}

protectedinherited

A pointer to the Track's Perigee parameters.

This will be null if the track does not contain a Perigee or MeasuredPerigee parameter

Definition at line 288 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_refitD0

bool Trk::AlignTrack::m_refitD0

private

Definition at line 251 of file AlignTrack.h.

◆ m_refitPhi

bool Trk::AlignTrack::m_refitPhi

private

Definition at line 253 of file AlignTrack.h.

◆ m_refitQovP

bool Trk::AlignTrack::m_refitQovP

private

Definition at line 255 of file AlignTrack.h.

◆ m_refitTheta

bool Trk::AlignTrack::m_refitTheta

private

Definition at line 254 of file AlignTrack.h.

◆ m_refitZ0

bool Trk::AlignTrack::m_refitZ0

private

Definition at line 252 of file AlignTrack.h.

◆ m_residuals

Amg::VectorX* Trk::AlignTrack::m_residuals

private

Definition at line 241 of file AlignTrack.h.

◆ m_trackAlignParamQuality

double* Trk::AlignTrack::m_trackAlignParamQuality

private

Definition at line 257 of file AlignTrack.h.

◆ m_trackInfo

Trk::TrackInfo Trk::Track::m_trackInfo {}

protectedinherited

This is a class which stores the identity of where the track was created, fitted, which properties the reconstruction had.

Definition at line 304 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_trackStateVector

std::unique_ptr<TrackStates> Trk::Track::m_trackStateVector = nullptr

protectedinherited

TrackStateOnSurface.

These objects link the various parameters related to a surface, for example, TrackParameter, RIO_OnTrack and FitQualityOnSurface

Definition at line 248 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_trackSummary

std::unique_ptr<Trk::TrackSummary> Trk::Track::m_trackSummary {nullptr}

protectedinherited

Datamember to cache the TrackSummary.

Definition at line 298 of file Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h.

◆ m_trackWithoutScattering

CxxUtils::CachedUniquePtr<Trk::Track> Trk::AlignTrack::m_trackWithoutScattering

private

describes the quality of contribution to track to alignment parameter (used by shifting derivatives)

Definition at line 259 of file AlignTrack.h.

◆ m_type

AlignTrackType Trk::AlignTrack::m_type

private

refit type of this alignTrack

Definition at line 217 of file AlignTrack.h.

◆ m_vtx

AlignVertex* Trk::AlignTrack::m_vtx

private

pointer to the associated vertex

Definition at line 249 of file AlignTrack.h.

◆ m_weights

Amg::MatrixX* Trk::AlignTrack::m_weights

private

Definition at line 242 of file AlignTrack.h.

◆ m_weightsFirstDeriv

Amg::MatrixX* Trk::AlignTrack::m_weightsFirstDeriv

private

Definition at line 243 of file AlignTrack.h.

◆ s_numberOfInstantiations

std::atomic_size_t Trk::ObjectCounter< Trk::Track >::s_numberOfInstantiations

inlinestaticinherited

Definition at line 22 of file TrkObjectCounter.h.

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

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ TSoS_iterator

Member Enumeration Documentation

◆ AlignTrackType

Constructor & Destructor Documentation

◆ AlignTrack() [1/4]

◆ AlignTrack() [2/4]

◆ AlignTrack() [3/4]

◆ ~AlignTrack()

◆ AlignTrack() [4/4]

Member Function Documentation

◆ actualSecondDerivatives()

◆ alignTSOSCollection()

◆ chi2()

◆ chi2dof()

◆ copyHelper()

◆ derivativeErr()

◆ derivativeMatrix()

◆ derivatives()

◆ doFindPerigee()

◆ dump()

◆ dumpLessTrackInfo()

◆ dumpTrackInfo()

◆ findPerigee()

◆ findPerigeeImpl()

◆ firstAtsos() [1/2]

◆ firstAtsos() [2/2]

◆ fitQuality() [1/2]

◆ fitQuality() [2/2]

◆ fullCovarianceMatrix()

◆ getVtx() [1/2]

◆ getVtx() [2/2]

◆ info() [1/2]

◆ info() [2/2]

◆ isSLTrack() [1/2]

◆ isSLTrack() [2/2]

◆ isValid()

◆ lastAtsos() [1/2]

◆ lastAtsos() [2/2]

◆ localErrorMatrix()

◆ localErrorMatrixInv()

◆ measurementsOnTrack()

◆ nAlignTSOSMeas()

◆ numberOfInstantiations()

◆ operator=()

◆ originalTrack()

◆ outliersOnTrack()

◆ perigeeParameters()

◆ refitD0()

◆ refitPhi()

◆ refitQovP()

◆ refitTheta()

◆ refitZ0()

◆ resetCaches()

◆ residualVector()

◆ setActualSecondDerivatives()

◆ setAlignParameters()

◆ setAlignTSOSCollection()

◆ setChi2()

◆ setChi2dof()

◆ setDerivativeErr()

◆ setDerivativeMatrix()

◆ setDerivatives()

◆ setFitQuality()

◆ setFullCovarianceMatrix()

◆ setInfo()

◆ setOriginalTrack()

◆ setRefitD0()

◆ setRefitPhi()

◆ setRefitQovP()