Belle II Software development
PXDUnpackerOTModule.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#include <pxd/modules/pxdUnpacking/PXDUnpackerOTModule.h>
10#include <pxd/unpacking/PXDRawDataDefinitions.h>
11#include <pxd/unpacking/PXDRawDataStructs.h>
12#include <pxd/unpacking/PXDMappingLookup.h>
13#include <pxd/dataobjects/PXDRawHit.h>
14#include <pxd/dataobjects/PXDRawAdc.h>
15#include <pxd/dataobjects/PXDRawROIs.h>
16#include <vxd/dataobjects/VxdID.h>
17#include <rawdata/dataobjects/RawPXD.h>
18#include <framework/dataobjects/EventMetaData.h>
19#include <framework/logging/Logger.h>
20
21using namespace std;
22using namespace Belle2;
23using namespace Belle2::PXD;
24using namespace Belle2::PXD::PXDError;
25
26//-----------------------------------------------------------------
27// Register the Module
28//-----------------------------------------------------------------
29REG_MODULE(PXDUnpackerOT);
30
31//-----------------------------------------------------------------
32// Implementation
33//-----------------------------------------------------------------
34
38
40 Module(),
44{
45 //Set module properties
46 setDescription("Unpack Raw PXD Hits from ONSEN data stream");
48
49 addParam("RawPXDsName", m_RawPXDsName, "The name of the StoreArray of RawPXDs to be processed", std::string(""));
50 addParam("PXDRawHitsName", m_PXDRawHitsName, "The name of the StoreArray of generated PXDRawHits", std::string(""));
51 addParam("PXDDAQEvtStatsName", m_PXDDAQEvtStatsName, "The name of the StoreObjPtr of generated PXDDAQEvtStats", std::string(""));
52 addParam("PXDRawAdcsName", m_PXDRawAdcsName, "The name of the StoreArray of generated PXDRawAdcs", std::string(""));
53 addParam("PXDRawROIsName", m_PXDRawROIsName, "The name of the StoreArray of generated PXDRawROIs", std::string(""));
54 addParam("DoNotStore", m_doNotStore, "only unpack and check, but do not store", false);
55// addParam("CriticalErrorMask", m_criticalErrorMask, "Set error mask which stops processing by returning false by task", 0);
56 // m_criticalErrorMask not longer a aparameter
57 // the constructor default to zero anyway
58// addParam("SuppressErrorMask", m_suppressErrorMask, "Set mask for errors msgs which are not printed", getSilenceMask());
59 // m_suppressErrorMask not longer a aparameter
60 m_suppressErrorMask = getSilenceMask();
61 addParam("ForceMapping", m_forceMapping, "Force Mapping even if DHH bit is NOT requesting it", false);
62 addParam("ForceNoMapping", m_forceNoMapping, "Force NO Mapping even if DHH bit is requesting it", false);
63 addParam("CheckPaddingCRC", m_checkPaddingCRC, "Check for susp. padding (debug option, many false positive)", false);
64 addParam("MaxDHPFrameDiff", m_maxDHPFrameDiff, "Maximum DHP Frame Nr Difference w/o reporting error", 2u);
65 addParam("FormatBonnDAQ", m_formatBonnDAQ, "ONSEN or BonnDAQ format", false);
66 addParam("Verbose", m_verbose, "Turn on extra verbosity for log-level debug", false);
67 addParam("ContinueOnError", m_continueOnError, "Continue package depacking on error (for debugging)", false);
68// (
69// /*EPXDErrFlag::c_DHC_END | EPXDErrFlag::c_DHE_START | EPXDErrFlag::c_DATA_OUTSIDE |*/
70// EPXDErrFlag::c_FIX_SIZE | EPXDErrFlag::c_DHE_CRC | EPXDErrFlag::c_DHC_UNKNOWN | /*EPXDErrFlag::c_MERGER_CRC |*/
71// EPXDErrFlag::c_DHP_SIZE | /*EPXDErrFlag::c_DHP_PIX_WO_ROW | EPXDErrFlag::c_DHE_START_END_ID | EPXDErrFlag::c_DHE_START_ID |*/
72// EPXDErrFlag::c_DHE_START_WO_END | EPXDErrFlag::c_DHP_NOT_CONT
73// ));
74
75 // this is not really a parameter, it should be fixed.
76 m_errorSkipPacketMask[c_nrDHE_CRC] = true;
77 m_errorSkipPacketMask[c_nrFIX_SIZE] = true;
78}
79
81{
82 // Required input
83 m_eventMetaData.isRequired();
84 // Optional input
85 m_storeRawPXD.isOptional(m_RawPXDsName);
86
87 //Register output collections
93
94 B2DEBUG(29, "ForceMapping: " << m_forceMapping);
95 B2DEBUG(29, "ForceNoMapping: " << m_forceNoMapping);
96 B2DEBUG(29, "CheckPaddingCRC: " << m_checkPaddingCRC);
97 B2DEBUG(29, "MaxDHPFrameDiff: " << m_maxDHPFrameDiff);
98
100 m_sendrois = 0;
101 m_notaccepted = 0;
103 for (int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) m_errorCounter[i] = 0;
104
105}
106
108{
109 int flag = 0;
110 string errstr = "Statistic ( ;";
111 errstr += to_string(m_unpackedEventsCount) + ";";
112 for (int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) { errstr += to_string(m_errorCounter[i]) + ";"; flag |= m_errorCounter[i];}
113 if (flag != 0) {
114 B2RESULT("PXD Unpacker --> Error Statistics (counted once per event!) in Events: " << m_unpackedEventsCount);
115 B2RESULT(errstr + " )");
116 for (int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) {
117 if (m_errorCounter[i]) {
118 B2RESULT(getPXDBitErrorName(i) << ": " << m_errorCounter[i]);
119 }
120 }
121 } else {
122 B2RESULT("PXD Unpacker --> No Error found in Events: " << m_unpackedEventsCount);
123 }
124 B2RESULT("Statistic 2: !Accepted: " << m_notaccepted << " SendROIs: " << m_sendrois << " Unfiltered: " << m_sendunfiltered);
125}
126
128{
129 m_storeDAQEvtStats.create();
130
131 m_errorMask = 0;
133
134 m_meta_event_nr = m_eventMetaData->getEvent();// used for error output below
135
136 if (!m_storeRawPXD) {// if no input, nothing to do
137 m_errorMask[c_nrNO_PXD] = true;
138 } else {
139 int nRaws = m_storeRawPXD.getEntries();
140 if (m_verbose) {
141 B2DEBUG(29, "PXD Unpacker --> RawPXD Objects in event: " << LogVar("Objects", nRaws));
142 };
143
144 m_meta_run_nr = m_eventMetaData->getRun();
145 m_meta_subrun_nr = m_eventMetaData->getSubrun();
146 m_meta_experiment = m_eventMetaData->getExperiment();
147 m_meta_time = m_eventMetaData->getTime();
148 m_meta_ticks = (unsigned int)std::round((m_meta_time % 1000000000ull) * 0.127216); // calculate ticks in 127MHz RF clock
149 m_meta_sec = (unsigned int)(m_meta_time / 1000000000ull) & 0x1FFFF;
150
151 int inx = 0; // count index for output objects
152 for (auto& it : m_storeRawPXD) {
153 if (m_verbose) {
154 B2DEBUG(29, "PXD Unpacker --> Unpack Objects: ");
155 };
156 unpack_rawpxd(it, inx++);
157 }
158
159 if (nRaws == 0) m_errorMask[c_nrNO_PXD] = true;
160 }
163
165 {
166 for (unsigned int i = 0; i < ONSEN_MAX_TYPE_ERR; i++) {
167 if (m_errorMaskEvent[i]) m_errorCounter[i]++;
168 }
169 }
170
171 if ((PXDErrorFlags(m_criticalErrorMask) & m_errorMaskEvent) != PXDErrorFlags(0)) B2ERROR("Error in PXD unpacking" <<
172 LogVar("event nr", m_meta_event_nr));
173 setReturnValue(PXDErrorFlags(0) == (PXDErrorFlags(m_criticalErrorMask) & m_errorMaskEvent));
174}
175
177{
178 int Frames_in_event;
179 int fullsize;
180 int datafullsize;
181
182 m_errorMaskDHE = 0;
183 m_errorMaskDHC = 0;
185 PXDDAQPacketStatus& daqpktstat = m_storeDAQEvtStats->newPacket(inx);
186
187 if (px.size() <= 0 || px.size() > 16 * 1024 * 1024) {
188 if (!(m_suppressErrorMask[c_nrPACKET_SIZE])) {
189 B2WARNING("PXD Unpacker --> invalid packet size" <<
190 LogVar("size [32bit words] $", static_cast < std::ostringstream && >(std::ostringstream() << hex << px.size()).str()));
191 }
192 m_errorMask[c_nrPACKET_SIZE] = true;
193 return;
194 }
195 std::vector<unsigned int> data(px.size());
196 fullsize = px.size() * 4;
197 std::copy_n(px.data(), px.size(), data.begin());
198
199 if (fullsize < 8) {
200 if (!(m_suppressErrorMask[c_nrPACKET_SIZE])) {
201 B2WARNING("Data is to small to hold a valid Header! Will not unpack anything." << LogVar("size [32bit words] $",
202 static_cast < std::ostringstream && >(std::ostringstream() << hex << fullsize).str()));
203 }
204 m_errorMask[c_nrPACKET_SIZE] = true;
205 return;
206 }
207
208 if (data[0] != 0xCAFEBABE && data[0] != 0xBEBAFECA) {
209 if (!(m_suppressErrorMask[c_nrMAGIC])) {
210 B2WARNING("Magic invalid: Will not unpack anything. Header corrupted." <<
211 LogVar("Header Magic $", static_cast < std::ostringstream && >(std::ostringstream() << hex << data[0]).str()));
212 }
213 m_errorMask[c_nrMAGIC] = true;
214 return;
215 }
216
217
218 Frames_in_event = ((ubig32_t*)data.data())[1];
219 if (Frames_in_event < 0 || Frames_in_event > 256) {
220 if (!(m_suppressErrorMask[c_nrFRAME_NR])) {
221 B2WARNING("Number of Frames invalid: Will not unpack anything. Header corrupted!" << LogVar("Frames in event", Frames_in_event));
222 }
223 m_errorMask[c_nrFRAME_NR] = true;
224 return;
225 }
226 if (Frames_in_event < 3) {
227 if (!(m_suppressErrorMask[c_nrNR_FRAMES_TO_SMALL])) {
228 B2WARNING("Number of Frames too small: It cannot contain anything useful." << LogVar("Frames in event", Frames_in_event));
229 }
230 m_errorMask[c_nrNR_FRAMES_TO_SMALL] = true;
231 }
232
234 if (m_verbose) {
235 B2DEBUG(29, "PXD Unpacker --> data[0]: <-- Magic $" << hex << data[0]);
236 B2DEBUG(29, "PXD Unpacker --> data[1]: <-- #Frames $" << hex << data[1]);
237 if (data[1] >= 1 && fullsize < 12) B2DEBUG(29, "PXD Unpacker --> data[2]: <-- Frame 1 len $" << hex << data[2]);
238 if (data[1] >= 2 && fullsize < 16) B2DEBUG(29, "PXD Unpacker --> data[3]: <-- Frame 2 len $" << hex << data[3]);
239 if (data[1] >= 3 && fullsize < 20) B2DEBUG(29, "PXD Unpacker --> data[4]: <-- Frame 3 len $" << hex << data[4]);
240 if (data[1] >= 4 && fullsize < 24) B2DEBUG(29, "PXD Unpacker --> data[5]: <-- Frame 4 len $" << hex << data[5]);
241 };
242
243 unsigned int* tableptr;
244 tableptr = &data[2]; // skip header!!!
245
246 unsigned int* dataptr;
247 dataptr = &tableptr[Frames_in_event];
248 datafullsize = fullsize - 2 * 4 - Frames_in_event * 4; // Size is fullsize minus header minus table
249
250 int ll = 0; // Offset in dataptr in bytes
251 for (int j = 0; j < Frames_in_event; j++) {
252 int lo;
253
254 lo = ((ubig32_t*)tableptr)[j];
255 if (lo <= 0) {
256 if (!(m_suppressErrorMask[c_nrFRAME_SIZE])) {
257 B2WARNING("size of frame invalid");
258 B2DEBUG(29, "size of frame invalid: " << j << "size " << lo << " at byte offset in dataptr " << ll);
259 }
260 m_errorMask[c_nrFRAME_SIZE] = true;
261 return;
262 }
263 if (ll + lo > datafullsize) {
264 if (!(m_suppressErrorMask[c_nrFRAME_SIZE])) {
265 B2WARNING("Frames exceed packet size");
266 B2DEBUG(29, "Frames exceed packet size: " << j << " size " << lo << " at byte offset in dataptr " << ll << " of datafullsize " <<
267 datafullsize << " of fullsize " << fullsize);
268 }
269 m_errorMask[c_nrFRAME_SIZE] = true;
270 return;
271 }
272 if (lo & 0x3) {
273 if (!(m_suppressErrorMask[c_nrFRAME_SIZE])) {
274 B2WARNING("SKIP Frame with Data with not MOD 4 length");
275 B2DEBUG(29, "SKIP Frame with Data with not MOD 4 length " << " ( $" << hex << lo << " ) ");
276 }
277 ll += (lo + 3) & 0xFFFFFFFC;
278 m_errorMask[c_nrFRAME_SIZE] = true;
279 } else {
280 B2DEBUG(29, "unpack DHE(C) frame: " << j << " with size " << lo << " at byte offset in dataptr " << ll);
281 unpack_dhc_frame(ll + (char*)dataptr, lo, j, Frames_in_event, daqpktstat);
282 ll += lo;
283 }
288 m_errorMask = 0;
289
290 if (!m_continueOnError && (m_errorMaskPacket & PXDErrorFlags(m_errorSkipPacketMask)) != PXDErrorFlags(0)) {
291 // skip full package on error, recovery to next DHC/DHE Start might be possible in some cases
292 // But that's to hard to implement
293 // Remark: PXD data for broken events is removed in next PXDPostChecker module, thus skipping the
294 // unpacking is not strictly necessary here.
295 break;
296 }
297 }
299}
300
301void PXDUnpackerOTModule::unpack_dhp_raw(void* data, unsigned int frame_len, unsigned int dhe_ID, unsigned dhe_DHPport,
302 VxdID vxd_id)
303{
304// unsigned int nr_words = frame_len / 2; // frame_len in bytes (excl. CRC)!!!
305 ubig16_t* dhp_pix = (ubig16_t*)data;
306
313
314 // Size: 64*768 + 8 bytes for a full frame readout
315 if (frame_len != 0xC008) {
316 if (!(m_suppressErrorMask[c_nrFIX_SIZE])) B2WARNING("Frame size unsupported for RAW ADC frame! $" <<
317 LogVar("size [bytes] $", static_cast < std::ostringstream && >(std::ostringstream() << hex << frame_len).str())
318 << LogVar("DHE", dhe_ID) << LogVar("DHP", dhe_DHPport));
319 m_errorMask[c_nrFIX_SIZE] = true;
320 return;
321 }
322 unsigned int dhp_header_type = 0;
323// unsigned int dhp_reserved = 0;
324 unsigned int dhp_dhe_id = 0;
325 unsigned int dhp_dhp_id = 0;
326
327 dhp_header_type = (dhp_pix[2] & 0xE000) >> 13;
328// dhp_reserved = (dhp_pix[2] >> 8) & 0x1F;
329 dhp_dhe_id = (dhp_pix[2] & 0x00FC) >> 2;
330 dhp_dhp_id = dhp_pix[2] & 0x0003;
331
332 if (dhe_ID != dhp_dhe_id) {
333 if (!(m_suppressErrorMask[c_nrDHE_DHP_DHEID])) {
334 B2WARNING("DHE ID in DHE and DHP header differ");
335 B2DEBUG(29, "DHE ID in DHE and DHP header differ $" << hex << dhe_ID << " != $" << dhp_dhe_id);
336 }
337 m_errorMask[c_nrDHE_DHP_DHEID] = true;
338 }
339 if (dhe_DHPport != dhp_dhp_id) {
340 if (!(m_suppressErrorMask[c_nrDHE_DHP_PORT])) {
341 B2WARNING("DHP ID (Chip/Port) in DHE and DHP header differ");
342 B2DEBUG(29, "DHP ID (Chip/Port) in DHE and DHP header differ $" << hex << dhe_DHPport << " != $" << dhp_dhp_id);
343 }
344 m_errorMask[c_nrDHE_DHP_PORT] = true;
345 }
346
347 if (dhp_header_type != EDHPFrameHeaderDataType::c_RAW) {
348 if (!(m_suppressErrorMask[c_nrHEADERTYPE_INV])) {
349 B2WARNING("Header type invalid for this kind of DHE frame");
350 B2DEBUG(29, "Header type invalid for this kind of DHE frame: $" << hex << dhp_header_type);
351 }
352 m_errorMask[c_nrHEADERTYPE_INV] = true;
353 return;
354 }
355
357 B2DEBUG(29, "Raw ADC Data");
358 // size checked already above
359 m_storeRawAdc.appendNew(vxd_id, data, frame_len);
360};
361
362void PXDUnpackerOTModule::unpack_fce([[maybe_unused]] unsigned short* data, [[maybe_unused]] unsigned int length,
363 [[maybe_unused]] VxdID vxd_id)
364{
371
372 B2WARNING("FCE (Cluster) Packet have not yet been tested with real HW clusters. Dont assume that this code is working!");
373 return;
374
375 // implement the unpacking here and not as a separate module ... when it is available in HW
376// ubig16_t* cluster = (ubig16_t*)data;
377// int nr_words; //words in dhp frame
378// unsigned int words_in_cluster = 0; //counts 16bit words in cluster
379// nr_words = length / 2;
380// ubig16_t sor;
381// sor = 0x0000;
382//
383// for (int i = 2 ; i < nr_words ; i++) {
384// if (i != 2) { //skip header
385// if ((((cluster[i] & 0x8000) == 0)
386// && ((cluster[i] & 0x4000) >> 14) == 1)) { //searches for start of row frame with start of cluster flag = 1 => new cluster
387// if (!m_doNotStore) m_storeRawCluster.appendNew(&data[i - words_in_cluster], words_in_cluster, vxd_id);
388// words_in_cluster = 0;
389// }
390// }
391// if ((cluster[i] & 0x8000) == 0) {
392// sor = cluster[i];
393// }
394// words_in_cluster++;
395//
396// if ((cluster[nr_words - 1] & 0xFFFF) == (sor &
397// 0xFFFF)) {//if frame is not 32bit aligned last word will be the last start of row word
398// cluster[nr_words - 1] = 0x0000;//overwrites the last redundant word with zero to make checking easier in PXDHardwareClusterUnpacker
399// }
400//
401// if (i == nr_words - 1) {
402// if (!m_doNotStore) m_storeRawCluster.appendNew(&data[i - words_in_cluster + 1], words_in_cluster, vxd_id);
403// }
404// }
405}
406
407void PXDUnpackerOTModule::dump_dhp(void* data, unsigned int frame_len)
408{
409 // called only for debugging purpose, will never be called in normal running
410 unsigned int w = frame_len / 2;
411 ubig16_t* d = (ubig16_t*)data;
412
413 B2WARNING("HEADER -- $" << hex << d[0] << ",$" << hex << d[1] << ",$" << hex << d[2] << ",$" << hex << d[3] << " -- ");
414
415 auto dhp_header_type = (d[2] & 0xE000) >> 13;
416 auto dhp_reserved = (d[2] & 0x1F00) >> 8;
417 auto dhp_dhe_id = (d[2] & 0x00FC) >> 2;
418 auto dhp_dhp_id = d[2] & 0x0003;
419
420 B2WARNING("DHP type | $" << hex << dhp_header_type << " ( " << dec << dhp_header_type << " ) ");
421 B2WARNING("DHP reserved | $" << hex << dhp_reserved << " ( " << dec << dhp_reserved << " ) ");
422 B2WARNING("DHP DHE ID | $" << hex << dhp_dhe_id << " ( " << dec << dhp_dhe_id << " ) ");
423 B2WARNING("DHP DHP ID | $" << hex << dhp_dhp_id << " ( " << dec << dhp_dhp_id << " ) ");
424 for (unsigned int i = 4; i < w; i++) {
425 B2WARNING("DHP DATA $" << hex << d[i]);
426 }
427 B2WARNING("DHP CRC $" << hex << d[w] << ",$" << hex << d[w + 1]);
428}
429
430void PXDUnpackerOTModule::dump_roi(void* data, unsigned int frame_len)
431{
432 // called only for debugging purpose, will never be called in normal running
433 unsigned int w = frame_len / 4;
434 ubig32_t* d = (ubig32_t*)data;
435
436 B2WARNING("HEADER -- $" << hex << d[0] << ",$" << hex << d[1] << ",$" << hex << d[2] << ",$" << hex << d[3] << " -- Len $" << hex
437 << frame_len);
438
439 for (unsigned int i = 0; i < w; i++) {
440 B2WARNING("ROI DATA $" << hex << d[i]);
441 }
442 B2WARNING("ROI CRC $" << hex << d[w]);
443}
444
445void PXDUnpackerOTModule::unpack_dhp(void* data, unsigned int frame_len, unsigned int dhe_first_readout_frame_id_lo,
446 unsigned int dhe_ID, unsigned dhe_DHPport, unsigned dhe_reformat, VxdID vxd_id,
447 PXDDAQPacketStatus& daqpktstat)
448{
449 unsigned int nr_words = frame_len / 2; // frame_len in bytes (excl. CRC)!!!
450 ubig16_t* dhp_pix = (ubig16_t*)data;
451
452 unsigned int dhp_readout_frame_lo = 0;
453 unsigned int dhp_header_type = 0;
454 unsigned int dhp_reserved = 0;
455 unsigned int dhp_dhe_id = 0;
456 unsigned int dhp_dhp_id = 0;
457 unsigned int wrap = 0; // workaround to recalc a relative frame number
458 int last_gate = -1; // workaround to recalc a relative frame number
459
460 // cppcheck-suppress unreadVariable
461 unsigned int dhp_row = 0, dhp_col = 0, dhp_cm = 0;
462// unsigned int dhp_offset = 0;
463 bool rowflag = false;
464 bool pixelflag = true; // just for first row start
465
466 if (nr_words < 4) {
467 if (!(m_suppressErrorMask[c_nrDHP_SIZE])) B2WARNING("DHP frame size error (too small)" << LogVar("Nr words", nr_words));
468 m_errorMask[c_nrDHP_SIZE] = true;
469 return;
470 }
471
472 B2DEBUG(29, "HEADER -- $" << hex << dhp_pix[0] << hex << dhp_pix[1] << hex << dhp_pix[2] << hex << dhp_pix[3] << " -- ");
473
474 B2DEBUG(29, "DHP Header | $" << hex << dhp_pix[2] << " ( " << dec << dhp_pix[2] << " ) ");
475 dhp_header_type = (dhp_pix[2] & 0xE000) >> 13;
476 dhp_reserved = (dhp_pix[2] & 0x1F00) >> 8;
477 dhp_dhe_id = (dhp_pix[2] & 0x00FC) >> 2;
478 dhp_dhp_id = dhp_pix[2] & 0x0003;
479
480 B2DEBUG(29, "DHP type | $" << hex << dhp_header_type << " ( " << dec << dhp_header_type << " ) ");
481 B2DEBUG(29, "DHP reserved | $" << hex << dhp_reserved << " ( " << dec << dhp_reserved << " ) ");
482 B2DEBUG(29, "DHP DHE ID | $" << hex << dhp_dhe_id << " ( " << dec << dhp_dhe_id << " ) ");
483 B2DEBUG(29, "DHP DHP ID | $" << hex << dhp_dhp_id << " ( " << dec << dhp_dhp_id << " ) ");
484
485 if (dhe_ID != dhp_dhe_id) {
486 if (!(m_suppressErrorMask[c_nrDHE_DHP_DHEID])) {
487 B2WARNING("DHE ID in DHE and DHP header differ");
488 B2DEBUG(29, "DHE ID in DHE and DHP header differ $" << hex << dhe_ID << " != $" << dhp_dhe_id);
489 }
490 m_errorMask[c_nrDHE_DHP_DHEID] = true;
491 }
492 if (dhe_DHPport != dhp_dhp_id) {
493 if (!(m_suppressErrorMask[c_nrDHE_DHP_PORT])) {
494 B2WARNING("DHP ID (Chip/Port) in DHE and DHP header differ");
495 B2DEBUG(29, "DHP ID (Chip/Port) in DHE and DHP header differ $" << hex << dhe_DHPport << " != $" << dhp_dhp_id);
496 }
497 m_errorMask[c_nrDHE_DHP_PORT] = true;
498 }
499
500 if (dhp_header_type != EDHPFrameHeaderDataType::c_ZSD) {
501 if (!(m_suppressErrorMask[c_nrHEADERTYPE_INV])) {
502 B2WARNING("Header type invalid for this kind of DHE frame");
503 B2DEBUG(29, "Header type invalid for this kind of DHE frame: $" << hex << dhp_header_type);
504 }
505 m_errorMask[c_nrHEADERTYPE_INV] = true;
506 return;
507 }
508
509// static int offtab[4] = {0, 64, 128, 192};
510// dhp_offset = offtab[dhp_dhp_id];
511
512 dhp_readout_frame_lo = dhp_pix[3] & 0xFFFF;
513 B2DEBUG(29, "DHP Frame Nr | $" << hex << dhp_readout_frame_lo << " ( " << dec << dhp_readout_frame_lo << " ) ");
514
515 /* // TODO removed because data format error is not to be fixed soon
516 if (((dhp_readout_frame_lo - dhe_first_readout_frame_id_lo) & 0x3F) > m_maxDHPFrameDiff) {
517 if(!m_suppressErrorMask&c_DHP_DHE_FRAME_DIFFER ) B2WARNING("DHP Frame Nr differ from DHE Frame Nr by >1 DHE " << dhe_first_readout_frame_id_lo << " != DHP " << (dhp_readout_frame_lo & 0x3F) << " delta "<< ((dhp_readout_frame_lo - dhe_first_readout_frame_id_lo) & 0x3F) );
518 m_errorMask[c_nrDHP_DHE_FRAME_DIFFER] = true;
519 }
520 */
521 /* // TODO removed because data format error is not to be fixed soon
522 if (m_last_dhp_readout_frame_lo[dhp_dhp_id] != -1) {
523 if (((dhp_readout_frame_lo - m_last_dhp_readout_frame_lo[dhp_dhp_id]) & 0xFFFF) > m_maxDHPFrameDiff) {
524 if(!m_suppressErrorMask&c_DHP_NOT_CONT ) B2WARNING("Two DHP Frames per sensor which frame number differ more than one! " << m_last_dhp_readout_frame_lo[dhp_dhp_id] << ", " <<
525 dhp_readout_frame_lo);
526 m_errorMask[c_nrDHP_NOT_CONT] = true;
527 }
528 }
529 */
530
531 if (daqpktstat.dhc_size() > 0) {
532 if (daqpktstat.dhc_back().dhe_size() > 0) {
533 // only is we have a DHC and DHE object... or back() is undefined
534 // Remark, if we have a broken data (DHE_START/END) structure, we might fill the
535 // previous DHE object ... but then the data is junk anyway
536 daqpktstat.dhc_back().dhe_back().newDHP(dhp_dhp_id, dhp_readout_frame_lo);
537 }
538 }
539
540 /* // TODO removed because the data is not ordered as expected in current firmware
541 for (auto j = 0; j < 4; j++) {
542 if (m_last_dhp_readout_frame_lo[j] != -1) {
543 if (((dhp_readout_frame_lo - m_last_dhp_readout_frame_lo[j]) & 0xFFFF) > m_maxDHPFrameDiff) {
544 if(!m_suppressErrorMask&c_DHP_DHP_FRAME_DIFFER ) B2WARNING("Two DHP Frames (different DHP) per sensor which frame number differ more than one! " << m_last_dhp_readout_frame_lo[j] <<
545 ", " <<
546 dhp_readout_frame_lo);
547 m_errorMask[c_nrDHP_DHP_FRAME_DIFFER] = true;
548 break;// give msg only once
549 }
550 }
551 }
552 */
553 m_last_dhp_readout_frame_lo[dhp_dhp_id] = dhp_readout_frame_lo;
554
555// TODO Please check if this can happen by accident with valid data!
556 if (dhp_pix[2] == dhp_pix[4] && dhp_pix[3] + 1 == dhp_pix[5]) {
557 // We see a second "header" with framenr+1 ...
558 if (!(m_suppressErrorMask[c_nrDHP_DBL_HEADER])) {
559 B2WARNING("DHP data: seems to be double header! skipping.");
560 B2DEBUG(29, "DHP data: seems to be double header! skipping." << LogVar("Length",
561 frame_len));
562 }
563 m_errorMask[c_nrDHP_DBL_HEADER] = true;
564 // dump_dhp(data, frame_len); print out guilty dhp packet
565 return;
566 }
567
568 // Start with offset 4, thus skipping header words
569 for (unsigned int i = 4; i < nr_words ; i++) {
570
571 B2DEBUG(29, "-- $" << hex << dhp_pix[i] << " -- " << dec << i);
572 {
573 if (((dhp_pix[i] & 0x8000) >> 15) == 0) {
574 rowflag = true;
575 if (!pixelflag) {
576 if (!(m_suppressErrorMask[c_nrDHP_ROW_WO_PIX])) B2WARNING("DHP Unpacking: Row w/o Pix");
577 m_errorMask[c_nrDHP_ROW_WO_PIX] = true;
578 }
579 pixelflag = false;
580 dhp_row = (dhp_pix[i] & 0xFFC0) >> 5;
581 dhp_cm = dhp_pix[i] & 0x3F;
582 if (last_gate != -1 && (int)dhp_row / 4 < last_gate) {
583 // B2DEBUF(29,"Wrap " << LogVar("last", last_gate) << LogVar("curr", dhp_row / 4) << LogVar("DHE", dhe_ID) << LogVar("DHP", dhp_dhp_id));
584 wrap++; // relies on the order of data before mapping and the fact that OT firmware delivers only one frame
585 }
586 last_gate = dhp_row / 4;
587
588 if (dhp_cm == 63) { // fifo overflow
589 B2WARNING("DHP data loss (CM=63) in " << LogVar("DHE", dhe_ID) << LogVar("DHP", dhp_dhp_id));
591 m_errorMask[c_nrDHH_MISC_ERROR] = true;
592 }
593 if (daqpktstat.dhc_size() > 0) {
594 if (daqpktstat.dhc_back().dhe_size() > 0) {
595 PXDDAQDHPComMode cm(dhp_dhp_id, dhp_row, dhp_cm);
596 // only is we have a DHC and DHE object... or back() is undefined
597 // Remark, if we have a broken data (DHE_START/END) structure, we might fill the
598 // previous DHE object ... but then the data is junk anyway
599 daqpktstat.dhc_back().dhe_back().addCM(cm);
600 }
601 }
602 B2DEBUG(29, "SetRow: $" << hex << dhp_row << " CM $" << hex << dhp_cm);
603 } else {
604 if (!rowflag) {
605 if (!(m_suppressErrorMask[c_nrDHP_PIX_WO_ROW])) B2WARNING("DHP Unpacking: Pix without Row!!! skip dhp data ");
606 m_errorMask[c_nrDHP_PIX_WO_ROW] = true;
607 // dump_dhp(data, frame_len);// print out faulty dhp frame
608 return;
609 } else {
610 pixelflag = true;
611 dhp_row = (dhp_row & 0xFFE) | ((dhp_pix[i] & 0x4000) >> 14);
612 dhp_col = ((dhp_pix[i] & 0x3F00) >> 8);
613 unsigned int v_cellID, u_cellID;
614 v_cellID = dhp_row;// defaults for no mapping
615 if (dhp_row >= 768) {
616 if (!(m_suppressErrorMask[c_nrROW_OVERFLOW])) B2WARNING("DHP ROW Overflow " << LogVar("Row", dhp_row));
617 m_errorMask[c_nrROW_OVERFLOW] = true;
618 }
619 // we cannot do col overflow check before mapping :-(
620
621 if ((dhe_reformat == 0 && !m_forceNoMapping) || m_forceMapping) {
622 u_cellID = dhp_col;// defaults for no mapping
623 // data has not been pre-processed by DHH, thus we have to do the mapping ourselves
624 if ((dhe_ID & 0x21) == 0x00 || (dhe_ID & 0x21) == 0x21) {
625 // if IFOB
626 PXDMappingLookup::map_rc_to_uv_IF_OB(v_cellID, u_cellID, dhp_dhp_id, dhe_ID);
627 } else { // else OFIB
628 PXDMappingLookup::map_rc_to_uv_IB_OF(v_cellID, u_cellID, dhp_dhp_id, dhe_ID);
629 }
630 } else {
631 u_cellID = dhp_col + 64 * dhp_dhp_id; // defaults for already mapped
632 }
633 if (u_cellID >= 250) {
634 if (!(m_suppressErrorMask[c_nrCOL_OVERFLOW])) {
635 B2WARNING("DHP COL Overflow (unconnected drain lines)");
636 B2DEBUG(29, "DHP COL Overflow (unconnected drain lines) " << u_cellID << ", reformat " << dhe_reformat << ", dhpcol " << dhp_col <<
637 ", id " << dhp_dhp_id);
638 }
639 m_errorMask[c_nrCOL_OVERFLOW] = true;
640 }
641 auto dhp_adc = dhp_pix[i] & 0xFF;
642 B2DEBUG(29, "SetPix: Row $" << hex << dhp_row << " Col $" << hex << dhp_col << " ADC $" << hex << dhp_adc
643 << " CM $" << hex << dhp_cm);
644
645 if (dhp_adc == 0) {
646 // if !suppress error flag
647 B2WARNING("DHE Event truncation in DHE " << dhe_ID << " DHP " << dhp_dhp_id);
648 // m_errorMask |= c_DHE_EVENT_TRUNC;
649 daqpktstat.dhc_back().dhe_back().dhp_back().setTruncated();
650 } else {
651 // in new firmware, (dhp_readout_frame_lo - dhe_first_readout_frame_id_lo) would always been 0
652 if (!m_doNotStore) m_storeRawHits.appendNew(vxd_id, v_cellID, u_cellID, dhp_adc,
653 (dhp_readout_frame_lo - dhe_first_readout_frame_id_lo + wrap) & 0x3F);
654 }
655 }
656 }
657 }
658 }
659
660 B2DEBUG(29, "(DHE) DHE_ID $" << hex << dhe_ID << " (DHE) DHP ID $" << hex << dhe_DHPport << " (DHP) DHE_ID $" << hex << dhp_dhe_id
661 << " (DHP) DHP ID $" << hex << dhp_dhp_id);
662 /*for (int i = 0; i < raw_nr_words ; i++) {
663 B2DEBUG(29, "RAW | " << hex << p_pix[i]);
664 printf("raw %08X | ", p_pix[i]);
665 B2DEBUG(29, "row " << hex << ((p_pix[i] >> 20) & 0xFFF) << dec << " ( " << ((p_pix[i] >> 20) & 0xFFF) << " ) " << " col " << hex << ((p_pix[i] >> 8) & 0xFFF)
666 << " ( " << dec << ((p_pix[i] >> 8) & 0xFFF) << " ) " << " adc " << hex << (p_pix[i] & 0xFF) << " ( " << (p_pix[i] & 0xFF) << " ) "
667 );
668 }*/
669}
670
671void PXDUnpackerOTModule::unpack_dhc_frame(void* data, const int len, const int Frame_Number, const int Frames_in_event,
672 PXDDAQPacketStatus& daqpktstat)
673{
677 static unsigned int eventNrOfOnsenTrgFrame = 0;
678 static int countedBytesInDHC = 0;
679 static bool cancheck_countedBytesInDHC = false;
680 static int countedBytesInDHE = 0;
681 static bool cancheck_countedBytesInDHE = false;
682 static int countedDHEStartFrames = 0;
683 static int countedDHEEndFrames = 0;
684 static int mask_active_dhe = 0;// DHE mask (5 bit)
685 static int nr_active_dhe =
686 0;// TODO just count the active DHEs. Until now, it is not possible to check for the bit mask. we would need the info on which DHE connects to which DHC at which port from gearbox/geometry?
687 static int mask_active_dhp = 0;// DHP active mask, 4 bit, per current DHE
688 static int found_mask_active_dhp = 0;// mask which DHP send data and check on DHE END frame if it matches
689 static int found_good_mask_active_dhp = 0;// mask which DHP send useful data
690 static unsigned int dhe_first_readout_frame_id_lo = 0;
691 // cppcheck-suppress variableScope
692 static unsigned int dhe_first_triggergate = 0;
693 static unsigned int currentDHCID = 0xFFFFFFFF;
694 static unsigned int currentDHEID = 0xFFFFFFFF;
695 static unsigned int currentVxdId = 0;
696 static bool isFakedData_event = false;
697 static bool isUnfiltered_event = false;
698
699
700 if (Frame_Number == 0) {
701 // We reset the counters on the first event
702 // we do this before any other check is done
703 eventNrOfOnsenTrgFrame = 0;
704 countedDHEStartFrames = 0;
705 countedDHEEndFrames = 0;
706 countedBytesInDHC = 0;
707 cancheck_countedBytesInDHC = false;
708 countedBytesInDHE = 0;
709 cancheck_countedBytesInDHE = false;
710 currentDHCID = 0xFFFFFFFF;
711 currentDHEID = 0xFFFFFFFF;
712 currentVxdId = 0;
713 isUnfiltered_event = false;
714 isFakedData_event = false;
715 mask_active_dhe = 0;
716 nr_active_dhe = 0;
717 mask_active_dhp = 0;
718 found_mask_active_dhp = 0;
719 found_good_mask_active_dhp = 0;
720 }
721
723
724 dhc_frames dhc;
725 dhc.set(data, hw->getFrameType(), len);
726
727 {
728 // if a fixed size frame has a different length, how can we rely on its content???
729 // AND we could by typecasting access memory beyond end of data (but very unlikely)
730 // for that reason this we have to check before any CRC and stop unpacking the frame
731 int s = dhc.getFixedSize();
732 if (len != s && s != 0) {
733 if (!(m_suppressErrorMask[c_nrFIX_SIZE])) {
734 B2WARNING("Fixed frame type size does not match specs" << LogVar("expected length",
735 len) << LogVar("length in data", s));
736 }
737 m_errorMask[c_nrFIX_SIZE] = true;
738 if (!m_continueOnError) return;
739 }
740 }
741
742 // What do we do with wrong checksum frames? As we do not know WHAT is wrong, we have to skip them altogether.
743 // As they might contain HEADER Info, we might better skip the processing of the full package, too.
744 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
745 if (!m_continueOnError && m_errorMask[c_nrDHE_CRC]) {
746 // if CRC is wrong, we cannot rely on the content of the frame, thus skipping is the best option
747 return;
748 }
749
750 unsigned int eventNrOfThisFrame = dhc.getEventNrLo();
751 int frame_type = dhc.getFrameType();
752
753 if (Frame_Number == 0) {
754 if (m_formatBonnDAQ) {
755 if (frame_type != EDHCFrameHeaderDataType::c_DHC_START) {
756 if (!(m_suppressErrorMask[c_nrEVENT_STRUCT])) B2WARNING("This looks not like BonnDAQ format.");
757 m_errorMask[c_nrEVENT_STRUCT] = true;
758// if (!m_continueOnError) return; // requires more testing
759 }
760 } else {
761 if (frame_type == EDHCFrameHeaderDataType::c_DHC_START) {
762 if (!(m_suppressErrorMask[c_nrEVENT_STRUCT]))
763 B2WARNING("This looks like BonnDAQ or old Desy 2013/14 testbeam format. Please use formatBonnDAQ or the pxdUnpackerDesy1314 module.");
764 m_errorMask[c_nrEVENT_STRUCT] = true;
765// if (!m_continueOnError) return; // requires more testing
766 }
767 }
768 }
769
770 if (!m_formatBonnDAQ) {
771 if (Frame_Number == 1) {
772 if (frame_type == EDHCFrameHeaderDataType::c_DHC_START) {
773 isFakedData_event = dhc.data_dhc_start_frame->isFakedData();
774 }
775 }
776
777 // please check if this mask is suitable. At least here we are limited by the 16 bit trigger number in the DHH packet header.
778 // we can use more bits in the DHC and DHE START Frame
779 if ((eventNrOfThisFrame & 0xFFFF) != (m_meta_event_nr & 0xFFFF)) {
780 if (!isFakedData_event) {
781 if (!(m_suppressErrorMask[c_nrMETA_MM])) {
782 B2WARNING("Event Numbers do not match for this frame");
783 B2DEBUG(29, "Event Numbers do not match for this frame" <<
784 LogVar("Event nr in frame $", static_cast < std::ostringstream
785 && >(std::ostringstream() << hex << eventNrOfThisFrame).str()) <<
786 LogVar("Event nr in MetaInfo (bits masked) $",
787 static_cast < std::ostringstream && >(std::ostringstream() << hex << m_meta_event_nr).str()));
788 }
789 m_errorMask[c_nrMETA_MM] = true;
790// if (!m_continueOnError) return; // requires more testing
791 }
792 }
793
794 if (Frame_Number > 1 && Frame_Number < Frames_in_event - 1) {
795 if (countedDHEStartFrames != countedDHEEndFrames + 1)
796 if (frame_type != EDHCFrameHeaderDataType::c_ONSEN_ROI && frame_type != EDHCFrameHeaderDataType::c_DHE_START) {
797 if (!(m_suppressErrorMask[c_nrDATA_OUTSIDE])) B2WARNING("Data Frame outside a DHE START/END");
798 m_errorMask[c_nrDATA_OUTSIDE] = true;
799// if (!m_continueOnError) return; // requires more testing
800 }
801 }
802 }
803
804 // TODO How do we handle Frames where Error Bit is set in header?
805 // Currently there is no documentation what it actually means... only an error bit is set (below)
806 // the following errors must be "accepted", as all firmware sets it wrong from Ghost frames.
807 if (hw->getErrorFlag()) {
808 if (frame_type != EDHCFrameHeaderDataType::c_GHOST) {
809 // We get ERROR bits in header even if only one module or DHP link is missing... thus
810 // we better filter a bit more ... but how?
811 if (!(m_suppressErrorMask[c_nrHEADER_ERR])) B2WARNING("Error Bit set in DHE Header");
812 m_errorMask[c_nrHEADER_ERR] = true;// TODO this should have some effect ... when does it mean something? documentation missing
813 }
814 } else {
815 if (frame_type == EDHCFrameHeaderDataType::c_GHOST) {
816 m_errorMask[c_nrHEADER_ERR_GHOST] = true;
817 }
818 }
819
820 switch (frame_type) {
821 case EDHCFrameHeaderDataType::c_DHP_RAW: {
822
824 if (currentDHEID != dhc.data_direct_readout_frame_raw->getDHEId()) {
825 if (!(m_suppressErrorMask[c_nrDHE_START_ID])) {
826 B2WARNING("DHE ID from DHE Start and this frame do not match");
827 B2DEBUG(29, "DHE ID from DHE Start and this frame do not match" <<
828 LogVar("DHEID in this frame $", static_cast < std::ostringstream
829 && >(std::ostringstream() << hex << dhc.data_direct_readout_frame_raw->getDHEId()).str()) <<
830 LogVar("DHEID expected $", static_cast < std::ostringstream && >(std::ostringstream() << hex << currentDHEID).str()));
831 }
832 m_errorMask[c_nrDHE_START_ID] = true;
833 }
834 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
835 if ((found_mask_active_dhp & (1 << dhc.data_direct_readout_frame->getDHPPort())) != 0) {
836 B2ERROR("Second DHP data packet (MEMDUMP) for " << LogVar("DHE", currentDHEID) << LogVar("DHP",
838 }
839
840 found_mask_active_dhp |= 1 << dhc.data_direct_readout_frame->getDHPPort();
841
842 unpack_dhp_raw(data, len - 4,
845 currentVxdId);
846
847 break;
848 };
849 case EDHCFrameHeaderDataType::c_ONSEN_DHP:
850 // Set the counted size invalid if negative, needs a large negative value because we are adding up to that
851 cancheck_countedBytesInDHC = false;
852 cancheck_countedBytesInDHE = false;
853 [[fallthrough]];
854 case EDHCFrameHeaderDataType::c_DHP_ZSD: {
855
857 if (isUnfiltered_event) {
858 if (frame_type == EDHCFrameHeaderDataType::c_ONSEN_DHP) m_errorMask[c_nrSENDALL_TYPE] = true;
859 } else {
860 if (frame_type == EDHCFrameHeaderDataType::c_DHP_ZSD) m_errorMask[c_nrNOTSENDALL_TYPE] = true;
861 }
862
863 //m_errorMask |= dhc.data_direct_readout_frame->check_error();
864
865 if (currentDHEID != dhc.data_direct_readout_frame_raw->getDHEId()) {
866 if (!(m_suppressErrorMask[c_nrDHE_START_ID])) {
867 B2WARNING("DHE ID from DHE Start and this frame do not match");
868 B2DEBUG(29, "DHE ID from DHE Start and this frame do not match" <<
869 LogVar("DHEID in this frame $", static_cast < std::ostringstream
870 && >(std::ostringstream() << hex << dhc.data_direct_readout_frame_raw->getDHEId()).str()) <<
871 LogVar("DHEID expected $", static_cast < std::ostringstream && >(std::ostringstream() << hex << currentDHEID).str()));
872 }
873 m_errorMask[c_nrDHE_START_ID] = true;
874 }
875 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
876 if ((found_mask_active_dhp & (1 << dhc.data_direct_readout_frame->getDHPPort())) != 0) {
877 B2ERROR("Second DHP data packet for " << LogVar("DHE", currentDHEID) << LogVar("DHP", dhc.data_direct_readout_frame->getDHPPort()));
878 }
879 found_mask_active_dhp |= 1 << dhc.data_direct_readout_frame->getDHPPort();
880 found_good_mask_active_dhp |= 1 << dhc.data_direct_readout_frame->getDHPPort();// only this frametype has useful data
881 if (m_checkPaddingCRC) dhc.check_padding(m_errorMask); // isUnfiltered_event
882
883
884 unpack_dhp(data, len - 4,
885 dhe_first_readout_frame_id_lo,
889 currentVxdId, daqpktstat);
890
891 break;
892 };
893 case EDHCFrameHeaderDataType::c_ONSEN_FCE:
894 // Set the counted size invalid if negative, needs a large negative value because we are adding up to that
895 cancheck_countedBytesInDHC = false;
896 cancheck_countedBytesInDHE = false;
897 [[fallthrough]];
898 case EDHCFrameHeaderDataType::c_FCE_RAW: {
899 if (!(m_suppressErrorMask[c_nrUNEXPECTED_FRAME_TYPE])) B2WARNING("Unexpected Frame Type (Clustering FCE)");
900 m_errorMask[c_nrUNEXPECTED_FRAME_TYPE] = true;
901 if (m_verbose) hw->print();
902 if (isUnfiltered_event) {
903 if (frame_type == EDHCFrameHeaderDataType::c_ONSEN_FCE) {
904 // TODO add error message
905 m_errorMask[c_nrSENDALL_TYPE] = true;
906 }
907 } else {
908 if (frame_type == EDHCFrameHeaderDataType::c_FCE_RAW) {
909 // TODO add error message
910 m_errorMask[c_nrNOTSENDALL_TYPE] = true;
911 }
912 }
913
914 if (currentDHEID != dhc.data_direct_readout_frame_raw->getDHEId()) {
915 if (!(m_suppressErrorMask[c_nrDHE_START_ID])) {
916 B2WARNING("DHE ID from DHE Start and this frame do not match");
917 B2DEBUG(29, "DHE ID from DHE Start and this frame do not match" <<
918 LogVar("DHEID in this frame $", static_cast < std::ostringstream
919 && >(std::ostringstream() << hex << dhc.data_direct_readout_frame_raw->getDHEId()).str()) <<
920 LogVar("DHEID expected $", static_cast < std::ostringstream && >(std::ostringstream() << hex << currentDHEID).str()));
921 }
922 m_errorMask[c_nrDHE_START_ID] = true;
923 }
924 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
925 if ((found_mask_active_dhp & (1 << dhc.data_direct_readout_frame->getDHPPort())) != 0) {
926 B2ERROR("Second DHP data packet (FCE) for " << LogVar("DHE", currentDHEID) << LogVar("DHP",
928 }
929 found_mask_active_dhp |= 1 << dhc.data_direct_readout_frame->getDHPPort();
930
931 B2DEBUG(29, "UNPACK FCE FRAME with len $" << hex << len);
932 unpack_fce((unsigned short*) data, len - 4, currentVxdId);
933
934 break;
935 };
936 case EDHCFrameHeaderDataType::c_COMMODE: {
937 // this frame type has up to now not been well defined, we do not expect it until
938 // the firmware supports clustering in hardware
939 if (!(m_suppressErrorMask[c_nrUNEXPECTED_FRAME_TYPE])) B2WARNING("Unexpected Frame Type (COMMODE)");
940 m_errorMask[c_nrUNEXPECTED_FRAME_TYPE] = true;
941
942 if (m_verbose) hw->print();
943 if (currentDHEID != dhc.data_commode_frame->getDHEId()) {
944 if (!(m_suppressErrorMask[c_nrDHE_START_ID])) {
945 B2WARNING("DHE ID from DHE Start and this frame do not match");
946 B2DEBUG(29, "DHE ID from DHE Start and this frame do not match" <<
947 LogVar("DHEID in this frame $", static_cast < std::ostringstream
948 && >(std::ostringstream() << hex << dhc.data_commode_frame->getDHEId()).str()) <<
949 LogVar("DHEID expected $", static_cast < std::ostringstream && >(std::ostringstream() << hex << currentDHEID).str()));
950 }
951 m_errorMask[c_nrDHE_START_ID] = true;
952 }
953 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
954 break;
955 };
956 case EDHCFrameHeaderDataType::c_DHC_START: {
957 countedBytesInDHC = 0;
958 cancheck_countedBytesInDHC = true;
959 if (isFakedData_event != dhc.data_dhc_start_frame->isFakedData()) {
960 if (!(m_suppressErrorMask[c_nrFAKE_NO_FAKE_DATA])) B2WARNING("DHC START mixed Fake/no Fake event.");
961 m_errorMask[c_nrFAKE_NO_FAKE_DATA] = true;
962 }
964 if (!(m_suppressErrorMask[c_nrFAKE_NO_DATA_TRIG])) B2WARNING("Faked DHC START Data -> trigger without Data!");
965 m_errorMask[c_nrFAKE_NO_DATA_TRIG] = true;
966 } else {
968 }
969
970// eventNrOfOnsenTrgFrame = eventNrOfThisFrame;
971 currentDHEID = 0xFFFFFFFF;
972 currentVxdId = 0;
973 currentDHCID = dhc.data_dhc_start_frame->get_dhc_id();
974 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
975
976 if (m_formatBonnDAQ) eventNrOfOnsenTrgFrame = eventNrOfThisFrame;
977
978 if (!isFakedData_event) {
982 if (!(m_suppressErrorMask[c_nrMETA_MM_DHC_ERS])) {
983 B2WARNING("DHC-Meta Experiment number mismatch");
984 B2DEBUG(29, "DHC-Meta Experiment number mismatch" <<
985 LogVar("DHC exp nr",
987 LogVar("META exp nr", m_meta_experiment));
988 }
989 m_errorMask[c_nrMETA_MM_DHC_ERS] = true;
990 }
992 if (!(m_suppressErrorMask[c_nrMETA_MM_DHC_ERS])) {
993 B2WARNING("DHC-Meta Run number mismatch");
994 B2DEBUG(29, "DHC-Meta Run number mismatch" <<
995 LogVar("DHC Run nr",
997 LogVar("META run nr", m_meta_run_nr));
998 }
999 m_errorMask[c_nrMETA_MM_DHC_ERS] = true;
1000 }
1002 if (!(m_suppressErrorMask[c_nrMETA_MM_DHC_ERS])) {
1003 B2WARNING("DHC-Meta Sub-Run number mismatch");
1004 B2DEBUG(29, "DHC-Meta Sub-Run number mismatch" <<
1005 LogVar("DHC subrun nr",
1007 LogVar("META subrun nr", m_meta_subrun_nr));
1008 }
1009 m_errorMask[c_nrMETA_MM_DHC_ERS] = true;
1010 }
1011 if ((((unsigned int)dhc.data_dhc_start_frame->getEventNrHi() << 16) | dhc.data_dhc_start_frame->getEventNrLo()) !=
1012 (m_meta_event_nr & 0xFFFFFFFF)) {
1013 if (!(m_suppressErrorMask[c_nrMETA_MM_DHC])) {
1014 B2WARNING("DHC-Meta 32 bit event number mismatch");
1015 B2DEBUG(29, "DHC-Meta 32 bit event number mismatch" <<
1016 LogVar("DHC trigger nr", (((unsigned int) dhc.data_dhc_start_frame->getEventNrHi() << 16) |
1018 LogVar("META trigger nr", (unsigned int)(m_meta_event_nr & 0xFFFFFFFF)));
1019 }
1020 m_errorMask[c_nrMETA_MM_DHC] = true;
1021 }
1022 uint32_t trig_ticks = (((unsigned int)dhc.data_dhc_start_frame->time_tag_mid & 0x7FFF) << 12) | ((unsigned int)
1024 uint32_t trig_sec = (dhc.data_dhc_start_frame->time_tag_hi * 2) ;
1025 if (dhc.data_dhc_start_frame->time_tag_mid & 0x8000) trig_sec++;
1026
1027 if ((trig_ticks - m_meta_ticks) != 0 || (trig_sec - m_meta_sec) != 0) {
1028 m_errorMask[c_nrMETA_MM_DHC_TT] = true;
1029 if (!(m_suppressErrorMask[c_nrMETA_MM_DHC_TT])) {
1030 B2WARNING("DHC-Meta TimeTag mismatch");
1031 B2DEBUG(29, "DHC-Meta TimeTag mismatch" <<
1032 LogVar("Header Time $", static_cast < std::ostringstream && >(std::ostringstream() <<
1033 hex << dhc.data_dhc_start_frame->time_tag_hi << "." <<
1034 dhc.data_dhc_start_frame->time_tag_mid << "." <<
1036 LogVar("Meta Time $", static_cast < std::ostringstream && >(std::ostringstream() << hex << m_meta_time).str()) <<
1037 LogVar("Trigger Type", static_cast < std::ostringstream
1038 && >(std::ostringstream() << hex << (dhc.data_dhc_start_frame->time_tag_lo_and_type & 0xF)).str()) <<
1039 LogVar("Meta seconds: $", static_cast < std::ostringstream && >(std::ostringstream() << hex << m_meta_sec).str()) <<
1040 LogVar("DHC seconds $", static_cast < std::ostringstream && >(std::ostringstream() << hex << trig_sec).str()) <<
1041 LogVar("Seconds difference $", static_cast < std::ostringstream
1042 && >(std::ostringstream() << hex << (trig_sec - m_meta_sec)).str()) <<
1043 LogVar("Meta ticks from 127MHz $", static_cast < std::ostringstream && >(std::ostringstream() << hex << m_meta_ticks).str()) <<
1044 LogVar("DHC ticks from 127MHz $", static_cast < std::ostringstream && >(std::ostringstream() << hex << trig_ticks).str()) <<
1045 LogVar("Tick difference $", static_cast < std::ostringstream
1046 && >(std::ostringstream() << hex << (trig_ticks - m_meta_ticks)).str()));
1047 }
1048 } else {
1049 B2DEBUG(29, "DHC TT: $" << hex << dhc.data_dhc_start_frame->time_tag_hi << "." << dhc.data_dhc_start_frame->time_tag_mid << "." <<
1050 dhc.data_dhc_start_frame->time_tag_lo_and_type << " META " << m_meta_time << " TRG Type " <<
1052 }
1053 }
1054 mask_active_dhe = dhc.data_dhc_start_frame->get_active_dhe_mask();
1055 nr_active_dhe = nr5bits(mask_active_dhe);
1056
1057 m_errorMaskDHC = m_errorMask; // forget about anything before this frame
1058 daqpktstat.newDHC(currentDHCID, m_errorMask);
1061
1062 break;
1063 };
1064 case EDHCFrameHeaderDataType::c_DHE_START: {
1065 countedBytesInDHE = 0;
1066 cancheck_countedBytesInDHE = true;
1072 dhe_first_readout_frame_id_lo = dhc.data_dhe_start_frame->getStartFrameNr();
1073 dhe_first_triggergate = dhc.data_dhe_start_frame->getTriggerGate();
1074 if (currentDHEID != 0xFFFFFFFF && (currentDHEID & 0xFFFF) >= dhc.data_dhe_start_frame->getDHEId()) {
1075 if (!(m_suppressErrorMask[c_nrDHE_WRONG_ID_SEQ])) {
1076 B2WARNING("DHH IDs are not in expected order");
1077 B2DEBUG(29, "DHH IDs are not in expected order" <<
1078 LogVar("Previous ID", (currentDHEID & 0xFFFF)) <<
1079 LogVar("Current ID", dhc.data_dhe_start_frame->getDHEId()));
1080 }
1081 m_errorMask[c_nrDHE_WRONG_ID_SEQ] = true;
1082 }
1083 currentDHEID = dhc.data_dhe_start_frame->getDHEId();
1084 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
1085
1086 if (countedDHEStartFrames > countedDHEEndFrames) {
1087 if (!(m_suppressErrorMask[c_nrDHE_START_WO_END])) B2WARNING("DHE_START without DHE_END");
1088 m_errorMask[c_nrDHE_START_WO_END] = true;
1089 }
1090 countedDHEStartFrames++;
1091
1092 found_mask_active_dhp = 0;
1093 found_good_mask_active_dhp = 0;
1094 mask_active_dhp = dhc.data_dhe_start_frame->getActiveDHPMask();
1095
1096 if ((((unsigned int)dhc.data_dhe_start_frame->getEventNrHi() << 16) | dhc.data_dhe_start_frame->getEventNrLo()) != (unsigned int)(
1097 m_meta_event_nr & 0xFFFFFFFF)) {
1098 if (!(m_suppressErrorMask[c_nrMETA_MM_DHE])) {
1099 B2WARNING("DHE START trigger mismatch in EVT32b/HI WORD");
1100 B2DEBUG(29, "DHE START trigger mismatch in EVT32b/HI WORD" <<
1101 LogVar("DHE Start trigger nr", (dhc.data_dhe_start_frame->getEventNrHi() << 16) | dhc.data_dhe_start_frame->getEventNrLo()) <<
1102 LogVar("Meta trigger nr", (m_meta_event_nr & 0xFFFFFFFF)));
1103 }
1104 m_errorMask[c_nrMETA_MM_DHE] = true;
1105 }
1106// B2WARNING("DHE TT: $" << hex << dhc.data_dhe_start_frame->dhe_time_tag_hi << "." << dhc.data_dhe_start_frame->dhe_time_tag_lo <<
1107// " META " << m_meta_time);
1108
1109 if (currentDHEID == 0) {
1110 if (!(m_suppressErrorMask[c_nrDHE_ID_INVALID])) B2WARNING("DHE ID is invalid=0 (not initialized)");
1111 m_errorMask[c_nrDHE_ID_INVALID] = true;
1112 }
1113 // calculate the VXDID for DHE and save them for DHP unpacking
1114 {
1121 unsigned short sensor, ladder, layer;
1122 sensor = (currentDHEID & 0x1) + 1;
1123 ladder = (currentDHEID & 0x1E) >> 1; // no +1
1124 layer = ((currentDHEID & 0x20) >> 5) + 1;
1125 currentVxdId = VxdID(layer, ladder, sensor);
1126 if (ladder == 0 || (layer == 1 && ladder > 8) || (layer == 2 && ladder > 12)) {
1127 if (!(m_suppressErrorMask[c_nrDHE_ID_INVALID])) {
1128 B2WARNING("DHE ID is invalid");
1129 B2DEBUG(29, "DHE ID is invalid" <<
1130 LogVar("DHE ID", currentDHEID) <<
1131 LogVar("Layer", layer) <<
1132 LogVar("Ladder", ladder) <<
1133 LogVar("Sensor", sensor));
1134 }
1135 m_errorMask[c_nrDHE_ID_INVALID] = true;
1136 }
1137 }
1138
1139 m_errorMaskDHE = m_errorMask; // forget about anything before this frame
1140 if (daqpktstat.dhc_size() > 0) {
1141 // if no DHC has been defined yet, do nothing!
1142 daqpktstat.dhc_back().newDHE(currentVxdId, currentDHEID, m_errorMask, dhe_first_triggergate, dhe_first_readout_frame_id_lo);
1143 }
1144 break;
1145 };
1146 case EDHCFrameHeaderDataType::c_GHOST:
1147 if (m_verbose) dhc.data_ghost_frame->print();
1148 if (currentDHEID != dhc.data_ghost_frame->getDHEId()) {
1149 if (!(m_suppressErrorMask[c_nrDHE_START_ID])) {
1150 B2WARNING("DHE ID from DHE Start and this frame do not match");
1151 B2DEBUG(29, "Start ID $" << hex << currentDHEID << " != $" << dhc.data_ghost_frame->getDHEId());
1152 }
1153 m_errorMask[c_nrDHE_START_ID] = true;
1154 }
1156 if ((found_mask_active_dhp & (1 << dhc.data_ghost_frame->getDHPPort())) != 0) {
1157 B2ERROR("Second DHP data packet (GHOST) for " << LogVar("DHE", currentDHEID) << LogVar("DHP", dhc.data_ghost_frame->getDHPPort()));
1158 }
1159 found_mask_active_dhp |= 1 << dhc.data_ghost_frame->getDHPPort();
1160
1161 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
1162
1163 break;
1164 case EDHCFrameHeaderDataType::c_DHC_END: {
1165 if (dhc.data_dhc_end_frame->isFakedData() != isFakedData_event) {
1166 if (!(m_suppressErrorMask[c_nrFAKE_NO_FAKE_DATA])) B2WARNING("DHC END mixed Fake/no Fake event.");
1167 m_errorMask[c_nrFAKE_NO_FAKE_DATA] = true;
1168 }
1169 if (dhc.data_dhc_end_frame->isFakedData()) {
1170 if (!(m_suppressErrorMask[c_nrFAKE_NO_DATA_TRIG])) B2WARNING("Faked DHC END Data -> trigger without Data!");
1171 m_errorMask[c_nrFAKE_NO_DATA_TRIG] = true;
1172 } else {
1174 }
1175
1176 if (!isFakedData_event) {
1177 if (dhc.data_dhc_end_frame->get_dhc_id() != currentDHCID) {
1178 if (!(m_suppressErrorMask[c_nrDHC_DHCID_START_END_MM])) {
1179 B2WARNING("DHC ID Mismatch between Start and End");
1180 B2DEBUG(29, "DHC ID Mismatch between Start and End $" << std::hex <<
1181 currentDHCID << "!=$" << dhc.data_dhc_end_frame->get_dhc_id());
1182 }
1183 m_errorMask[c_nrDHC_DHCID_START_END_MM] = true;
1184 }
1185 int w;
1186 w = dhc.data_dhc_end_frame->get_words() * 4;
1187 if (cancheck_countedBytesInDHC) {
1188 if (countedBytesInDHC != w) {
1189 if (!(m_suppressErrorMask[c_nrDHC_WIE])) {
1190 B2WARNING("Number of Words in DHC END does not match");
1191 B2DEBUG(29, "Number of Words in DHC END does not match: WIE $" << hex << countedBytesInDHC << " != DHC END $" << hex << w);
1192 }
1193 m_errorMask[c_nrDHC_WIE] = true;
1194 } else {
1195 if (m_verbose)
1196 B2DEBUG(29, "EVT END: WIE $" << hex << countedBytesInDHC << " == DHC END $" << hex << w);
1197 }
1198 // else ... processed data -> length invalid
1199 }
1200 }
1202 if (dhc.data_dhc_end_frame->getErrorInfo() != 0) {
1203 if (!(m_suppressErrorMask[c_nrDHH_END_ERRORBITS])) B2ERROR("DHC END Error Info set to $" << hex <<
1205 m_errorMask[c_nrDHH_END_ERRORBITS] = true;
1206 }
1207 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
1208 m_errorMaskDHC |= m_errorMask; // do latest updates
1209
1210 if (daqpktstat.dhc_size() > 0) {
1211 // only is we have a DHC object... or back() is undefined
1212 // Remark, if we have a broken data (DHC_START/END) structure, we might fill the
1213 // previous DHC object ... but then the data is junk anyway
1214 daqpktstat.dhc_back().setErrorMask(m_errorMaskDHC);
1215 //B2DEBUG(98,"** DHC "<<currentDHCID<<" Raw"<<dhc.data_dhc_end_frame->get_words() * 4 <<" Red"<<countedBytesInDHC);
1216 daqpktstat.dhc_back().setCounters(dhc.data_dhc_end_frame->get_words() * 4, countedBytesInDHC);
1218 }
1219 m_errorMaskDHC = 0;
1220 currentDHEID = 0xFFFFFFFF;
1221 currentDHCID = 0xFFFFFFFF;
1222 currentVxdId = 0;
1223 break;
1224 };
1225 case EDHCFrameHeaderDataType::c_DHE_END: {
1227 if (currentDHEID != dhc.data_dhe_end_frame->getDHEId()) {
1228 if (!(m_suppressErrorMask[c_nrDHE_START_END_ID])) {
1229 B2WARNING("DHE ID from DHE Start and this frame do not match");
1230 B2DEBUG(29, "DHE ID from DHE Start and this frame do not match $" << hex << currentDHEID << " != $" <<
1232 }
1233 m_errorMask[c_nrDHE_START_END_ID] = true;
1234 }
1236 if (dhc.data_dhe_end_frame->getErrorInfo() != 0) {
1237 if (!(m_suppressErrorMask[c_nrDHH_END_ERRORBITS])) {
1238 B2ERROR("DHE END Error Info set to $" << hex << dhc.data_dhe_end_frame->getErrorInfo());
1239 }
1240 m_errorMask[c_nrDHH_END_ERRORBITS] = true;
1241 }
1242 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
1243 if (found_mask_active_dhp != mask_active_dhp) {
1244 if (!(m_suppressErrorMask[c_nrDHP_ACTIVE])) {
1245 B2WARNING("DHE_END: DHP active mask differs from found data");
1246 B2DEBUG(29, "DHE_END: DHP active mask differs from found data $" << hex << mask_active_dhp << " != $" << hex <<
1247 found_mask_active_dhp
1248 << " mask of found dhp/ghost frames");
1249 }
1250 m_errorMask[c_nrDHP_ACTIVE] = true;
1251 }
1252 countedDHEEndFrames++;
1253 if (countedDHEStartFrames < countedDHEEndFrames) {
1254 // the other case is checked in Start
1255 if (!(m_suppressErrorMask[c_nrDHE_END_WO_START])) B2WARNING("DHE_END without DHE_START");
1256 m_errorMask[c_nrDHE_END_WO_START] = true;
1257 }
1258 {
1259 int w;
1260 w = dhc.data_dhe_end_frame->get_words() * 2;
1261 if (cancheck_countedBytesInDHE) {
1262 if (countedBytesInDHE != w) {
1263 if (!(m_suppressErrorMask[c_nrDHE_WIE])) {
1264 B2WARNING("Number of Words in DHE END does not match");
1265 B2DEBUG(29, "Number of Words in DHE END does not match: WIE $" << hex << countedBytesInDHE << " != DHE END $" << hex << w);
1266 }
1267 m_errorMask[c_nrDHE_WIE] = true;
1268 } else {
1269 if (m_verbose)
1270 B2DEBUG(29, "EVT END: WIE $" << hex << countedBytesInDHE << " == DHE END $" << hex << w);
1271 }
1272 // else ... processed data -> length invalid
1273 }
1274 }
1275 m_errorMaskDHE |= m_errorMask; // do latest updates
1276
1277 if (daqpktstat.dhc_size() > 0) {
1278 if (daqpktstat.dhc_back().dhe_size() > 0) {
1279 // only is we have a DHC and DHE object... or back() is undefined
1280 // Remark, if we have a broken data (DHE_START/END) structure, we might fill the
1281 // previous DHE object ... but then the data is junk anyway
1283 // B2DEBUG(98,"** DHC "<<currentDHEID<<" Raw "<<dhc.data_dhe_end_frame->get_words() * 2 <<" Red"<<countedBytesInDHE);
1284 daqpktstat.dhc_back().dhe_back().setCounters(dhc.data_dhe_end_frame->get_words() * 2, countedBytesInDHE);
1285 daqpktstat.dhc_back().dhe_back().setDHPFoundMask(found_good_mask_active_dhp);
1287 }
1288 }
1289 m_errorMaskDHE = 0;
1290 currentDHEID |= 0xFF000000;// differentiate from 0xFFFFFFFFF as initial value
1291 currentVxdId = 0;
1292 break;
1293 };
1294 case EDHCFrameHeaderDataType::c_ONSEN_ROI:
1296 dhc.data_onsen_roi_frame->check_error(m_errorMask, len, m_suppressErrorMask[c_nrROI_PACKET_INV_SIZE]);
1298 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
1299 if (!m_doNotStore) {
1300 //dhc.data_onsen_roi_frame->save(m_storeROIs, len, (unsigned int*) data);
1301 // void save(StoreArray<PXDRawROIs>& sa, unsigned int length, unsigned int* data) const
1302 // 4 byte header, ROIS (n*8), 4 byte copy of inner CRC, 4 byte outer CRC
1303 if (len >= dhc.data_onsen_roi_frame->getMinSize()) {
1304 //if ((len - dhc.data_onsen_roi_frame->getMinSize()) % 8 != 0) {
1305 // error checking in check_error() above, this is only for dump-ing
1306 // dump_roi(data, len - 4); // dump ROI payload, minus CRC
1307 //}
1308 unsigned int l;
1309 l = (len - dhc.data_onsen_roi_frame->getMinSize()) / 8;
1310 // Endian swapping is done in Constructor of RawRoi object
1311 m_storeROIs.appendNew(l, &((unsigned int*) data)[1]);
1312 }
1313 }
1314 break;
1315 case EDHCFrameHeaderDataType::c_ONSEN_TRG:
1316 eventNrOfOnsenTrgFrame = eventNrOfThisFrame;
1317 if (dhc.data_onsen_trigger_frame->get_trig_nr1() != (unsigned int)(m_meta_event_nr & 0xFFFFFFFF)) {
1318 if (!(m_suppressErrorMask[c_nrMETA_MM_ONS_HLT])) {
1319 B2WARNING("Trigger Frame HLT Trigger Nr mismatch");
1320 B2DEBUG(29, "Trigger Frame HLT Trigger Nr mismatch: HLT $" <<
1321 dhc.data_onsen_trigger_frame->get_trig_nr1() << " META " << (m_meta_event_nr & 0xFFFFFFFF));
1322 }
1323 m_errorMask[c_nrMETA_MM_ONS_HLT] = true;
1324 }
1328 if (!(m_suppressErrorMask[c_nrMETA_MM_ONS_HLT])) {
1329 B2WARNING("Trigger Frame HLT Exp/Run/Subrun Nr mismatch");
1330 B2DEBUG(29, "Trigger Frame HLT Exp/Run/Subrun Nr mismatch: Exp HLT $" <<
1332 " Run HLT $" << dhc.data_onsen_trigger_frame->get_run1() << " META " << m_meta_run_nr <<
1333 " Subrun HLT $" << dhc.data_onsen_trigger_frame->get_subrun1() << " META " << m_meta_subrun_nr);
1334 }
1335 m_errorMask[c_nrMETA_MM_ONS_HLT] = true;
1336 }
1337
1339 if (dhc.data_onsen_trigger_frame->get_trig_nr2() != (unsigned int)(m_meta_event_nr & 0xFFFFFFFF)) {
1340 if (!(m_suppressErrorMask[c_nrMETA_MM_ONS_DC])) {
1341 B2WARNING("Trigger Frame DATCON Trigger Nr mismatch");
1342 B2DEBUG(29, "Trigger Frame DATCON Trigger Nr mismatch: DC $" <<
1343 dhc.data_onsen_trigger_frame->get_trig_nr2() << " META " << (m_meta_event_nr & 0xFFFFFFFF));
1344 }
1345 m_errorMask[c_nrMETA_MM_ONS_DC] = true;
1346 }
1350 if (!(m_suppressErrorMask[c_nrMETA_MM_ONS_DC])) {
1351 B2WARNING("Trigger Frame DATCON Exp/Run/Subrun Nr mismatch");
1352 B2DEBUG(29, "Trigger Frame DATCON Exp/Run/Subrun Nr mismatch: Exp DC $" <<
1354 " Run DC $" << dhc.data_onsen_trigger_frame->get_run2() << " META " << m_meta_run_nr <<
1355 " Subrun DC $" << dhc.data_onsen_trigger_frame->get_subrun2() << " META " << m_meta_subrun_nr);
1356 }
1357 m_errorMask[c_nrMETA_MM_ONS_DC] = true;
1358 }
1359 }
1360
1361// B2WARNING("TRG TAG HLT: $" << hex << dhc.data_onsen_trigger_frame->get_trig_tag1() << " DATCON $" << dhc.data_onsen_trigger_frame->get_trig_tag2() << " META " << m_meta_time);
1362
1365 m_suppressErrorMask[c_nrMERGER_TRIGNR]);
1366 dhc.check_crc(m_errorMask, m_suppressErrorMask[c_nrDHE_CRC]);
1367 if (Frame_Number != 0) {
1368 if (!(m_suppressErrorMask[c_nrEVENT_STRUCT])) B2WARNING("ONSEN TRG Frame must be the first one.");
1369 m_errorMask[c_nrEVENT_STRUCT] = true;
1370 }
1371 isUnfiltered_event = dhc.data_onsen_trigger_frame->is_SendUnfiltered();
1372 if (isUnfiltered_event) m_sendunfiltered++;
1375 break;
1376 default:
1377 if (!(m_suppressErrorMask[c_nrDHC_UNKNOWN])) B2WARNING("UNKNOWN DHC frame type");
1378 m_errorMask[c_nrDHC_UNKNOWN] = true;
1379 if (m_verbose) hw->print();
1380 break;
1381 }
1382
1383 if (eventNrOfThisFrame != eventNrOfOnsenTrgFrame && !isFakedData_event) {
1384 if (!(m_suppressErrorMask[c_nrFRAME_TNR_MM])) {
1385 B2WARNING("Frame TrigNr != ONSEN Trig Nr");
1386 B2DEBUG(29, "Frame TrigNr != ONSEN Trig Nr $" << hex << eventNrOfThisFrame << " != $" << eventNrOfOnsenTrgFrame);
1387 }
1388 m_errorMask[c_nrFRAME_TNR_MM] = true;
1389 }
1390
1391 if (Frame_Number == 0) {
1393 if (frame_type != EDHCFrameHeaderDataType::c_ONSEN_TRG) {
1394 if (!m_formatBonnDAQ) {
1395 if (!(m_suppressErrorMask[c_nrONSEN_TRG_FIRST])) B2WARNING("First frame is not a ONSEN Trigger frame");
1396 m_errorMask[c_nrONSEN_TRG_FIRST] = true;
1397 }
1398 }
1399 } else { // (Frame_Number != 0 &&
1401 if (frame_type == EDHCFrameHeaderDataType::c_ONSEN_TRG) {
1402 if (!(m_suppressErrorMask[c_nrONSEN_TRG_FIRST])) B2WARNING("More than one ONSEN Trigger frame");
1403 m_errorMask[c_nrONSEN_TRG_FIRST] = true;
1404 }
1405 }
1406
1407 if (!m_formatBonnDAQ) {
1408 if (Frame_Number == 1) {
1410 if (frame_type != EDHCFrameHeaderDataType::c_DHC_START) {
1411 if (!(m_suppressErrorMask[c_nrDHC_START_SECOND])) B2WARNING("Second frame is not a DHC start of subevent frame");
1412 m_errorMask[c_nrDHC_START_SECOND] = true;
1413 }
1414 } else { // (Frame_Number != 0 &&
1416 if (frame_type == EDHCFrameHeaderDataType::c_DHC_START) {
1417 if (!(m_suppressErrorMask[c_nrDHC_START_SECOND])) B2WARNING("More than one DHC start of subevent frame");
1418 m_errorMask[c_nrDHC_START_SECOND] = true;
1419 }
1420 }
1421 }
1422
1423 if (Frame_Number == Frames_in_event - 1) {
1425 if (frame_type != EDHCFrameHeaderDataType::c_DHC_END) {
1426 if (!(m_suppressErrorMask[c_nrDHC_END_MISS])) B2WARNING("Last frame is not a DHC end of subevent frame");
1427 m_errorMask[c_nrDHC_END_MISS] = true;
1428 }
1429
1431 if (countedDHEStartFrames != countedDHEEndFrames || countedDHEStartFrames != nr_active_dhe) {
1432 if (!(m_suppressErrorMask[c_nrDHE_ACTIVE]) || !(m_suppressErrorMask[c_nrDHE_START_WO_END])
1433 || !(m_suppressErrorMask[c_nrDHE_END_WO_START])) {
1434 B2WARNING("The number of DHE Start/End does not match the number of active DHE in DHC Header!");
1435 B2DEBUG(29, "The number of DHE Start/End does not match the number of active DHE in DHC Header! Header: " << nr_active_dhe <<
1436 " Start: " << countedDHEStartFrames << " End: " << countedDHEEndFrames << " Mask: $" << hex << mask_active_dhe << " in Event Nr " <<
1437 eventNrOfThisFrame);
1438 }
1439 if (countedDHEStartFrames == countedDHEEndFrames) m_errorMask[c_nrDHE_ACTIVE] = true;
1440 if (countedDHEStartFrames > countedDHEEndFrames) m_errorMask[c_nrDHE_START_WO_END] = true;
1441 if (countedDHEStartFrames < countedDHEEndFrames) m_errorMask[c_nrDHE_END_WO_START] = true;
1442 }
1443
1444 } else { // (Frame_Number != Frames_in_event - 1 &&
1446 if (frame_type == EDHCFrameHeaderDataType::c_DHC_END) {
1447 if (!(m_suppressErrorMask[c_nrDHC_END_DBL])) B2WARNING("More than one DHC end of subevent frame");
1448 m_errorMask[c_nrDHC_END_DBL] = true;
1449 }
1450 }
1451
1452 if (!m_formatBonnDAQ) {
1454 if (Frame_Number == 2 && nr_active_dhe != 0 && frame_type != EDHCFrameHeaderDataType::c_DHE_START) {
1455 if (!(m_suppressErrorMask[c_nrDHE_START_THIRD])) B2WARNING("Third frame is not a DHE start frame");
1456 m_errorMask[c_nrDHE_START_THIRD] = true;
1457 }
1458 }
1459
1460 if (frame_type != EDHCFrameHeaderDataType::c_ONSEN_ROI && frame_type != EDHCFrameHeaderDataType::c_ONSEN_TRG) {
1461 // actually, they should not be within Start and End, but better be sure.
1462 countedBytesInDHC += len;
1463 countedBytesInDHE += len;
1464 }
1465 B2DEBUG(29, "DHC/DHE $" << hex << countedBytesInDHC << ", $" << hex << countedBytesInDHE);
1466}
1467
1469{
1471 const int lut[32] = {
1472 0, 1, 1, 2, 1, 2, 2, 3, 1, 2, 2, 3, 2, 3, 3, 4,
1473 1, 2, 2, 3, 2, 3, 3, 4, 2, 3, 3, 4, 3, 4, 4, 5
1474 };
1475 return lut[i & 0x1F];
1476}
1477
@ c_ErrorIfAlreadyRegistered
If the object/array was already registered, produce an error (aborting initialisation).
Definition DataStore.h:72
void setDescription(const std::string &description)
Sets the description of the module.
Definition Module.cc:214
void setPropertyFlags(unsigned int propertyFlags)
Sets the flags for the module properties.
Definition Module.cc:208
Module()
Constructor.
Definition Module.cc:30
void setReturnValue(int value)
Sets the return value for this module as integer.
Definition Module.cc:220
@ c_ParallelProcessingCertified
This module can be run in parallel processing mode safely (All I/O must be done through the data stor...
Definition Module.h:80
void setErrorMask(const PXDErrorFlags &mask)
Set Error bit mask This should be the OR of error masks of all sub-objects (DHC, DHE)
void setGatedFlag(bool gm)
set gating info from the DHC END
PXDDAQDHEStatus & newDHE(Args &&... params)
Add new DHE information.
void setGatedHER(bool isher)
set HER/LER gating info from the DHC END
PXDDAQDHEStatus & dhe_back()
Returns PXDDAQDHEStatus for last DHE.
void setEndErrorInfo(uint32_t e)
set errorinfo from the DHC END
size_t dhe_size() const
Returns number of DHEs.
void setCounters(uint32_t raw, uint32_t red)
Set Data counters for reduction calculation.
PXDDAQDHPStatus & newDHP(Args &&... params)
New DHP information.
void setDHPFoundMask(unsigned short dhpmask)
set Mask for found DHPs with valid data
void setErrorMask(const PXDErrorFlags &mask)
Set Error bit mask.
auto addCM(PXDDAQDHPComMode &daqcm)
Add Common Mode information.
PXDDAQDHPStatus & dhp_back()
Returns PXDDAQDHPStatus for the last DHP.
void setEndErrorInfo(uint32_t e)
set erroinfo from the DHE END
void setCounters(uint32_t raw, uint32_t red)
Set Data counters for reduction calculation.
void setTruncated(void)
set Truncation
The PXD DAQ Packet Status class.
void setErrorMask(const PXDErrorFlags &mask)
Set Error bit mask This should be the OR of error masks of all sub-objects (DHC, DHE)
size_t dhc_size() const
Returns number of DHCs.
PXDDAQDHCStatus & dhc_back()
Returns PXDDAQDHCStatus for last DHC.
PXDDAQDHCStatus & newDHC(Args &&... params)
Add new DHC information.
static void map_rc_to_uv_IF_OB(unsigned int &row_u, unsigned int &col_v, const unsigned int dhp_id, const unsigned int dhe_ID)
Maps row/col of inner forward (IF) and outer backward (OB) modules of the PXD to U/V cell.
static void map_rc_to_uv_IB_OF(unsigned int &row_u, unsigned int &col_v, const unsigned int dhp_id, const unsigned int dhe_ID)
Maps row/cols of inner backward (IB) and outer forward (OF) modules of the PXD to U/V cell.
std::string m_PXDRawHitsName
The name of the StoreArray of PXDRawHits to be generated.
void initialize() override final
Initialize the module.
unsigned long m_meta_experiment
Experiment from MetaInfo.
PXDError::PXDErrorFlags m_errorMaskEvent
Error Mask set per packet / event.
StoreObjPtr< PXDDAQStatus > m_storeDAQEvtStats
Output array for DAQ Status.
bool m_doNotStore
Only unpack, but Do Not Store anything to file.
StoreArray< RawPXD > m_storeRawPXD
Input array for PXD Raw.
bool m_forceNoMapping
Force No Mapping even if DHH bit is requesting it.
std::string m_PXDDAQEvtStatsName
The name of the StoreObjPtr of PXDDAQStatus to be generated.
static void dump_roi(void *data, unsigned int frame_len)
dump to a file, helper function for debugging.
PXDError::PXDErrorFlags m_suppressErrorMask
Mask for suppressing selected error messages.
bool m_formatBonnDAQ
flag ONSEN or BonnDAQ format
unsigned int m_errorCounter[PXDError::ONSEN_MAX_TYPE_ERR]
Error counters.
PXDUnpackerOTModule()
Constructor defining the parameters.
StoreArray< PXDRawROIs > m_storeROIs
Output array for Raw ROIs.
unsigned int m_sendunfiltered
counter for send unfiltered
PXDError::PXDErrorFlags m_errorMaskDHC
Error Mask set per packet / DHC.
PXDError::PXDErrorFlags m_errorMaskPacket
Error Mask set per packet / packet.
std::string m_RawPXDsName
The name of the StoreArray of processed RawPXDs.
void unpack_dhp(void *data, unsigned int len, unsigned int dhe_first_readout_frame_lo, unsigned int dhe_ID, unsigned dhe_DHPport, unsigned dhe_reformat, VxdID vxd_id, PXDDAQPacketStatus &daqpktstat)
Unpack DHP data within one DHE frame.
void terminate() override final
Terminate the module.
StoreObjPtr< EventMetaData > m_eventMetaData
Input ptr for EventMetaData.
int m_last_dhp_readout_frame_lo[4]
some workaround check for continouous frame ids
unsigned long m_meta_subrun_nr
Subrun Number from MetaInfo.
void event() override final
do the unpacking
StoreArray< PXDRawAdc > m_storeRawAdc
Output array for Raw Adcs.
unsigned long m_meta_event_nr
Event Number from MetaInfo.
PXDError::PXDErrorFlags m_errorSkipPacketMask
Mask for error which stop package unpacking directly.
static int nr5bits(int i)
helper function to "count" nr of set bits within lower 5 bits.
void unpack_fce(unsigned short *data, unsigned int length, VxdID vxd_id)
Unpack DHP/FCE data within one DHE frame Not fully implemented as cluster format not 100% fixed.
std::string m_PXDRawAdcsName
The name of the StoreArray of PXDRawAdcs to be generated.
unsigned long m_meta_run_nr
Run Number from MetaInfo.
bool m_continueOnError
flag continue unpacking of frames even after error (for debugging)
unsigned int m_notaccepted
counter for not accepted events... should not happen TODO discussion ongoing with DAQ group
unsigned int m_meta_sec
Time(Tag) from MetaInfo, seconds (masked to lower bits)
void unpack_dhp_raw(void *data, unsigned int len, unsigned int dhe_ID, unsigned dhe_DHPport, VxdID vxd_id)
Unpack DHP RAW data within one DHE frame (pedestals, etc)
PXDError::PXDErrorFlags m_criticalErrorMask
Critical error mask which defines return value of task.
unsigned int m_sendrois
counter for send debug rois
static void dump_dhp(void *data, unsigned int frame_len)
dump to a file, helper function for debugging.
void unpack_rawpxd(RawPXD &px, int inx)
Unpack one event (several frames) stored in RawPXD object.
unsigned long long int m_meta_time
Time(Tag) from MetaInfo.
unsigned int m_maxDHPFrameDiff
Maximum DHP frame difference until error is reported.
PXDError::PXDErrorFlags m_errorMask
Error Mask set per packet / frame.
unsigned int m_unpackedEventsCount
Event counter.
StoreArray< PXDRawHit > m_storeRawHits
Output array for Raw Hits.
void unpack_dhc_frame(void *data, const int len, const int Frame_Number, const int Frames_in_event, PXDDAQPacketStatus &daqpktstat)
Unpack one frame (within an event).
bool m_forceMapping
Force Mapping even if DHH bit is not requesting it.
unsigned int m_meta_ticks
Time(Tag) from MetaInfo, Ticks of 127MHz.
bool m_verbose
give verbose unpacking information
PXDError::PXDErrorFlags m_errorMaskDHE
Error Mask set per packet / DHE.
std::string m_PXDRawROIsName
The name of the StoreArray of PXDRawROIs to be generated.
DHC frame wrapper class.
const dhc_ghost_frame * data_ghost_frame
data_ghost_frame
unsigned int getEventNrLo(void) const
get event nr lo (from data)
const dhc_end_frame * data_dhc_end_frame
data_dhc_end_frame
unsigned int getFixedSize(void)
get fixed size
int getFrameType(void)
get type of frame
const dhc_dhe_start_frame * data_dhe_start_frame
data_dhe_start_frame
void set(const void *d, unsigned int t)
set data and type (and length to 0)
const dhc_start_frame * data_dhc_start_frame
data_dhc_start_frame
const dhc_direct_readout_frame * data_direct_readout_frame
data_direct_readout_frame
const dhc_dhe_end_frame * data_dhe_end_frame
data_dhe_end_frame
void check_padding(PXDErrorFlags &errormask)
check padding and return it
const dhc_onsen_roi_frame * data_onsen_roi_frame
data_onsen_roi_frame
const dhc_direct_readout_frame_raw * data_direct_readout_frame_raw
data_direct_readout_frame_raw
const dhc_commode_frame * data_commode_frame
data_commode_frame
void check_crc(PXDErrorFlags &errormask, bool ignore_crc_flag=false)
check crc and return it
const dhc_onsen_trigger_frame * data_onsen_trigger_frame
data_onsen_trigger_frame
The Raw PXD class.
Definition RawPXD.h:27
virtual int * data(void)
get pointer to data
Definition RawPXD.cc:81
virtual int size() const
get size of buffer in 32 Bit words
Definition RawPXD.cc:76
Class to uniquely identify a any structure of the PXD and SVD.
Definition VxdID.h:32
Class to store variables with their name which were sent to the logging service.
void addParam(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.
Definition Module.h:559
#define REG_MODULE(moduleName)
Register the given module (without 'Module' suffix) with the framework.
Definition Module.h:649
std::tuple< uint8_t, uint16_t, uint8_t > PXDDAQDHPComMode
tuple of Chip ID (2 bit), Row (10 bit), Common Mode (6 bit)
Namespace to encapsulate code needed for simulation and reconstrucion of the PXD.
boost::endian::big_uint16_t ubig16_t
define alias ubig16_t
boost::endian::big_uint32_t ubig32_t
define alias ubig32_t
Abstract base class for different kinds of events.
STL namespace.
unsigned int getDHEId(void) const
get DHE Id (from word0)
unsigned int get_words(void) const
get words
unsigned int getDHEId(void) const
get DHE Id
unsigned int getErrorInfo(void) const
get error info
unsigned int getDHEId(void) const
get DHE Id (from word0)
unsigned short getEventNrLo(void) const
get trigger_nr_lo
unsigned short getTriggerGate(void) const
trigger gate (updated to 8 bit, before 10!)
unsigned short getStartFrameNr(void) const
last DHP frame before trigger
unsigned short getEventNrHi(void) const
get trigger_nr_hi
unsigned int getActiveDHPMask(void) const
get Active DHP Mask (from word0)
bool getDataReformattedFlag(void) const
get DataReformattedFlag (from word0)
unsigned short getDHEId(void) const
get DHE Id (from word0)
unsigned short getDHPPort(void) const
get DHP Port (from word0)
void print(void) const
print
unsigned int get_words(void) const
get words
bool isFakedData(void) const
is faked data
unsigned int get_dhc_id(void) const
get dhc id (from word0)
unsigned int getErrorInfo(void) const
get error info
DHC frame header word data struct.
unsigned short getFrameType(void) const
get type of frame
unsigned short getErrorFlag(void) const
get error flag
unsigned short getDHEId(void) const
get DHE Id (from word0)
unsigned short getDHPPort(void) const
get DDHP port (from word0)
unsigned int check_inner_crc(PXDErrorFlags &, unsigned int) const
check inner crc (currently not implemented/needed)
void check_error(PXDErrorFlags &errormask, int length, bool ignore_inv_size_flag=false) const
check error and return error mask
int getMinSize(void) const
4 byte header, ROIS (n*8), 4 byte copy of inner CRC, 4 byte outer CRC
unsigned short get_subrun1(void) const
get subrun1 (from trigtag1)
unsigned int get_trig_nr1(void) const
get trignr1
bool is_SendUnfiltered(void) const
is sendUnfiltered
bool is_SendROIs(void) const
is sendROIs
void check_error(PXDErrorFlags &errormask, bool ignore_datcon_flag=false, bool ignore_hltroi_magic_flag=false, bool ignore_merger_mm_flag=false) const
check error and return error mask
unsigned short get_run2(void) const
get run2 (from trigtag2)
bool is_fake_datcon(void) const
is fake datcon
unsigned short get_experiment1(void) const
get experiment1 (from trigtag1)
unsigned int get_trig_nr2(void) const
get trignr2
unsigned short get_experiment2(void) const
get experiment2
unsigned short get_subrun2(void) const
get subrun2 (from trigtag2)
bool is_Accepted(void) const
is accepted
unsigned short get_run1(void) const
get run1 (from trigtag1)
unsigned short get_gated_isher(void) const
get gated_isher (from word0)
unsigned short get_subrun(void) const
get subrun (from run_subrun)
unsigned short get_run(void) const
get run (from run_subrun)
unsigned short get_dhc_id(void) const
get dhc_id (from word0)
const ubig16_t time_tag_hi
time_tag_hi
const ubig16_t time_tag_mid
time_tag_mid
bool isFakedData(void) const
isFakedData
unsigned short getEventNrLo(void) const
get trigger_nr_lo
unsigned short get_experiment(void) const
get experiment (from exp_run)
unsigned short get_gated_flag(void) const
get gated_flag (from word0)
unsigned short get_active_dhe_mask(void) const
get active_dhe_mask (from word0)
const ubig16_t time_tag_lo_and_type
time_tag_lo_and_type
unsigned short getEventNrHi(void) const
get trigger_nr_hi