#include <MomentumKinematicConstraint.h>

Inheritance diagram for MomentumKinematicConstraint:

Public Member Functions
MomentumKinematicConstraint *	clone () const override

std::pair< AlgebraicMatrix, AlgebraicVector >	derivative (const AlgebraicVector &exPoint) const override

std::pair< AlgebraicMatrix, AlgebraicVector >	derivative (const std::vector< RefCountedKinematicParticle > &par) const override

AlgebraicVector	deviations (int nStates) const override

	MomentumKinematicConstraint (const AlgebraicVector &momentum, const AlgebraicVector &dev)

int	numberOfEquations () const override

std::pair< AlgebraicVector, AlgebraicVector >	value (const AlgebraicVector &exPoint) const override

std::pair< AlgebraicVector, AlgebraicVector >	value (const std::vector< RefCountedKinematicParticle > &par) const override

Public Member Functions inherited from KinematicConstraint
	KinematicConstraint ()

virtual	~KinematicConstraint ()

Private Attributes
AlgebraicVector	dd

AlgebraicVector	mm

Detailed Description

Class constraining total 3-momentum of the particle (p_x,p_y,p_z) This constraint can not be implemented on multiple particles, without fitting the vertex. Current version supports one state refit only.

Kirill Prokofiev, October 2003 MultiState version: July 2004

Definition at line 15 of file MomentumKinematicConstraint.h.

Constructor & Destructor Documentation

◆ MomentumKinematicConstraint()

MomentumKinematicConstraint::MomentumKinematicConstraint	(	const AlgebraicVector &	momentum,
		const AlgebraicVector &	dev
	)

Constructor with the 4-momentum vector as an argument

Definition at line 4 of file MomentumKinematicConstraint.cc.

References dd, and mm.

Referenced by clone().

                                                                                                                     {
   if ((momentum.num_row() != 3) || (dev.num_row() != 3))
     throw VertexException("MomentumKinemticConstraint::Momentum or Deviation vector passed is not 3-dimensional");
   mm = momentum;
   AlgebraicVector dev_l(7, 0);
   dev_l(4) = dev(1) * dev(1);
   dev_l(5) = dev(2) * dev(2);
   dev_l(6) = dev(3) * dev(3);
   dd = dev_l;
 }

Member Function Documentation

◆ clone()

MomentumKinematicConstraint* MomentumKinematicConstraint::clone ( ) const

inlineoverridevirtual

Clone method

Implements KinematicConstraint.

Definition at line 48 of file MomentumKinematicConstraint.h.

References MomentumKinematicConstraint().

48 { return new MomentumKinematicConstraint(*this); }

MomentumKinematicConstraint::MomentumKinematicConstraint

MomentumKinematicConstraint(const AlgebraicVector &momentum, const AlgebraicVector &dev)

Definition: MomentumKinematicConstraint.cc:4

◆ derivative() [1/2]

std::pair< AlgebraicMatrix, AlgebraicVector > MomentumKinematicConstraint::derivative ( const AlgebraicVector & exPoint ) const

overridevirtual

Implements KinematicConstraint.

Definition at line 34 of file MomentumKinematicConstraint.cc.

References flavorHistoryFilter_cfi::dr.

                                           {
   if (exPoint.num_row() == 0)
     throw VertexException("MomentumKinematicConstraint::derivative requested for zero Linearization point");
 
   //security check for extended cartesian parametrization
   int inSize = exPoint.num_row();
   if ((inSize % 7) != 0)
     throw VertexException("MomentumKinematicConstraint::linearization point has a wrong dimension");
   int nStates = inSize / 7;
   if (nStates != 1)
     throw VertexException("MomentumKinematicConstraint::Multistate refit is not foreseen for this constraint");
 
   AlgebraicVector pr = exPoint;
   AlgebraicMatrix dr(3, 7, 0);
   dr(1, 4) = 1.;
   dr(2, 5) = 1.;
   dr(3, 6) = 1.;
   return std::pair<AlgebraicMatrix, AlgebraicVector>(dr, pr);
 }

◆ derivative() [2/2]

std::pair< AlgebraicMatrix, AlgebraicVector > MomentumKinematicConstraint::derivative ( const std::vector< RefCountedKinematicParticle > & par ) const

overridevirtual

Vector of values and matrix of derivatives calculated using current state parameters as expansion point

Implements KinematicConstraint.

Definition at line 70 of file MomentumKinematicConstraint.cc.

References flavorHistoryFilter_cfi::dr, and point.

                                                              {
   int nStates = par.size();
   if (nStates == 0)
     throw VertexException("MomentumKinematicConstraint::Empty vector of particles passed");
   if (nStates != 1)
     throw VertexException("MomentumKinematicConstraint::Multistate refit is not foreseen for this constraint");
   AlgebraicVector point = asHepVector<7>(par.front()->currentState().kinematicParameters().vector());
   AlgebraicMatrix dr(3, 7, 0);
   dr(1, 4) = 1.;
   dr(2, 5) = 1.;
   dr(3, 6) = 1.;
   return std::pair<AlgebraicMatrix, AlgebraicVector>(dr, point);
 }

◆ deviations()

AlgebraicVector MomentumKinematicConstraint::deviations ( int nStates ) const

overridevirtual

Returns vector of sigma squared associated to the KinematicParameters of refitted particles Initial deviations are given by user for the constraining parameters (mass, momentum components etc). In case of multiple states exactly the same values are added to every particle parameters

Implements KinematicConstraint.

Definition at line 85 of file MomentumKinematicConstraint.cc.

References dd.

                                                                          {
   if (nStates == 0)
     throw VertexException("MomentumKinematicConstraint::Empty vector of particles passed");
   if (nStates != 1)
     throw VertexException("MomentumKinematicConstraint::Multistate refit is not foreseen for this constraint");
   AlgebraicVector res = dd;
   return res;
 }

◆ numberOfEquations()

int MomentumKinematicConstraint::numberOfEquations ( ) const

overridevirtual

Returns number of constraint equations used for fitting. Method is relevant for proper NDF calculations.

Implements KinematicConstraint.

Definition at line 94 of file MomentumKinematicConstraint.cc.

94 { return 3; }

◆ value() [1/2]

std::pair< AlgebraicVector, AlgebraicVector > MomentumKinematicConstraint::value ( const AlgebraicVector & exPoint ) const

overridevirtual

Vector of values and matrix of derivatives calculated at given expansion 7xNumberOfStates point

Implements KinematicConstraint.

Definition at line 15 of file MomentumKinematicConstraint.cc.

References mm.

                                                                                                                  {
   if (exPoint.num_row() == 0)
     throw VertexException("MomentumKinematicConstraint::value requested for zero Linearization point");
 
   //security check for extended cartesian parametrization
   int inSize = exPoint.num_row();
   if ((inSize % 7) != 0)
     throw VertexException("MomentumKinematicConstraint::linearization point has a wrong dimension");
   int nStates = inSize / 7;
   if (nStates != 1)
     throw VertexException("MomentumKinematicConstraint::Multistate refit is not foreseen for this constraint");
   AlgebraicVector pr = exPoint;
   AlgebraicVector vl(3, 0);
   vl(1) = pr(4) - mm(1);
   vl(2) = pr(5) - mm(2);
   vl(3) = pr(6) - mm(3);
   return std::pair<AlgebraicVector, AlgebraicVector>(vl, pr);
 }

◆ value() [2/2]

std::pair< AlgebraicVector, AlgebraicVector > MomentumKinematicConstraint::value ( const std::vector< RefCountedKinematicParticle > & par ) const

overridevirtual

Methods making value and derivative matrix using current state parameters as expansion 7-point. Constraint can be made equaly for single and multiple states

Implements KinematicConstraint.

Definition at line 55 of file MomentumKinematicConstraint.cc.

References mm, and point.

                                                              {
   int nStates = par.size();
   if (nStates == 0)
     throw VertexException("MomentumKinematicConstraint::Empty vector of particles passed");
   if (nStates != 1)
     throw VertexException("MomentumKinematicConstraint::Multistate refit is not foreseen for this constraint");
   AlgebraicVector point = asHepVector<7>(par.front()->currentState().kinematicParameters().vector());
   AlgebraicVector vl(3, 0);
   vl(1) = point(4) - mm(1);
   vl(2) = point(5) - mm(2);
   vl(3) = point(6) - mm(3);
   return std::pair<AlgebraicVector, AlgebraicVector>(vl, point);
 }