|
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
|
/// \file
///
/// $Date: 2007/12/06 01:55:18 $
/// $Revision: 1.2 $
///
/// $Author: ratnik $
/// \author Frederic Ronga - CERN-PH-CMG
#include <string>
#include <iostream>
#include <sstream>
// Framework
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
// Alignment
#include "Alignment/TrackerAlignment/interface/TrackerScenarioBuilder.h"
#include "Alignment/TrackerAlignment/interface/AlignableTracker.h"
//__________________________________________________________________________________________________
TrackerScenarioBuilder::TrackerScenarioBuilder(AlignableTracker* alignable)
: MisalignmentScenarioBuilder(alignable->objectIdProvider().geometry()), theAlignableTracker(alignable) {
if (!theAlignableTracker) {
throw cms::Exception("TypeMismatch") << "Pointer to AlignableTracker is empty.\n";
}
// Fill what is needed for possiblyPartOf(..):
theSubdets.push_back(stripOffModule(align::TPBModule)); // Take care, order matters: 1st pixel, 2nd strip.
theSubdets.push_back(stripOffModule(align::TPEModule));
theFirstStripIndex = theSubdets.size();
theSubdets.push_back(stripOffModule(align::TIBModule));
theSubdets.push_back(stripOffModule(align::TIDModule));
theSubdets.push_back(stripOffModule(align::TOBModule));
theSubdets.push_back(stripOffModule(align::TECModule));
}
//__________________________________________________________________________________________________
void TrackerScenarioBuilder::applyScenario(const edm::ParameterSet& scenario) {
// Apply the scenario to all main components of tracker.
theModifierCounter = 0;
// Seed is set at top-level, and is mandatory
if (this->hasParameter_("seed", scenario))
theModifier.setSeed(static_cast<long>(scenario.getParameter<int>("seed")));
else
throw cms::Exception("BadConfig") << "No generator seed defined!";
// misalignment applied recursively ('subStructures("Tracker")' contains only tracker itself)
this->decodeMovements_(scenario, theAlignableTracker->subStructures("Tracker"));
edm::LogInfo("TrackerScenarioBuilder") << "Applied modifications to " << theModifierCounter << " alignables";
}
//__________________________________________________________________________________________________
bool TrackerScenarioBuilder::isTopLevel_(const std::string& parameterSetName) const {
// Get root name (strip last character [s])
std::string root = this->rootName_(parameterSetName);
if (root == "Tracker")
return true;
return false;
}
//__________________________________________________________________________________________________
bool TrackerScenarioBuilder::possiblyPartOf(const std::string& subStruct, const std::string& largeStr) const {
// string::find(s) != nPos => 's' is contained in string!
const std::string::size_type nPos = std::string::npos;
// First check whether anything from pixel in strip.
if (largeStr.find("Strip") != nPos) {
if (subStruct.find("Pixel") != nPos)
return false;
for (unsigned int iPix = 0; iPix < theFirstStripIndex; ++iPix) {
if (subStruct.find(theSubdets[iPix]) != nPos)
return false;
}
}
// Now check whether anything from strip in pixel.
if (largeStr.find("Pixel") != nPos) {
if (subStruct.find("Strip") != nPos)
return false;
for (unsigned int iStrip = theFirstStripIndex; iStrip < theSubdets.size(); ++iStrip) {
if (subStruct.find(theSubdets[iStrip]) != nPos)
return false;
}
}
// Finally check for any different detector parts, e.g. TIDEndcap/TIBString gives false.
for (unsigned int iSub = 0; iSub < theSubdets.size(); ++iSub) {
for (unsigned int iLarge = 0; iLarge < theSubdets.size(); ++iLarge) {
if (iLarge == iSub)
continue;
if (largeStr.find(theSubdets[iLarge]) != nPos && subStruct.find(theSubdets[iSub]) != nPos) {
return false;
}
}
}
// It seems like a possible combination:
return true;
}
//__________________________________________________________________________________________________
std::string TrackerScenarioBuilder::stripOffModule(const align::StructureType& type) const {
const std::string module{"Module"};
std::string name{theAlignableTracker->objectIdProvider().typeToName(type)};
auto start = name.find(module);
if (start == std::string::npos) {
throw cms::Exception("LogicError") << "[TrackerScenarioBuilder] '" << name << "' is not a module type";
}
name.replace(start, module.length(), "");
return name;
}
|