|
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
|
#ifndef Alignment_CommonAlignmentAlgorithm_AlignmentCorrelationsStore_h
#define Alignment_CommonAlignmentAlgorithm_AlignmentCorrelationsStore_h
/// This class manages the storage and retrieval of correlations between Alignables
/// for the AlignmentParameterStore. This basic implementation simply stores the
/// off diagonal entries of the "big covariance matrix". ATTENTION: Definition of
/// data structure "Correlations" differs from definition in AlignmentParameterStore,
/// but is used only internally.
#include <map>
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
class Alignable;
class AlignmentCorrelationsStore {
public:
typedef std::map<Alignable*, AlgebraicMatrix> CorrelationsTable;
typedef std::map<Alignable*, CorrelationsTable*> Correlations;
AlignmentCorrelationsStore(void);
virtual ~AlignmentCorrelationsStore(void) {}
/// Write correlations directly to the covariance matrix starting at the
/// given position. Indices are assumed to start from 0.
virtual void correlations(Alignable* ap1, Alignable* ap2, AlgebraicSymMatrix& cov, int row, int col) const;
/// Get correlations directly from the given position of the covariance
/// matrix and store them. Indices are assumed to start from 0.
virtual void setCorrelations(Alignable* ap1, Alignable* ap2, const AlgebraicSymMatrix& cov, int row, int col);
/// Set correlations.
virtual void setCorrelations(Alignable* ap1, Alignable* ap2, AlgebraicMatrix& mat);
/// Check whether correlations are stored for a given pair of alignables.
virtual bool correlationsAvailable(Alignable* ap1, Alignable* ap2) const;
/// Reset correlations.
virtual void resetCorrelations(void);
/// Get number of stored correlations.
virtual unsigned int size(void) const;
private:
void fillCorrelationsTable(Alignable* ap1,
Alignable* ap2,
CorrelationsTable* table,
const AlgebraicSymMatrix& cov,
int row,
int col,
bool transpose);
void fillCovariance(
Alignable* ap1, Alignable* ap2, const AlgebraicMatrix& entry, AlgebraicSymMatrix& cov, int row, int col) const;
void fillCovarianceT(
Alignable* ap1, Alignable* ap2, const AlgebraicMatrix& entry, AlgebraicSymMatrix& cov, int row, int col) const;
void readFromCovariance(
Alignable* ap1, Alignable* ap2, AlgebraicMatrix& entry, const AlgebraicSymMatrix& cov, int row, int col);
void readFromCovarianceT(
Alignable* ap1, Alignable* ap2, AlgebraicMatrix& entry, const AlgebraicSymMatrix& cov, int row, int col);
Correlations theCorrelations;
};
#endif
|