|
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
|
#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerOTLayerBuilder.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 "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerOTRingBuilder.h"
#include "Geometry/TrackerNumberingBuilder/plugins/CmsTrackerLadderBuilder.h"
#include "Geometry/TrackerNumberingBuilder/interface/trackerStablePhiSort.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <vector>
#include <bitset>
template <class FilteredView>
void CmsTrackerOTLayerBuilder<FilteredView>::buildComponent(FilteredView& fv, GeometricDet* g, const std::string& s) {
LogTrace("DetConstruction") << " CmsTrackerOTLayerBuilder::buildComponent ";
CmsTrackerLadderBuilder<FilteredView> theCmsTrackerLadderBuilder;
CmsTrackerOTRingBuilder<FilteredView> theCmsTrackerOTRingBuilder;
GeometricDet* subdet = new GeometricDet(&fv,
CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
ExtractStringFromDDD<FilteredView>::getString(s, &fv)));
switch (CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
ExtractStringFromDDD<FilteredView>::getString(s, &fv))) {
case GeometricDet::ladder:
theCmsTrackerLadderBuilder.build(fv, subdet, s);
break;
case GeometricDet::panel:
theCmsTrackerOTRingBuilder.build(fv, subdet, s);
break;
default:
edm::LogError("CmsTrackerOTLayerBuilder") << " ERROR - I was expecting a ladder or a panel, I got a "
<< ExtractStringFromDDD<FilteredView>::getString(s, &fv);
}
g->addComponent(subdet);
}
template <class FilteredView>
void CmsTrackerOTLayerBuilder<FilteredView>::sortNS(FilteredView& fv, GeometricDet* det) {
GeometricDet::ConstGeometricDetContainer comp = det->components();
//order ladder and rings together
GeometricDet::GeometricDetContainer rods;
GeometricDet::GeometricDetContainer ringsNeg;
GeometricDet::GeometricDetContainer ringsPos;
rods.clear();
ringsNeg.clear();
ringsPos.clear();
for (uint32_t i = 0; i < comp.size(); i++) {
auto component = det->component(i);
if (component->type() == GeometricDet::ladder) {
rods.emplace_back(component);
} else if (component->type() == GeometricDet::panel) {
if (component->translation().z() < 0.) {
ringsNeg.emplace_back(component);
} else if (component->translation().z() > 0.) {
ringsPos.emplace_back(component);
}
} else {
edm::LogError("CmsTrackerOTLayerBuilder")
<< "ERROR - wrong SubDet to sort..... " << det->components().front()->type();
}
}
// negative rings
if (!ringsNeg.empty()) {
std::sort(ringsNeg.begin(), ringsNeg.end(), CmsTrackerLevelBuilderHelper::isLessZ);
uint32_t totalringsNeg = ringsNeg.size();
LogTrace("DetConstruction") << " Neg rings ordered by z: ";
for (uint32_t ring = 0; ring < totalringsNeg; ring++) {
uint32_t temp = ring + 1;
temp |= (1 << 8);
ringsNeg[ring]->setGeographicalID(temp);
LogTrace("BuildingTrackerDetId") << "\t\t\t DetId >> " << temp
<< "(r: " << sqrt(ringsNeg[ring]->translation().Perp2())
<< ", phi: " << ringsNeg[ring]->phi()
<< ", z: " << ringsNeg[ring]->translation().z() << ")";
}
}
// rods
if (!rods.empty()) {
trackerStablePhiSort(rods.begin(), rods.end(), CmsTrackerLevelBuilderHelper::getPhi);
uint32_t totalrods = rods.size();
LogTrace("DetConstruction") << " Rods ordered by phi: ";
for (uint32_t rod = 0; rod < totalrods; rod++) {
uint32_t temp = rod + 1;
temp |= (3 << 8);
rods[rod]->setGeographicalID(DetId(temp));
LogTrace("BuildingTrackerDetId") << "\t\t\t DetId >> " << temp << "(r: " << sqrt(rods[rod]->translation().Perp2())
<< ", phi: " << rods[rod]->phi() << ", z: " << rods[rod]->translation().z()
<< ")";
}
}
// positive rings
if (!ringsPos.empty()) {
std::sort(ringsPos.begin(), ringsPos.end(), CmsTrackerLevelBuilderHelper::isLessZ);
uint32_t totalringsPos = ringsPos.size();
LogTrace("DetConstruction") << " Pos rings ordered by z: ";
for (uint32_t ring = 0; ring < totalringsPos; ring++) {
uint32_t temp = ring + 1;
temp |= (2 << 8);
ringsPos[ring]->setGeographicalID(temp);
LogTrace("BuildingTrackerDetId") << "\t\t\t DetId >> " << temp
<< "(r: " << sqrt(ringsPos[ring]->translation().Perp2())
<< ", phi: " << ringsPos[ring]->phi()
<< ", z: " << ringsPos[ring]->translation().z() << ")";
}
}
det->clearComponents();
det->addComponents(ringsNeg);
det->addComponents(rods);
det->addComponents(ringsPos);
}
template class CmsTrackerOTLayerBuilder<DDFilteredView>;
template class CmsTrackerOTLayerBuilder<cms::DDFilteredView>;
|