SegmentAlignmentDerivatives4D Class Reference

#include <SegmentAlignmentDerivatives4D.h>

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

Detailed Description

Calculates derivatives for 4D Segments

Date:

2007/03/12 21:28:48

Revision:

1.4

(last update by

Author:

P. Martinez

)

Definition at line 16 of file SegmentAlignmentDerivatives4D.h.

Member Function Documentation

AlgebraicMatrix SegmentAlignmentDerivatives4D::operator()

(

const TrajectoryStateOnSurface &

tsos

)

const

Returns 6x4 jacobian matrix.

Definition at line 13 of file SegmentAlignmentDerivatives4D.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,4);

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

}