#include <AlignVertex.h>

Collaboration diagram for Trk::AlignVertex:

Public Types
enum	AlignVertexType { Unknown = 0, Primary = 1, Refitted = 2, Accumulated = 3, NVertexTypes }

Public Member Functions
	AlignVertex ()
	default constructor More...

	AlignVertex (const xAOD::Vertex *originalVertex)
	constructor takes the original vertex candidate More...

	AlignVertex (const AlignVertex &rhs)
	copy constructor More...

AlignVertex &	operator= (const AlignVertex &rhs)
	assignment More...

	~AlignVertex ()
	destructor More...

const AlignTrackCollection *	alignTrackCollection () const
	returns collection of alignTracks More...

AlignTrackCollection::const_iterator	firstATrack () const
	retrieve iterator pointer to first element in collection More...

AlignTrackCollection::const_iterator	lastATrack () const
	returns iterator pointer to last element in collection More...

const xAOD::Vertex *	originalVertex () const
	retrieve pointer to original vertex and its position More...

const Amg::Vector3D *	originalPosition () const

void	setOriginalVertex (const xAOD::Vertex *vertex)
	set pointer to original vertex More...

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

void	setType (AlignVertexType type)

const Amg::Vector3D *	position () const
	get the vertex position and covariance More...

const	AmgSymMatrix (3) *covariance() const

const std::vector< AlignModuleVertexDerivatives > *	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< AlignModuleVertexDerivatives > *vec)

void	addDerivatives (std::vector< AlignModuleVertexDerivatives > *vec)

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

void	incrementVector (const Amg::Vector3D &vtxV)
	increment algebra objects for this verterx: More...

void	incrementMatrix (const AmgSymMatrix(3) vtxM)

void	fitVertex ()
	fit the vertex internally More...

int	Ntracks () const
	get the number of contributing tracks More...

void	setConstraint (AmgSymMatrix(3) , Amg::Vector3D )
	set and get the constraint on VTX position More...

const	AmgSymMatrix (3) *Qmatrix() const

const Amg::Vector3D *	Vvector () const

void	setUnconstrained ()

bool	Constrained () const

Private Member Functions
	AmgSymMatrix (3) *m_matrix

	AmgSymMatrix (3) *m_covariance

	AmgSymMatrix (3) *m_qMatrix

Private Attributes
int	m_nTracks

Amg::Vector3D *	m_vector

const xAOD::Vertex *	m_original

Amg::Vector3D *	m_originalPosition

AlignTrackCollection	m_alignTracks

Amg::Vector3D *	m_position

std::vector< AlignModuleVertexDerivatives > *	m_derivatives

Amg::Vector3D *	m_vVector

bool	m_constrained

AlignVertexType	m_type

Detailed Description

Definition at line 39 of file AlignVertex.h.

Member Enumeration Documentation

◆ AlignVertexType

enum Trk::AlignVertex::AlignVertexType

Enumerator
Unknown	default type
Primary	not refitted, just copy constructed from original Track
Refitted	normally refitted, without adding any pseudo-measurement
Accumulated	accumulated by the GX algorithm.
NVertexTypes	total number of track types

Definition at line 43 of file AlignVertex.h.

                         {
         Unknown = 0,               
         Primary = 1,               
         Refitted = 2,              
         Accumulated = 3,           
  
         NVertexTypes               
     };

Constructor & Destructor Documentation

◆ AlignVertex() [1/3]

AlignVertex::AlignVertex ( )

default constructor

Definition at line 11 of file AlignVertex.cxx.

     : m_matrix(nullptr)
     , m_vector(nullptr)
     , m_original(nullptr)
     , m_originalPosition(nullptr)
     , m_position(nullptr)
     , m_covariance(nullptr)
     , m_derivatives(nullptr)
     , m_qMatrix(nullptr)
     , m_vVector(nullptr)
     , m_constrained(false)
     , m_type(Unknown)
   {
     m_vector=new Amg::Vector3D(0.,0.,0.); 
     m_matrix=new AmgSymMatrix(3);
     m_derivatives=new std::vector<AlignModuleVertexDerivatives >(0);
     m_qMatrix= new AmgSymMatrix(3);
     m_vVector=new Amg::Vector3D(0.,0.,0.); 
     m_nTracks=0;
   }

◆ AlignVertex() [2/3]

AlignVertex::AlignVertex ( const xAOD::Vertex * originalVertex )

constructor takes the original vertex candidate

Definition at line 73 of file AlignVertex.cxx.

     : m_matrix(nullptr)
     , m_vector(nullptr)
     , m_original(vertex)
     , m_originalPosition(nullptr)
     , m_position(nullptr)
     , m_covariance(nullptr)
     , m_derivatives(nullptr)
     , m_qMatrix(nullptr)
     , m_vVector(nullptr)
     , m_constrained(false)
     , m_type(Unknown)
   {
     m_vector=new Amg::Vector3D(0.,0.,0.); 
     m_matrix=new AmgSymMatrix(3);
     m_original = vertex;
     m_originalPosition = new Amg::Vector3D( vertex->position() );
     m_derivatives=new std::vector<AlignModuleVertexDerivatives >(0);
     m_qMatrix=new AmgSymMatrix(3);
     m_vVector=new Amg::Vector3D(0.,0.,0.); 
     m_nTracks=0;
   }

◆ AlignVertex() [3/3]

AlignVertex::AlignVertex ( const AlignVertex & rhs )

copy constructor

Definition at line 33 of file AlignVertex.cxx.

     : m_nTracks(rhs.m_nTracks)
     , m_matrix(new AmgSymMatrix(3)(*(rhs.m_matrix)))
     , m_vector(new Amg::Vector3D(*(rhs.m_vector)))
     , m_original(new xAOD::Vertex(*(rhs.m_original)))
     , m_originalPosition(new Amg::Vector3D( *(rhs.m_originalPosition )))
     , m_position(new Amg::Vector3D(*(rhs.m_position)))
     , m_covariance(new AmgSymMatrix(3) (*(rhs.m_covariance)) )
     , m_derivatives(new auto(*(rhs.m_derivatives)))
     , m_qMatrix(new AmgSymMatrix(3) (*(rhs.m_qMatrix)) )
     , m_vVector(new Amg::Vector3D(*(rhs.m_vVector)))
     , m_constrained(rhs.m_constrained)
     , m_type(rhs.m_type)
   {
    
   }

◆ ~AlignVertex()

AlignVertex::~AlignVertex ( )

destructor

Definition at line 98 of file AlignVertex.cxx.

   {
     if( m_matrix )           delete m_matrix;
     if( m_vector )           delete m_vector;
  
     if( m_originalPosition ) delete m_originalPosition;
  
     if( m_position )         delete m_position;
     if( m_covariance )       delete m_covariance;
     if( m_derivatives )      delete m_derivatives;
  
     if( m_qMatrix )          delete m_qMatrix;
     if( m_vVector )          delete m_vVector;
   }

Member Function Documentation

◆ addDerivatives()

void AlignVertex::addDerivatives ( std::vector< AlignModuleVertexDerivatives > * vec )

Definition at line 132 of file AlignVertex.cxx.

   {    
     std::vector<AlignModuleVertexDerivatives>::iterator derivIt     = vec->begin();
     std::vector<AlignModuleVertexDerivatives>::iterator derivIt_end = vec->end();
  
     for ( ; derivIt!=derivIt_end ; ++derivIt) {
       m_derivatives->push_back(*derivIt);
     }
  
      }

◆ alignTrackCollection()

const AlignTrackCollection* Trk::AlignVertex::alignTrackCollection ( ) const

returns collection of alignTracks

◆ AmgSymMatrix() [1/5]

const Trk::AlignVertex::AmgSymMatrix ( 3 ) const

inline

Definition at line 90 of file AlignVertex.h.

90 {return m_covariance;}

◆ AmgSymMatrix() [2/5]

Trk::AlignVertex::AmgSymMatrix ( 3 )

private

◆ AmgSymMatrix() [3/5]

Trk::AlignVertex::AmgSymMatrix ( 3 )

private

◆ AmgSymMatrix() [4/5]

Trk::AlignVertex::AmgSymMatrix ( 3 )

private

◆ AmgSymMatrix() [5/5]

const Trk::AlignVertex::AmgSymMatrix ( 3 ) const

inline

Definition at line 117 of file AlignVertex.h.

117 { return m_qMatrix; }

◆ Constrained()

bool Trk::AlignVertex::Constrained ( ) const

inline

Definition at line 123 of file AlignVertex.h.

123 { return m_constrained; }

◆ derivatives()

const std::vector<AlignModuleVertexDerivatives>* Trk::AlignVertex::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 times W*F. There is one Amg::VectorX in the vector for each alignment parameter of each align module. Set by accumulateVTX.

Definition at line 95 of file AlignVertex.h.

95 { return m_derivatives; }

◆ dump()

void AlignVertex::dump ( MsgStream & msg )

dump align vertex information

Definition at line 218 of file AlignVertex.cxx.

   {
     msg<<"dumping AlignVertex:  "<<endmsg;
  
     msg<<" vertex Position:    "<< (*m_position) << endmsg;
     msg<<" vertex Covariance:  "<< (*m_covariance) << endmsg;
     msg<< endmsg;
  
   }

◆ firstATrack()

AlignTrackCollection::const_iterator Trk::AlignVertex::firstATrack ( ) const

retrieve iterator pointer to first element in collection

◆ fitVertex()

void AlignVertex::fitVertex ( )

fit the vertex internally

Definition at line 144 of file AlignVertex.cxx.

   {
     
     m_type = Unknown;
  
     if( Ntracks() < 2 )   return;    // this is not a good vertex!
  
  
     if( !m_matrix ) {
       std::cout <<"NULL pointer to the matrix! Bailing out..."<< std::endl;
       return;
     }
  
     Amg::Vector3D  vec = (*m_vector);
     AmgSymMatrix(3)   cov = (*m_matrix);
     AmgSymMatrix(3)   covcons = (*m_matrix);
  
     AmgSymMatrix(3)  Qinv;
     if( m_constrained && m_qMatrix->determinant() > 1.0e-24 )   {     // just my guess sigma>0.1 micron ???
       
       bool invertible;
       m_qMatrix->computeInverseWithCheck(Qinv,invertible);
   
       if(!invertible) {
         std::cout <<"fitVertex: Q inversion failed. " << std::endl;
         return;
       }
  
   
       Amg::Vector3D vtemp(3, 0);
       
       vtemp  = *m_originalPosition;
       vtemp -= *m_vVector;
  
       covcons += 2.0*Qinv;
       vec     += 2.0*Qinv*vtemp;
     }
  
     // invert the constrained covariance:
     AmgSymMatrix(3)   invcovcons;
     
     bool invertible;
     covcons.computeInverseWithCheck(invcovcons,invertible);
  
     if(!invertible) {
       std::cout <<"fitVertex: covcons inversion failed. " << std::endl;
       return;
     }
     
     // invert the covariance matrix:
     AmgSymMatrix(3)   invcov;
     cov.computeInverseWithCheck(invcov,invertible);
  
     if(!invertible) {
       std::cout <<"fitVertex: cov inversion failed. " << std::endl;
       return;
     }
  
  
  
     // calculate corrections    (mind the sign!)
     Amg::Vector3D delta(invcovcons * vec);
  
     if ( !m_position ) m_position=new Amg::Vector3D(0.,0.,0.);   
     *m_position = *m_originalPosition;
     *m_position -= delta;
     
  
     if ( !m_covariance ) m_covariance= new AmgSymMatrix(3);   
     (*m_covariance)=invcov;                       // this one is unconstrained!
  
     m_type = Refitted;
  }

◆ incrementMatrix()

void Trk::AlignVertex::incrementMatrix ( const AmgSymMatrix(3) vtxM )

inline

Definition at line 105 of file AlignVertex.h.

105 { (*m_matrix)+=vtxM; }

◆ incrementVector()

void Trk::AlignVertex::incrementVector ( const Amg::Vector3D & vtxV )

inline

increment algebra objects for this verterx:

Definition at line 104 of file AlignVertex.h.

104 { (*m_vVector)+=vtxV; }

◆ lastATrack()

AlignTrackCollection::const_iterator Trk::AlignVertex::lastATrack ( ) const

returns iterator pointer to last element in collection

◆ Ntracks()

int Trk::AlignVertex::Ntracks ( ) const

inline

get the number of contributing tracks

Definition at line 112 of file AlignVertex.h.

112 { return m_nTracks; }

◆ operator=()

AlignVertex & AlignVertex::operator= ( const AlignVertex & rhs )

assignment

Definition at line 51 of file AlignVertex.cxx.

                                                                {
     if (&rhs!=this){
       m_nTracks = rhs.m_nTracks;
       delete m_matrix; m_matrix = new AmgSymMatrix(3)(*(rhs.m_matrix));
       delete m_vector; m_vector = new Amg::Vector3D(*(rhs.m_vector));
       delete m_original; m_original = new xAOD::Vertex(*(rhs.m_original));
         delete m_originalPosition; m_originalPosition = new Amg::Vector3D( *(rhs.m_originalPosition ));
         delete m_position; m_position = new Amg::Vector3D(*(rhs.m_position));
         delete m_covariance; m_covariance = new AmgSymMatrix(3) (*(rhs.m_covariance));
       delete m_derivatives; m_derivatives = new auto(*(rhs.m_derivatives));
       delete m_qMatrix; m_qMatrix = new AmgSymMatrix(3) (*(rhs.m_qMatrix));
       delete m_vVector; m_vVector = new Amg::Vector3D(*(rhs.m_vVector));
       m_constrained = rhs.m_constrained;
       m_type = rhs.m_type;
     }
     return *this;
   }

◆ originalPosition()

const Amg::Vector3D* Trk::AlignVertex::originalPosition ( ) const

inline

Definition at line 79 of file AlignVertex.h.

79 {return m_originalPosition;}

◆ originalVertex()

const xAOD::Vertex* Trk::AlignVertex::originalVertex ( ) const

inline

retrieve pointer to original vertex and its position

Definition at line 78 of file AlignVertex.h.

78 { return m_original; }

◆ position()

const Amg::Vector3D* Trk::AlignVertex::position ( ) const

inline

get the vertex position and covariance

Definition at line 89 of file AlignVertex.h.

89 {return m_position;}

◆ setConstraint()

void AlignVertex::setConstraint	(	AmgSymMatrix(3) *	Q,
		Amg::Vector3D *	V
	)

set and get the constraint on VTX position

Definition at line 114 of file AlignVertex.cxx.

   {
     if( Q->rows()!=3 || V->rows()!=3 ) {
       std::cout <<"Unexpected size of the constraint:"<< std::endl;
       std::cout <<"  Q size = "<< Q->rows() << ",   V size = "<< V->rows() << std::endl;
       return;
     }
     if( Q->determinant() < 1.0e-24 ) {
       std::cout <<" Suspicious Q determinant:  detQ = "<< Q->determinant() << std::endl;
       return;
     }
     (*m_qMatrix) = (*Q);
     (*m_vVector) = (*V);
     m_constrained = true;
   }

◆ setDerivatives()

void Trk::AlignVertex::setDerivatives ( std::vector< AlignModuleVertexDerivatives > * vec )

inline

Definition at line 96 of file AlignVertex.h.

96 { m_derivatives=vec; }

◆ setOriginalVertex()

void Trk::AlignVertex::setOriginalVertex ( const xAOD::Vertex * vertex )

inline

set pointer to original vertex

Definition at line 82 of file AlignVertex.h.

82 { m_original = vertex; }

◆ setType()

void Trk::AlignVertex::setType ( AlignVertexType type )

inline

Definition at line 86 of file AlignVertex.h.

86 {m_type = type;}

◆ setUnconstrained()

void Trk::AlignVertex::setUnconstrained ( )

inline

Definition at line 121 of file AlignVertex.h.

121 { m_constrained=false; }

◆ type()

AlignVertexType Trk::AlignVertex::type ( ) const

inline

get and set the refit type

Definition at line 85 of file AlignVertex.h.

85 {return m_type;}

◆ Vvector()

const Amg::Vector3D* Trk::AlignVertex::Vvector ( ) const

inline

Definition at line 119 of file AlignVertex.h.

119 { return m_vVector; }

Member Data Documentation

◆ m_alignTracks

AlignTrackCollection Trk::AlignVertex::m_alignTracks

private

Definition at line 135 of file AlignVertex.h.

◆ m_constrained

bool Trk::AlignVertex::m_constrained

private

Definition at line 144 of file AlignVertex.h.

◆ m_derivatives

std::vector<AlignModuleVertexDerivatives >* Trk::AlignVertex::m_derivatives

private

Definition at line 140 of file AlignVertex.h.

◆ m_nTracks

int Trk::AlignVertex::m_nTracks

private

Definition at line 127 of file AlignVertex.h.

◆ m_original

const xAOD::Vertex* Trk::AlignVertex::m_original

private

Definition at line 132 of file AlignVertex.h.

◆ m_originalPosition

Amg::Vector3D* Trk::AlignVertex::m_originalPosition

private

Definition at line 133 of file AlignVertex.h.

◆ m_position

Amg::Vector3D* Trk::AlignVertex::m_position

private

Definition at line 137 of file AlignVertex.h.

◆ m_type

AlignVertexType Trk::AlignVertex::m_type

private

Definition at line 147 of file AlignVertex.h.

◆ m_vector

Amg::Vector3D* Trk::AlignVertex::m_vector

private

Definition at line 130 of file AlignVertex.h.

◆ m_vVector

Amg::Vector3D* Trk::AlignVertex::m_vVector

private

Definition at line 143 of file AlignVertex.h.

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

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ AlignVertexType

Constructor & Destructor Documentation

◆ AlignVertex() [1/3]

◆ AlignVertex() [2/3]

◆ AlignVertex() [3/3]

◆ ~AlignVertex()

Member Function Documentation

◆ addDerivatives()

◆ alignTrackCollection()

◆ AmgSymMatrix() [1/5]

◆ AmgSymMatrix() [2/5]

◆ AmgSymMatrix() [3/5]

◆ AmgSymMatrix() [4/5]

◆ AmgSymMatrix() [5/5]

◆ Constrained()

◆ derivatives()

◆ dump()

◆ firstATrack()

◆ fitVertex()

◆ incrementMatrix()

◆ incrementVector()

◆ lastATrack()

◆ Ntracks()

◆ operator=()

◆ originalPosition()

◆ originalVertex()

◆ position()

◆ setConstraint()

◆ setDerivatives()

◆ setOriginalVertex()

◆ setType()

◆ setUnconstrained()

◆ type()

◆ Vvector()

Member Data Documentation

◆ m_alignTracks

◆ m_constrained

◆ m_derivatives

◆ m_nTracks

◆ m_original

◆ m_originalPosition

◆ m_position

◆ m_type

◆ m_vector

◆ m_vVector