|
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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
|
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingTree.h"
#include <algorithm>
#include <sstream>
#include <iostream>
using namespace std;
using namespace sipixelobjects;
typedef std::unordered_map<int, SiPixelFedCablingTree::PixelFEDCabling>::const_iterator IMAP;
std::vector<sipixelobjects::CablingPathToDetUnit> SiPixelFedCablingTree::pathToDetUnit(uint32_t rawDetId) const {
std::vector<sipixelobjects::CablingPathToDetUnit> result;
for (auto im = theFedCablings.begin(); im != theFedCablings.end(); ++im) {
const PixelFEDCabling& aFed = im->second;
for (unsigned int idxLink = 1; idxLink <= aFed.numberOfLinks(); idxLink++) {
const PixelFEDLink* link = aFed.link(idxLink);
if (!link)
continue;
unsigned int numberOfRocs = link->numberOfROCs();
for (unsigned int idxRoc = 1; idxRoc <= numberOfRocs; idxRoc++) {
const PixelROC* roc = link->roc(idxRoc);
if (rawDetId == roc->rawId()) {
CablingPathToDetUnit path = {aFed.id(), link->id(), roc->idInLink()};
result.push_back(path);
}
}
}
}
return result;
}
bool SiPixelFedCablingTree::pathToDetUnitHasDetUnit(uint32_t rawDetId, unsigned int fedId) const {
for (auto im = theFedCablings.begin(); im != theFedCablings.end(); ++im) {
const PixelFEDCabling& aFed = im->second;
if (aFed.id() == fedId) {
for (unsigned int idxLink = 1; idxLink <= aFed.numberOfLinks(); idxLink++) {
const PixelFEDLink* link = aFed.link(idxLink);
if (!link)
continue;
unsigned int numberOfRocs = link->numberOfROCs();
for (unsigned int idxRoc = 1; idxRoc <= numberOfRocs; idxRoc++) {
const PixelROC* roc = link->roc(idxRoc);
if (rawDetId == roc->rawId()) {
return true;
}
}
}
}
}
return false;
}
std::unordered_map<uint32_t, unsigned int> SiPixelFedCablingTree::det2fedMap() const {
std::unordered_map<uint32_t, unsigned int> result;
for (auto im = theFedCablings.begin(); im != theFedCablings.end(); ++im) {
auto const& aFed = im->second;
for (unsigned int idxLink = 1; idxLink <= aFed.numberOfLinks(); idxLink++) {
auto link = aFed.link(idxLink);
if (!link)
continue;
unsigned int numberOfRocs = link->numberOfROCs();
for (unsigned int idxRoc = 1; idxRoc <= numberOfRocs; idxRoc++) {
auto roc = link->roc(idxRoc);
result[roc->rawId()] = aFed.id(); // we know that a det is in one fed only...
}
}
}
return result;
}
std::map<uint32_t, std::vector<sipixelobjects::CablingPathToDetUnit> > SiPixelFedCablingTree::det2PathMap() const {
std::map<uint32_t, std::vector<sipixelobjects::CablingPathToDetUnit> > result;
for (auto im = theFedCablings.begin(); im != theFedCablings.end(); ++im) {
auto const& aFed = im->second;
for (unsigned int idxLink = 1; idxLink <= aFed.numberOfLinks(); idxLink++) {
auto link = aFed.link(idxLink);
if (!link)
continue;
unsigned int numberOfRocs = link->numberOfROCs();
for (unsigned int idxRoc = 1; idxRoc <= numberOfRocs; idxRoc++) {
auto roc = link->roc(idxRoc);
CablingPathToDetUnit path = {aFed.id(), link->id(), roc->idInLink()};
result[roc->rawId()].push_back(path);
}
}
}
return result;
}
void SiPixelFedCablingTree::addFed(const PixelFEDCabling& f) {
int id = f.id();
theFedCablings[id] = f;
}
const PixelFEDCabling* SiPixelFedCablingTree::fed(unsigned int id) const {
auto it = theFedCablings.find(id);
return (it == theFedCablings.end()) ? nullptr : &(*it).second;
}
string SiPixelFedCablingTree::print(int depth) const {
ostringstream out;
if (depth-- >= 0) {
out << theVersion << endl;
for (IMAP it = theFedCablings.begin(); it != theFedCablings.end(); it++) {
out << (*it).second.print(depth);
}
}
out << endl;
return out.str();
}
std::vector<const PixelFEDCabling*> SiPixelFedCablingTree::fedList() const {
std::vector<const PixelFEDCabling*> result;
for (IMAP im = theFedCablings.begin(); im != theFedCablings.end(); im++) {
result.push_back(&(im->second));
}
std::sort(result.begin(), result.end(), [](const PixelFEDCabling* a, const PixelFEDCabling* b) {
return a->id() < b->id();
});
return result;
}
void SiPixelFedCablingTree::addItem(unsigned int fedId, unsigned int linkId, const PixelROC& roc) {
PixelFEDCabling& cabling = theFedCablings[fedId];
if (cabling.id() != fedId)
cabling = PixelFEDCabling(fedId);
cabling.addItem(linkId, roc);
}
const sipixelobjects::PixelROC* SiPixelFedCablingTree::findItem(const CablingPathToDetUnit& path) const {
const PixelROC* roc = nullptr;
const PixelFEDCabling* aFed = fed(path.fed);
if (aFed) {
const PixelFEDLink* aLink = aFed->link(path.link);
if (aLink)
roc = aLink->roc(path.roc);
}
return roc;
}
const sipixelobjects::PixelROC* SiPixelFedCablingTree::findItemInFed(const CablingPathToDetUnit& path,
const PixelFEDCabling* aFed) const {
const PixelROC* roc = nullptr;
const PixelFEDLink* aLink = aFed->link(path.link);
if (aLink)
roc = aLink->roc(path.roc);
return roc;
}
int SiPixelFedCablingTree::checkNumbering() const {
int status = 0;
for (auto im = theFedCablings.begin(); im != theFedCablings.end(); ++im) {
if (im->first != static_cast<int>(im->second.id())) {
status = 1;
std::cout << "PROBLEM WITH FED ID!!" << im->first << " vs: " << im->second.id() << std::endl;
}
im->second.checkLinkNumbering();
}
return status;
}
|