CMS 3D CMS Logo

List of all members | Public Types | Static Public Member Functions
ErrorFrameTransformer Struct Reference

#include <ErrorFrameTransformer.h>

Public Types

typedef Surface::Scalar Scalar
 

Static Public Member Functions

static GlobalError transform (const LocalError &le, const Surface &surf)
 
static LocalError transform (const GlobalError &ge, const Surface &surf)
 
static LocalErrorExtended transform46 (const GlobalErrorExtended &ge, const AlgebraicVector &positions, const AlgebraicVector &directions)
 

Detailed Description

Definition at line 12 of file ErrorFrameTransformer.h.

Member Typedef Documentation

◆ Scalar

Definition at line 13 of file ErrorFrameTransformer.h.

Member Function Documentation

◆ transform() [1/2]

static GlobalError ErrorFrameTransformer::transform ( const LocalError le,
const Surface surf 
)
inlinestatic

Definition at line 16 of file ErrorFrameTransformer.h.

References alignCSCRings::r, GloballyPositioned< T >::rotation(), LocalError::xx(), LocalError::xy(), and LocalError::yy().

Referenced by TestOutliers::analyze(), GEMEfficiencyAnalyzer::analyze(), GEMEfficiencyAnalyzer::checkBounds(), GsfConstraintAtVertex::constrainAtPoint(), DTSegmentUpdator::fit(), MuonTransientTrackingRecHit::globalDirectionError(), TValidTrackingRecHit::globalPositionError(), TRecHit5DParamConstraint::globalPositionError(), TRecHit2DPosConstraint::globalPositionError(), BaseTrackerRecHit::globalPositionError(), TValidTrackingRecHit::globalState(), BaseTrackerRecHit::globalState(), TkMSParameterizationBuilder::produce(), TrackerGeomDet::setAlignmentPositionError(), DynamicTruncation::testCSCstation(), and DynamicTruncation::testDTstation().

16  {
17  // the GlobalError is a sym matrix, initialisation takes only
18  // 6 T because GlobalError is stored as a lower triangular matrix.
19  Scalar cxx = le.xx();
20  Scalar cxy = le.xy();
21  Scalar cyy = le.yy();
22 
24 
25  return GlobalError(r.xx() * (r.xx() * cxx + r.yx() * cxy) + r.yx() * (r.xx() * cxy + r.yx() * cyy),
26  r.xx() * (r.xy() * cxx + r.yy() * cxy) + r.yx() * (r.xy() * cxy + r.yy() * cyy),
27  r.xy() * (r.xy() * cxx + r.yy() * cxy) + r.yy() * (r.xy() * cxy + r.yy() * cyy),
28  r.xx() * (r.xz() * cxx + r.yz() * cxy) + r.yx() * (r.xz() * cxy + r.yz() * cyy),
29  r.xy() * (r.xz() * cxx + r.yz() * cxy) + r.yy() * (r.xz() * cxy + r.yz() * cyy),
30  r.xz() * (r.xz() * cxx + r.yz() * cxy) + r.yz() * (r.xz() * cxy + r.yz() * cyy));
31  }
GlobalErrorBase< double, ErrorMatrixTag > GlobalError
Definition: GlobalError.h:13
float yy() const
Definition: LocalError.h:24
float xy() const
Definition: LocalError.h:23
const RotationType & rotation() const
float xx() const
Definition: LocalError.h:22

◆ transform() [2/2]

static LocalError ErrorFrameTransformer::transform ( const GlobalError ge,
const Surface surf 
)
inlinestatic

Definition at line 33 of file ErrorFrameTransformer.h.

References GlobalErrorBase< T, ErrorWeightType >::cxx(), GlobalErrorBase< T, ErrorWeightType >::cyx(), GlobalErrorBase< T, ErrorWeightType >::cyy(), GlobalErrorBase< T, ErrorWeightType >::czx(), GlobalErrorBase< T, ErrorWeightType >::czy(), GlobalErrorBase< T, ErrorWeightType >::czz(), alignCSCRings::r, and GloballyPositioned< T >::rotation().

33  {
34  Scalar cxx = ge.cxx();
35  Scalar cyx = ge.cyx();
36  Scalar cyy = ge.cyy();
37  Scalar czx = ge.czx();
38  Scalar czy = ge.czy();
39  Scalar czz = ge.czz();
40 
42 
43  Scalar l11 = r.xx() * (r.xx() * cxx + r.xy() * cyx + r.xz() * czx) +
44  r.xy() * (r.xx() * cyx + r.xy() * cyy + r.xz() * czy) +
45  r.xz() * (r.xx() * czx + r.xy() * czy + r.xz() * czz);
46  Scalar l12 = r.yx() * (r.xx() * cxx + r.xy() * cyx + r.xz() * czx) +
47  r.yy() * (r.xx() * cyx + r.xy() * cyy + r.xz() * czy) +
48  r.yz() * (r.xx() * czx + r.xy() * czy + r.xz() * czz);
49  Scalar l22 = r.yx() * (r.yx() * cxx + r.yy() * cyx + r.yz() * czx) +
50  r.yy() * (r.yx() * cyx + r.yy() * cyy + r.yz() * czy) +
51  r.yz() * (r.yx() * czx + r.yy() * czy + r.yz() * czz);
52 
53  return LocalError(l11, l12, l22);
54  }
const RotationType & rotation() const

◆ transform46()

static LocalErrorExtended ErrorFrameTransformer::transform46 ( const GlobalErrorExtended ge,
const AlgebraicVector positions,
const AlgebraicVector directions 
)
inlinestatic

Definition at line 57 of file ErrorFrameTransformer.h.

References GlobalErrorBaseExtended< T, ErrorWeightType >::matrix(), and MillePedeFileConverter_cfg::out.

Referenced by MuonTransientTrackingRecHit::parametersError().

59  {
60  const AlgebraicSymMatrix66& as(ge.matrix());
61 
62  AlgebraicMatrix46 jacobian46;
63  jacobian46[0][0] = 1.;
64  jacobian46[0][1] = 0.;
65  jacobian46[0][2] = -directions[0];
66  jacobian46[0][3] = -positions[1] * directions[0];
67  jacobian46[0][4] = positions[0] * directions[0];
68  jacobian46[0][5] = -positions[1];
69 
70  jacobian46[1][0] = 0.;
71  jacobian46[1][1] = 1.;
72  jacobian46[1][2] = -directions[1];
73  jacobian46[1][3] = -positions[1] * directions[1];
74  jacobian46[1][4] = positions[0] * directions[1];
75  jacobian46[1][5] = positions[0];
76 
77  jacobian46[2][0] = 0.;
78  jacobian46[2][1] = 0.;
79  jacobian46[2][2] = 0.;
80  jacobian46[2][3] = -directions[1] * directions[0];
81  jacobian46[2][4] = 1. + directions[0] * directions[0];
82  jacobian46[2][5] = -directions[1];
83 
84  jacobian46[3][0] = 0.;
85  jacobian46[3][1] = 0.;
86  jacobian46[3][2] = 0.;
87  jacobian46[3][3] = -1. - directions[1] * directions[1];
88  jacobian46[3][4] = directions[0] * directions[1];
89  jacobian46[3][5] = directions[0];
90 
91  AlgebraicSymMatrix44 out = ROOT::Math::Similarity(jacobian46, as);
92 
93  LocalErrorExtended newError(out);
94 
95  return newError;
96  }
const AlgebraicSymMatrix66 & matrix() const
ROOT::Math::SMatrix< double, 4, 6, ROOT::Math::MatRepStd< double, 4, 6 > > AlgebraicMatrix46
ROOT::Math::SMatrix< double, 4, 4, ROOT::Math::MatRepSym< double, 4 > > AlgebraicSymMatrix44
ROOT::Math::SMatrix< double, 6, 6, ROOT::Math::MatRepSym< double, 6 > > AlgebraicSymMatrix66