#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
	AlignVertex (const xAOD::Vertex *originalVertex)
	constructor takes the original vertex candidate but doesn't own it
	AlignVertex (const AlignVertex &rhs)
	copy constructor
AlignVertex &	operator= (const AlignVertex &rhs)
	assignment
	~AlignVertex ()
	destructor
const xAOD::Vertex *	originalVertex () const
	retrieve pointer to original vertex and its position
const Amg::Vector3D *	originalPosition () const
void	setOriginalVertex (const xAOD::Vertex *vertex)
	set pointer to original vertex
AlignVertexType	type () const
	get and set the refit type
void	setType (AlignVertexType type)
const Amg::Vector3D *	position () const
	get the vertex position and covariance
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.
void	setDerivatives (std::vector< AlignModuleVertexDerivatives > *vecvec)
void	addDerivatives (std::vector< AlignModuleVertexDerivatives > *vecvec)
void	dump (MsgStream &msg)
	dump align vertex information
void	incrementVector (const Amg::Vector3D &vtxV)
	increment algebra objects for this verterx:
void	incrementMatrix (const AmgSymMatrix(3) vtxM)
void	fitVertex ()
	fit the vertex internally
int	Ntracks () const
	get the number of contributing tracks
void	setConstraint (AmgSymMatrix(3) , Amg::Vector3D )
	set and get the constraint on VTX position
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
Amg::Vector3D *	m_position
std::vector< AlignModuleVertexDerivatives > *	m_derivatives
Amg::Vector3D *	m_vVector
bool	m_constrained
AlignVertexType	m_type

Detailed Description

Definition at line 34 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 38 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 12 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 but doesn't own it

Definition at line 74 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 34 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(rhs.m_original) //non-owning
    , 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 99 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 133 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);
    }
 
     }

◆ AmgSymMatrix() [1/5]

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

inline

Definition at line 77 of file AlignVertex.h.

77{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 103 of file AlignVertex.h.

103{ return m_qMatrix; }

◆ Constrained()

bool Trk::AlignVertex::Constrained ( ) const

inline

Definition at line 109 of file AlignVertex.h.

109{ 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 82 of file AlignVertex.h.

82{ return m_derivatives; }

◆ dump()

void AlignVertex::dump ( MsgStream & msg )

dump align vertex information

Definition at line 219 of file AlignVertex.cxx.

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

◆ fitVertex()

void AlignVertex::fitVertex ( )

fit the vertex internally

Definition at line 145 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 91 of file AlignVertex.h.

91{ (*m_matrix)+=vtxM; }

◆ incrementVector()

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

inline

increment algebra objects for this verterx:

Definition at line 90 of file AlignVertex.h.

90{ (*m_vVector)+=vtxV; }

◆ Ntracks()

int Trk::AlignVertex::Ntracks ( ) const

inline

get the number of contributing tracks

Definition at line 98 of file AlignVertex.h.

98{ return m_nTracks; }

◆ operator=()

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

assignment

Definition at line 52 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));
      m_original = rhs.m_original; //non-owning
        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 66 of file AlignVertex.h.

66{return m_originalPosition;}

◆ originalVertex()

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

inline

retrieve pointer to original vertex and its position

Definition at line 65 of file AlignVertex.h.

65{ return m_original; }

◆ position()

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

inline

get the vertex position and covariance

Definition at line 76 of file AlignVertex.h.

76{return m_position;}

◆ setConstraint()

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

set and get the constraint on VTX position

Definition at line 115 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 83 of file AlignVertex.h.

83{ m_derivatives=vec; }

◆ setOriginalVertex()

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

inline

set pointer to original vertex

Definition at line 69 of file AlignVertex.h.

69{ m_original = vertex; }

◆ setType()

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

inline

Definition at line 73 of file AlignVertex.h.

73{m_type = type;}

Trk::AlignVertex::type

AlignVertexType type() const

get and set the refit type

Definition AlignVertex.h:72

◆ setUnconstrained()

void Trk::AlignVertex::setUnconstrained ( )

inline

Definition at line 107 of file AlignVertex.h.

107{ m_constrained=false; }

◆ type()

AlignVertexType Trk::AlignVertex::type ( ) const

inline

get and set the refit type

Definition at line 72 of file AlignVertex.h.

72{return m_type;}

◆ Vvector()

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

inline

Definition at line 105 of file AlignVertex.h.

105{ return m_vVector; }

Member Data Documentation

◆ m_constrained

bool Trk::AlignVertex::m_constrained

private

Definition at line 129 of file AlignVertex.h.

◆ m_derivatives

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

private

Definition at line 125 of file AlignVertex.h.

◆ m_nTracks

int Trk::AlignVertex::m_nTracks

private

Definition at line 113 of file AlignVertex.h.

◆ m_original

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

private

Definition at line 118 of file AlignVertex.h.

◆ m_originalPosition

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

private

Definition at line 119 of file AlignVertex.h.

◆ m_position

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

private

Definition at line 122 of file AlignVertex.h.

◆ m_type

AlignVertexType Trk::AlignVertex::m_type

private

Definition at line 132 of file AlignVertex.h.

◆ m_vector

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

private

Definition at line 116 of file AlignVertex.h.

◆ m_vVector

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

private

Definition at line 128 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()

◆ AmgSymMatrix() [1/5]

◆ AmgSymMatrix() [2/5]

◆ AmgSymMatrix() [3/5]

◆ AmgSymMatrix() [4/5]

◆ AmgSymMatrix() [5/5]

◆ Constrained()

◆ derivatives()

◆ dump()

◆ fitVertex()

◆ incrementMatrix()

◆ incrementVector()

◆ Ntracks()

◆ operator=()

◆ originalPosition()

◆ originalVertex()

◆ position()

◆ setConstraint()

◆ setDerivatives()

◆ setOriginalVertex()

◆ setType()

◆ setUnconstrained()

◆ type()

◆ Vvector()

Member Data Documentation

◆ m_constrained

◆ m_derivatives

◆ m_nTracks

◆ m_original

◆ m_originalPosition

◆ m_position

◆ m_type

◆ m_vector

◆ m_vVector