CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions
KarimakiAlignmentDerivatives Class Reference

#include <KarimakiAlignmentDerivatives.h>

Public Member Functions

AlgebraicMatrix operator() (const TrajectoryStateOnSurface &tsos) const
 

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 of derivatives of residuals in x and y with respect to rigid body aligment parameters:

/ dr_x/du dr_y/du | | dr_x/dv dr_y/dv | | dr_x/dw dr_y/dw | | dr_x/da dr_y/da | | dr_x/db dr_y/db | \ dr_x/dg dr_y/dg /

Definition at line 13 of file KarimakiAlignmentDerivatives.cc.

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

14 {
15 
16  // Get track parameters on surface
18 
19  // [0] q/p : charged: charge (+ or - one) divided by magnitude of momentum
20  // neutral : inverse magnitude of momentum
21  // [1] dxdz : direction tangent in local xz-plane
22  // [2] dydz : direction tangent in local yz-plane
23  // [3] x : local x-coordinate
24  // [4] y : local y-coordinate
25 
26  double tanpsi = alivec[1];
27  double tantheta = alivec[2];
28  double ux = alivec[3];
29  double vx = alivec[4];
30 
31  AlgebraicMatrix aliderivs(6,2);
32 
33  aliderivs[0][0]= -1.0;
34  aliderivs[0][1]= 0.0;
35  aliderivs[1][0]= 0.0;
36  aliderivs[1][1]= -1.0;
37  aliderivs[2][0]= tanpsi;
38  aliderivs[2][1]= tantheta;
39  aliderivs[3][0]= vx*tanpsi;
40  aliderivs[3][1]= vx*tantheta;
41  aliderivs[4][0]= -ux*tanpsi; // New beta sign convention
42  aliderivs[4][1]= -ux*tantheta; // New beta sign convention
43  aliderivs[5][0]= vx;
44  aliderivs[5][1]= -ux;
45 
46  return(aliderivs);
47 
48 }
const LocalTrajectoryParameters & localParameters() const
CLHEP::HepMatrix AlgebraicMatrix
ROOT::Math::SVector< double, 5 > AlgebraicVector5
AlgebraicVector5 mixedFormatVector() const