KarimakiAlignmentDerivatives Class Reference

#include <KarimakiAlignmentDerivatives.h>

Public Member Functions
AlgebraicMatrix	operator() (const TrajectoryStateOnSurface &tsos) const
	Returns 6x2 jacobian matrix. More...

Detailed Description

Calculates derivatives à la Karimaki (cf. CR-2003/022)

Date:

2007/03/12 21:28:48

Revision:

1.4

(last update by

Author:

cklae

)

Definition at line 16 of file KarimakiAlignmentDerivatives.h.

Member Function Documentation

AlgebraicMatrix KarimakiAlignmentDerivatives::operator()

(

const TrajectoryStateOnSurface &

tsos

)

const

Returns 6x2 jacobian matrix.

Definition at line 13 of file KarimakiAlignmentDerivatives.cc.

References TrajectoryStateOnSurface::localParameters(), and LocalTrajectoryParameters::mixedFormatVector().

{ 

  // Get track parameters on surface
  AlgebraicVector5 alivec = tsos.localParameters().mixedFormatVector();

  // [0] q/p  : charged: charge (+ or - one) divided by magnitude of momentum
  //            neutral : inverse magnitude of momentum
  // [1] dxdz : direction tangent in local xz-plane
  // [2] dydz : direction tangent in local yz-plane
  // [3] x    : local x-coordinate
  // [4] y    : local y-coordinate

  double tanpsi   = alivec[1];
  double tantheta = alivec[2];
  double ux       = alivec[3];
  double vx       = alivec[4];

  AlgebraicMatrix aliderivs(6,2);

  aliderivs[0][0]= -1.0;
  aliderivs[0][1]=  0.0;
  aliderivs[1][0]=  0.0;
  aliderivs[1][1]= -1.0;
  aliderivs[2][0]=  tanpsi;
  aliderivs[2][1]=  tantheta;
  aliderivs[3][0]=  vx*tanpsi;
  aliderivs[3][1]=  vx*tantheta;
  aliderivs[4][0]= -ux*tanpsi;   // New beta sign convention
  aliderivs[4][1]= -ux*tantheta; // New beta sign convention
  aliderivs[5][0]=  vx;
  aliderivs[5][1]= -ux;
   
  return(aliderivs);

}