|
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
|
#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerRingBuilder.h"
#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
#include "Geometry/TrackerNumberingBuilder/plugins/ExtractStringFromDDD.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "Geometry/TrackerNumberingBuilder/plugins/CmsDetConstruction.h"
#include "Geometry/TrackerNumberingBuilder/interface/trackerStablePhiSort.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <vector>
#include <bitset>
template <class FilteredView>
void CmsTrackerRingBuilder<FilteredView>::buildComponent(FilteredView& fv, GeometricDet* g, const std::string& s) {
moduleName_ = ExtractStringFromDDD<FilteredView>::getString("TkDDDStructure", &fv);
CmsDetConstruction<FilteredView> theCmsDetConstruction;
theCmsDetConstruction.buildComponent(fv, g, s);
}
template <class FilteredView>
void CmsTrackerRingBuilder<FilteredView>::sortNS(FilteredView& fv, GeometricDet* det) {
GeometricDet::ConstGeometricDetContainer& comp = det->components();
GeometricDet::GeometricDetContainer compfw;
GeometricDet::GeometricDetContainer compbw;
switch (comp.front()->type()) {
case GeometricDet::mergedDet:
trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhiGluedModule);
break;
case GeometricDet::DetUnit:
trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhi);
break;
default:
edm::LogError("CmsTrackerRingBuilder")
<< "ERROR - wrong SubDet to sort..... " << det->components().front()->type();
}
static std::string const TECGluedDet("TECGluedDet");
static std::string const TECDet("TECDet");
// TEC
// Module Number: 3 bits [1,...,5 at most]
if (moduleName_ == TECGluedDet || moduleName_ == TECDet) {
// TEC-
if (det->translation().z() < 0 && moduleName_ == TECDet) {
trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhiMirror);
}
if (det->translation().z() < 0 && moduleName_ == TECGluedDet) {
trackerStablePhiSort(comp.begin(), comp.end(), CmsTrackerLevelBuilderHelper::getPhiGluedModuleMirror);
}
for (uint32_t i = 0; i < comp.size(); i++)
det->component(i)->setGeographicalID(i + 1);
} else {
// TID
// Ring Side: 2 bits [back:1 front:2]
// Module Number: 5 bits [1,...,20 at most]
//
for (uint32_t i = 0; i < comp.size(); i++) {
if (std::abs(comp[i]->translation().z()) < std::abs(det->translation().z())) {
compfw.emplace_back(det->component(i));
} else {
compbw.emplace_back(det->component(i));
}
}
for (uint32_t i = 0; i < compbw.size(); i++) {
uint32_t temp = i + 1;
temp |= (1 << 5);
compbw[i]->setGeographicalID(temp);
}
for (uint32_t i = 0; i < compfw.size(); i++) {
uint32_t temp = i + 1;
temp |= (2 << 5);
compfw[i]->setGeographicalID(temp);
}
det->clearComponents();
det->addComponents(compfw);
det->addComponents(compbw);
}
}
template class CmsTrackerRingBuilder<DDFilteredView>;
template class CmsTrackerRingBuilder<cms::DDFilteredView>;
|