Belle II Software development
KLMElectronicsMapImporter.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9/* Own header. */
10#include <klm/calibration/KLMElectronicsMapImporter.h>
11
12/* KLM headers. */
13#include <klm/dataobjects/KLMChannelIndex.h>
14#include <klm/dbobjects/KLMElectronicsMap.h>
15
16/* Basf2 headers. */
17#include <framework/database/DBImportObjPtr.h>
18#include <framework/database/IntervalOfValidity.h>
19#include <framework/logging/Logger.h>
20#include <rawdata/dataobjects/RawCOPPERFormat.h>
21
22using namespace Belle2;
23
28
32
33void KLMElectronicsMapImporter::setIOV(int experimentLow, int runLow,
34 int experimentHigh, int runHigh)
35{
36 m_ExperimentLow = experimentLow;
37 m_RunLow = runLow;
38 m_ExperimentHigh = experimentHigh;
39 m_RunHigh = runHigh;
40}
41
46
47void KLMElectronicsMapImporter::loadBKLMElectronicsMap(int version)
48{
49 const int minimalVersion = 1;
50 const int maximalVersion = 2;
51 if (version < minimalVersion || version > maximalVersion) {
52 B2FATAL("Incorrect version (" << version << ") of BKLM electronics map. "
53 "It must be from " << minimalVersion << " to " << maximalVersion);
54 }
55 int copperId = 0;
56 int slotId = 0;
57 int laneId;
58 int axisId = 0;
59 KLMChannelIndex bklmPlanes(KLMChannelIndex::c_IndexLevelPlane);
60 for (KLMChannelIndex bklmPlane = bklmPlanes.beginBKLM();
61 bklmPlane != bklmPlanes.endBKLM(); ++bklmPlane) {
62 int section = bklmPlane.getSection();
63 int sector = bklmPlane.getSector();
64 int layer = bklmPlane.getLayer();
65 int plane = bklmPlane.getPlane();
66
67 if (section == BKLMElementNumbers::c_ForwardSection) {
68 if (sector == 3 || sector == 4 || sector == 5 || sector == 6)
69 copperId = 1 + BKLM_ID;
70 if (sector == 1 || sector == 2 || sector == 7 || sector == 8)
71 copperId = 2 + BKLM_ID;
72 }
73 if (section == BKLMElementNumbers::c_BackwardSection) {
74 if (sector == 3 || sector == 4 || sector == 5 || sector == 6)
75 copperId = 3 + BKLM_ID;
76 if (sector == 1 || sector == 2 || sector == 7 || sector == 8)
77 copperId = 4 + BKLM_ID;
78 }
79
80 if (sector == 3 || sector == 4 || sector == 5 || sector == 6)
81 slotId = sector - 2;
82 if (sector == 1 || sector == 2)
83 slotId = sector + 2;
84 if (sector == 7 || sector == 8)
85 slotId = sector - 6;
86
87 if (layer >= BKLMElementNumbers::c_FirstRPCLayer) {
88 laneId = layer + 5;
89 axisId = plane;
90 } else {
91 laneId = layer;
92 if (plane == BKLMElementNumbers::c_ZPlane)
93 axisId = 1;
94 if (plane == BKLMElementNumbers::c_PhiPlane)
95 axisId = 0;
96 }
97
98 int MaxiChannel = BKLMElementNumbers::getNStrips(
99 section, sector, layer, plane);
100
101 bool dontFlip = false;
102 if (section == BKLMElementNumbers::c_ForwardSection &&
103 (sector == 7 || sector == 8 || sector == 1 || sector == 2))
104 dontFlip = true;
105 if (section == BKLMElementNumbers::c_BackwardSection &&
106 (sector == 4 || sector == 5 || sector == 6 || sector == 7))
107 dontFlip = true;
108
109 for (int iStrip = 1; iStrip <= MaxiChannel; iStrip++) {
110 int channelId = iStrip;
111
112 if (!(dontFlip && layer >= BKLMElementNumbers::c_FirstRPCLayer && plane == BKLMElementNumbers::c_PhiPlane))
113 channelId = MaxiChannel - iStrip + 1;
114
115 if (plane == BKLMElementNumbers::c_PhiPlane) {
116 if (layer == 1)
117 channelId += 4;
118 if (layer == 2)
119 channelId += 2;
120 }
121
122 if (plane == BKLMElementNumbers::c_ZPlane) {
123 if (layer < BKLMElementNumbers::c_FirstRPCLayer) {
124 int channelCheck = channelId;
125 if (section == BKLMElementNumbers::c_BackwardSection
126 && sector == BKLMElementNumbers::c_ChimneySector) {
127 if (layer == 1) {
128 if (version == 1) {
129 if (channelCheck > 0 && channelCheck < 9)
130 channelId = 9 - channelId;
131 if (channelCheck > 8 && channelCheck < 24)
132 channelId = 54 - channelId;
133 if (channelCheck > 23 && channelCheck < 39)
134 channelId = 54 - channelId;
135 } else {
136 if (channelCheck > 0 && channelCheck < 9)
137 channelId = 9 - channelId; // 8 : 1
138 if (channelCheck > 8 && channelCheck < 24)
139 channelId = 39 - channelId; // 30 : 16
140 if (channelCheck > 23 && channelCheck < 39)
141 channelId = 69 - channelId; // 45 : 31
142 }
143 }
144 if (layer == 2) {
145 if (channelCheck > 0 && channelCheck < 10)
146 channelId = 10 - channelId;
147 if (channelCheck > 9 && channelCheck < 24)
148 channelId = 40 - channelId;
149 if (channelCheck > 23 && channelCheck < 39)
150 channelId = 69 - channelId;
151 }
152 } else { // All the sectors except the chimney one
153 if (channelCheck > 0 && channelCheck < 10)
154 channelId = 10 - channelId;
155 if (channelCheck > 9 && channelCheck < 25)
156 channelId = 40 - channelId;
157 if (channelCheck > 24 && channelCheck < 40)
158 channelId = 70 - channelId;
159 if (channelCheck > 39 && channelCheck < 55)
160 channelId = 100 - channelId;
161 }
162 }
163 }
164
165 KLMChannelNumber detectorChannel = m_ElementNumbers->channelNumberBKLM(
166 section, sector, layer, plane, iStrip);
167 m_ChannelMap.insert(
168 std::pair<KLMChannelNumber, KLMElectronicsChannel>(
169 detectorChannel,
170 KLMElectronicsChannel(copperId, slotId, laneId, axisId, channelId)));
171 }
172 }
173}
174
175int KLMElectronicsMapImporter::getEKLMStripFirmwareBySoftware(int stripSoftware) const
176{
177 int segment, strip;
178 /* Segment is 0-based here. */
179 segment = (stripSoftware - 1) / EKLMElementNumbers::getNStripsSegment();
180 /* Order of segment readout boards in the firmware is opposite. */
181 segment = 4 - segment;
182 strip = segment * EKLMElementNumbers::getNStripsSegment() +
183 (stripSoftware - 1) % EKLMElementNumbers::getNStripsSegment() + 1;
184 return strip;
185}
186
187void KLMElectronicsMapImporter::addEKLMLane(
188 int section, int sector, int layer, int copper, int slot, int lane)
189{
190 for (int plane = 1; plane <= EKLMElementNumbers::getMaximalPlaneNumber(); ++plane) {
191 for (int strip = 1; strip <= EKLMElementNumbers::getMaximalStripNumber(); ++strip) {
192 int axis = plane - 1;
193 int channel = getEKLMStripFirmwareBySoftware(strip);
194 KLMChannelNumber detectorChannel = m_ElementNumbers->channelNumberEKLM(
195 section, sector, layer, plane, strip);
196 m_ChannelMap.insert(
197 std::pair<KLMChannelNumber, KLMElectronicsChannel>(
198 detectorChannel,
199 KLMElectronicsChannel(EKLM_ID + copper, slot, lane, axis, channel)));
200 }
201 }
202}
203
204void KLMElectronicsMapImporter::setChannelsEKLMSegment(
205 int section, int sector, int layer, int plane, int segment,
206 int firmwareSegment)
207{
208 std::map<KLMChannelNumber, KLMElectronicsChannel>:: iterator it;
209 for (int strip = 1; strip <= EKLMElementNumbers::getNStripsSegment(); ++strip) {
210 int stripPlane = strip + EKLMElementNumbers::getNStripsSegment() * (segment - 1);
211 KLMChannelNumber detectorChannel = m_ElementNumbers->channelNumberEKLM(
212 section, sector, layer, plane, stripPlane);
213 it = m_ChannelMap.find(detectorChannel);
214 if (it == m_ChannelMap.end())
215 B2FATAL("The KLM electronics map is not loaded or incomplete.");
216 int channel = (firmwareSegment - 1) * EKLMElementNumbers::getNStripsSegment() +
217 strip;
218 it->second.setChannel(channel);
219 }
220}
221
222void KLMElectronicsMapImporter::loadEKLMElectronicsMap(int version, bool mc)
223{
224 const int minimalVersion = 1;
225 const int maximalVersion = 2;
226 if (version < minimalVersion || version > maximalVersion) {
227 B2FATAL("Incorrect version (" << version << ") of EKLM electronics map. "
228 "It must be from " << minimalVersion << " to " << maximalVersion);
229 }
230 /* Backward section. */
231 /* Sector 1. */
232 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 1, 3, 1, 1);
233 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 2, 3, 1, 2);
234 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 3, 3, 1, 3);
235 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 4, 3, 1, 4);
236 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 5, 3, 1, 5);
237 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 6, 3, 1, 6);
238 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 7, 3, 2, 1);
239 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 8, 3, 2, 2);
240 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 9, 3, 2, 3);
241 /* Wrong connection was fixed between phase 2 and phase 3. */
242 if (mc || (version >= 2)) {
243 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 10, 3, 2, 4);
244 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 11, 3, 2, 5);
245 } else {
246 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 10, 3, 2, 5);
247 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 11, 3, 2, 4);
248 }
249 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 1, 12, 3, 2, 6);
250 /* Sector 2. */
251 if (version == 1) {
252 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 1, 3, 3, 1);
253 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 2, 3, 3, 2);
254 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 3, 3, 3, 3);
255 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 4, 3, 3, 4);
256 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 5, 3, 3, 5);
257 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 6, 3, 3, 6);
258 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 7, 3, 4, 1);
259 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 8, 3, 4, 2);
260 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 9, 3, 4, 3);
261 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 10, 3, 4, 4);
262 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 11, 3, 4, 5);
263 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 12, 3, 4, 6);
264 } else {
265 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 1, 3, 4, 6);
266 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 2, 3, 4, 5);
267 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 3, 3, 4, 4);
268 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 4, 3, 4, 3);
269 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 5, 3, 4, 2);
270 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 6, 3, 4, 1);
271 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 7, 3, 3, 6);
272 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 8, 3, 3, 5);
273 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 9, 3, 3, 4);
274 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 10, 3, 3, 3);
275 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 11, 3, 3, 2);
276 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 2, 12, 3, 3, 1);
277 }
278 /* Sector 3. */
279 if (version == 1) {
280 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 1, 4, 1, 1);
281 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 2, 4, 1, 2);
282 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 3, 4, 1, 3);
283 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 4, 4, 1, 4);
284 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 5, 4, 1, 5);
285 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 6, 4, 1, 6);
286 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 7, 4, 2, 1);
287 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 8, 4, 2, 2);
288 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 9, 4, 2, 3);
289 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 10, 4, 2, 4);
290 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 11, 4, 2, 5);
291 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 12, 4, 2, 6);
292 } else {
293 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 1, 4, 2, 6);
294 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 2, 4, 2, 5);
295 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 3, 4, 2, 4);
296 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 4, 4, 2, 3);
297 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 5, 4, 2, 2);
298 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 6, 4, 2, 1);
299 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 7, 4, 1, 6);
300 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 8, 4, 1, 5);
301 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 9, 4, 1, 4);
302 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 10, 4, 1, 3);
303 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 11, 4, 1, 2);
304 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 3, 12, 4, 1, 1);
305 }
306 /* Sector 4. */
307 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 1, 4, 3, 1);
308 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 2, 4, 3, 2);
309 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 3, 4, 3, 3);
310 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 4, 4, 3, 4);
311 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 5, 4, 3, 5);
312 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 6, 4, 3, 6);
313 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 7, 4, 4, 1);
314 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 8, 4, 4, 2);
315 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 9, 4, 4, 3);
316 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 10, 4, 4, 4);
317 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 11, 4, 4, 5);
318 addEKLMLane(EKLMElementNumbers::c_BackwardSection, 4, 12, 4, 4, 6);
319 /* Forward section. */
320 /* Sector 1. */
321 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 1, 1, 3, 1);
322 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 2, 1, 3, 2);
323 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 3, 1, 3, 3);
324 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 4, 1, 3, 4);
325 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 5, 1, 3, 5);
326 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 6, 1, 3, 6);
327 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 7, 1, 3, 7);
328 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 8, 1, 4, 1);
329 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 9, 1, 4, 2);
330 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 10, 1, 4, 3);
331 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 11, 1, 4, 4);
332 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 12, 1, 4, 5);
333 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 13, 1, 4, 6);
334 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 1, 14, 1, 4, 7);
335 /* Sector 2. */
336 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 1, 1, 1, 1);
337 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 2, 1, 1, 2);
338 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 3, 1, 1, 3);
339 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 4, 1, 1, 4);
340 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 5, 1, 1, 5);
341 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 6, 1, 1, 6);
342 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 7, 1, 1, 7);
343 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 8, 1, 2, 1);
344 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 9, 1, 2, 2);
345 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 10, 1, 2, 3);
346 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 11, 1, 2, 4);
347 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 12, 1, 2, 5);
348 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 13, 1, 2, 6);
349 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 2, 14, 1, 2, 7);
350 /* Sector 3. */
351 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 1, 2, 3, 1);
352 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 2, 2, 3, 2);
353 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 3, 2, 3, 3);
354 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 4, 2, 3, 4);
355 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 5, 2, 3, 5);
356 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 6, 2, 3, 6);
357 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 7, 2, 3, 7);
358 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 8, 2, 4, 1);
359 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 9, 2, 4, 2);
360 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 10, 2, 4, 3);
361 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 11, 2, 4, 4);
362 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 12, 2, 4, 5);
363 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 13, 2, 4, 6);
364 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 3, 14, 2, 4, 7);
365 /* Sector 4. */
366 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 1, 2, 1, 1);
367 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 2, 2, 1, 2);
368 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 3, 2, 1, 3);
369 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 4, 2, 1, 4);
370 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 5, 2, 1, 5);
371 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 6, 2, 1, 6);
372 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 7, 2, 1, 7);
373 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 8, 2, 2, 1);
374 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 9, 2, 2, 2);
375 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 10, 2, 2, 3);
376 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 11, 2, 2, 4);
377 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 12, 2, 2, 5);
378 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 13, 2, 2, 6);
379 addEKLMLane(EKLMElementNumbers::c_ForwardSection, 4, 14, 2, 2, 7);
380 /* Switch of internal cables. */
381 if (!mc) {
382 setChannelsEKLMSegment(1, 1, 5, 1, 1, 4);
383 setChannelsEKLMSegment(1, 1, 5, 1, 2, 5);
384 setChannelsEKLMSegment(1, 1, 5, 2, 1, 4);
385 setChannelsEKLMSegment(1, 1, 5, 2, 2, 5);
386 setChannelsEKLMSegment(2, 1, 10, 2, 3, 2);
387 setChannelsEKLMSegment(2, 1, 10, 2, 4, 3);
388 }
389}
390
391void KLMElectronicsMapImporter::setLane(
392 int subdetector, int section, int sector, int layer, int lane)
393{
394 std::map<KLMChannelNumber, KLMElectronicsChannel>::iterator it;
395 int minimalPlane = m_ElementNumbers->getMinimalPlaneNumber(subdetector);
396 KLMChannelIndex klmChannel(subdetector, section, sector, layer, minimalPlane, 1);
397 KLMChannelIndex klmModule(klmChannel);
398 klmModule.setIndexLevel(KLMChannelIndex::c_IndexLevelLayer);
399 KLMChannelIndex klmNextModule(klmModule);
400 ++klmNextModule;
401 for (; klmChannel != klmNextModule; ++klmChannel) {
402 KLMChannelNumber channel = klmChannel.getKLMChannelNumber();
403 it = m_ChannelMap.find(channel);
404 if (it == m_ChannelMap.end())
405 B2FATAL("The KLM electronics map is not loaded or incomplete.");
406 it->second.setLane(lane);
407 }
408}
409
410void KLMElectronicsMapImporter::setLane(
411 int subdetector, int section, int sector, int layer, int plane, int lane)
412{
413 std::map<KLMChannelNumber, KLMElectronicsChannel>::iterator it;
414 KLMChannelIndex klmChannel(subdetector, section, sector, layer, plane, 1);
415 KLMChannelIndex klmPlane(klmChannel);
416 klmPlane.setIndexLevel(KLMChannelIndex::c_IndexLevelPlane);
417 KLMChannelIndex klmNextPlane(klmPlane);
418 ++klmNextPlane;
419 for (; klmChannel != klmNextPlane; ++klmChannel) {
420 KLMChannelNumber channel = klmChannel.getKLMChannelNumber();
421 it = m_ChannelMap.find(channel);
422 if (it == m_ChannelMap.end())
423 B2FATAL("The KLM electronics map is not loaded or incomplete.");
424 it->second.setLane(lane);
425 }
426}
427
428void KLMElectronicsMapImporter::importElectronicsMap()
429{
430 DBImportObjPtr<KLMElectronicsMap> electronicsMap;
431 electronicsMap.construct();
432 std::map<KLMChannelNumber, KLMElectronicsChannel>::iterator it;
433 for (it = m_ChannelMap.begin(); it != m_ChannelMap.end(); ++it) {
434 electronicsMap->addChannel(
435 it->first,
436 it->second.getCopper(),
437 it->second.getSlot(),
438 it->second.getLane(),
439 it->second.getAxis(),
440 it->second.getChannel());
441 }
442 IntervalOfValidity iov(m_ExperimentLow, m_RunLow,
443 m_ExperimentHigh, m_RunHigh);
444 electronicsMap.import(iov);
445}
Belle2::BKLMElementNumbers::c_ChimneySector
@ c_ChimneySector
Chimney sector: BB3 in 1-based notation; BB2 in 0-based notation.
Definition BKLMElementNumbers.h:54
Belle2::BKLMElementNumbers::c_FirstRPCLayer
@ c_FirstRPCLayer
First RPC layer.
Definition BKLMElementNumbers.h:64
Belle2::BKLMElementNumbers::getNStrips
static int getNStrips(int section, int sector, int layer, int plane)
Get number of strips.
Definition BKLMElementNumbers.cc:107
Belle2::DBImportBase::import
bool import(const IntervalOfValidity &iov)
Import the object to database.
Definition DBImportBase.cc:36
Belle2::DBImportObjPtr
Class for importing a single object to the database.
Definition DBImportObjPtr.h:23
Belle2::DBImportObjPtr::construct
void construct(Args &&... params)
Construct an object of type T in this DBImportObjPtr using the provided constructor arguments.
Definition DBImportObjPtr.h:47
Belle2::EKLMElementNumbers::getMaximalStripNumber
static constexpr int getMaximalStripNumber()
Get maximal strip number.
Definition EKLMElementNumbers.h:319
Belle2::EKLMElementNumbers::getNStripsSegment
static constexpr int getNStripsSegment()
Get number of strips in a segment.
Definition EKLMElementNumbers.h:391
Belle2::EKLMElementNumbers::getMaximalPlaneNumber
static constexpr int getMaximalPlaneNumber()
Get maximal plane number.
Definition EKLMElementNumbers.h:303
Belle2::IntervalOfValidity
A class that describes the interval of experiments/runs for which an object in the database is valid.
Definition IntervalOfValidity.h:25
Belle2::KLMChannelIndex
KLM channel index.
Definition KLMChannelIndex.h:23
Belle2::KLMChannelIndex::beginBKLM
KLMChannelIndex beginBKLM()
First channel for BKLM.
Definition KLMChannelIndex.cc:187
Belle2::KLMChannelIndex::getKLMChannelNumber
KLMChannelNumber getKLMChannelNumber() const
Get KLM channel number.
Definition KLMChannelIndex.cc:143
Belle2::KLMChannelIndex::setIndexLevel
void setIndexLevel(enum IndexLevel indexLevel)
Set index level.
Definition KLMChannelIndex.cc:65
Belle2::KLMChannelIndex::endBKLM
KLMChannelIndex & endBKLM()
Last channel for BKLM.
Definition KLMChannelIndex.cc:193
Belle2::KLMElectronicsMapImporter::m_ElementNumbers
const KLMElementNumbers * m_ElementNumbers
Element numbers.
Definition KLMElectronicsMapImporter.h:138
Belle2::KLMElectronicsMapImporter::setChannelsEKLMSegment
void setChannelsEKLMSegment(int section, int sector, int layer, int plane, int segment, int firmwareSegment)
Set channels for EKLM segment.
Definition KLMElectronicsMapImporter.cc:204
Belle2::KLMElectronicsMapImporter::setIOV
void setIOV(int experimentLow, int runLow, int experimentHigh, int runHigh)
Set interval of validity.
Definition KLMElectronicsMapImporter.cc:33
Belle2::KLMElectronicsMapImporter::clearElectronicsMap
void clearElectronicsMap()
Clear electronics map (to be able to import its multiple versions).
Definition KLMElectronicsMapImporter.cc:42
Belle2::KLMElectronicsMapImporter::addEKLMLane
void addEKLMLane(int section, int sector, int layer, int copper, int slot, int lane)
Add EKLM electronics map lane.
Definition KLMElectronicsMapImporter.cc:187
Belle2::KLMElectronicsMapImporter::setLane
void setLane(int subdetector, int section, int sector, int layer, int lane)
Set non-default lane for all channels in a module.
Definition KLMElectronicsMapImporter.cc:391
Belle2::KLMElectronicsMapImporter::loadEKLMElectronicsMap
void loadEKLMElectronicsMap(int version, bool mc)
Load EKLM electronics map.
Definition KLMElectronicsMapImporter.cc:222
Belle2::KLMElectronicsMapImporter::loadBKLMElectronicsMap
void loadBKLMElectronicsMap(int version)
Load BKLM electronics map.
Definition KLMElectronicsMapImporter.cc:47
Belle2::KLMElectronicsMapImporter::m_ChannelMap
std::map< KLMChannelNumber, KLMElectronicsChannel > m_ChannelMap
Data for creation of the electronics map.
Definition KLMElectronicsMapImporter.h:141
Belle2::KLMElectronicsMapImporter::importElectronicsMap
void importElectronicsMap()
Import electronics map.
Definition KLMElectronicsMapImporter.cc:428
Belle2::KLMElectronicsMapImporter::getEKLMStripFirmwareBySoftware
int getEKLMStripFirmwareBySoftware(int stripSoftware) const
Get EKLM firmware strip number by software strip number.
Definition KLMElectronicsMapImporter.cc:175
Belle2::KLMChannelNumber
uint16_t KLMChannelNumber
Channel number.
Definition KLMElementNumberDefinitions.h:36
Belle2
Abstract base class for different kinds of events.
Definition MillepedeAlgorithm.h:17