Trk::JacobianCartesianToPerigee Class Reference

Jacobian from cartesian parameters \( (p_x, p_y, p_z, E, x, y, z) \) to perigee parameters \( (d_0, z_0, \phi_0, \theta, q/p) \). More...

#include <JacobianCartesianToPerigee.h>

Inheritance diagram for Trk::JacobianCartesianToPerigee:

Collaboration diagram for Trk::JacobianCartesianToPerigee:

Public Member Functions
	JacobianCartesianToPerigee (const double px, const double py, const double pz, const double E, const double x, const double y, const double z, const double charge)

Detailed Description

Jacobian from cartesian parameters \( (p_x, p_y, p_z, E, x, y, z) \) to perigee parameters \( (d_0, z_0, \phi_0, \theta, q/p) \).

It needs the particle charge to provide the correct perigee \( q/p \) sign.

Author

Maaik.nosp@m.e.Li.nosp@m.mper@.nosp@m.cern.nosp@m..ch

Definition at line 28 of file JacobianCartesianToPerigee.h.

Constructor & Destructor Documentation

JacobianCartesianToPerigee()

Trk::JacobianCartesianToPerigee::JacobianCartesianToPerigee	(	const double	px,
		const double	py,
		const double	pz,
		const double	E,
		const double	x,
		const double	y,
		const double	z,
		const double	charge
	)

Definition at line 12 of file JacobianCartesianToPerigee.cxx.

                                                                                                                                                                                          :
  AmgMatrix(5,7)()
{
  // initialize to zero  
  this->setZero();
      
  double ptot = std::sqrt(px*px+py*py+pz*pz);
  double pt =   std::sqrt(px*px+py*py);
 
  double tan_theta = pt/pz;
  
  double cos_phi = px/pt;
  double sin_phi = py/pt;
 
  (*this)(0,4) =  -sin_phi ; // dd0/dx
  (*this)(0,5) =  cos_phi ; // dd0/dy 
  (*this)(1,6) =  1.; // dz0/dz
  (*this)(2,0) =  -py/(px*px+py*py) ; // dphi/dpx 
  (*this)(2,1) =  px/(px*px+py*py) ; // dphi/dpy 
  (*this)(3,0) =  (cos_phi/pz)*(1./(1+tan_theta*tan_theta)); // dtheta/dpx
  (*this)(3,1) =  (sin_phi/pz)*(1./(1+tan_theta*tan_theta)); // dtheta/dpy
  (*this)(3,2) =  (-tan_theta/pz)*(1./(1+tan_theta*tan_theta)); // dtheta/dpz
  (*this)(4,3) =  -charge*E/(ptot*ptot*ptot); // dqOverP/dE
 
}

---

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

• JacobianCartesianToPerigee.h
• JacobianCartesianToPerigee.cxx