#include <VertexKinematicConstraintT.h>

Inheritance diagram for VertexKinematicConstraintT:

Public Member Functions
VertexKinematicConstraintT *	clone () const override

void	init (const std::vector< KinematicState > &states, const GlobalPoint &point, const GlobalVector &mf) override

int	numberOfEquations () const override

	VertexKinematicConstraintT ()

	~VertexKinematicConstraintT () override

Public Member Functions inherited from MultiTrackKinematicConstraintT< 2, 4 >
parametersDerivativeType const &	parametersDerivative () const

positionDerivativeType const &	positionDerivative () const

valueType const &	value () const

	~MultiTrackKinematicConstraintT () override

Public Member Functions inherited from MultiTrackKinematicConstraintBaseT
virtual	~MultiTrackKinematicConstraintBaseT ()

Private Types
typedef MultiTrackKinematicConstraintT< 2, 4 >	super

Private Member Functions
void	fillParametersDerivative () const override

void	fillPositionDerivative () const override

void	fillValue () const override

Private Attributes
double	a_i [2]

double	delta [2]

GlobalVector	dpos [2]

double	k [2]

double	m [2]

GlobalVector	mom [2]

double	n [2]

double	novera [2]

Additional Inherited Members
Public Types inherited from MultiTrackKinematicConstraintT< 2, 4 >
enum

typedef ROOT::Math::SMatrix< double, DIM, 7 *NTRK >	parametersDerivativeType

typedef ROOT::Math::SMatrix< double, DIM, 3 >	positionDerivativeType

typedef MultiTrackKinematicConstraintT< NTRK, DIM >	self

typedef ROOT::Math::SVector< double, DIM >	valueType

Protected Member Functions inherited from MultiTrackKinematicConstraintT< 2, 4 >
parametersDerivativeType &	jac_d () const

double &	jac_d (size_t i, size_t j) const

positionDerivativeType &	jac_e () const

double &	jac_e (size_t i, size_t j) const

valueType &	vl () const

double &	vl (size_t i) const

Detailed Description

Class implementing the vertexing constraint for extended cartesian parametrization (x,y,z,p_x,p_y,p_z,m). The equations and derivatives in general follow the P.Avery's "Applied Fitting Theory-VI" CBX 98-37

Definition at line 14 of file VertexKinematicConstraintT.h.

Member Typedef Documentation

◆ super

typedef MultiTrackKinematicConstraintT<2, 4> VertexKinematicConstraintT::super

private

Definition at line 16 of file VertexKinematicConstraintT.h.

Constructor & Destructor Documentation

◆ VertexKinematicConstraintT()

VertexKinematicConstraintT::VertexKinematicConstraintT ( )

Definition at line 5 of file VertexKinematicConstraintT.cc.

5 {}

Referenced by clone().

◆ ~VertexKinematicConstraintT()

VertexKinematicConstraintT::~VertexKinematicConstraintT ( )

override

Definition at line 7 of file VertexKinematicConstraintT.cc.

7 {}

Member Function Documentation

◆ clone()

VertexKinematicConstraintT* VertexKinematicConstraintT::clone ( void ) const

inlineoverridevirtual

Implements MultiTrackKinematicConstraintBaseT.

Definition at line 38 of file VertexKinematicConstraintT.h.

38 { return new VertexKinematicConstraintT(*this); }

References VertexKinematicConstraintT().

◆ fillParametersDerivative()

void VertexKinematicConstraintT::fillParametersDerivative ( ) const

overrideprivatevirtual

fills a matrix of derivatives of constraint equations w.r.t. particle parameters

Implements MultiTrackKinematicConstraintT< 2, 4 >.

Definition at line 54 of file VertexKinematicConstraintT.cc.

                                                                 {
   ROOT::Math::SMatrix<double, 2, 7> el_part_d;
   for (int j = 0; j != 2; ++j) {
     if (a_i[j] != 0) {
       //charged particle
  
       //D Jacobian matrix
       el_part_d(0, 0) = mom[j].y() + a_i[j] * dpos[j].x();
       el_part_d(0, 1) = -mom[j].x() + a_i[j] * dpos[j].y();
       el_part_d(1, 0) = -k[j] * (m[j] * mom[j].x() - n[j] * mom[j].y());
       el_part_d(1, 1) = -k[j] * (m[j] * mom[j].y() + n[j] * mom[j].x());
       el_part_d(1, 2) = -1.;
       el_part_d(0, 3) = dpos[j].y();
       el_part_d(0, 4) = -dpos[j].x();
       el_part_d(1, 3) = k[j] * (m[j] * dpos[j].x() - novera[j] * (2 * mom[j].x() - a_i[j] * dpos[j].y()));
       el_part_d(1, 4) = k[j] * (m[j] * dpos[j].y() - novera[j] * (2 * mom[j].y() + a_i[j] * dpos[j].x()));
       el_part_d(1, 5) = -delta[j] / a_i[j];
       super::jac_d().Place_at(el_part_d, j * 2, j * 7);
     } else {
       //neutral particle
       double pt2Inverse = 1. / mom[j].perp2();
       el_part_d(0, 0) = mom[j].y();
       el_part_d(0, 1) = -mom[j].x();
       el_part_d(1, 0) = mom[j].x() * (mom[j].z() * pt2Inverse);
       el_part_d(1, 1) = mom[j].y() * (mom[j].z() * pt2Inverse);
       el_part_d(1, 2) = -1.;
       el_part_d(0, 3) = dpos[j].y();
       el_part_d(0, 4) = -dpos[j].x();
       el_part_d(1, 3) = 2 * (dpos[j].x() * mom[j].x() + dpos[j].y() * mom[j].y()) * pt2Inverse * mom[j].x() *
                             (mom[j].z() * pt2Inverse) -
                         dpos[j].x() * (mom[j].z() * pt2Inverse);
       el_part_d(1, 4) = 2 * (dpos[j].x() * mom[j].x() + dpos[j].y() * mom[j].y()) * pt2Inverse * mom[j].y() *
                             (mom[j].z() * pt2Inverse) -
                         dpos[j].x() * (mom[j].z() * pt2Inverse);
       el_part_d(1, 5) = -(dpos[j].x() * mom[j].x() + dpos[j].y() * mom[j].y()) * pt2Inverse;
       super::jac_d().Place_at(el_part_d, j * 2, j * 7);
     }
   }
 }

References a_i, delta, dpos, dqmiolumiharvest::j, MultiTrackKinematicConstraintT< 2, 4 >::jac_d(), k, m, mom, n, novera, PV3DBase< T, PVType, FrameType >::perp2(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ fillPositionDerivative()

void VertexKinematicConstraintT::fillPositionDerivative ( ) const

overrideprivatevirtual

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

Implements MultiTrackKinematicConstraintT< 2, 4 >.

Definition at line 94 of file VertexKinematicConstraintT.cc.

                                                               {
   ROOT::Math::SMatrix<double, 2, 3> el_part_e;
   for (int j = 0; j != 2; ++j) {
     if (a_i[j] != 0) {
       //charged particle
  
       //E jacobian matrix
       el_part_e(0, 0) = -(mom[j].y() + a_i[j] * dpos[j].x());
       el_part_e(0, 1) = mom[j].x() - a_i[j] * dpos[j].y();
       el_part_e(1, 0) = k[j] * (m[j] * mom[j].x() - n[j] * mom[j].y());
       el_part_e(1, 1) = k[j] * (m[j] * mom[j].y() + n[j] * mom[j].x());
       el_part_e(1, 2) = 1;
       super::jac_e().Place_at(el_part_e, 2 * j, 0);
     } else {
       //neutral particle
       double pt2Inverse = 1. / mom[j].perp2();
       el_part_e(0, 0) = -mom[j].y();
       el_part_e(0, 1) = mom[j].x();
       el_part_e(1, 0) = -mom[j].x() * mom[j].z() * pt2Inverse;
       el_part_e(1, 1) = -mom[j].y() * mom[j].z() * pt2Inverse;
       el_part_e(1, 2) = 1;
       super::jac_e().Place_at(el_part_e, 2 * j, 0);
     }
   }
 }

References a_i, dpos, dqmiolumiharvest::j, MultiTrackKinematicConstraintT< 2, 4 >::jac_e(), k, m, mom, n, PV3DBase< T, PVType, FrameType >::perp2(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ fillValue()

void VertexKinematicConstraintT::fillValue ( ) const

overrideprivatevirtual

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

Implements MultiTrackKinematicConstraintT< 2, 4 >.

Definition at line 36 of file VertexKinematicConstraintT.cc.

                                                  {
   //it is 2 equations per track
   for (int j = 0; j != 2; ++j) {
     if (a_i[j] != 0) {
       //vector of values
       super::vl(j * 2) = dpos[j].y() * mom[j].x() - dpos[j].x() * mom[j].y() -
                          a_i[j] * (dpos[j].x() * dpos[j].x() + dpos[j].y() * dpos[j].y()) * 0.5;
       super::vl(j * 2 + 1) = dpos[j].z() - mom[j].z() * delta[j] / a_i[j];
     } else {
       //neutral particle
       double pt2Inverse = 1. / mom[j].perp2();
       super::vl(j * 2) = dpos[j].y() * mom[j].x() - dpos[j].x() * mom[j].y();
       super::vl(j * 2 + 1) =
           dpos[j].z() - mom[j].z() * (dpos[j].x() * mom[j].x() + dpos[j].y() * mom[j].y()) * pt2Inverse;
     }
   }
 }

References a_i, delta, dpos, dqmiolumiharvest::j, mom, PV3DBase< T, PVType, FrameType >::perp2(), MultiTrackKinematicConstraintT< 2, 4 >::vl(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ init()

void VertexKinematicConstraintT::init	(	const std::vector< KinematicState > &	states,
		const GlobalPoint &	point,
		const GlobalVector &	mf
	)

overridevirtual

Implements MultiTrackKinematicConstraintBaseT.

Definition at line 9 of file VertexKinematicConstraintT.cc.

                                                                       {
   int num = states.size();
   if (num != 2)
     throw VertexException("VertexKinematicConstraintT !=2 states passed");
   double mfz = fieldValue.z();
  
   int j = 0;
   for (std::vector<KinematicState>::const_iterator i = states.begin(); i != states.end(); i++) {
     mom[j] = i->globalMomentum();
     dpos[j] = ipoint - i->globalPosition();
     a_i[j] = -i->particleCharge() * mfz;
  
     double pvx = mom[j].x() - a_i[j] * dpos[j].y();
     double pvy = mom[j].y() + a_i[j] * dpos[j].x();
     double pvt2 = pvx * pvx + pvy * pvy;
     novera[j] = (dpos[j].x() * mom[j].x() + dpos[j].y() * mom[j].y());
     n[j] = a_i[j] * novera[j];
     m[j] = (pvx * mom[j].x() + pvy * mom[j].y());
     k[j] = -mom[j].z() / (mom[j].perp2() * pvt2);
     delta[j] = std::atan2(n[j], m[j]);
  
     ++j;
   }
 }