#include <CombinedKinematicConstraint.h>

Inheritance diagram for CombinedKinematicConstraint:

Public Member Functions
CombinedKinematicConstraint *	clone () const override

	CombinedKinematicConstraint (const std::vector< MultiTrackKinematicConstraint *> &constraintVector)

int	numberOfEquations () const override

AlgebraicMatrix	parametersDerivative (const std::vector< KinematicState > &states, const GlobalPoint &point) const override

AlgebraicMatrix	positionDerivative (const std::vector< KinematicState > &states, const GlobalPoint &point) const override

AlgebraicVector	value (const std::vector< KinematicState > &states, const GlobalPoint &point) const override

Public Member Functions inherited from MultiTrackKinematicConstraint
	MultiTrackKinematicConstraint ()

virtual	~MultiTrackKinematicConstraint ()

Private Attributes
std::vector< MultiTrackKinematicConstraint * >	constraints

Detailed Description

This class combines several user defined constraints (by expanding the vector d and the matrices D and E). Usage: Add each constraint to a std::vector<MultiTrackKinematicConstraint* >. This vector has to be used in the constructor: MultiTrackKinematicConstraint combiC = new CombinedKinematicConstraint(std::vector<MultiTrackKinematicConstraint >) The produced object can be used by KinematicConstrainedVertexFitter.fit()

Lars Perchalla, Philip Sauerland, Dec 2009

Definition at line 20 of file CombinedKinematicConstraint.h.

Constructor & Destructor Documentation

◆ CombinedKinematicConstraint()

CombinedKinematicConstraint::CombinedKinematicConstraint ( const std::vector< MultiTrackKinematicConstraint *> & constraintVector )

inline

Definition at line 22 of file CombinedKinematicConstraint.h.

References constraints.

Referenced by clone().

       : constraints(constraintVector) {
     if (constraints.empty())
       throw VertexException("CombinedKinematicConstraint::<1 constraints passed.");
   }

Member Function Documentation

◆ clone()

CombinedKinematicConstraint* CombinedKinematicConstraint::clone ( void ) const

inlineoverridevirtual

Implements MultiTrackKinematicConstraint.

Definition at line 56 of file CombinedKinematicConstraint.h.

References CombinedKinematicConstraint().

56 { return new CombinedKinematicConstraint(*this); }

CombinedKinematicConstraint::CombinedKinematicConstraint

CombinedKinematicConstraint(const std::vector< MultiTrackKinematicConstraint *> &constraintVector)

Definition: CombinedKinematicConstraint.h:22

◆ numberOfEquations()

int CombinedKinematicConstraint::numberOfEquations ( ) const

overridevirtual

Number of equations per track used for the combined fit

Implements MultiTrackKinematicConstraint.

Definition at line 66 of file CombinedKinematicConstraint.cc.

References constraints, and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

                                                          {
   int noEq = 0;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it)
     noEq += (*it)->numberOfEquations();
 
   return noEq;
 }

◆ parametersDerivative()

AlgebraicMatrix CombinedKinematicConstraint::parametersDerivative	(	const std::vector< KinematicState > &	states,
		const GlobalPoint &	point
	)		const

overridevirtual

Returns a matrix of derivatives of the combined constraint equations w.r.t. particle parameters

Implements MultiTrackKinematicConstraint.

Definition at line 24 of file CombinedKinematicConstraint.cc.

References constraints, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeMuonMisalignmentScenario::matrix, and point.

                                                                                                   {
   AlgebraicMatrix tmpMatrix;
   int row = 0;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it) {
     tmpMatrix = (*it)->parametersDerivative(states, point);
     row += tmpMatrix.num_row();
   }
   AlgebraicMatrix matrix(row, 7 * states.size(), 0);
   int posRow = 1;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it) {
     tmpMatrix = (*it)->parametersDerivative(states, point);
     matrix.sub(posRow, 1, tmpMatrix);
     posRow += tmpMatrix.num_row();
   }
 
   return matrix;
 }

◆ positionDerivative()

AlgebraicMatrix CombinedKinematicConstraint::positionDerivative	(	const std::vector< KinematicState > &	states,
		const GlobalPoint &	point
	)		const

overridevirtual

Returns a matrix of derivatives of constraint equations w.r.t. vertex position

Implements MultiTrackKinematicConstraint.

Definition at line 45 of file CombinedKinematicConstraint.cc.

References constraints, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, makeMuonMisalignmentScenario::matrix, and point.

                                                                                                 {
   AlgebraicMatrix tmpMatrix;
   int row = 0;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it) {
     tmpMatrix = (*it)->positionDerivative(states, point);
     row += tmpMatrix.num_row();
   }
   AlgebraicMatrix matrix(row, 3, 0);
   int posRow = 1;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it) {
     tmpMatrix = (*it)->positionDerivative(states, point);
     matrix.sub(posRow, 1, tmpMatrix);
     posRow += tmpMatrix.num_row();
   }
 
   return matrix;
 }

◆ value()

AlgebraicVector CombinedKinematicConstraint::value	(	const std::vector< KinematicState > &	states,
		const GlobalPoint &	point
	)		const

overridevirtual

Returns a vector of values of the combined constraint equations at the point where the input particles are defined.

Implements MultiTrackKinematicConstraint.

Definition at line 3 of file CombinedKinematicConstraint.cc.

References constraints, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, point, position, findQualityFiles::size, and contentValuesCheck::values.

                                                                                    {
   AlgebraicVector tmpValue;
   int size = 0;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it) {
     tmpValue = (*it)->value(states, point);
     size += tmpValue.num_row();
   }
   AlgebraicVector values(size);
   int position = 1;
   for (std::vector<MultiTrackKinematicConstraint*>::const_iterator it = constraints.begin(); it != constraints.end();
        ++it) {
     tmpValue = (*it)->value(states, point);
     values.sub(position, tmpValue);
     position += tmpValue.num_row();
   }
 
   return values;
 }