Trk::JacobianCotThetaPtToThetaP Class Reference

#include <JacobianCotThetaPtToThetaP.h>

Inheritance diagram for Trk::JacobianCotThetaPtToThetaP:

Collaboration diagram for Trk::JacobianCotThetaPtToThetaP:

Public Member Functions
	JacobianCotThetaPtToThetaP (const double cotTheta, const double qpT)

Detailed Description

This is the 5x5 jacobian for the transformation of track parameters and errors having a momentum representation of \( (\phi, cot\theta, q/p_{T}) \) to the new standard representation of \( (\phi, \theta, q/p) \), while the first two variables that express the local position representation keep unchanged.

Mathematical motivation:
assuming \( u_{\mu} = ( d_{0}', z_{0}', \phi_{0}', cot\theta, q/p_{T} ) \) and \( v_{\mu} = (d_{0}, z_{0}, \phi_{0}, \theta, q/p)\), the set of \( v_{\mu}\) parameters can be expressed by the set of \( u_{\mu} \) parameters by following functions:

• \( v_{1} = d_{0} = d_{0}' \)
• \( v_{2} = z_{0} = z_{0}' \)
• \( v_{3} = \phi_{0} = \phi_{0}' \)
• \( v_{4} = \theta = atan(\frac{1}{cot\theta}) \)
• \( v_{5} = q/p = \frac{q}{p_{T}} \cdot sin(atan(\frac{1}{cot\theta})) \)
  this results into three non-trivial terms of the corresponding jacobian:
• \( \frac{\partial \theta}{\partial cot\theta} = - \frac{1}{1+cot^{2}\theta} \)
• \( \frac{\partial (q/p)}{\partial cot\theta} = - \frac{cot\theta}{(1+cot^{2}\theta)^{3/2}} \frac{q}{p_{T}}\)
• \( \frac{\partial (q/p)}{\partial (q/p_{T})} = \frac{1}{sqrt{1+cot^2\theta}} \)

Author

Andre.nosp@m.as.S.nosp@m.alzbu.nosp@m.rger.nosp@m.@cern.nosp@m..ch

Definition at line 42 of file JacobianCotThetaPtToThetaP.h.

Constructor & Destructor Documentation

JacobianCotThetaPtToThetaP()

Trk::JacobianCotThetaPtToThetaP::JacobianCotThetaPtToThetaP	(	const double	cotTheta,
		const double	qpT
	)

Definition at line 12 of file JacobianCotThetaPtToThetaP.cxx.

                                                                                                :
  AmgMatrix(5,5)()
{
    
   // initialize to identify matrix  
   this->setIdentity();
        
   double onePlusCotTSq = 1 + cotTheta*cotTheta;
   (*this)(3,3) = - 1./onePlusCotTSq;
   (*this)(4,3) = -cotTheta/(onePlusCotTSq*std::sqrt(onePlusCotTSq)) * qpT;
   (*this)(4,4) = 1./std::sqrt(onePlusCotTSq);
}

---

The documentation for this class was generated from the following files:

• JacobianCotThetaPtToThetaP.h
• JacobianCotThetaPtToThetaP.cxx