|
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
114
115
116
117
118
119
120
121
122
123
124
|
/// $Date: 2010/09/10 11:19:37 $
/// $Revision: 1.3 $
/// (last update by $Author: mussgill $)
#include "Alignment/CommonAlignmentParametrization/interface/AlignmentParametersFactory.h"
#include "Alignment/CommonAlignmentParametrization/interface/BeamSpotAlignmentParameters.h"
#include "Alignment/CommonAlignmentParametrization/interface/BowedSurfaceAlignmentParameters.h"
#include "Alignment/CommonAlignmentParametrization/interface/RigidBodyAlignmentParameters.h"
#include "Alignment/CommonAlignmentParametrization/interface/RigidBodyAlignmentParameters4D.h"
#include "Alignment/CommonAlignmentParametrization/interface/TwoBowedSurfacesAlignmentParameters.h"
//#include "Alignment/SurveyAnalysis/interface/SurveyParameters.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <string>
#include <vector>
namespace AlignmentParametersFactory {
//_______________________________________________________________________________________
ParametersType parametersType(const std::string &typeString) {
if (typeString == "RigidBody")
return kRigidBody;
else if (typeString == "Survey")
return kSurvey; // GF: do not belong here, so remove in the long term...
else if (typeString == "RigidBody4D")
return kRigidBody4D;
else if (typeString == "BeamSpot")
return kBeamSpot;
else if (typeString == "BowedSurface")
return kBowedSurface;
else if (typeString == "TwoBowedSurfaces")
return kTwoBowedSurfaces;
throw cms::Exception("BadConfig") << "AlignmentParametersFactory"
<< " No AlignmentParameters with name '" << typeString << "'.";
return kRigidBody; // to please compiler...
}
//_______________________________________________________________________________________
ParametersType parametersType(int typeInt) {
if (typeInt == kRigidBody)
return kRigidBody;
if (typeInt == kSurvey)
return kSurvey; // GF: do not belong here, so remove in the long term...
if (typeInt == kRigidBody4D)
return kRigidBody4D;
if (typeInt == kBeamSpot)
return kBeamSpot;
if (typeInt == kBowedSurface)
return kBowedSurface;
if (typeInt == kTwoBowedSurfaces)
return kTwoBowedSurfaces;
throw cms::Exception("BadConfig") << "AlignmentParametersFactory"
<< " No AlignmentParameters with number " << typeInt << ".";
return kRigidBody; // to please compiler...
}
//_______________________________________________________________________________________
std::string parametersTypeName(ParametersType parType) {
switch (parType) {
case kRigidBody:
return "RigiBody";
case kSurvey: // GF: do not belong here, so remove in the long term...
return "Survey";
case kRigidBody4D:
return "RigiBody4D";
case kBeamSpot:
return "BeamSpot";
case kBowedSurface:
return "BowedSurface";
case kTwoBowedSurfaces:
return "TwoBowedSurfaces";
}
return "unknown_should_never_reach"; // to please the compiler
}
//_______________________________________________________________________________________
AlignmentParameters *createParameters(Alignable *ali, ParametersType parType, const std::vector<bool> &sel) {
switch (parType) {
case kRigidBody: {
const AlgebraicVector par(RigidBodyAlignmentParameters::N_PARAM, 0);
const AlgebraicSymMatrix cov(RigidBodyAlignmentParameters::N_PARAM, 0);
return new RigidBodyAlignmentParameters(ali, par, cov, sel);
} break;
case kSurvey:
// creates some unwanted dependencies - and does not fit into
// AlignmentParameters anyway!
throw cms::Exception("BadConfig") << "AlignmentParametersFactory cannot create SurveyParameters.";
// edm::LogWarning("Alignment") << "@SUB=createParameters"
// << "Creating SurveyParameters of length 0!";
// return new SurveyParameters(ali, AlgebraicVector(),
// AlgebraicSymMatrix());
break;
case kRigidBody4D: {
const AlgebraicVector par(RigidBodyAlignmentParameters4D::N_PARAM, 0);
const AlgebraicSymMatrix cov(RigidBodyAlignmentParameters4D::N_PARAM, 0);
return new RigidBodyAlignmentParameters4D(ali, par, cov, sel);
} break;
case kBeamSpot: {
const AlgebraicVector par(BeamSpotAlignmentParameters::N_PARAM, 0);
const AlgebraicSymMatrix cov(BeamSpotAlignmentParameters::N_PARAM, 0);
return new BeamSpotAlignmentParameters(ali, par, cov, sel);
} break;
case kBowedSurface: {
const AlgebraicVector par(BowedSurfaceAlignmentParameters::N_PARAM, 0);
const AlgebraicSymMatrix cov(BowedSurfaceAlignmentParameters::N_PARAM, 0);
return new BowedSurfaceAlignmentParameters(ali, par, cov, sel);
} break;
case kTwoBowedSurfaces: {
const AlgebraicVector par(TwoBowedSurfacesAlignmentParameters::N_PARAM, 0);
const AlgebraicSymMatrix cov(TwoBowedSurfacesAlignmentParameters::N_PARAM, 0);
return new TwoBowedSurfacesAlignmentParameters(ali, par, cov, sel);
} break;
}
return nullptr; // unreached (all ParametersType appear in switch), to please
// the compiler
}
} // namespace AlignmentParametersFactory
|