---

JacobianLocalToCurvilinear Class Reference

Class which calculates the Jacobian matrix of the transformation from the local to the curvilinear frame. More...

#include <TrackingTools/AnalyticalJacobians/interface/JacobianLocalToCurvilinear.h>

List of all members.

Public Member Functions
const AlgebraicMatrix55 &	jacobian () const
	Access to Jacobian.
const AlgebraicMatrix	jacobian_old () const
	JacobianLocalToCurvilinear (const Surface &surface, const LocalTrajectoryParameters &localParameters, const MagneticField &magField)
	Constructor from local trajectory parameters and surface defining the local frame.
Private Attributes
AlgebraicMatrix55	theJacobian

---

Detailed Description

Class which calculates the Jacobian matrix of the transformation from the local to the curvilinear frame.

The Jacobian is calculated during construction and thereafter cached, enabling reuse of the same Jacobian without calculating it again.

Definition at line 16 of file JacobianLocalToCurvilinear.h.

---

Constructor & Destructor Documentation

JacobianLocalToCurvilinear::JacobianLocalToCurvilinear	(	const Surface &	surface,
		const LocalTrajectoryParameters &	localParameters,
		const MagneticField &	magField
	)

Constructor from local trajectory parameters and surface defining the local frame.

NB!! No default constructor exists!

Definition at line 7 of file JacobianLocalToCurvilinear.cc.

References Vector3DBase< T, FrameTag >::dot(), e, h, MagneticField::inTesla(), PV3DBase< T, PVType, FrameType >::mag(), LocalTrajectoryParameters::momentum(), p, PV3DBase< T, PVType, FrameType >::perp(), LocalTrajectoryParameters::position(), LocalTrajectoryParameters::signedInverseMomentum(), theJacobian, Surface::toGlobal(), Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                          : theJacobian() {
  
    // Origin: TRSDSC
  GlobalPoint  x = surface.toGlobal(localParameters.position());
  LocalVector tnl = localParameters.momentum().unit();
  GlobalVector dj = surface.toGlobal(LocalVector(1., 0., 0.));
  GlobalVector dk = surface.toGlobal(LocalVector(0., 1., 0.));
  GlobalVector tn = surface.toGlobal(tnl);

  GlobalVector p = surface.toGlobal(localParameters.momentum());
  GlobalVector pt(p.x(), p.y(), 0.);
  pt = pt.unit();
  //  GlobalVector di = tsos.surface().toGlobal(LocalVector(0., 0., 1.));

  // rotate coordinates because of wrong coordinate system in orca
  LocalVector tvw(tnl.z(), tnl.x(), tnl.y());
  double cosl = tn.perp(); if (cosl < 1.e-30) cosl = 1.e-30;
  double cosl1 = 1./cosl;
  GlobalVector un(-tn.y()*cosl1, tn.x()*cosl1, 0.);
  GlobalVector vn(-tn.z()*un.y(), tn.z()*un.x(), cosl);
  double uj = un.dot(dj);
  double uk = un.dot(dk);
  double vj = vn.dot(dj);
  double vk = vn.dot(dk);
  theJacobian(0,0) = 1.;
  theJacobian(1,1) = tvw.x()*vj;
  theJacobian(1,2) = tvw.x()*vk;
  theJacobian(2,1) = tvw.x()*uj*cosl1;
  theJacobian(2,2) = tvw.x()*uk*cosl1;
  theJacobian(3,3) = uj;
  theJacobian(3,4) = uk;
  theJacobian(4,3) = vj;
  theJacobian(4,4) = vk;
  // GlobalVector h = MagneticField::inInverseGeV(x);
  GlobalVector h  = magField.inTesla(x) * 2.99792458e-3;
  double q = -h.mag() * localParameters.signedInverseMomentum();
  double sinz =-un.dot(h.unit());
  double cosz = vn.dot(h.unit());
  theJacobian(1,3) = -q*tvw.y()*sinz;
  theJacobian(1,4) = -q*tvw.z()*sinz;
  theJacobian(2,3) = -q*tvw.y()*cosz*cosl1;
  theJacobian(2,4) = -q*tvw.z()*cosz*cosl1;
  // end of TRSDSC

  //dbg::dbg_trace(1,"Loc2Cu", localParameters.vector(),x,dj,dk,theJacobian);
}

---

Member Function Documentation

const AlgebraicMatrix55 & JacobianLocalToCurvilinear::jacobian

(

)

const

Access to Jacobian.

Definition at line 58 of file JacobianLocalToCurvilinear.cc.

References theJacobian.

Referenced by BasicSingleTrajectoryState::checkCurvilinError().

                                                                   {
  return theJacobian;
}

const AlgebraicMatrix JacobianLocalToCurvilinear::jacobian_old

(

)

const

Definition at line 55 of file JacobianLocalToCurvilinear.cc.

References asHepMatrix(), and theJacobian.

                                                                    {
  return asHepMatrix(theJacobian);
}

---

Member Data Documentation

AlgebraicMatrix55 JacobianLocalToCurvilinear::theJacobian [private]

Definition at line 37 of file JacobianLocalToCurvilinear.h.

Referenced by jacobian(), jacobian_old(), and JacobianLocalToCurvilinear().

---

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

• TrackingTools/AnalyticalJacobians/interface/JacobianLocalToCurvilinear.h
• TrackingTools/AnalyticalJacobians/src/JacobianLocalToCurvilinear.cc

---