|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
/** \file KarimakiAlignmentDerivatives.cc
*
* $Date: 2007/05/02 21:01:53 $
* $Revision: 1.7 $
*/
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "Alignment/CommonAlignmentParametrization/interface/SegmentAlignmentDerivatives4D.h"
AlgebraicMatrix SegmentAlignmentDerivatives4D::operator()(const TrajectoryStateOnSurface &tsos) const {
// 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);
}
|