|
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
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
|
#ifndef Alignment_CommonAlignment_TwoBowedSurfacesAlignmentParameters_h
#define Alignment_CommonAlignment_TwoBowedSurfacesAlignmentParameters_h
#include "Alignment/CommonAlignment/interface/AlignmentParameters.h"
#include "CondFormats/Alignment/interface/Definitions.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "Alignment/CommonAlignmentParametrization/interface/BowedSurfaceAlignmentDerivatives.h"
/// \class TwoBowedSurfacesAlignmentParameters
///
/// Concrete class for alignment parameters and associated quantities
/// [derived from AlignmentParameters].
/// The alignable is assumed to have in fact two surfaces devided at
/// a specific local ySplit.
/// The number of parameters N_PARAM is 2x9, i.e. one set of
/// - 3 translations
/// - 3 rotations/slopes
/// - 3 bows/sagittas (u, v and mixed term)
/// for the two surfaces (similar as in BowedSurfaceAlignmentParameters).
/// Storage is done differently:
/// The movements and translations are averaged and applied to the
/// alignable in the 'classical' way.
/// In SurfaceDeformation we store:
/// - the half differences of the movements and rotations, to be taken into
/// account,
/// with positive/negative sign, for corrections for the surface at
/// lower/higher y,
/// - the mean values 'm' and half differences 'hd' of the bows/sagittas of each
/// surface,
/// where these corrections have to be applied as m+hd for the one and m-hd
/// for the other sensor
///
///
/// $Date: 2010/10/26 20:41:08 $
/// $Revision: 1.1 $
/// (last update by $Author: flucke $)
class Alignable;
class AlignableDetOrUnitPtr;
class TrajectoryStateOnSurface;
class TwoBowedSurfacesAlignmentParameters : public AlignmentParameters {
public:
/// Give parameters a name (do not change order, see derivatives(..)!)
typedef BowedSurfaceAlignmentDerivatives BowedDerivs;
enum AlignmentParameterName {
// 1st surface
dx1 = BowedDerivs::dx,
dy1 = BowedDerivs::dy,
dz1 = BowedDerivs::dz,
dslopeX1 = BowedDerivs::dslopeX, // NOTE: slope(u) -> halfWidth*tan(beta),
dslopeY1 = BowedDerivs::dslopeY, // slope(v) -> halfLength*tan(alpha)
drotZ1 = BowedDerivs::drotZ, // rot(w) -> g-scale*gamma
dsagittaX1 = BowedDerivs::dsagittaX,
dsagittaXY1 = BowedDerivs::dsagittaXY,
dsagittaY1 = BowedDerivs::dsagittaY,
// 2nd surface
dx2 = BowedDerivs::dx + BowedDerivs::N_PARAM,
dy2 = BowedDerivs::dy + BowedDerivs::N_PARAM,
dz2 = BowedDerivs::dz + BowedDerivs::N_PARAM,
dslopeX2 = BowedDerivs::dslopeX + BowedDerivs::N_PARAM, // NOTE: slope(u) -> k*tan(beta),
dslopeY2 = BowedDerivs::dslopeY + BowedDerivs::N_PARAM, // slope(v) -> k*tan(alpha)
drotZ2 = BowedDerivs::drotZ + BowedDerivs::N_PARAM, // rot(w) -> m*gamma
dsagittaX2 = BowedDerivs::dsagittaX + BowedDerivs::N_PARAM,
dsagittaXY2 = BowedDerivs::dsagittaXY + BowedDerivs::N_PARAM,
dsagittaY2 = BowedDerivs::dsagittaY + BowedDerivs::N_PARAM,
// number of parameters
N_PARAM = BowedDerivs::N_PARAM + BowedDerivs::N_PARAM
};
/// Constructor with empty parameters/covariance
TwoBowedSurfacesAlignmentParameters(Alignable *alignable);
/// Constructor for full set of parameters
TwoBowedSurfacesAlignmentParameters(Alignable *alignable,
const AlgebraicVector ¶meters,
const AlgebraicSymMatrix &covMatrix);
/// Constructor for selection
TwoBowedSurfacesAlignmentParameters(Alignable *alignable,
const AlgebraicVector ¶meters,
const AlgebraicSymMatrix &covMatrix,
const std::vector<bool> &selection);
/// Destructor
~TwoBowedSurfacesAlignmentParameters() override {}
void apply() override;
int type() const override;
/// Clone all parameters (for update of parameters)
TwoBowedSurfacesAlignmentParameters *clone(const AlgebraicVector ¶meters,
const AlgebraicSymMatrix &covMatrix) const override;
/// Clone selected parameters (for update of parameters)
TwoBowedSurfacesAlignmentParameters *cloneFromSelected(const AlgebraicVector ¶meters,
const AlgebraicSymMatrix &covMatrix) const override;
/// Get all derivatives
AlgebraicMatrix derivatives(const TrajectoryStateOnSurface &tsos, const AlignableDetOrUnitPtr &aliDet) const override;
/// print parameters to screen
virtual void print() const;
double ySplit() const { return ySplit_; }
private:
double ySplitFromAlignable(const Alignable *ali) const;
double ySplit_;
};
#endif
|