|
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
|
#include "EventFilter/L1TRawToDigi/interface/OmtfLinkMappingRpc.h"
#include <fstream>
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CondFormats/RPCObjects/interface/RPCReadOutMapping.h"
#include "CondFormats/RPCObjects/interface/RPCEMap.h"
#include "CondFormats/DataRecord/interface/RPCEMapRcd.h"
#include "CondFormats/RPCObjects/interface/DccSpec.h"
#include "CondFormats/RPCObjects/interface/LinkBoardElectronicIndex.h"
#include "CondFormats/RPCObjects/interface/LinkBoardPackedStrip.h"
#include "CondFormats/RPCObjects/interface/TriggerBoardSpec.h"
#include "CondFormats/RPCObjects/interface/LinkBoardSpec.h"
#include "CondFormats/RPCObjects/interface/RPCAMCLinkMap.h"
namespace omtf {
bool lessLinkBoardElectronicIndex::operator()(const LinkBoardElectronicIndex& o1,
const LinkBoardElectronicIndex& o2) const {
if (o1.dccId < o2.dccId)
return true;
else if (o1.dccId == o2.dccId && o1.dccInputChannelNum < o2.dccInputChannelNum)
return true;
else if (o1.dccId == o2.dccId && o1.dccInputChannelNum == o2.dccInputChannelNum &&
o1.tbLinkInputNum < o2.tbLinkInputNum)
return true;
else
return false;
}
MapEleIndex2LBIndex translateOmtf2Pact(const RpcLinkMap& omtfLink2Ele, const RPCReadOutMapping* pactCabling) {
MapEleIndex2LBIndex omtf2rpc;
std::vector<const DccSpec*> dccs = pactCabling->dccList();
for (auto it1 : dccs) {
const std::vector<TriggerBoardSpec>& rmbs = it1->triggerBoards();
for (auto const& it2 : rmbs) {
const std::vector<LinkConnSpec>& links = it2.linkConns();
for (auto const& it3 : links) {
const std::vector<LinkBoardSpec>& lbs = it3.linkBoards();
for (std::vector<LinkBoardSpec>::const_iterator it4 = lbs.begin(); it4 != lbs.end(); ++it4) {
std::string lbNameCH = it4->linkBoardName();
std::string lbName = lbNameCH.substr(0, lbNameCH.size() - 4);
std::vector<EleIndex> omtfEles = omtfLink2Ele.omtfEleIndex(lbName);
// if (!omtfEles.empty()) std::cout <<" isOK ! " << it4->linkBoardName() <<" has: " << omtfEles.size() << " first: "<< omtfEles[0] << std::endl;
LinkBoardElectronicIndex rpcEle = {
it1->id(), it2.dccInputChannelNum(), it3.triggerBoardInputNumber(), it4->linkBoardNumInLink()};
for (const auto& omtfEle : omtfEles)
omtf2rpc[omtfEle] = rpcEle;
}
}
}
}
LogTrace(" ") << " SIZE OF OMTF to RPC map TRANSLATION is: " << omtf2rpc.size() << std::endl;
return omtf2rpc;
}
MapLBIndex2EleIndex translatePact2Omtf(const RpcLinkMap& omtfLink2Ele, const RPCReadOutMapping* pactCabling) {
MapLBIndex2EleIndex pact2omtfs;
MapEleIndex2LBIndex omtf2rpcs = translateOmtf2Pact(omtfLink2Ele, pactCabling);
for (const auto& omtf2rpc : omtf2rpcs) {
std::pair<EleIndex, EleIndex>& omtfs = pact2omtfs[omtf2rpc.second];
if (omtfs.first.fed() == 0)
omtfs.first = omtf2rpc.first;
else if (omtfs.second.fed() == 0)
omtfs.second = omtf2rpc.first;
else
edm::LogError(" translatePact2Omtf ") << " PROBLEM LinkBoardElectronicIndex already USED!!!! ";
}
return pact2omtfs;
}
void RpcLinkMap::init(const RPCAMCLinkMap& amcMapping) {
const RPCAMCLinkMap::map_type& amcMap = amcMapping.getMap();
for (const auto& item : amcMap) {
unsigned int fedId = item.first.getFED();
unsigned int amcSlot = item.first.getAMCNumber();
unsigned int link = item.first.getAMCInput();
std::string lbName = item.second.getName();
std::string processorNameStr = "OMTF";
;
if (fedId == 1380)
processorNameStr = "OMTFn";
else
processorNameStr = "OMTFp";
processorNameStr += std::to_string(amcSlot / 2 + 1);
std::string processorName = processorNameStr;
std::map<unsigned int, std::string>& li2lb = link2lbName[processorName];
std::map<std::string, unsigned int>& lb2li = lbName2link[processorName];
li2lb[link] = lbName;
lb2li[lbName] = link;
EleIndex ele(processorName, link);
lbName2OmtfIndex[lbName].push_back(ele);
}
}
void RpcLinkMap::init(const std::string& fName) {
std::ifstream inFile;
inFile.open(fName);
if (inFile) {
LogTrace("") << " reading OmtfRpcLinksMap from: " << fName;
} else {
LogTrace("") << " Unable to open file " << fName;
throw std::runtime_error("Unable to open OmtfRpcLinksMap file " + fName);
}
std::string line;
while (std::getline(inFile, line)) {
line.erase(0, line.find_first_not_of(" \t\r\n")); //cut first character
if (line.empty() || !line.compare(0, 2, "--"))
continue; // empty or comment line
std::stringstream ss(line);
std::string processorName, lbName;
unsigned int link, dbId;
if (ss >> processorName >> link >> lbName >> dbId) {
std::map<unsigned int, std::string>& li2lb = link2lbName[processorName];
std::map<std::string, unsigned int>& lb2li = lbName2link[processorName];
li2lb[link] = lbName;
lb2li[lbName] = link;
EleIndex ele(processorName, link);
lbName2OmtfIndex[lbName].push_back(ele);
}
}
inFile.close();
}
std::vector<EleIndex> RpcLinkMap::omtfEleIndex(const std::string& lbName) const {
const auto pos = lbName2OmtfIndex.find(lbName);
return pos != lbName2OmtfIndex.end() ? pos->second : std::vector<EleIndex>();
}
} // namespace omtf
|