#include <ColinearityKinematicConstraintT.h>

Inheritance diagram for ColinearityKinematicConstraintT< Dim >:

Public Types
typedef MultiTrackKinematicConstraintT< 2, int(Dim)>	super

Public Types inherited from MultiTrackKinematicConstraintT< 2, int(Dim)>
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

Public Member Functions
ColinearityKinematicConstraintT< Dim > *	clone () const override

	ColinearityKinematicConstraintT ()

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

int	numberOfEquations () const override

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

positionDerivativeType const &	positionDerivative () const

valueType const &	value () const

	~MultiTrackKinematicConstraintT () override

Public Member Functions inherited from MultiTrackKinematicConstraintBaseT
virtual	~MultiTrackKinematicConstraintBaseT ()

Private Member Functions
void	fillParametersDerivative () const override

void	fillPositionDerivative () const override

void	fillValue () const override

Private Attributes
double	a_1

double	a_2

AlgebraicVector7	p1

AlgebraicVector7	p2

GlobalPoint	point

Additional Inherited Members
Protected Member Functions inherited from MultiTrackKinematicConstraintT< 2, int(Dim)>
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

template<enum colinearityKinematic::ConstraintDim Dim>
class ColinearityKinematicConstraintT< Dim >

Definition at line 23 of file ColinearityKinematicConstraintT.h.

Member Typedef Documentation

◆ super

template<enum colinearityKinematic::ConstraintDim Dim>

typedef MultiTrackKinematicConstraintT<2, int(Dim)> ColinearityKinematicConstraintT< Dim >::super

Definition at line 34 of file ColinearityKinematicConstraintT.h.

Constructor & Destructor Documentation

◆ ColinearityKinematicConstraintT()

template<enum colinearityKinematic::ConstraintDim Dim>

ColinearityKinematicConstraintT< Dim >::ColinearityKinematicConstraintT ( )

inline

Definition at line 36 of file ColinearityKinematicConstraintT.h.

36 {}

Member Function Documentation

◆ clone()

template<enum colinearityKinematic::ConstraintDim Dim>

ColinearityKinematicConstraintT<Dim>* ColinearityKinematicConstraintT< Dim >::clone ( void ) const

inlineoverridevirtual

Implements MultiTrackKinematicConstraintBaseT.

Definition at line 59 of file ColinearityKinematicConstraintT.h.

                                                                {
     return new ColinearityKinematicConstraintT<Dim>(*this);
   }

◆ fillParametersDerivative()

template<enum colinearityKinematic::ConstraintDim Dim>

void ColinearityKinematicConstraintT< Dim >::fillParametersDerivative ( ) const

overrideprivatevirtual

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

Implements MultiTrackKinematicConstraintT< 2, int(Dim)>.

Definition at line 115 of file ColinearityKinematicConstraintT.h.

                                                                           {
   typename super::parametersDerivativeType& jac_d = super::jac_d();
  
   double p1vx = p1(3) - a_1 * (point.y() - p1(1));
   double p1vy = p1(4) + a_1 * (point.x() - p1(0));
   double k1 = 1.0 / (p1vx * p1vx + p1vy * p1vy);
  
   double p2vx = p2(3) - a_2 * (point.y() - p2(1));
   double p2vy = p2(4) + a_2 * (point.x() - p2(0));
   double k2 = 1.0 / (p2vx * p2vx + p2vy * p2vy);
  
   // H_phi:
  
   //x1 and x2 derivatives: 1st and 8th elements
   jac_d(0, 0) = -k1 * a_1 * p1vx;
   jac_d(0, 7) = k2 * a_2 * p2vx;
  
   //y1 and y2 derivatives: 2nd and 9th elements:
   jac_d(0, 1) = -k1 * a_1 * p1vy;
   jac_d(0, 8) = k2 * a_2 * p2vy;
  
   //z1 and z2 components: 3d and 10th elmnets stay 0:
   jac_d(0, 2) = 0.;
   jac_d(0, 9) = 0.;
  
   //px1 and px2 components: 4th and 11th elements:
   jac_d(0, 3) = -k1 * p1vy;
   jac_d(0, 10) = k2 * p2vy;
  
   //py1 and py2 components: 5th and 12 elements:
   jac_d(0, 4) = k1 * p1vx;
   jac_d(0, 11) = -k2 * p2vx;
  
   //pz1 and pz2 components: 6th and 13 elements:
   jac_d(0, 5) = 0.;
   jac_d(0, 12) = 0.;
   //mass components: 7th and 14th elements:
   jac_d(0, 6) = 0.;
   jac_d(0, 13) = 0.;
  
   if (Dim == colinearityKinematic::PhiTheta) {
     double pt1 = sqrt(p1(3) * p1(3) + p1(4) * p1(4));
     double pTot1 = p1(3) * p1(3) + p1(4) * p1(4) + p1(5) * p1(5);
     double pt2 = sqrt(p2(3) * p2(3) + p2(4) * p2(4));
     double pTot2 = p2(3) * p2(3) + p2(4) * p2(4) + p2(5) * p2(5);
  
     // H_theta:
     //x1 and x2 derivatives: 1st and 8th elements
     jac_d(1, 0) = 0.;
     jac_d(1, 7) = 0.;
  
     //y1 and y2 derivatives: 2nd and 9th elements:
     jac_d(1, 1) = 0.;
     jac_d(1, 8) = 0.;
  
     //z1 and z2 components: 3d and 10th elmnets stay 0:
     jac_d(1, 2) = 0.;
     jac_d(1, 9) = 0.;
  
     jac_d(1, 3) = p1(3) * (p1(5) / (pTot1 * pt1));
     jac_d(1, 10) = p2(3) * (-p2(5) / (pTot2 * pt2));
  
     //py1 and py2 components: 5th and 12 elements:
     jac_d(1, 4) = p1(4) * (p1(5) / (pTot1 * pt1));
     jac_d(1, 11) = p2(4) * (-p2(5) / (pTot2 * pt2));
  
     //pz1 and pz2 components: 6th and 13 elements:
     jac_d(1, 5) = -pt1 / pTot1;
     jac_d(1, 12) = pt2 / pTot2;
  
     //mass components: 7th and 14th elements:
     jac_d(1, 6) = 0.;
     jac_d(1, 13) = 0.;
   }
 }

References p1, p2, colinearityKinematic::PhiTheta, point, HLT_FULL_cff::pt1, HLT_FULL_cff::pt2, and mathSSE::sqrt().

◆ fillPositionDerivative()

template<enum colinearityKinematic::ConstraintDim Dim>

void ColinearityKinematicConstraintT< Dim >::fillPositionDerivative ( ) const

overrideprivatevirtual

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

Implements MultiTrackKinematicConstraintT< 2, int(Dim)>.

Definition at line 192 of file ColinearityKinematicConstraintT.h.

                                                                         {
   typename super::positionDerivativeType& jac_e = super::jac_e();
  
   double p1vx = p1(3) - a_1 * (point.y() - p1(1));
   double p1vy = p1(4) + a_1 * (point.x() - p1(0));
   double k1 = 1.0 / (p1vx * p1vx + p1vy * p1vy);
  
   double p2vx = p2(3) - a_2 * (point.y() - p2(1));
   double p2vy = p2(4) + a_2 * (point.x() - p2(0));
   double k2 = 1.0 / (p2vx * p2vx + p2vy * p2vy);
  
   // H_phi:
  
   // xv component
   jac_e(0, 0) = k1 * a_1 * p1vx - k2 * a_2 * p2vx;
  
   //yv component
   jac_e(0, 1) = k1 * a_1 * p1vy - k2 * a_2 * p2vy;
  
   //zv component
   jac_e(0, 2) = 0.;
  
   // H_theta: no correlation with vertex position
   if (Dim == colinearityKinematic::PhiTheta) {
     jac_e(1, 0) = 0.;
     jac_e(1, 1) = 0.;
     jac_e(1, 2) = 0.;
   }
 }

References p1, p2, colinearityKinematic::PhiTheta, and point.

◆ fillValue()

template<enum colinearityKinematic::ConstraintDim Dim>

void ColinearityKinematicConstraintT< Dim >::fillValue ( ) const

overrideprivatevirtual

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

Implements MultiTrackKinematicConstraintT< 2, int(Dim)>.

Definition at line 87 of file ColinearityKinematicConstraintT.h.

                                                            {
   typename super::valueType& vl = super::vl();
  
   double p1vx = p1(3) - a_1 * (point.y() - p1(1));
   double p1vy = p1(4) + a_1 * (point.x() - p1(0));
  
   double p2vx = p2(3) - a_2 * (point.y() - p2(1));
   double p2vy = p2(4) + a_2 * (point.x() - p2(0));
  
   // H_phi:
   vl(0) = atan2(p1vy, p1vx) - atan2(p2vy, p2vx);
   if (vl(0) > M_PI)
     vl(0) -= 2.0 * M_PI;
   if (vl(0) <= -M_PI)
     vl(0) += 2.0 * M_PI;
   // H_theta:
   if (Dim == colinearityKinematic::PhiTheta) {
     double pt1 = sqrt(p1(3) * p1(3) + p1(4) * p1(4));
     double pt2 = sqrt(p2(3) * p2(3) + p2(4) * p2(4));
     vl(1) = atan2(pt1, p1(5)) - atan2(pt2, p2(5));
     if (vl(1) > M_PI)
       vl(1) -= 2.0 * M_PI;
     if (vl(1) <= -M_PI)
       vl(1) += 2.0 * M_PI;
   }
 }

References M_PI, p1, p2, colinearityKinematic::PhiTheta, point, HLT_FULL_cff::pt1, HLT_FULL_cff::pt2, and mathSSE::sqrt().

◆ init()

template<enum colinearityKinematic::ConstraintDim Dim>

void ColinearityKinematicConstraintT< Dim >::init	(	const std::vector< KinematicState > &	states,
		const GlobalPoint &	ipoint,
		const GlobalVector &	fieldValue
	)

inlineoverridevirtual

Implements MultiTrackKinematicConstraintBaseT.

Definition at line 39 of file ColinearityKinematicConstraintT.h.

                                                      {
     if (states.size() != 2)
       throw VertexException("ColinearityKinematicConstraint::<2 states passed");
  
     point = ipoint;
  
     a_1 = -states[0].particleCharge() * fieldValue.z();
     a_2 = -states[1].particleCharge() * fieldValue.z();
  
     p1 = states[0].kinematicParameters().vector();
     p2 = states[1].kinematicParameters().vector();
   }