|
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
|
#ifndef AlignTransformError_H
#define AlignTransformError_H
#include "CondFormats/Serialization/interface/Serializable.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Vector/RotationInterfaces.h"
#include "CondFormats/Alignment/interface/Definitions.h"
/// Class holding error due to an Alignment transformation
/// It contains the raw detector id and the symmetrical error matrix.
/// It is optimized for storage (error matrix is stored as C-array)
class AlignTransformError {
public:
typedef CLHEP::HepSymMatrix SymMatrix;
AlignTransformError() {}
AlignTransformError(const SymMatrix& symMatrix, align::ID irawId) : m_rawId(irawId) {
for (unsigned int i = 0; i < m_nPars; ++i)
for (unsigned int j = 0; j <= i; ++j)
m_Parameters[i * (i + 1) / 2 + j] = symMatrix[i][j];
}
SymMatrix matrix() const {
SymMatrix result(m_nPars);
for (unsigned int i = 0; i < m_nPars; ++i)
for (unsigned int j = 0; j <= i; ++j)
result[i][j] = m_Parameters[i * (i + 1) / 2 + j];
return result;
}
align::ID rawId() const { return m_rawId; }
private:
static const unsigned int m_nPars = 3;
double m_Parameters[m_nPars * (m_nPars + 1) / 2];
align::ID m_rawId;
COND_SERIALIZABLE;
};
#endif //AlignTransformError_H
|