|
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
|
#ifndef MTDTOPOLOGY_H
#define MTDTOPOLOGY_H
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/ForwardDetId/interface/ForwardSubdetector.h"
#include "DataFormats/ForwardDetId/interface/BTLDetId.h"
#include "DataFormats/ForwardDetId/interface/ETLDetId.h"
#include <Geometry/CommonDetUnit/interface/GeomDet.h>
#include <vector>
#include <string>
class MTDTopology {
public:
struct BTLLayout {
// number of logical rods, i.e. rows of sensor modules along eta/z in phi, and of modules per rod
static constexpr uint32_t nBTLphi_ = BTLDetId::HALF_ROD * BTLDetId::kModulesPerTrkV2;
static constexpr uint32_t nBTLeta_ =
2 * BTLDetId::kRUPerTypeV2 * BTLDetId::kCrystalTypes * BTLDetId::kModulesPerRUV2 / BTLDetId::kModulesPerTrkV2;
static constexpr uint32_t nBTLmodules_ = nBTLphi_ * nBTLeta_;
std::array<uint32_t, nBTLmodules_> btlDetId_;
std::array<uint32_t, nBTLmodules_> btlPhi_;
std::array<uint32_t, nBTLmodules_> btlEta_;
};
using BTLValues = BTLLayout;
struct ETLfaceLayout {
uint32_t idDiscSide_; // disc face identifier: 0 disc1 F, 1 disc1 B, 2 disc2 F, 3 disc2 B
uint32_t idDetType1_; // module type id identifier for first row
std::array<std::vector<int>, 2> start_copy_; // start copy per row, first of type idDetType1_
std::array<std::vector<int>, 2> offset_; // offset per row, first of type idDetType1_
};
using ETLValues = std::vector<ETLfaceLayout>;
MTDTopology(const int& topologyMode, const BTLValues& btl, const ETLValues& etl);
int getMTDTopologyMode() const { return mtdTopologyMode_; }
uint32_t btlRods() const { return btlVals_.nBTLphi_; }
uint32_t btlModulesPerRod() const { return btlVals_.nBTLeta_; }
uint32_t btlModules() const { return btlVals_.nBTLmodules_; }
// BTL topology navigation is based on a predefined order of dets in MTDGeometry, mapped onto phi/eta grid
std::pair<uint32_t, uint32_t> btlIndex(const uint32_t detId) const;
uint32_t btlidFromIndex(const uint32_t iphi, const uint32_t ieta) const;
// BTL topology navigation methods, find index of closest module along eta or phi
uint32_t phishiftBTL(const uint32_t detid, const int phiShift) const;
uint32_t etashiftBTL(const uint32_t detid, const int etaShift) const;
// ETL topology navigation is based on a predefined order of dets in sector
static bool orderETLSector(const GeomDet*& gd1, const GeomDet*& gd2);
// navigation methods in ETL topology, provide the index of the det next to DetId for
// horizontal and vertical shifts in both directions, assuming the predefined order in a sector
size_t hshiftETL(const uint32_t detid, const int horizontalShift) const;
size_t vshiftETL(const uint32_t detid, const int verticalShift, size_t& closest) const;
private:
const int mtdTopologyMode_;
const BTLValues btlVals_;
const ETLValues etlVals_;
static constexpr size_t failIndex_ =
std::numeric_limits<unsigned int>::max(); // return out-of-range value for any failure
};
#endif
|