|
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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
/// \author : Gero Flucke
/// date : October 2010
#include "Geometry/CommonTopologies/interface/BowedSurfaceDeformation.h"
#include "Geometry/CommonTopologies/interface/SurfaceDeformationFactory.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
// already included via header:
// #include <vector>
//------------------------------------------------------------------------------
BowedSurfaceDeformation::BowedSurfaceDeformation(const std::vector<double> &pars)
: theSagittaX(!pars.empty() ? pars[0] : 0.),
theSagittaY(pars.size() > 2 ? pars[2] : 0.),
theSagittaXY(pars.size() > 1 ? pars[1] : 0.) {
if (pars.size() != minParameterSize()) {
edm::LogError("BadSetup") << "@SUB=BowedSurfaceDeformation"
<< "Input vector of wrong size " << pars.size() << " instead of " << minParameterSize()
<< ", filled up with zeros!";
}
}
//------------------------------------------------------------------------------
BowedSurfaceDeformation *BowedSurfaceDeformation::clone() const {
return new BowedSurfaceDeformation(theSagittaX, theSagittaXY, theSagittaY);
}
//------------------------------------------------------------------------------
int BowedSurfaceDeformation::type() const { return SurfaceDeformationFactory::kBowedSurface; }
//------------------------------------------------------------------------------
SurfaceDeformation::Local2DVector BowedSurfaceDeformation::positionCorrection(const Local2DPoint &localPos,
const LocalTrackAngles &localAngles,
double length,
double width) const {
// different widthes at high/low y could somehow be treated by theRelWidthLowY
// if (widthLowY > 0. && widthHighY != widthLowY) {
// // TEC would always create a warning...
// edm::LogWarning("UnusableData") << "@SUB=BowedSurfaceDeformation::positionCorrection"
// << "Cannot yet deal with different widthes, take "
// << widthHighY << " not " << widthLowY;
// }
// const double width = widthHighY;
double uRel = (width ? 2. * localPos.x() / width : 0.); // relative u (-1 .. +1)
double vRel = (length ? 2. * localPos.y() / length : 0.); // relative v (-1 .. +1)
// 'range check':
const double cutOff = 1.5;
if (uRel < -cutOff) {
uRel = -cutOff;
} else if (uRel > cutOff) {
uRel = cutOff;
}
if (vRel < -cutOff) {
vRel = -cutOff;
} else if (vRel > cutOff) {
vRel = cutOff;
}
// apply coefficients to Legendre polynomials
// to get local height relative to 'average'
const double dw =
(uRel * uRel - 1. / 3.) * theSagittaX + uRel * vRel * theSagittaXY + (vRel * vRel - 1. / 3.) * theSagittaY;
// positive dxdz/dydz and positive dw mean negative shift in x/y:
return Local2DVector(-dw * localAngles);
}
//------------------------------------------------------------------------------
bool BowedSurfaceDeformation::add(const SurfaceDeformation &other) {
if (other.type() == this->type()) {
const std::vector<double> otherParams(other.parameters());
if (otherParams.size() == parameterSize()) { // double check!
theSagittaX += otherParams[0]; // bows can simply be added up
theSagittaXY += otherParams[1];
theSagittaY += otherParams[2];
return true;
}
}
return false;
}
//------------------------------------------------------------------------------
std::vector<double> BowedSurfaceDeformation::parameters() const {
std::vector<double> result(parameterSize());
result[0] = theSagittaX;
result[1] = theSagittaXY;
result[2] = theSagittaY;
return result;
}
|