|
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
|
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
CaloTopology::CaloTopology() {}
CaloTopology::~CaloTopology() = default;
int CaloTopology::makeIndex(DetId::Detector det, int subdet) const { return (int(det) << 4) | (subdet & 0xF); }
void CaloTopology::setSubdetTopology(DetId::Detector det,
int subdet,
std::unique_ptr<const CaloSubdetectorTopology> geom) {
int index = makeIndex(det, subdet);
theTopologies_[index] = std::move(geom);
}
const CaloSubdetectorTopology* CaloTopology::getSubdetectorTopology(const DetId& id) const {
auto i = theTopologies_.find(makeIndex(id.det(), id.subdetId()));
return (i == theTopologies_.end()) ? (nullptr) : (i->second.get());
}
const CaloSubdetectorTopology* CaloTopology::getSubdetectorTopology(DetId::Detector det, int subdet) const {
auto i = theTopologies_.find(makeIndex(det, subdet));
return (i == theTopologies_.end()) ? (nullptr) : (i->second.get());
}
static const std::vector<DetId> emptyDetIdVector;
std::vector<DetId> CaloTopology::east(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->east(id));
}
std::vector<DetId> CaloTopology::west(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->west(id));
}
std::vector<DetId> CaloTopology::north(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->north(id));
}
std::vector<DetId> CaloTopology::south(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->south(id));
}
std::vector<DetId> CaloTopology::up(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->up(id));
}
std::vector<DetId> CaloTopology::down(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->down(id));
}
std::vector<DetId> CaloTopology::getNeighbours(const DetId& id, const CaloDirection& dir) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->getNeighbours(id, dir));
}
std::vector<DetId> CaloTopology::getWindow(const DetId& id, const int& northSouthSize, const int& eastWestSize) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->getWindow(id, northSouthSize, eastWestSize));
}
std::vector<DetId> CaloTopology::getAllNeighbours(const DetId& id) const {
const CaloSubdetectorTopology* topology = getSubdetectorTopology(id);
return (topology == nullptr) ? (emptyDetIdVector) : (topology->getAllNeighbours(id));
}
bool CaloTopology::valid(const DetId& id) const {
const CaloSubdetectorTopology* geom = getSubdetectorTopology(id);
return (geom == nullptr) ? (false) : (geom->valid(id));
}
|