]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HLT/MUON/utils/dHLTdumpraw.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Update output information for new HLT-CDH bit.
[u/mrichter/AliRoot.git] / HLT / MUON / utils / dHLTdumpraw.cxx
CommitLineData
3dad52b0 1/**************************************************************************
2 * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
1d8ae082 16// $Id$
3dad52b0 17
18/**
19 * @file dHLTdumpraw.cxx
20 * @author Artur Szostak <artursz@iafrica.com>,
21 * Seforo Mohlalisi <seforomohlalisi@yahoo.co.uk>
887a669c 22 * @date 1 July 2007
3dad52b0 23 * @brief Command line utility to dump dHLT's internal raw data blocks.
24 */
25
26// We define NDEBUG for the AliHLTMUONDataBlockReader.h header file since this
27// program by definition handles corrupt data. So we do not need the assertions
28// in the AliHLTMUONDataBlockReader class to be checked.
9c234d7b 29#ifndef NDEBUG
3dad52b0 30#define NDEBUG
9c234d7b 31#endif
3dad52b0 32#include "AliHLTMUONDataBlockReader.h"
9c234d7b 33#if defined(DEBUG) && defined(NDEBUG)
3dad52b0 34#undef NDEBUG
9c234d7b 35#endif
36
3dad52b0 37#include "AliHLTMUONUtils.h"
887a669c 38#include "AliHLTMUONConstants.h"
39#include "Rtypes.h"
40#include "AliRawDataHeader.h"
41#include "AliMUONTrackerDDLDecoder.h"
42#include "AliMUONTriggerDDLDecoder.h"
43#include "AliHLTSystem.h"
44#include "AliHLTConfiguration.h"
45#include "AliLog.h"
46#include "TClassTable.h"
47#include "TString.h"
48#include "TRegexp.h"
3dad52b0 49
b09247a2 50#include <cstring>
fdee2dad 51#include <cstdlib>
887a669c 52#include <cstdio>
3dad52b0 53#include <cassert>
54#include <new>
55#include <fstream>
56
57#include <iostream>
58using std::cout;
59using std::cerr;
60using std::endl;
61using std::showbase;
62using std::noshowbase;
63using std::hex;
64using std::dec;
65
66#include <iomanip>
67using std::setw;
68using std::left;
887a669c 69using std::right;
3dad52b0 70using std::internal;
71
72
73#define CMDLINE_ERROR 1
74#define PARSE_ERROR 2
75#define SYSTEM_ERROR 3
76#define FATAL_ERROR 4
887a669c 77#define HLTSYSTEM_ERROR 5
78
79
80// Adding enum types for extending AliHLTMUONDataBlockType with the
81// raw DDL data types.
82enum AliHLTMUONDDLRawDataType
83{
84 kTrackerDDLRawData = 10,
85 kTriggerDDLRawData = 11
86};
3dad52b0 87
88
89void PrintRubbishData(AliHLTUInt32_t offset, const char* padByte, AliHLTUInt32_t padCount)
90{
91 if (padCount == 0) return;
92
93 cerr << "ERROR: Found the following unexpected rubbish data at the"
94 " end of the data block:" << endl;
95 cerr << "Byte #\tValue\tCharacter" << endl;
96 for (AliHLTUInt32_t i = 0; i < padCount; i++)
97 {
98 short value = short(padByte[i]) & 0xFF;
99 char character = padByte[i];
100 cerr << offset + i + 1 << "\t"
101 << noshowbase << hex << "0x" << value << dec << "\t"
102 << character << endl;
103 }
104}
105
106
a090ff22 107void PrintBitPattern(AliHLTUInt32_t value, int bitcount = 32)
108{
109 // Prints a bit pattern to cout.
110
111 for (int i = bitcount-1; i >= 0; i--)
112 {
113 if ( ((value >> i) & 0x1) == 1 )
114 cout << "1";
115 else
116 cout << "0";
117 }
118}
119
120
3dad52b0 121template <typename FieldType>
122int CheckHeaderField(
123 FieldType& field, const char* buffer, unsigned long bufferSize,
124 bool continueParse
125 )
126{
127 const char* fieldptr = reinterpret_cast<const char*>(&field);
128 const char* endptr = buffer + bufferSize;
129 AliHLTUInt32_t bufferRemaining = endptr > fieldptr ? endptr - fieldptr : 0;
130
131 if (bufferRemaining < sizeof(field))
132 {
133 cout << "..." << endl; // We may be half way through printing a line so end it.
134 cerr << "ERROR: The data block is too short. The header is corrupt." << endl;
135 if (continueParse)
136 {
137 AliHLTUInt32_t offset = fieldptr - buffer;
138 PrintRubbishData(offset, fieldptr, bufferRemaining);
139 }
140 return PARSE_ERROR;
141 }
142 return EXIT_SUCCESS;
143}
144
145
146template <typename FieldType>
147int CheckField(
148 FieldType& field, const char* buffer, unsigned long bufferSize,
149 bool continueParse
150 )
151{
152 const char* fieldptr = reinterpret_cast<const char*>(&field);
153 const char* endptr = buffer + bufferSize;
154 AliHLTUInt32_t bufferRemaining = endptr > fieldptr ? endptr - fieldptr : 0;
155
156 if (bufferRemaining < sizeof(field))
157 {
158 cout << "..." << endl; // We may be half way through printing a line so end it.
159 cerr << "ERROR: The data block is too short. The data is corrupt." << endl;
160 if (continueParse)
161 {
162 AliHLTUInt32_t offset = fieldptr - buffer;
163 PrintRubbishData(offset, fieldptr, bufferRemaining);
164 }
165 return PARSE_ERROR;
166 }
167 return EXIT_SUCCESS;
168}
169
170
171template <typename BlockType>
172int CheckCommonHeader(
63ae97e0 173 BlockType& block, const char* /*buffer*/, unsigned long bufferSize,
3dad52b0 174 bool continueParse
175 )
176{
177 int result = EXIT_SUCCESS;
178
179 // Check the fRecordWidth field in the common header.
180 if (block.CommonBlockHeader().fRecordWidth !=
181 sizeof(typename BlockType::ElementType))
182 {
183 cerr << "ERROR: The record width found in the header is incorrect."
184 " Found a record width of "
185 << block.CommonBlockHeader().fRecordWidth << " bytes, but expected"
186 " a value of " << sizeof(typename BlockType::ElementType)
187 << " bytes." << endl;
188 result = PARSE_ERROR;
189 if (not continueParse) return result;
190 }
191
192 if (not block.BufferSizeOk())
193 {
194 cerr << "ERROR: The size of the file is incorrect. It is "
195 << bufferSize << " bytes big, but according"
196 " to the data block header it should be " << block.BytesUsed()
197 << " bytes." << endl;
198 result = PARSE_ERROR;
199 if (not continueParse) return result;
200 }
201
202 return result;
203}
204
205
206template <typename BlockType>
207AliHLTUInt32_t CalculateNEntries(BlockType& block, unsigned long bufferSize)
208{
209 // Calculate how many entries we can display. If the buffer size is correct
210 // we just use the number of entries the block specifies. Otherwise we need
211 // to calculate it from the buffer size.
212 AliHLTUInt32_t nentries;
213 if (block.BytesUsed() == bufferSize)
214 {
215 nentries = block.Nentries();
216 }
217 else
218 {
219 AliHLTInt32_t dataSize = bufferSize
220 - sizeof(typename BlockType::HeaderType);
221 nentries = dataSize / sizeof(typename BlockType::ElementType);
222 if (dataSize % sizeof(typename BlockType::ElementType) > 0)
223 nentries++;
224 }
225 return nentries;
226}
227
228
887a669c 229namespace
230{
231 /**
232 * Common methods for DDL decoder event handlers.
233 */
234 class AliDecoderHandler
235 {
236 public:
237 AliDecoderHandler() :
238 fBufferStart(NULL),
239 fDumpStart(NULL),
240 fDumpData(false)
241 {
242 }
243
244 virtual ~AliDecoderHandler() {}
245
246 protected:
247
248 // Do not allow copying of this class.
249 AliDecoderHandler(const AliDecoderHandler& obj);
250 AliDecoderHandler& operator = (const AliDecoderHandler& obj);
251
252 void HandleError(
253 const char* errorMessage, int errorCode,
254 const char* errorCodeString, const void* location
a090ff22 255 );
887a669c 256
a090ff22 257 void TryDumpCorruptData(const void* dumpEnd);
887a669c 258
259 const void* fBufferStart; ///< Start location of buffer.
260 const void* fDumpStart; ///< Start location of corrupt data to dump.
261 bool fDumpData; ///< Flag indicating if fDumpStart points to corrupt data and should be dumped.
262 };
a090ff22 263
264
265 void AliDecoderHandler::HandleError(
266 const char* errorMessage, int errorCode,
267 const char* errorCodeString, const void* location
268 )
269 {
270 unsigned long offset = (unsigned long)location - (unsigned long)fBufferStart
271 + sizeof(AliRawDataHeader);
272
273 cerr << "ERROR: " << errorMessage
274 << " [Error code = " << errorCode << " ("
275 << errorCodeString << "), at byte "
276 << offset << " (" << noshowbase << hex << "0x"
277 << offset << dec << ")]" << endl;
278
279 if (fDumpStart == NULL) fDumpStart = location;
280 fDumpData = true;
281 }
282
283
284 void AliDecoderHandler::TryDumpCorruptData(const void* dumpEnd)
285 {
286 if (dumpEnd < fDumpStart) return;
287 if (not fDumpData) return;
288
289 unsigned long startOffset = (unsigned long)fDumpStart - (unsigned long)fBufferStart
290 + sizeof(AliRawDataHeader);
291 unsigned long endOffset = (unsigned long)dumpEnd - (unsigned long)fBufferStart
292 + sizeof(AliRawDataHeader);
293 if (endOffset - startOffset > 264)
294 {
295 endOffset = startOffset + 264;
296 dumpEnd = reinterpret_cast<const char*>(fBufferStart) + endOffset;
297 }
298 cerr << "Dumping corrupt data words from byte " << startOffset
299 << " (" << noshowbase << hex << "0x" << startOffset
300 << dec << "), to byte " << endOffset << " (" << noshowbase
301 << hex << "0x" << endOffset << dec << "):" << endl;
302 const UInt_t* start = reinterpret_cast<const UInt_t*>(fDumpStart);
303 const UInt_t* end = reinterpret_cast<const UInt_t*>(dumpEnd);
304 cerr << " Start byte | Data words" << endl;
305 for (const UInt_t* current = start; current < end; current++)
306 {
307 unsigned long currentByte = (unsigned long)current
308 - (unsigned long)fBufferStart + sizeof(AliRawDataHeader);
309 cerr << right << setw(9) << dec << currentByte << setw(0)
310 << " 0x" << left << setw(7) << noshowbase << hex
311 << currentByte << setw(0) << right << " | ";
312 char fillChar = cerr.fill();
313 cerr.fill('0');
314 for (int i = 0; i < 4 and current < end; i++, current++)
315 {
316 cerr << noshowbase << hex << "0x" << setw(8)
317 << (*current) << setw(0) << dec << " ";
318 }
319 cerr.fill(fillChar);
320 cerr << endl;
321 }
322 fDumpStart = NULL;
323 fDumpData = false;
324 }
887a669c 325
326 /**
327 * Event handler for the tracker DDL decoder.
328 * It simply prints the structure to standard output.
329 */
330 class AliTrackerDecoderHandler :
331 public AliMUONTrackerDDLDecoderEventHandler, public AliDecoderHandler
332 {
333 public:
334 AliTrackerDecoderHandler() :
335 AliMUONTrackerDDLDecoderEventHandler(),
a090ff22 336 AliDecoderHandler(),
337 fBlockNum(0),
338 fDSPNum(0),
339 fBusPatchNum(0)
887a669c 340 {}
341
342 virtual ~AliTrackerDecoderHandler() {}
343
344 void OnNewBuffer(const void* buffer, UInt_t /*bufferSize*/)
345 {
346 fBufferStart = buffer;
a090ff22 347 fBlockNum = fDSPNum = fBusPatchNum = 0;
887a669c 348 }
349
350 void OnEndOfBuffer(const void* buffer, UInt_t bufferSize)
351 {
352 const char* bufferEnd =
353 reinterpret_cast<const char*>(buffer) + bufferSize;
354 TryDumpCorruptData(bufferEnd);
355 }
356
a090ff22 357 void OnNewBlock(const AliMUONBlockHeaderStruct* header, const void* /*data*/);
887a669c 358
a090ff22 359 void OnNewDSP(const AliMUONDSPHeaderStruct* header, const void* /*data*/);
360
361 void OnNewBusPatch(const AliMUONBusPatchHeaderStruct* header, const void* /*data*/);
362
363 void OnData(UInt_t data, bool parityError);
364
365 void OnError(ErrorCode error, const void* location);
887a669c 366
a090ff22 367 private:
368
369 UInt_t fBlockNum; ///< Number of block being processed [1..maxBlock].
370 UInt_t fDSPNum; ///< Number of DSP being processed [1..maxDSP].
371 UInt_t fBusPatchNum; ///< Number of bus patch being processed [1..maxBusPatch].
372 };
373
374
375 void AliTrackerDecoderHandler::OnNewBlock(const AliMUONBlockHeaderStruct* header, const void* /*data*/)
376 {
377 TryDumpCorruptData(header);
378 fBlockNum++;
379 char fillChar = cout.fill(); // store current fill char to restore it.
380 cout << "================================ Block header "
381 << setw(3) << fBlockNum << setw(0)
382 << " ================================" << endl;
383 cout << " Data key word for block : 0x" << noshowbase << setw(8) << setfill('0')
384 << hex << header->fDataKey << dec << setfill(fillChar) << setw(0) << endl;
385 cout << " Total Length of block : " << header->fTotalLength << endl;
386 cout << " Length of raw data : " << header->fLength << endl;
387 cout << " DSP ID : " << header->fDSPId << endl;
388 cout << " L0 trigger word : " << header->fL0Trigger << " (0x"
389 << noshowbase << setw(8) << setfill('0') << hex << header->fL0Trigger
390 << dec << setw(0) << setfill(fillChar) << ")" << endl;
391 cout << " Mini event ID : " << header->fMiniEventId << " (0x"
392 << noshowbase << setw(8) << setfill('0') << hex << header->fMiniEventId
393 << dec << setw(0) << setfill(fillChar) << ")" << endl;
394 cout << "Event ID In Bunch Crossing : " << header->fEventId1 << " (0x"
395 << noshowbase << setw(8) << setfill('0') << hex << header->fEventId1
396 << dec << setw(0) << setfill(fillChar) << ")" << endl;
397 cout << " Event ID In Orbit Number : " << header->fEventId2 << " (0x"
398 << noshowbase << setw(8) << setfill('0') << hex << header->fEventId2
399 << dec << setw(0) << setfill(fillChar) << ")" << endl;
400 }
401
402
403 void AliTrackerDecoderHandler::OnNewDSP(const AliMUONDSPHeaderStruct* header, const void* /*data*/)
404 {
405 TryDumpCorruptData(header);
406 fDSPNum++;
407 char fillChar = cout.fill(); // store current fill char to restore it.
408 cout << "================================= DSP header "
409 << setw(3) << fDSPNum << setw(0)
410 << " =================================" << endl;
411 cout << " Data key word for FRT : 0x" << noshowbase
412 << setw(8) << setfill('0') << hex << header->fDataKey << dec
413 << setfill(fillChar) << setw(0) << endl;
414 cout << " Total length of structure : " << header->fTotalLength << endl;
415 cout << " Length of raw data : " << header->fLength << endl;
416 cout << " DSP ID : " << header->fDSPId << endl;
417 cout << " L1 accept in block Structure (CRT) : " << header->fBlkL1ATrigger
418 << " (0x" << noshowbase << setw(8) << setfill('0') << hex
419 << header->fBlkL1ATrigger << dec << setw(0) << setfill(fillChar) << ")" << endl;
420 cout << " Mini event ID in bunch crossing : " << header->fMiniEventId
421 << " (0x" << noshowbase << setw(8) << setfill('0') << hex
422 << header->fMiniEventId << dec << setw(0) << setfill(fillChar) << ")" << endl;
423 cout << "Number of L1 accept in DSP Structure (FRT) : " << header->fL1ATrigger << endl;
424 cout << "Number of L1 reject in DSP Structure (FRT) : " << header->fL1RTrigger << endl;
425 const char* paddingWordValue = ") (INVALID)";
426 if (header->fPaddingWord == 0) paddingWordValue = ") (false)";
427 if (header->fPaddingWord == 1) paddingWordValue = ") (true)";
428 cout << " Padding word set for 64 bits transfer : " << header->fPaddingWord
429 << " (0x" << noshowbase << setw(8) << setfill('0') << hex
430 << header->fPaddingWord << setw(0) << setfill(fillChar) << dec
431 << paddingWordValue << endl;
432 cout << " Error word : " << header->fErrorWord
433 << " (" << noshowbase << setw(8) << setfill('0') << hex
434 << header->fErrorWord << setw(0) << setfill(fillChar)
435 << dec << ")" << endl;
436 }
437
438
439 void AliTrackerDecoderHandler::OnNewBusPatch(const AliMUONBusPatchHeaderStruct* header, const void* /*data*/)
440 {
441 TryDumpCorruptData(header);
442 fBusPatchNum++;
443 char fillChar = cout.fill(); // store current fill char to restore it.
444 cout << "============================== Bus patch header "
445 << setw(3) << fBusPatchNum << setw(0)
446 << " ==============================" << endl;
447 cout << "Data key word for bus patch : 0x" << noshowbase << setw(8)
448 << setfill('0') << hex << header->fDataKey << dec
449 << setfill(fillChar) << setw(0) << endl;
450 cout << " Total length of structure : " << header->fTotalLength << endl;
451 cout << " Length of raw data : " << header->fLength << endl;
452 cout << " Bus patch ID : " << header->fBusPatchId << endl;
453 if (header->fLength > 0)
887a669c 454 {
a090ff22 455 cout << " Raw bits | Manu ID | Manu channel | ADC Signal" << endl;
456 cout << "----------------------------------------------------------" << endl;
887a669c 457 }
a090ff22 458 }
459
460
461 void AliTrackerDecoderHandler::OnData(UInt_t data, bool parityError)
462 {
463 UShort_t manuId, adc;
464 UChar_t manuChannel;
465 UnpackADC(data, manuId, manuChannel, adc);
466 char fillChar = cout.fill(); // store current fill char to restore it.
467 cout << noshowbase << " 0x" << setfill('0') << hex << setw(8) << data
468 << setw(0) << dec << setfill(fillChar) << " "
469 << setw(12) << manuId << setw(0) << " "
470 << setw(12) << (UInt_t)manuChannel << setw(0) << " "
471 << setw(12) << adc << setw(0);
472 if (parityError)
887a669c 473 {
a090ff22 474 cout << " <= WARNING: Raw data word with parity error!" << endl;
887a669c 475 }
a090ff22 476 else
887a669c 477 {
a090ff22 478 cout << endl;
887a669c 479 }
a090ff22 480 }
481
482
483 void AliTrackerDecoderHandler::OnError(ErrorCode error, const void* location)
484 {
485 TryDumpCorruptData(location);
486 HandleError(
487 ErrorCodeToMessage(error), error,
488 ErrorCodeToString(error), location
489 );
490 }
887a669c 491
492 /**
493 * Event handler for the trigger DDL decoder.
494 * It simply prints the structure to standard output.
495 */
496 class AliTriggerDecoderHandler :
497 public AliMUONTriggerDDLDecoderEventHandler, public AliDecoderHandler
498 {
499 public:
500 AliTriggerDecoderHandler() :
501 AliMUONTriggerDDLDecoderEventHandler(),
a090ff22 502 AliDecoderHandler(),
503 fRegNum(0),
504 fLocNum(0)
887a669c 505 {}
506
507 virtual ~AliTriggerDecoderHandler() {}
508
509 void OnNewBuffer(const void* buffer, UInt_t /*bufferSize*/)
510 {
511 fBufferStart = buffer;
a090ff22 512 fRegNum = fLocNum = 0;
887a669c 513 }
514
515 void OnEndOfBuffer(const void* buffer, UInt_t bufferSize)
516 {
517 const char* bufferEnd =
518 reinterpret_cast<const char*>(buffer) + bufferSize;
519 TryDumpCorruptData(bufferEnd);
520 }
521
a090ff22 522 const char* EventTypeToString(UInt_t type);
523
524 const char* DarcTypeToString(UInt_t type);
525
887a669c 526 void OnDarcHeader(
527 UInt_t header,
528 const AliMUONDarcScalarsStruct* scalars,
529 const void* data
a090ff22 530 );
887a669c 531
532 void OnGlobalHeader(
533 const AliMUONGlobalHeaderStruct* header,
a090ff22 534 const AliMUONGlobalScalarsStruct* scalars,
887a669c 535 const void* /*data*/
a090ff22 536 );
887a669c 537
538 void OnNewRegionalStruct(
539 const AliMUONRegionalHeaderStruct* regionalStruct,
a090ff22 540 const AliMUONRegionalScalarsStruct* scalars,
887a669c 541 const void* /*data*/
a090ff22 542 );
887a669c 543
544 void OnLocalStruct(
545 const AliMUONLocalInfoStruct* localStruct,
a090ff22 546 const AliMUONLocalScalarsStruct* scalars
547 );
548
549 void OnError(ErrorCode error, const void* location);
550
551 private:
552
553 UInt_t fRegNum; ///< Number of block being processed [1..maxReg].
554 UInt_t fLocNum; ///< Number of DSP being processed [1..maxLoc].
555 };
556
557
558 const char* AliTriggerDecoderHandler::EventTypeToString(UInt_t type)
559 {
560 switch (type)
887a669c 561 {
a090ff22 562 case 0: return "Other software trigger";
563 case 1: return "Physics";
564 case 2: return "Start of run";
565 case 3: return "End of run";
566 default: return "UNKNOWN";
887a669c 567 }
a090ff22 568 }
569
570
571 const char* AliTriggerDecoderHandler::DarcTypeToString(UInt_t type)
572 {
573 typedef AliMUONTriggerDDLDecoder<AliMUONTriggerDDLDecoderEventHandler> AliDec;
574 if (type == AliDec::DarcDefaultType()) return "Default";
575 if (type == AliDec::DarcVadorhType()) return "Vadorh";
576 return "UNKNOWN";
577 }
578
579
580 void AliTriggerDecoderHandler::OnDarcHeader(
581 UInt_t header,
582 const AliMUONDarcScalarsStruct* scalars,
583 const void* data
584 )
585 {
586 if (scalars != NULL)
587 TryDumpCorruptData(scalars);
588 else
589 TryDumpCorruptData(data);
887a669c 590
a090ff22 591 cout << "================================ DARC header =====================================" << endl;
592 char fillChar = cout.fill(); // store current fill char to restore it.
593 cout << noshowbase << " Header bit pattern : 0x" << setfill('0') << hex << setw(8)
594 << header << setw(0) << dec << setfill(fillChar) << endl;
595 UInt_t eventType = GetDarcEventType(header);
596 cout << " Event type : " << eventType << " ("
597 << EventTypeToString(eventType) << ")" << endl;
598 cout << " Application type : " << ((header >> 27) & 0x7) << endl;
599 UInt_t darcType = GetDarcType(header);
600 cout << " DARC type : " << darcType << " ("
601 << DarcTypeToString(darcType) << ")" << endl;
602 cout << " Serial number : " << (UInt_t)GetDarcSerialNb(header) << endl;
603 cout << " Version : " << (UInt_t)GetDarcVersion(header) << endl;
604 cout << " VME trigger used : " << boolalpha << GetDarcVMETrig(header) << endl;
605 cout << "Global card data occurrence : " << boolalpha << GetDarcGlobalFlag(header) << endl;
606 cout << " CTP or LTU interfaced : " << boolalpha << GetDarcCTPTrig(header) << endl;
607 cout << " DAQ interfaced : " << boolalpha << GetDarcDAQFlag(header) << endl;
608 cout << " Regional cards occurrence : 0x" << noshowbase << hex << setw(2) << setfill('0')
609 << (UInt_t)GetDarcRegPattern(header) << dec << setfill(fillChar) << setw(0)
610 << " [bits: ";
611 PrintBitPattern(((UInt_t)GetDarcRegPattern(header) >> 4) & 0xF, 4); cout << " ";
612 PrintBitPattern(((UInt_t)GetDarcRegPattern(header) >> 0) & 0xF, 4);
613 cout << "]" << endl;
614
615 cout << "================================ DARC scalars ====================================" << endl;
616 if (scalars != NULL)
887a669c 617 {
a090ff22 618 cout << " Trigger | Received | Used" << endl;
619 cout << "----------------------------------" << endl;
620 cout << " L0 " << setw(8) << ((scalars->fL0R >> 16) & 0xFFFF) << setw(0)
621 << " " << setw(8) << ((scalars->fL0R >> 0) & 0xFFFF) << setw(0) << endl;
622 cout << " L1 physics " << setw(8) << ((scalars->fL1P >> 16) & 0xFFFF) << setw(0)
623 << " " << setw(8) << ((scalars->fL1P >> 0) & 0xFFFF) << setw(0) << endl;
624 cout << "L1 software " << setw(8) << ((scalars->fL1S >> 16) & 0xFFFF) << setw(0)
625 << " " << setw(8) << ((scalars->fL1S >> 0) & 0xFFFF) << setw(0) << endl;
626 cout << " L2 accept " << setw(8) << ((scalars->fL2A >> 16) & 0xFFFF) << setw(0)
627 << " " << setw(8) << ((scalars->fL2A >> 0) & 0xFFFF) << setw(0) << endl;
628 cout << " L2 reject " << setw(8) << ((scalars->fL2R >> 16) & 0xFFFF) << setw(0)
629 << " " << setw(8) << ((scalars->fL2R >> 0) & 0xFFFF) << setw(0) << endl;
630 cout << " Clock : " << scalars->fClk << endl;
631 cout << "Hold (dead time) : " << scalars->fHold << endl;
632 cout << " Spare word : " << scalars->fSpare << " (0x"
633 << setw(8) << setfill('0') << hex << scalars->fSpare <<
634 setw(0) << setfill(fillChar) << dec << ")" << endl;
887a669c 635 }
a090ff22 636 else
637 {
638 cout << "(None found)" << endl;
639 }
640 }
641
642
643 void AliTriggerDecoderHandler::OnGlobalHeader(
644 const AliMUONGlobalHeaderStruct* header,
645 const AliMUONGlobalScalarsStruct* scalars,
646 const void* /*data*/
647 )
648 {
649 TryDumpCorruptData(header);
650
651 cout << "=============================== Global header ====================================" << endl;
652 cout << "Global input from regional controllers:" << endl;
653 cout << " | | High pT | Low pT " << endl;
654 cout << " | Bit pattern | Single mu | mu pair | Single mu | mu pair " << endl;
655 cout << "Input | Hex | Binary | -ve | +ve | unlike | like | -ve | +ve | unlike | like" << endl;
656 cout << "--------------------------------------------------------------------------------" << endl;
657 char fillChar = cout.fill(); // store current fill char to restore it.
658 for (int i = 0; i < 4; i++)
659 {
660 union
661 {
662 UInt_t fWord;
663 UChar_t fByte[4];
664 } input;
665 input.fWord = header->fInput[i];
666 for (int j = 0; j < 4; j++)
667 {
668 cout << setw(5) << (i*4+j) << " 0x" << noshowbase << setw(2)
669 << setfill('0') << hex << (UInt_t)input.fByte[j] << setw(0)
670 << setfill(fillChar) << dec << " ";
671 PrintBitPattern(input.fByte[j], 8);
672 cout << ((((input.fByte[j] >> 7) & 0x1) == 1) ? " yes" : " no");
673 cout << ((((input.fByte[j] >> 6) & 0x1) == 1) ? " yes" : " no");
674 cout << ((((input.fByte[j] >> 5) & 0x1) == 1) ? " yes " : " no ");
675 cout << ((((input.fByte[j] >> 4) & 0x1) == 1) ? " yes" : " no");
676 cout << ((((input.fByte[j] >> 3) & 0x1) == 1) ? " yes" : " no");
677 cout << ((((input.fByte[j] >> 2) & 0x1) == 1) ? " yes" : " no");
678 cout << ((((input.fByte[j] >> 1) & 0x1) == 1) ? " yes " : " no ");
679 cout << ((((input.fByte[j] >> 0) & 0x1) == 1) ? " yes" : " no");
680 cout << endl;
681 }
682 }
683 cout << " Global ouput : 0x" << noshowbase << setw(2) << setfill('0') << hex
684 << (UInt_t)GetGlobalOutput(header) << setw(0) << setfill(fillChar) << dec << " [Bits: ";
685 PrintBitPattern(((UInt_t)GetGlobalOutput(header) >> 4) & 0xF, 4); cout << " ";
686 PrintBitPattern(((UInt_t)GetGlobalOutput(header) >> 0) & 0xF, 4);
687 cout << "]" << endl;
688 cout << " [ unlike sign | like sign | single muon ]" << endl;
689 cout << " [ high | low | high | low | high | low ]" << endl;
690 cout << " [-----------------------------------------]" << endl;
691 cout << " [ ";
692 cout << ((((GetGlobalOutput(header) >> 5) & 0x1) == 1) ? " yes" : " no");
693 cout << ((((GetGlobalOutput(header) >> 4) & 0x1) == 1) ? " yes" : " no");
694 cout << ((((GetGlobalOutput(header) >> 3) & 0x1) == 1) ? " yes" : " no");
695 cout << ((((GetGlobalOutput(header) >> 2) & 0x1) == 1) ? " yes" : " no");
696 cout << ((((GetGlobalOutput(header) >> 1) & 0x1) == 1) ? " yes" : " no");
697 cout << ((((GetGlobalOutput(header) >> 0) & 0x1) == 1) ? " yes" : " no");
698 cout << " ]" << endl;
699 cout << "Global configuration : 0x" << noshowbase << setw(4) << setfill('0') << hex
700 << GetGlobalConfig(header) << setw(0) << setfill(fillChar) << dec << " [Bits: ";
701 PrintBitPattern(((UInt_t)GetGlobalConfig(header) >> 12) & 0xF, 4); cout << " ";
702 PrintBitPattern(((UInt_t)GetGlobalConfig(header) >> 8) & 0xF, 4); cout << " ";
703 PrintBitPattern(((UInt_t)GetGlobalConfig(header) >> 4) & 0xF, 4); cout << " ";
704 PrintBitPattern(((UInt_t)GetGlobalConfig(header) >> 0) & 0xF, 4);
705 cout << "]" << endl;
706
707 cout << "=============================== Global scalars ===================================" << endl;
708 if (scalars != NULL)
709 {
710 cout << " Number of L0 triggers : " << scalars->fL0 << endl;
711 cout << " Number of clock cycles : " << scalars->fClk << endl;
712 cout << " Number of unlike mu pair low pT : " << scalars->fScaler[0] << endl;
713 cout << "Number of unlike mu pair high pT : " << scalars->fScaler[1] << endl;
714 cout << " Number of like mu pair low pT : " << scalars->fScaler[2] << endl;
715 cout << " Number of like mu pair high pT : " << scalars->fScaler[3] << endl;
716 cout << " Number of single mu low pT : " << scalars->fScaler[4] << endl;
717 cout << " Number of single mu high pT : " << scalars->fScaler[5] << endl;
718 cout << " Hold (dead time) : " << scalars->fHold << endl;
719 cout << " Spare word : " << scalars->fSpare << " (0x"
720 << setw(8) << setfill('0') << hex << scalars->fSpare <<
721 setw(0) << setfill(fillChar) << dec << ")" << endl;
722 }
723 else
724 {
725 cout << "(None found)" << endl;
726 }
727 }
728
729
730 void AliTriggerDecoderHandler::OnNewRegionalStruct(
731 const AliMUONRegionalHeaderStruct* regionalStruct,
732 const AliMUONRegionalScalarsStruct* scalars,
733 const void* /*data*/
734 )
735 {
736 TryDumpCorruptData(regionalStruct);
737
738 fRegNum++;
739 cout << "========================= Regional structure header "
740 << setw(3) << fRegNum << setw(0)
741 << " ==========================" << endl;
742 char fillChar = cout.fill(); // store current fill char to restore it.
743 cout << "Darc word bit pattern : 0x" << noshowbase << setw(8) << setfill('0')
744 << hex << regionalStruct->fDarcWord << setw(0)
745 << setfill(fillChar) << dec << endl;
746 UShort_t errBits = GetRegionalErrorBits(regionalStruct);
747 cout << " [ Error type : "
748 << (((errBits >> 7) & 0x1) == 1 ? "1 (Physics) " : "0 (Software)")
749 << " ]" << endl;
750 cout << " [ Regional error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
751 cout << " [ Full error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
752 cout << " [ Empty error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
753 cout << " [ DARC L2 reject error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
754 cout << " [ DARC L2 error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
755 cout << " [ DARC L1 error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
756 cout << " [ DARC L0 error : " << ((errBits >> 6) & 0x1) << " ]" << endl;
757
758 cout << " [ FPGA number in card : "
759 << (UInt_t)GetRegionalFPGANumber(regionalStruct) << " (";
760 PrintBitPattern((UInt_t)GetRegionalFPGANumber(regionalStruct), 3);
761 cout << "b) ]" << endl;
762
763 cout << " [ Physics trigger : " << setw(13) << boolalpha << left
764 << GetRegionalDarcPhysFlag(regionalStruct) << setw(0)
765 << right << "]" << endl;
766 cout << " [ Regional present : " << setw(13) << boolalpha << left
767 << GetRegionalPresentFlag(regionalStruct) << setw(0)
768 << right << "]" << endl;
769 cout << " [ RAM not full : " << setw(13) << boolalpha << left
770 << GetRegionalRamNotFullFlag(regionalStruct) << setw(0)
771 << right << "]" << endl;
772 cout << " [ RAM not empty : " << setw(13) << boolalpha << left
773 << GetRegionalRamNotEmptyFlag(regionalStruct) << setw(0)
774 << right << "]" << endl;
775 cout << " [ L2 rejected : " << setw(13) << boolalpha << left
776 << GetRegionalL2RejStatus(regionalStruct) << setw(0)
777 << right << "]" << endl;
778 cout << " [ L2 accepted : " << setw(13) << boolalpha << left
779 << GetRegionalL2AccStatus(regionalStruct) << setw(0)
780 << right << "]" << endl;
781 cout << " [ L1 found : " << setw(13) << boolalpha << left
782 << GetRegionalL1Status(regionalStruct) << setw(0)
783 << right << "]" << endl;
784 cout << " [ L0 found : " << setw(13) << boolalpha << left
785 << GetRegionalL0Status(regionalStruct) << setw(0)
786 << right << "]" << endl;
787 cout << " [ No. of events in RAM : " << setw(13) << left
788 << (UInt_t)GetRegionalEventInRam(regionalStruct) << setw(0)
789 << right << "]" << endl;
790 cout << " [ Busy word : 0x"
791 << (UInt_t)GetRegionalBusy(regionalStruct) << " (";
792 PrintBitPattern((UInt_t)GetRegionalBusy(regionalStruct), 4);
793 cout << "b) ]" << endl;
794
795 cout << "Regional word bit pattern : 0x" << noshowbase << setw(8) << setfill('0')
796 << hex << regionalStruct->fWord << setw(0)
797 << setfill(fillChar) << dec << endl;
798 cout << " [ Physics trigger occurred : " << setw(27) << boolalpha << left
799 << (UInt_t)GetRegionalPhysFlag(regionalStruct) << setw(0)
800 << right << "]" << endl;
801 cout << " [ Number of resets : " << setw(27) << left
802 << (UInt_t)GetRegionalResetNb(regionalStruct) << setw(0)
803 << right << "]" << endl;
804 cout << " [ Controller Serial number : " << setw(27) << left
805 << (UInt_t)GetRegionalSerialNb(regionalStruct) << setw(0)
806 << right << "]" << endl;
807 cout << " [ Regional crate ID : " << setw(27) << left
808 << (UInt_t)GetRegionalId(regionalStruct) << setw(0)
809 << right << "]" << endl;
810 cout << " [ FPGA software version : " << setw(27) << left
811 << (UInt_t)GetRegionalVersion(regionalStruct) << setw(0)
812 << right << "]" << endl;
813 UInt_t output = GetRegionalOutput(regionalStruct);
814 cout << " [ Regional output : 0x" << setw(2) << setfill('0') << hex
815 << output << dec << setw(0) << setfill(fillChar) << setw(0) << right << " (";
816 PrintBitPattern(output, 8);
817 cout << "b) ]" << endl;
818
819 cout << " [ High pT | Low pT ]" << endl;
820 cout << " [ Single mu | mu pair | Single mu | mu pair ]" << endl;
821 cout << " [ -ve | +ve | unlike | like | -ve | +ve | unlike | like ]" << endl;
822 cout << " [-------------------------------------------------------]" << endl;
823 cout << ((((output >> 7) & 0x1) == 1) ? " [ yes" : " [ no");
824 cout << ((((output >> 6) & 0x1) == 1) ? " yes" : " no");
825 cout << ((((output >> 5) & 0x1) == 1) ? " yes " : " no ");
826 cout << ((((output >> 4) & 0x1) == 1) ? " yes" : " no");
827 cout << ((((output >> 3) & 0x1) == 1) ? " yes" : " no");
828 cout << ((((output >> 2) & 0x1) == 1) ? " yes" : " no");
829 cout << ((((output >> 1) & 0x1) == 1) ? " yes " : " no ");
830 cout << ((((output >> 0) & 0x1) == 1) ? " yes ]" : " no ]");
831 cout << endl;
832
833 cout << "Regional input (low pT) bit pattern : 0x" << noshowbase << setw(8)
834 << setfill('0') << hex << regionalStruct->fInput[0]
835 << setw(0) << setfill(fillChar) << dec << endl;
836 cout << "Regional input (high pT) bit pattern : 0x" << noshowbase << setw(8)
837 << setfill('0') << hex << regionalStruct->fInput[1]
838 << setw(0) << setfill(fillChar) << dec << endl;
839 cout << "Regional input logical dump: " << endl;
840 cout << " Local | Local | low pT mu | high pT mu" << endl;
841 cout << "structure no. | board | -ve | +ve | -ve | +ve" << endl;
842 cout << "-----------------------------------------------" << endl;
843 for (int i = 15; i >= 0 ; i--)
844 {
845 cout << setw(13) << (fLocNum + 16 - i) << setw(0) << " ";
846 cout << setw(5) << (15 - i) << setw(0) << " ";
847 cout << ((((regionalStruct->fInput[0] >> (i*2+1)) & 0x1) == 1) ? "yes" : " no");
848 cout << ((((regionalStruct->fInput[0] >> (i*2+0)) & 0x1) == 1) ? " yes" : " no");
849 cout << ((((regionalStruct->fInput[1] >> (i*2+1)) & 0x1) == 1) ? " yes" : " no");
850 cout << ((((regionalStruct->fInput[1] >> (i*2+0)) & 0x1) == 1) ? " yes" : " no");
851 cout << endl;
852 }
853
854 UInt_t mask = GetRegionalMask(regionalStruct);
855 cout << "Local mask : 0x" << noshowbase << setw(4) << setfill('0') << hex
856 << mask << dec << setw(0) << setfill(fillChar) << " [Bits: ";
857 PrintBitPattern(mask >> 12, 4); cout << " ";
858 PrintBitPattern(mask >> 8, 4); cout << " ";
859 PrintBitPattern(mask >> 4, 4); cout << " ";
860 PrintBitPattern(mask >> 0, 4); cout << "]" << endl;
861
862 cout << "L0 counter : " << GetRegionalL0(regionalStruct) << endl;
863
864 cout << "========================= Regional structure scalars "
865 << setw(3) << fRegNum << setw(0)
866 << " =========================" << endl;
867 if (scalars != NULL)
868 {
869 cout << " Number of clock cycles : " << scalars->fClk << endl;
870 cout << " Number of high pT +ve single mu : " << scalars->fScaler[0] << endl;
871 cout << " Number of high pT -ve single mu : " << scalars->fScaler[1] << endl;
872 cout << "Number of high pT unlike sign pair : " << scalars->fScaler[2] << endl;
873 cout << " Number of high pT like sign pair : " << scalars->fScaler[3] << endl;
874 cout << " Number of low pT +ve single mu : " << scalars->fScaler[4] << endl;
875 cout << " Number of low pT -ve single mu : " << scalars->fScaler[5] << endl;
876 cout << " Number of low pT unlike sign pair : " << scalars->fScaler[6] << endl;
877 cout << " Number of low pT like sign pair : " << scalars->fScaler[7] << endl;
878 cout << " Hold (dead time) : " << scalars->fHold << endl;
879 }
880 else
881 {
882 cout << "(None found)" << endl;
883 }
884 }
885
886
887 void AliTriggerDecoderHandler::OnLocalStruct(
888 const AliMUONLocalInfoStruct* localStruct,
889 const AliMUONLocalScalarsStruct* scalars
890 )
891 {
892 TryDumpCorruptData(localStruct);
893
894 fLocNum++;
895 cout << "=========================== Local structure header "
896 << setw(3) << fLocNum << setw(0)
897 << " ===========================" << endl;
898 char fillChar = cout.fill(); // store current fill char to restore it.
899 cout << "L0 strip patterns:" << endl;
900 cout << "Chamber | X | Y " << endl;
901 cout << "----------------------------------------------" << endl;
902 cout << " 11 ";
903 PrintBitPattern(GetLocalX1(localStruct), 16);
904 cout << " ";
905 PrintBitPattern(GetLocalY1(localStruct), 16);
906 cout << endl;
907 cout << " 12 ";
908 PrintBitPattern(GetLocalX2(localStruct), 16);
909 cout << " ";
910 PrintBitPattern(GetLocalY2(localStruct), 16);
911 cout << endl;
912 cout << " 13 ";
913 PrintBitPattern(GetLocalX3(localStruct), 16);
914 cout << " ";
915 PrintBitPattern(GetLocalY3(localStruct), 16);
916 cout << endl;
7f62a2ea 917 cout << " 14 ";
a090ff22 918 PrintBitPattern(GetLocalX4(localStruct), 16);
919 cout << " ";
920 PrintBitPattern(GetLocalY4(localStruct), 16);
921 cout << endl;
922
923 cout << "L0 trigger bits: (word = ";
924 cout << showbase << hex << localStruct->fTriggerBits
925 << noshowbase << dec << ")" << endl;
926 cout << " ID | Dec | TrigY | YPos | Sign XDev | XDev | XPos " << endl;
927 cout << "----------------------------------------------------------------" << endl;
928 cout << "Decimal : " << setw(4) << (UInt_t)GetLocalId(localStruct) << setw(0) << " ";
929 cout << setw(4) << (UInt_t)GetLocalDec(localStruct) << setw(0) << " ";
930 cout << setw(3) << (UInt_t)GetLocalTrigY(localStruct) << setw(0) << " ";
931 cout << setw(4) << (UInt_t)GetLocalYPos(localStruct) << setw(0) << " ";
932 cout << setw(5) << (UInt_t)GetLocalSXDev(localStruct) << setw(0) << " ";
933 cout << setw(4) << (UInt_t)GetLocalXDev(localStruct) << setw(0) << " ";
934 cout << setw(4) << (UInt_t)GetLocalXPos(localStruct) << setw(0) << endl;
935 cout << " Binary : ";
936 PrintBitPattern(AliHLTUInt32_t(GetLocalId(localStruct)), 4);
937 cout << " ";
938 PrintBitPattern(AliHLTUInt32_t(GetLocalDec(localStruct)), 4);
939 cout << " ";
940 PrintBitPattern(AliHLTUInt32_t(GetLocalTrigY(localStruct)), 1);
941 cout << " ";
942 PrintBitPattern(AliHLTUInt32_t(GetLocalYPos(localStruct)), 4);
943 cout << " ";
944 PrintBitPattern(AliHLTUInt32_t(GetLocalSXDev(localStruct)), 1);
945 cout << " ";
946 PrintBitPattern(AliHLTUInt32_t(GetLocalXDev(localStruct)), 4);
947 cout << " ";
948 PrintBitPattern(AliHLTUInt32_t(GetLocalXPos(localStruct)), 5);
949 cout << endl;
950
951 cout << "L0 decision (Dec): [high pT: ";
952 PrintBitPattern((UInt_t)GetLocalHpt(localStruct), 2);
953 cout << "b (";
954 switch ((UInt_t)GetLocalHpt(localStruct))
955 {
956 case 0: cout << "No trigger"; break;
957 case 1: cout << "-ve particle"; break;
958 case 2: cout << "+ve particle"; break;
959 case 3: cout << "No deviation trigger"; break;
960 default: cout << "UNKNOWN"; break;
961 }
962 cout << "), low pT: ";
963 PrintBitPattern((UInt_t)GetLocalLpt(localStruct), 2);
964 cout << "b (";
965 switch ((UInt_t)GetLocalLpt(localStruct))
966 {
967 case 0: cout << "No trigger"; break;
968 case 1: cout << "-ve particle"; break;
969 case 2: cout << "+ve particle"; break;
970 case 3: cout << "No deviation trigger"; break;
971 default: cout << "UNKNOWN"; break;
972 }
973 cout << ")]" << endl;
974
975 cout << "=========================== Local structure scalars "
976 << setw(3) << fLocNum << setw(0)
977 << " ==========================" << endl;
978 if (scalars != NULL)
979 {
980 cout << " Number of L0 triggers : " << scalars->fL0 << endl;
981 cout << " Hold (dead time) : " << scalars->fHold << endl;
982 cout << " Number of clock cycles : " << scalars->fClk << endl;
983 cout << " Number of low pT no triggers : " << scalars->fLPtNTrig << endl;
984 cout << " Number of high pT no triggers : " << scalars->fHPtNTrig << endl;
985 cout << " Number of low pT right triggers : " << scalars->fLPtRTrig << endl;
986 cout << " Number of high pT right triggers : " << scalars->fHPtRTrig << endl;
987 cout << " Number of low pT left triggers : " << scalars->fLPtLTrig << endl;
988 cout << " Number of high pT left triggers : " << scalars->fHPtLTrig << endl;
989 cout << " Number of low pT straight triggers : " << scalars->fLPtSTrig << endl;
990 cout << "Number of high pT straight triggers : " << scalars->fHPtSTrig << endl;
991
992 UInt_t xoryflag = GetLocalComptXY(scalars);
993 if (xoryflag == 1)
994 {
995 cout << "Y strip counts:" << endl;
996 }
997 else
998 {
999 cout << "X strip counts:" << endl;
1000 }
1001 cout << " | Chamber " << endl;
1002 cout << "Strip | 11 | 12 | 13 | 14" << endl;
1003 cout << "------------------------------------------" << endl;
1004 for (int i = 0; i < 16; i++)
1005 {
1006 cout << setw(5) << i << setw(0) << " "
1007 << setw(6) << (UInt_t)GetLocalXY1(scalars, i) << setw(0) << " "
1008 << setw(6) << (UInt_t)GetLocalXY2(scalars, i) << setw(0) << " "
1009 << setw(6) << (UInt_t)GetLocalXY3(scalars, i) << setw(0) << " "
1010 << setw(6) << (UInt_t)GetLocalXY4(scalars, i) << setw(0) << endl;
1011 }
1012
1013 cout << " EOS word : 0x" << setw(8) << setfill('0')
1014 << hex << scalars->fEOS << setw(0) << setfill(fillChar) << dec << endl;
1015 cout << " [ Switches : 0x" << setw(3)
1016 << setfill('0') << hex << (UInt_t)GetLocalSwitch(scalars)
1017 << setw(0) << setfill(fillChar) << dec << " (";
1018 PrintBitPattern((UInt_t)GetLocalSwitch(scalars) >> 8, 2); cout << " ";
1019 PrintBitPattern((UInt_t)GetLocalSwitch(scalars) >> 4, 4); cout << " ";
1020 PrintBitPattern((UInt_t)GetLocalSwitch(scalars) >> 0, 4);
1021 cout << "b) ]" << endl;
1022 cout << " [ X or Y : " << xoryflag
1023 << ((xoryflag == 1) ? " (scalars for Y strips) ]" : " (scalars for X strips) ]")
1024 << endl;
1025
1026 cout << "Reset signal : " << scalars->fReset << endl;
1027 }
1028 else
1029 {
1030 cout << "(None found)" << endl;
1031 }
1032 }
1033
1034
1035 void AliTriggerDecoderHandler::OnError(ErrorCode error, const void* location)
1036 {
1037 TryDumpCorruptData(location);
1038 HandleError(
1039 ErrorCodeToMessage(error), error,
1040 ErrorCodeToString(error), location
1041 );
1042 }
887a669c 1043
1044} // end of namespace
1045
1046
a090ff22 1047int DumpRawDataHeader(
1048 const char* buffer, unsigned long bufferSize, const AliRawDataHeader* header,
1049 bool continueParse
1050 )
1051{
a8982f78 1052 // Dumps the common DDL raw data block header.
1053
a090ff22 1054 cout << "*************************** Common DDL data header *******************************" << endl;
1055 char fillChar = cout.fill(); // remember fill char to set back to original later.
1056 int result = CheckHeaderField(header->fSize, buffer, bufferSize, continueParse);
1057 if (result != EXIT_SUCCESS) return result;
1058 cout << "Size of the raw data in bytes : ";
1059 if (header->fSize != 0xFFFFFFFF)
1060 cout << header->fSize;
1061 else
1062 cout << "0xFFFFFFFF (unknown)";
1063 cout << endl;
1064
1065 result = CheckHeaderField(header->fWord2, buffer, bufferSize, continueParse);
1066 if (result != EXIT_SUCCESS) return result;
1067 cout << " Format version : " << UInt_t(header->GetVersion()) << endl;
1068 UInt_t l1msg = header->GetL1TriggerMessage();
1069 cout << " L1 trigger message : 0x"
1070 << noshowbase << hex << setfill('0') << setw(2) << l1msg
1071 << setw(0) << setfill(fillChar) << dec
1072 << " [Spare: " << ((l1msg >> 7) & 0x1)
1073 << ", ClT: " << ((l1msg >> 6) & 0x1)
1074 << ", RoC: 0x" << noshowbase << hex << ((l1msg >> 2) & 0x4) << dec
1075 << ", ESR: " << ((l1msg >> 1) & 0x1)
1076 << ", L1SwC: " << ((l1msg >> 0) & 0x1) << "]" << endl;
1077 cout << " Bunch crossing (Event ID1) : " << header->GetEventID1() << " (0x"
1078 << noshowbase << hex << setfill('0') << setw(3) << header->GetEventID1()
1079 << dec << setw(0) << setfill(fillChar) << ")" << endl;
1080
1081 result = CheckHeaderField(header->fEventID2, buffer, bufferSize, continueParse);
1082 if (result != EXIT_SUCCESS) return result;
1083 cout << " Orbit number (Event ID2) : " << header->fEventID2 << " (0x"
1084 << noshowbase << hex << setfill('0') << setw(6) << header->fEventID2
1085 << dec << setw(0) << setfill(fillChar) << ")" << endl;
1086
1087 result = CheckHeaderField(header->fAttributesSubDetectors, buffer, bufferSize, continueParse);
1088 if (result != EXIT_SUCCESS) return result;
1089 cout << " Block attributes : " << UInt_t(header->GetAttributes()) << " (0x"
1090 << noshowbase << hex << setfill('0') << setw(2) << UInt_t(header->GetAttributes())
1091 << dec << setw(0) << setfill(fillChar) << ")" << endl;
1092 cout << " Participating sub-detectors : 0x" << noshowbase << hex
1093 << setfill('0') << setw(6) << header->GetSubDetectors() << dec
1094 << setw(0) << setfill(fillChar) << " [Bits: ";
1095 for (int i = 5; i >= 0; i--)
1096 {
1097 PrintBitPattern(header->GetSubDetectors() >> (i*4), 4);
1098 if (i != 0)
1099 cout << " ";
1100 else
1101 cout << "]" << endl;
1102 }
1103
1104 result = CheckHeaderField(header->fStatusMiniEventID, buffer, bufferSize, continueParse);
1105 if (result != EXIT_SUCCESS) return result;
1106 UInt_t statusBits = header->GetStatus();
1107 cout << " Status & error bits : 0x" << noshowbase << hex
1108 << setfill('0') << setw(4) << statusBits << setw(0) << setfill(fillChar)
1109 << dec << endl;
37d9b8cd 1110 cout << " [ Original data bit : " << ((statusBits >> 15) & 0x1) << " ]" << endl;
a090ff22 1111 cout << " [ Multi-event buffer error : " << ((statusBits >> 14) & 0x1) << " ]" << endl;
1112 cout << " [ Trigger L1 missing error : " << ((statusBits >> 13) & 0x1) << " ]" << endl;
1113 cout << " [ Trigger error (other) : " << ((statusBits >> 12) & 0x1) << " ]" << endl;
1114 cout << " [ Pre-pulse error : " << ((statusBits >> 11) & 0x1) << " ]" << endl;
1115 cout << " [ Trigger L2 time violation : " << ((statusBits >> 10) & 0x1) << " ]" << endl;
1116 cout << " [ Trigger L1 time violation : " << ((statusBits >> 9) & 0x1) << " ]" << endl;
1117 cout << " [ DDG payload flag : " << ((statusBits >> 8) & 0x1) << " ]" << endl;
1118 cout << " [ HLT payload flag : " << ((statusBits >> 7) & 0x1) << " ]" << endl;
1119 cout << " [ HLT decision flag : " << ((statusBits >> 6) & 0x1) << " ]" << endl;
1120 cout << " [ FEE error : " << ((statusBits >> 5) & 0x1) << " ]" << endl;
1121 cout << " [ Trigger information unavailable : " << ((statusBits >> 4) & 0x1) << " ]" << endl;
1122 cout << " [ Control parity error : " << ((statusBits >> 3) & 0x1) << " ]" << endl;
1123 cout << " [ Data parity error : " << ((statusBits >> 2) & 0x1) << " ]" << endl;
1124 cout << " [ Trigger missing error : " << ((statusBits >> 1) & 0x1) << " ]" << endl;
1125 cout << " [ Trigger overlap error : " << ((statusBits >> 0) & 0x1) << " ]" << endl;
1126 cout << " Mini event ID : " << header->GetMiniEventID() << " (0x"
1127 << noshowbase << hex << setfill('0') << setw(3) << header->GetMiniEventID()
1128 << dec << setw(0) << setfill(fillChar) << ")" << endl;
1129
1130 result = CheckHeaderField(header->fTriggerClassLow, buffer, bufferSize, continueParse);
1131 if (result != EXIT_SUCCESS) return result;
1132 result = CheckHeaderField(header->fROILowTriggerClassHigh, buffer, bufferSize, continueParse);
1133 if (result != EXIT_SUCCESS) return result;
1134 ULong64_t triggerClasses = header->GetTriggerClasses();
1135 cout << " Trigger classes : 0x" << noshowbase << hex << setw(13)
1136 << setfill('0') << triggerClasses << setw(0) << setfill(fillChar)
1137 << dec << endl;
1138 cout << " [Bits: ";
1139 PrintBitPattern(triggerClasses >> (12*4), 2);
1140 cout << " ";
1141 for (int i = 11; i >= 0; i--)
1142 {
1143 PrintBitPattern(triggerClasses >> (i*4), 4);
1144 if (i != 0)
1145 cout << " ";
1146 else
1147 cout << "]" << endl;
1148 }
1149
1150 result = CheckHeaderField(header->fROIHigh, buffer, bufferSize, continueParse);
1151 if (result != EXIT_SUCCESS) return result;
1152 ULong64_t roiBits = header->GetROI();
1153 cout << " Region of interest : 0x" << noshowbase << hex << setw(9)
1154 << setfill('0') << roiBits << setw(0) << setfill(fillChar)
1155 << dec << endl;
1156 cout << " [Bits: ";
1157 for (int i = 8; i >= 0; i--)
1158 {
1159 PrintBitPattern(roiBits >> (i*4), 4);
1160 if (i != 0)
1161 cout << " ";
1162 else
1163 cout << "]" << endl;
1164 }
1165 cout << "**********************************************************************************" <<endl;
1166 return EXIT_SUCCESS;
1167}
1168
1169
887a669c 1170int DumpTrackerDDLRawStream(
daab896b 1171 char* buffer, unsigned long bufferSize,
1172 bool continueParse, bool tryrecover
887a669c 1173 )
1174{
a8982f78 1175 // Dumps a tracker DDL raw stream data.
1176
daab896b 1177 AliRawDataHeader* header = reinterpret_cast<AliRawDataHeader*>(buffer);
a090ff22 1178 int result = DumpRawDataHeader(buffer, bufferSize, header, continueParse);
1179 if (result != EXIT_SUCCESS) return result;
887a669c 1180
1181 // Setup the decoder for the DDL payload.
1182 AliMUONTrackerDDLDecoder<AliTrackerDecoderHandler> decoder;
1183 decoder.ExitOnError(not continueParse);
a090ff22 1184 decoder.SendDataOnParityError(true);
daab896b 1185 decoder.TryRecover(tryrecover);
887a669c 1186 decoder.AutoDetectTrailer(true);
1187 decoder.CheckForTrailer(true);
daab896b 1188 char* payload = buffer + sizeof(AliRawDataHeader);
887a669c 1189 UInt_t payloadSize = bufferSize - sizeof(AliRawDataHeader);
1190 if (decoder.Decode(payload, payloadSize))
1191 {
1192 return EXIT_SUCCESS;
1193 }
1194 else
1195 {
1196 return PARSE_ERROR;
1197 }
1198}
1199
1200
1201int DumpTriggerDDLRawStream(
1202 const char* buffer, unsigned long bufferSize,
daab896b 1203 bool continueParse, bool tryrecover
887a669c 1204 )
1205{
a8982f78 1206 // Dumps a trigger DDL raw stream data.
1207
887a669c 1208 const AliRawDataHeader* header =
1209 reinterpret_cast<const AliRawDataHeader*>(buffer);
a090ff22 1210 int result = DumpRawDataHeader(buffer, bufferSize, header, continueParse);
1211 if(result != EXIT_SUCCESS) return result;
1212 bool scalarEvent = ((header->GetL1TriggerMessage() & 0x1) == 0x1);
887a669c 1213
1214 AliMUONTriggerDDLDecoder<AliTriggerDecoderHandler> decoder;
1215 decoder.ExitOnError(not continueParse);
daab896b 1216 decoder.TryRecover(tryrecover);
887a669c 1217 decoder.AutoDetectScalars(false);
1218 const char* payload = buffer + sizeof(AliRawDataHeader);
1219 UInt_t payloadSize = bufferSize - sizeof(AliRawDataHeader);
1220 if (decoder.Decode(payload, payloadSize, scalarEvent))
1221 {
1222 return EXIT_SUCCESS;
1223 }
1224 else
1225 {
1226 return PARSE_ERROR;
1227 }
1228}
1229
1230
a68bc1e6 1231int DumpRecHitStruct(
3dad52b0 1232 const char* buffer, unsigned long bufferSize,
a68bc1e6 1233 const AliHLTMUONRecHitStruct* hit,
3dad52b0 1234 bool continueParse
1235 )
1236{
a68bc1e6 1237 // Step through the fields trying to print them.
1238 // At each step check if we have not overflowed the buffer. If we have
1239 // not, then we can print the field, otherwise we print the left over
1240 // bytes assumed to be corrupted rubbish.
a090ff22 1241
1242 int result = CheckField(hit->fFlags, buffer, bufferSize, continueParse);
1243 if (result != EXIT_SUCCESS) return result;
1244 AliHLTUInt8_t chamber = 0xFF;
1245 AliHLTUInt16_t detElemId = 0xFFFF;
1246 AliHLTMUONUtils::UnpackRecHitFlags(hit->fFlags, chamber, detElemId);
1247 if (chamber == 0 and detElemId == 0)
1248 {
1249 cout << setw(10) << left << (int)(chamber) << setw(0);
1250 cout << setw(12) << left << (int)detElemId << setw(0);
1251 }
1252 else
1253 {
1254 cout << setw(10) << left << (int)(chamber+1) << setw(0);
1255 cout << setw(12) << left << (int)detElemId << setw(0);
1256 }
1257
1258 result = CheckField(hit->fX, buffer, bufferSize, continueParse);
a68bc1e6 1259 if (result != EXIT_SUCCESS) return result;
1260 cout << setw(13) << left << hit->fX << setw(0);
1261
1262 result = CheckField(hit->fY, buffer, bufferSize, continueParse);
1263 if (result != EXIT_SUCCESS) return result;
1264 cout << setw(13) << left << hit->fY << setw(0);
1265
1266 result = CheckField(hit->fZ, buffer, bufferSize, continueParse);
1267 if (result != EXIT_SUCCESS) return result;
1268 cout << hit->fZ << setw(0) << endl;
1269
1270 return result;
3dad52b0 1271}
1272
1273
a68bc1e6 1274int DumpRecHitsBlock(
3dad52b0 1275 const char* buffer, unsigned long bufferSize,
1276 bool continueParse
1277 )
1278{
a8982f78 1279 // Dumps a reconstructed hits data block.
1280
a68bc1e6 1281 int result = EXIT_SUCCESS;
1282 AliHLTMUONRecHitsBlockReader block(buffer, bufferSize);
1283
1284 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1285 if (result != EXIT_SUCCESS and not continueParse) return result;
1286
1287 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1288
1289 // Print the data block record entries.
a090ff22 1290 cout << "Chamber | DetElemID | X (cm) | Y (cm) | Z (cm)" << endl;
1291 cout << "-----------------------------------------------------------" << endl;
a68bc1e6 1292 const AliHLTMUONRecHitStruct* entry = block.GetArray();
1293 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1294 {
1295 int subResult = DumpRecHitStruct(buffer, bufferSize, entry++, continueParse);
1296 if (subResult != EXIT_SUCCESS) return subResult;
1297 }
1298
1299 return result;
3dad52b0 1300}
1301
1302
a68bc1e6 1303int DumpTriggerRecordStruct(
dba14d7d 1304 const char* buffer, unsigned long bufferSize,
1305 const AliHLTMUONTriggerRecordStruct* record,
1306 bool continueParse
1307 )
a68bc1e6 1308{
1309 // Step through the fields trying to print them.
1310 // At each step check if we have not overflowed the buffer. If we have
1311 // not, then we can print the field, otherwise we print the left over
1312 // bytes assumed to be corrupted rubbish.
1313 int result = CheckField(record->fId, buffer, bufferSize, continueParse);
1314 if (result != EXIT_SUCCESS) return result;
1315 cout << "Trigger Record ID: " << record->fId <<endl;
1316
1317 result = CheckField(record->fFlags, buffer, bufferSize, continueParse);
1318 if (result != EXIT_SUCCESS) return result;
710ddf4b 1319 ios::fmtflags oldflags = cout.flags();
a68bc1e6 1320 cout << "Flags: " << showbase << hex << record->fFlags << dec;
710ddf4b 1321 cout.flags(oldflags);
a68bc1e6 1322
1323 // Print the individual trigger bits.
1324 AliHLTMUONParticleSign sign;
1325 bool hitset[4];
1326 AliHLTMUONUtils::UnpackTriggerRecordFlags(record->fFlags, sign, hitset);
1327 cout << " [Sign: " << sign << ", Hits set on chambers: ";
1328 bool first = true;
1329 for (AliHLTUInt32_t i = 0; i < 4; i++)
1330 {
1331 if (hitset[i])
1332 {
1333 cout << (first ? "" : ", ") << i+11;
1334 first = false;
1335 }
1336 }
1337 cout << (first ? "none]" : "]") << endl;
1338
1339 result = CheckField(record->fPx, buffer, bufferSize, continueParse);
1340 if (result != EXIT_SUCCESS) return result;
1341 cout << "Momentum: (px = " << record->fPx << ", ";
1342
1343 result = CheckField(record->fPy, buffer, bufferSize, continueParse);
1344 if (result != EXIT_SUCCESS) return result;
1345 cout << "py = " << record->fPy << ", ";
1346
1347 result = CheckField(record->fPz, buffer, bufferSize, continueParse);
1348 if (result != EXIT_SUCCESS) return result;
1349 cout << "pz = " << record->fPz << ") GeV/c"<<endl;
1350
a090ff22 1351 cout << "Hits on chambers:" << endl;
1352 cout << "Chamber | DetElemID | X (cm) | Y (cm) | Z (cm)" << endl;
1353 cout << "-----------------------------------------------------------" << endl;
a68bc1e6 1354 const AliHLTMUONRecHitStruct* hit = &record->fHit[0];
1355 for(AliHLTUInt32_t ch = 0; ch < 4; ch++)
1356 {
a68bc1e6 1357 result = DumpRecHitStruct(buffer, bufferSize, hit++, continueParse);
1358 if (result != EXIT_SUCCESS) return result;
1359 }
1360
1361 return result;
1362
1363}
1364
1365
1366int DumpTriggerRecordsBlock(
3dad52b0 1367 const char* buffer, unsigned long bufferSize,
1368 bool continueParse
1369 )
1370{
a8982f78 1371 // Dumps a trigger records data block.
1372
a68bc1e6 1373 AliHLTMUONTriggerRecordsBlockReader block(buffer, bufferSize);
1374
1375 int result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1376 if (result != EXIT_SUCCESS and not continueParse) return result;
1377
1378 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1379
1380 // Print the data block record entries.
1381 const AliHLTMUONTriggerRecordStruct* entry = block.GetArray();
1382 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1383 {
1384 cout << "============================== Trigger Record number " << i+1
1385 << " of " << nentries << " ==============================" << endl;
1386 int subResult = DumpTriggerRecordStruct(buffer, bufferSize, entry++, continueParse);
1387 if (subResult != EXIT_SUCCESS) return subResult;
1388 }
1389
3dad52b0 1390 return EXIT_SUCCESS;
1391}
1392
1393
462e3880 1394int DumpLocalStruct(
1395 const char* buffer, unsigned long bufferSize,
1396 const AliMUONLocalInfoStruct* localStruct,
1397 bool continueParse,
1398 const char* title = ""
1399 )
1400{
1401 // Prints the fields of a L0 local structure as found in the buffer.
1402
1403 typedef AliMUONTriggerDDLDecoderEventHandler AliH;
1404
1405 cout << "L0 strip patterns" << title << ":" << endl;
1406 cout << "Chamber | X | Y " << endl;
1407 cout << "----------------------------------------------" << endl;
1408 int result = CheckField(localStruct->fX2X1, buffer, bufferSize, continueParse);
1409 if (result != EXIT_SUCCESS) return result;
1410 cout << " 11 ";
1411 PrintBitPattern(AliH::GetLocalX1(localStruct), 16);
1412 result = CheckField(localStruct->fY2Y1, buffer, bufferSize, continueParse);
1413 if (result != EXIT_SUCCESS) return result;
1414 cout << " ";
1415 PrintBitPattern(AliH::GetLocalY1(localStruct), 16);
1416 cout << endl;
1417 cout << " 12 ";
1418 PrintBitPattern(AliH::GetLocalX2(localStruct), 16);
1419 cout << " ";
1420 PrintBitPattern(AliH::GetLocalY2(localStruct), 16);
1421 cout << endl;
1422
1423 result = CheckField(localStruct->fX4X3, buffer, bufferSize, continueParse);
1424 if (result != EXIT_SUCCESS) return result;
1425 cout << " 13 ";
1426 PrintBitPattern(AliH::GetLocalX3(localStruct), 16);
1427 result = CheckField(localStruct->fY4Y3, buffer, bufferSize, continueParse);
1428 if (result != EXIT_SUCCESS) return result;
1429 cout << " ";
1430 PrintBitPattern(AliH::GetLocalY3(localStruct), 16);
1431 cout << endl;
1432 cout << " 12 ";
1433 PrintBitPattern(AliH::GetLocalX4(localStruct), 16);
1434 cout << " ";
1435 PrintBitPattern(AliH::GetLocalY4(localStruct), 16);
1436 cout << endl;
1437
1438 cout << "L0 trigger bits" << title << ": (word = ";
1439 result = CheckField(localStruct->fTriggerBits, buffer, bufferSize, continueParse);
1440 if (result != EXIT_SUCCESS) return result;
1441 cout << showbase << hex << localStruct->fTriggerBits
1442 << noshowbase << dec << ")" << endl;
1443 cout << " ID | Dec | TrigY | YPos | Sign XDev | XDev | XPos " << endl;
1444 cout << "------------------------------------------------------" << endl;
1445 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalId(localStruct)), 4);
1446 cout << " ";
1447 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalDec(localStruct)), 4);
1448 cout << " ";
1449 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalTrigY(localStruct)), 1);
1450 cout << " ";
1451 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalYPos(localStruct)), 4);
1452 cout << " ";
1453 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalSXDev(localStruct)), 1);
1454 cout << " ";
1455 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalXDev(localStruct)), 4);
1456 cout << " ";
1457 PrintBitPattern(AliHLTUInt32_t(AliH::GetLocalXPos(localStruct)), 5);
1458 cout << endl;
1459
1460 return result;
1461}
1462
1463
dba14d7d 1464int DumpTrigRecInfoStruct(
1465 const char* buffer, unsigned long bufferSize,
1466 const AliHLTMUONTrigRecInfoStruct* debuginfo,
1467 bool continueParse
1468 )
a68bc1e6 1469{
1470 // Step through the fields trying to print them.
1471 // At each step check if we have not overflowed the buffer. If we have
1472 // not, then we can print the field, otherwise we print the left over
1473 // bytes assumed to be corrupted rubbish.
a090ff22 1474
a68bc1e6 1475 int result = CheckField(debuginfo->fTrigRecId, buffer, bufferSize, continueParse);
1476 if (result != EXIT_SUCCESS) return result;
a090ff22 1477 cout << "Trigger Record ID: " << debuginfo->fTrigRecId << endl;
a68bc1e6 1478
a090ff22 1479 cout << "Detector element IDs:" << endl;
1480 cout << " Chamber : 11 | 12 | 13 | 14" << endl;
1481 cout << " ID : ";
1482 for (int i = 0; i < 4; i++)
1483 {
1484 result = CheckField(debuginfo->fDetElemId[i], buffer, bufferSize, continueParse);
1485 if (result != EXIT_SUCCESS) return result;
710ddf4b 1486 ios::fmtflags oldflags = cout.flags();
a090ff22 1487 cout << setw(12) << right << debuginfo->fDetElemId[i] << setw(0);
710ddf4b 1488 cout.flags(oldflags);
a090ff22 1489 if (i != 3) cout << " | ";
1490 }
1491 cout << endl;
1492
1493 cout << "Momentum estimation parameters:" << endl;
1494 cout << " Parameter : Z_middle coordinate (cm) | Integrated magnetic field (T.m)" << endl;
1495 cout << " Value : ";
a68bc1e6 1496 result = CheckField(debuginfo->fZmiddle, buffer, bufferSize, continueParse);
1497 if(result != EXIT_SUCCESS) return result;
a090ff22 1498 cout << setw(24) << right << debuginfo->fZmiddle << setw(0) << " | ";
1499
a68bc1e6 1500 result = CheckField(debuginfo->fBl, buffer, bufferSize, continueParse);
1501 if (result != EXIT_SUCCESS) return result;
a090ff22 1502 cout << setw(31) << right << debuginfo->fBl << setw(0) << endl;
1503
462e3880 1504 result = DumpLocalStruct(buffer, bufferSize, &debuginfo->fL0Struct, continueParse, " for central local structure");
a090ff22 1505 if (result != EXIT_SUCCESS) return result;
462e3880 1506 result = DumpLocalStruct(buffer, bufferSize, &debuginfo->fL0StructPrev, continueParse, " for previous local structure");
a090ff22 1507 if (result != EXIT_SUCCESS) return result;
462e3880 1508 result = DumpLocalStruct(buffer, bufferSize, &debuginfo->fL0StructNext, continueParse, " for next local structure");
a090ff22 1509 if (result != EXIT_SUCCESS) return result;
a090ff22 1510
a68bc1e6 1511 return result;
1512}
1513
1514
1515int DumpTrigRecsDebugBlock(
3dad52b0 1516 const char* buffer, unsigned long bufferSize,
3dad52b0 1517 bool continueParse
1518 )
a68bc1e6 1519{
a8982f78 1520 // Dumps the debugging information for trigger records.
1521
a68bc1e6 1522 AliHLTMUONTrigRecsDebugBlockReader block(buffer, bufferSize);
1523
1524 int result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1525 if (result != EXIT_SUCCESS and not continueParse) return result;
1526
1527 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1528
1529 // Print the data block record entries.
a68bc1e6 1530 const AliHLTMUONTrigRecInfoStruct* entry = block.GetArray();
1531 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1532 {
3d3c03a0 1533 cout << "======================= Trigger Record debug data " << i+1
1534 << " of " << nentries << " =======================" << endl;
a68bc1e6 1535 int subResult = DumpTrigRecInfoStruct(buffer, bufferSize, entry++, continueParse);
1536 if (subResult != EXIT_SUCCESS) return subResult;
1537 }
1538
1539 return EXIT_SUCCESS;
1540}
1541
1542
a68bc1e6 1543int DumpClusterStruct(
dba14d7d 1544 const char* buffer, unsigned long bufferSize,
1545 const AliHLTMUONClusterStruct* cluster,
1546 bool continueParse
1547 )
a68bc1e6 1548{
1549 // Step through the fields trying to print them.
1550 // At each step check if we have not overflowed the buffer. If we have
1551 // not, then we can print the field, otherwise we print the left over
1552 // bytes assumed to be corrupted rubbish.
1553 int result = CheckField(cluster->fId, buffer, bufferSize, continueParse);
1554 if (result != EXIT_SUCCESS) return result;
66622a82 1555 cout << " Cluster ID: " << cluster->fId << endl;
a68bc1e6 1556
1557 result = CheckField(cluster->fDetElemId, buffer, bufferSize, continueParse);
1558 if (result != EXIT_SUCCESS) return result;
66622a82 1559 cout << " Detector Element ID: " << cluster->fDetElemId << endl;
a68bc1e6 1560
66622a82 1561 result = CheckField(cluster->fNchannelsB, buffer, bufferSize, continueParse);
83d66053 1562 if (result != EXIT_SUCCESS) return result;
66622a82 1563 cout << " Number of channels in bending plane: " << cluster->fNchannelsB << endl;
1564
1565 result = CheckField(cluster->fNchannelsNB, buffer, bufferSize, continueParse);
1566 if (result != EXIT_SUCCESS) return result;
1567 cout << "Number of channels in non-bending plane: " << cluster->fNchannelsNB << endl;
83d66053 1568
66622a82 1569 result = CheckField(cluster->fChargeB, buffer, bufferSize, continueParse);
1570 if (result != EXIT_SUCCESS) return result;
1571 cout << " Charge on bending plane: " << cluster->fChargeB << endl;
1572
1573 result = CheckField(cluster->fChargeNB, buffer, bufferSize, continueParse);
83d66053 1574 if (result != EXIT_SUCCESS) return result;
66622a82 1575 cout << " Charge on non bending plane: " << cluster->fChargeNB << endl;
a68bc1e6 1576
a090ff22 1577 cout << "Corresponding Hit: "<< endl;
1578 cout << "Chamber | DetElemID | X (cm) | Y (cm) | Z (cm)" << endl;
1579 cout << "-----------------------------------------------------------" << endl;
1580 result = DumpRecHitStruct(buffer, bufferSize, &cluster->fHit, continueParse);
a68bc1e6 1581
3dad52b0 1582 return result;
1583}
1584
1585
1586int DumpClustersBlock(
1587 const char* buffer, unsigned long bufferSize,
1588 bool continueParse
1589 )
1590{
a8982f78 1591 // Dumps a clusters block structure.
1592
a68bc1e6 1593 int result = EXIT_SUCCESS;
3dad52b0 1594 AliHLTMUONClustersBlockReader block(buffer, bufferSize);
a68bc1e6 1595
1596 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1597 if (result != EXIT_SUCCESS and not continueParse) return result;
1598
1599 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1600
1601 // Print the data block record entries.
1602 const AliHLTMUONClusterStruct* entry = block.GetArray();
1603 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1604 {
3d3c03a0 1605 cout << "======================= Cluster Number "
1606 << i+1 << " =======================" << endl;
a68bc1e6 1607 int subResult = DumpClusterStruct(buffer, bufferSize, entry++, continueParse);
1608 if (subResult != EXIT_SUCCESS) return subResult;
1609 }
1610
1611 return result;
1612}
1613
1614
1615int DumpChannelStruct(
dba14d7d 1616 const char* buffer, unsigned long bufferSize,
1617 const AliHLTMUONChannelStruct* channel,
1618 bool continueParse
1619 )
a68bc1e6 1620{
1621 // Step through the fields trying to print them.
1622 // At each step check if we have not overflowed the buffer. If we have
1623 // not, then we can print the field, otherwise we print the left over
1624 // bytes assumed to be corrupted rubbish.
1625 int result = CheckField(channel->fClusterId, buffer, bufferSize, continueParse);
1626 if (result != EXIT_SUCCESS) return result;
710ddf4b 1627 ios::fmtflags oldflags = cout.flags();
a68bc1e6 1628 cout << setw(16) << left << channel->fClusterId << setw(0);
710ddf4b 1629 cout.flags(oldflags);
a68bc1e6 1630
dba14d7d 1631 result = CheckField(channel->fBusPatch, buffer, bufferSize, continueParse);
1632 if (result != EXIT_SUCCESS) return result;
710ddf4b 1633 oldflags = cout.flags();
dba14d7d 1634 cout << setw(16) << left << channel->fBusPatch << setw(0);
710ddf4b 1635 cout.flags(oldflags);
dba14d7d 1636
a68bc1e6 1637 result = CheckField(channel->fManu, buffer, bufferSize, continueParse);
1638 if (result != EXIT_SUCCESS) return result;
710ddf4b 1639 oldflags = cout.flags();
a68bc1e6 1640 cout << setw(16) << left << channel->fManu << setw(0);
710ddf4b 1641 cout.flags(oldflags);
a68bc1e6 1642
1643 result = CheckField(channel->fChannelAddress, buffer, bufferSize, continueParse);
1644 if (result != EXIT_SUCCESS) return result;
710ddf4b 1645 oldflags = cout.flags();
dba14d7d 1646 cout << setw(16) << left << channel->fChannelAddress << setw(0);
710ddf4b 1647 cout.flags(oldflags);
a68bc1e6 1648
1649 result = CheckField(channel->fSignal, buffer, bufferSize, continueParse);
1650 if(result != EXIT_SUCCESS) return result;
710ddf4b 1651 oldflags = cout.flags();
a68bc1e6 1652 cout << setw(16) << left << channel->fSignal << setw(0);
710ddf4b 1653 cout.flags(oldflags);
a68bc1e6 1654
1655 result = CheckField(channel->fRawDataWord, buffer, bufferSize, continueParse);
1656 if(result != EXIT_SUCCESS) return result;
710ddf4b 1657 oldflags = cout.flags();
1658 cout << showbase << hex << channel->fRawDataWord << dec << setw(0) << endl;
1659 cout.flags(oldflags);
a68bc1e6 1660
1661 return result;
3dad52b0 1662}
1663
1664
1665int DumpChannelsBlock(
1666 const char* buffer, unsigned long bufferSize,
1667 bool continueParse
1668 )
1669{
a8982f78 1670 // Dumps a channels block structure.
1671
a68bc1e6 1672 int result = EXIT_SUCCESS;
3dad52b0 1673 AliHLTMUONChannelsBlockReader block(buffer, bufferSize);
a68bc1e6 1674
1675 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1676 if (result != EXIT_SUCCESS and not continueParse) return result;
1677
1678 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1679
1680 // Print the data block record entries.
a090ff22 1681 cout << "Cluster ID | Bus Patch | Manu Address | Channel Addr | Signal Value | Raw Data Word " << endl;
dba14d7d 1682 cout << "----------------------------------------------------------------------------------------------" << endl;
a68bc1e6 1683 const AliHLTMUONChannelStruct* entry = block.GetArray();
1684 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1685 {
1686 int subResult = DumpChannelStruct(buffer, bufferSize, entry++, continueParse);
1687 if (subResult != EXIT_SUCCESS) return subResult;
1688 }
3dad52b0 1689 return EXIT_SUCCESS;
1690}
1691
1692
1693int DumpMansoTrackStruct(
1694 const char* buffer, unsigned long bufferSize,
1695 const AliHLTMUONMansoTrackStruct* track,
1696 bool continueParse
1697 )
1698{
1699 // Step through the fields trying to print them.
1700 // At each step check if we have not overflowed the buffer. If we have
1701 // not, then we can print the field, otherwise we print the left over
1702 // bytes assumed to be corrupted rubbish.
1703 int result = CheckField(track->fId, buffer, bufferSize, continueParse);
1704 if (result != EXIT_SUCCESS) return result;
1705 cout << "Track ID: " << track->fId << "\t";
1706
1707 result = CheckField(track->fTrigRec, buffer, bufferSize, continueParse);
1708 if (result != EXIT_SUCCESS) return result;
1709 cout << "Trigger Record ID: " << track->fTrigRec << endl;
1710
1711 result = CheckField(track->fFlags, buffer, bufferSize, continueParse);
1712 if (result != EXIT_SUCCESS) return result;
710ddf4b 1713 ios::fmtflags oldflags = cout.flags();
3dad52b0 1714 cout << "Flags: " << showbase << hex << track->fFlags << dec;
710ddf4b 1715 cout.flags(oldflags);
a68bc1e6 1716
3dad52b0 1717 // Print the individual trigger bits.
1718 AliHLTMUONParticleSign sign;
1719 bool hitset[4];
1720 AliHLTMUONUtils::UnpackMansoTrackFlags(track->fFlags, sign, hitset);
1721 cout << " [Sign: " << sign << ", Hits set on chambers: ";
1722 bool first = true;
1723 for (AliHLTUInt32_t i = 0; i < 4; i++)
1724 {
1725 if (hitset[i])
1726 {
1727 cout << (first ? "" : ", ") << i+7;
1728 first = false;
1729 }
1730 }
1731 cout << (first ? "none]" : "]") << endl;
1732
1733 result = CheckField(track->fPx, buffer, bufferSize, continueParse);
1734 if (result != EXIT_SUCCESS) return result;
1735 cout << "Momentum: (px = " << track->fPx << ", ";
1736
1737 result = CheckField(track->fPy, buffer, bufferSize, continueParse);
1738 if (result != EXIT_SUCCESS) return result;
1739 cout << "py = " << track->fPy << ", ";
1740
1741 result = CheckField(track->fPz, buffer, bufferSize, continueParse);
1742 if (result != EXIT_SUCCESS) return result;
1743 cout << "pz = " << track->fPz << ") GeV/c\t";
1744
1745 result = CheckField(track->fChi2, buffer, bufferSize, continueParse);
1746 if (result != EXIT_SUCCESS) return result;
1747 cout << "Chi squared fit: " << track->fChi2 << endl;
1748
1749 cout << "Track hits:" << endl;
a090ff22 1750 cout << "Chamber | DetElemID | X (cm) | Y (cm) | Z (cm)" << endl;
1751 cout << "-----------------------------------------------------------" << endl;
3dad52b0 1752 const AliHLTMUONRecHitStruct* hit = &track->fHit[0];
1753 for(AliHLTUInt32_t ch = 0; ch < 4; ch++)
1754 {
3dad52b0 1755 result = DumpRecHitStruct(buffer, bufferSize, hit++, continueParse);
1756 if (result != EXIT_SUCCESS) return result;
1757 }
1758
1759 return result;
1760}
1761
1762
1763int DumpMansoTracksBlock(
1764 const char* buffer, unsigned long bufferSize,
1765 bool continueParse
1766 )
1767{
a8982f78 1768 // Dumps the Manso tracks block structure.
1769
3dad52b0 1770 int result = EXIT_SUCCESS;
1771 AliHLTMUONMansoTracksBlockReader block(buffer, bufferSize);
1772
1773 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1774 if (result != EXIT_SUCCESS and not continueParse) return result;
1775
1776 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1777
1778 // Print the data block record entries.
1779 const AliHLTMUONMansoTrackStruct* entry = block.GetArray();
1780 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1781 {
1782 cout << "============================== Manso track number " << i+1
1783 << " of " << nentries << " ==============================" << endl;
1784 int subResult = DumpMansoTrackStruct(buffer, bufferSize, entry++, continueParse);
1785 if (subResult != EXIT_SUCCESS) return subResult;
1786 }
1787
1788 return result;
1789}
1790
1791
1792int DumpMansoRoIStruct(
1793 const char* buffer, unsigned long bufferSize,
1794 const AliHLTMUONMansoRoIStruct* roi,
1795 bool continueParse
1796 )
1797{
1798 // Step through the fields trying to print them.
1799 // At each step check if we have not overflowed the buffer. If we have
1800 // not, then we can print the field, otherwise we print the left over
1801 // bytes assumed to be corrupted rubbish.
1802 int result = CheckField(roi->fX, buffer, bufferSize, continueParse);
1803 if (result != EXIT_SUCCESS) return result;
710ddf4b 1804 ios::fmtflags oldflags = cout.flags();
3dad52b0 1805 cout << setw(13) << left << roi->fX << setw(0);
710ddf4b 1806 cout.flags(oldflags);
3dad52b0 1807
1808 result = CheckField(roi->fY, buffer, bufferSize, continueParse);
1809 if (result != EXIT_SUCCESS) return result;
710ddf4b 1810 oldflags = cout.flags();
3dad52b0 1811 cout << setw(13) << left << roi->fY << setw(0);
710ddf4b 1812 cout.flags(oldflags);
3dad52b0 1813
1814 result = CheckField(roi->fZ, buffer, bufferSize, continueParse);
1815 if (result != EXIT_SUCCESS) return result;
710ddf4b 1816 oldflags = cout.flags();
3dad52b0 1817 cout << setw(13) << left << roi->fZ << setw(0);
710ddf4b 1818 cout.flags(oldflags);
3dad52b0 1819
1820 result = CheckField(roi->fRadius, buffer, bufferSize, continueParse);
1821 if (result != EXIT_SUCCESS) return result;
710ddf4b 1822 oldflags = cout.flags();
3dad52b0 1823 cout << roi->fRadius << setw(0) << endl;
710ddf4b 1824 cout.flags(oldflags);
3dad52b0 1825
1826 return result;
1827}
1828
1829
1830int DumpMansoCandidateStruct(
1831 const char* buffer, unsigned long bufferSize,
1832 const AliHLTMUONMansoCandidateStruct* candidate,
1833 bool continueParse
1834 )
1835{
a8982f78 1836 // Dumps the manso candidate structure.
1837
3dad52b0 1838 int result = DumpMansoTrackStruct(buffer, bufferSize, &candidate->fTrack, continueParse);
1839 if (result != EXIT_SUCCESS) return result;
1840
1841 cout << "Regions of interest:" << endl;
1842 cout << "Chamber | X (cm) | Y (cm) | Z (cm) | Radius (cm)" << endl;
1843 cout << "-------------------------------------------------------------" << endl;
1844 const AliHLTMUONMansoRoIStruct* roi = &candidate->fRoI[0];
1845 for(AliHLTUInt32_t ch = 0; ch < 4; ch++)
1846 {
1847 cout << setw(10) << ch + 7;
1848 result = DumpMansoRoIStruct(buffer, bufferSize, roi++, continueParse);
1849 if (result != EXIT_SUCCESS) return result;
1850 }
a090ff22 1851
1852 cout << "Momentum estimation parameters:" << endl;
1853 cout << " Parameter : Z_middle coordinate (cm) | Integrated magnetic field (T.m)" << endl;
1854 cout << " Value : ";
1855 result = CheckField(candidate->fZmiddle, buffer, bufferSize, continueParse);
1856 if(result != EXIT_SUCCESS) return result;
710ddf4b 1857 ios::fmtflags oldflags = cout.flags();
a090ff22 1858 cout << setw(24) << right << candidate->fZmiddle << setw(0) << " | ";
710ddf4b 1859 cout.flags(oldflags);
a090ff22 1860
1861 result = CheckField(candidate->fBl, buffer, bufferSize, continueParse);
1862 if (result != EXIT_SUCCESS) return result;
710ddf4b 1863 oldflags = cout.flags();
a090ff22 1864 cout << setw(31) << right << candidate->fBl << setw(0) << endl;
710ddf4b 1865 cout.flags(oldflags);
a090ff22 1866
3dad52b0 1867 return result;
1868}
1869
1870
1871int DumpMansoCandidatesBlock(
1872 const char* buffer, unsigned long bufferSize,
1873 bool continueParse
1874 )
1875{
a8982f78 1876 // Dumps the manso candidates block structure.
1877
3dad52b0 1878 int result = EXIT_SUCCESS;
1879 AliHLTMUONMansoCandidatesBlockReader block(buffer, bufferSize);
1880
1881 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
1882 if (result != EXIT_SUCCESS and not continueParse) return result;
1883
1884 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
1885
1886 // Print the data block record entries.
1887 const AliHLTMUONMansoCandidateStruct* entry = block.GetArray();
1888 for(AliHLTUInt32_t i = 0; i < nentries; i++)
1889 {
1890 cout << "=========================== Manso track candidate number " << i+1
1891 << " of " << nentries << " ===========================" << endl;
1892 int subResult = DumpMansoCandidateStruct(buffer, bufferSize, entry++, continueParse);
1893 if (subResult != EXIT_SUCCESS) return subResult;
1894 }
1895
1896 return result;
1897}
1898
1899
a8982f78 1900int DumpTrackStruct(
1901 const char* buffer, unsigned long bufferSize,
1902 const AliHLTMUONTrackStruct* track,
1903 bool continueParse
1904 )
1905{
1906 // Step through the fields trying to print them.
1907 // At each step check if we have not overflowed the buffer. If we have
1908 // not, then we can print the field, otherwise we print the left over
1909 // bytes assumed to be corrupted rubbish.
1910 int result = CheckField(track->fId, buffer, bufferSize, continueParse);
1911 if (result != EXIT_SUCCESS) return result;
1912 cout << "Track ID: " << track->fId << "\t";
1913
1914 result = CheckField(track->fTrigRec, buffer, bufferSize, continueParse);
1915 if (result != EXIT_SUCCESS) return result;
1916 cout << "Trigger Record ID: " << track->fTrigRec << endl;
1917
1918 result = CheckField(track->fFlags, buffer, bufferSize, continueParse);
1919 if (result != EXIT_SUCCESS) return result;
710ddf4b 1920 ios::fmtflags oldflags = cout.flags();
a8982f78 1921 cout << "Flags: " << showbase << hex << track->fFlags << dec;
710ddf4b 1922 cout.flags(oldflags);
a8982f78 1923
1924 // Print the individual trigger bits.
1925 AliHLTMUONParticleSign sign;
1926 bool hitset[16];
1927 AliHLTMUONUtils::UnpackTrackFlags(track->fFlags, sign, hitset);
1928 cout << " [Sign: " << sign << ", Hits set: ";
1929 bool first = true;
1930 for (AliHLTUInt32_t i = 0; i < 16; i++)
1931 {
1932 if (hitset[i])
1933 {
1934 cout << (first ? "" : ", ") << i;
1935 first = false;
1936 }
1937 }
1938 cout << (first ? "none]" : "]") << endl;
1939
1940 result = CheckField(track->fPx, buffer, bufferSize, continueParse);
1941 if (result != EXIT_SUCCESS) return result;
1942 cout << "Momentum: (px = " << track->fPx << ", ";
1943
1944 result = CheckField(track->fPy, buffer, bufferSize, continueParse);
1945 if (result != EXIT_SUCCESS) return result;
1946 cout << "py = " << track->fPy << ", ";
1947
1948 result = CheckField(track->fPz, buffer, bufferSize, continueParse);
1949 if (result != EXIT_SUCCESS) return result;
1950 cout << "pz = " << track->fPz << ") GeV/c" << endl;
1951
1952 result = CheckField(track->fInverseBendingMomentum, buffer, bufferSize, continueParse);
1953 if (result != EXIT_SUCCESS) return result;
1954 cout << "Inverse bending momentum: " << track->fInverseBendingMomentum << " c/GeV" << endl;
1955
1956 result = CheckField(track->fThetaX, buffer, bufferSize, continueParse);
1957 if (result != EXIT_SUCCESS) return result;
1958 cout << "Slope: (non-bending plane = " << track->fThetaX << ", ";
1959
1960 result = CheckField(track->fThetaY, buffer, bufferSize, continueParse);
1961 if (result != EXIT_SUCCESS) return result;
1962 cout << "bending plane = " << track->fThetaY << ")" << endl;
1963
1964 result = CheckField(track->fX, buffer, bufferSize, continueParse);
1965 if (result != EXIT_SUCCESS) return result;
1966 cout << "DCA vertex: (x = " << track->fX << ", ";
1967
1968 result = CheckField(track->fY, buffer, bufferSize, continueParse);
1969 if (result != EXIT_SUCCESS) return result;
1970 cout << "y = " << track->fY << ", ";
1971
1972 result = CheckField(track->fZ, buffer, bufferSize, continueParse);
1973 if (result != EXIT_SUCCESS) return result;
1974 cout << "z = " << track->fZ << ") cm" << endl;
1975
1976 result = CheckField(track->fChi2, buffer, bufferSize, continueParse);
1977 if (result != EXIT_SUCCESS) return result;
1978 cout << "Chi squared fit: " << track->fChi2 << endl;
1979
1980 cout << "Track hits:" << endl;
1981 cout << "Chamber | DetElemID | X (cm) | Y (cm) | Z (cm)" << endl;
1982 cout << "-----------------------------------------------------------" << endl;
1983 const AliHLTMUONRecHitStruct* hit = &track->fHit[0];
1984 for(AliHLTUInt32_t ch = 0; ch < 16; ch++)
1985 {
1986 result = DumpRecHitStruct(buffer, bufferSize, hit++, continueParse);
1987 if (result != EXIT_SUCCESS) return result;
1988 }
1989
1990 return result;
1991}
1992
1993
1994int DumpTracksBlock(
1995 const char* buffer, unsigned long bufferSize,
1996 bool continueParse
1997 )
1998{
1999 // Dumps the full tracks block structure.
2000
2001 int result = EXIT_SUCCESS;
2002 AliHLTMUONTracksBlockReader block(buffer, bufferSize);
2003
2004 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
2005 if (result != EXIT_SUCCESS and not continueParse) return result;
2006
2007 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
2008
2009 // Print the data block record entries.
2010 const AliHLTMUONTrackStruct* entry = block.GetArray();
2011 for(AliHLTUInt32_t i = 0; i < nentries; i++)
2012 {
2013 cout << "================================ track number " << i+1
2014 << " of " << nentries << " ================================" << endl;
2015 int subResult = DumpTrackStruct(buffer, bufferSize, entry++, continueParse);
2016 if (subResult != EXIT_SUCCESS) return subResult;
2017 }
2018
2019 return result;
2020}
2021
2022
3dad52b0 2023int DumpSinglesDecisionBlockHeader(
2024 const char* buffer, unsigned long bufferSize,
2025 const AliHLTMUONSinglesDecisionBlockStruct* header,
2026 bool continueParse
2027 )
2028{
2029 // Step through the header fields trying to print them.
2030 // At each step check if we have not overflowed the buffer, if we have
2031 // not then we can print the field, otherwise we print the left over
2032 // bytes assumed to be corrupted rubbish.
2033 int result = CheckHeaderField(header->fNlowPt, buffer, bufferSize, continueParse);
2034 if (result != EXIT_SUCCESS) return result;
2035 cout << " Number of low pt triggers: " << header->fNlowPt << endl;
2036
2037 result = CheckHeaderField(header->fNhighPt, buffer, bufferSize, continueParse);
2038 if (result != EXIT_SUCCESS) return result;
2039 cout << "Number of high pt triggers: " << header->fNhighPt << endl;
2040
2041 return result;
2042}
2043
2044
2045int DumpTrackDecisionStruct(
2046 const char* buffer, unsigned long bufferSize,
2047 const AliHLTMUONTrackDecisionStruct* decision,
2048 bool continueParse
2049 )
2050{
2051 // Step through the fields trying to print them.
2052 // At each step check if we have not overflowed the buffer. If we have
2053 // not, then we can print the field, otherwise we print the left over
2054 // bytes assumed to be corrupted rubbish.
2055 int result = CheckField(decision->fTrackId, buffer, bufferSize, continueParse);
2056 if (result != EXIT_SUCCESS) return result;
710ddf4b 2057 ios::fmtflags oldflags = cout.flags();
3dad52b0 2058 cout << setw(13) << left << decision->fTrackId << setw(0);
710ddf4b 2059 cout.flags(oldflags);
3dad52b0 2060
2061 result = CheckField(decision->fTriggerBits, buffer, bufferSize, continueParse);
2062 if (result != EXIT_SUCCESS) return result;
710ddf4b 2063 oldflags = cout.flags();
3dad52b0 2064 cout << setw(12) << left << showbase << hex << decision->fTriggerBits
2065 << setw(0) << dec;
2066
2067 // Print the individual trigger bits.
2068 bool highPt, lowPt;
2069 AliHLTMUONUtils::UnpackTrackDecisionBits(decision->fTriggerBits, highPt, lowPt);
2070 cout << setw(7) << left << (highPt ? "yes" : "no");
2071 cout << setw(8) << left << (lowPt ? "yes" : "no");
710ddf4b 2072 cout.flags(oldflags);
c9537879 2073
2074 result = CheckField(decision->fPt, buffer, bufferSize, continueParse);
2075 if (result != EXIT_SUCCESS) return result;
2076 cout << setw(0) << decision->fPt << endl;
3dad52b0 2077
2078 return result;
2079}
2080
2081
2082int DumpSinglesDecisionBlock(
2083 const char* buffer, unsigned long bufferSize,
2084 bool continueParse
2085 )
2086{
2087 int result = EXIT_SUCCESS;
2088 AliHLTMUONSinglesDecisionBlockReader block(buffer, bufferSize);
2089
2090 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
2091 if (result != EXIT_SUCCESS and not continueParse) return result;
2092
2093 // Dump the rest of the block header.
2094 const AliHLTMUONSinglesDecisionBlockStruct* header = &block.BlockHeader();
2095 int subResult = DumpSinglesDecisionBlockHeader(buffer, bufferSize, header, continueParse);
2096 if (subResult != EXIT_SUCCESS) return subResult;
2097
2098 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
2099
2100 // Print the data block record entries.
c9537879 2101 cout << " | Trigger Bits |" << endl;
2102 cout << "Track ID | Raw HighPt LowPt | pT" << endl;
2103 cout << "----------------------------------------------------" << endl;
3dad52b0 2104 const AliHLTMUONTrackDecisionStruct* entry = block.GetArray();
2105 for(AliHLTUInt32_t i = 0; i < nentries; i++)
2106 {
2107 subResult = DumpTrackDecisionStruct(buffer, bufferSize, entry++, continueParse);
2108 if (subResult != EXIT_SUCCESS) return subResult;
2109 }
2110
2111 return result;
2112}
2113
2114
2115int DumpPairsDecisionBlockHeader(
2116 const char* buffer, unsigned long bufferSize,
2117 const AliHLTMUONPairsDecisionBlockStruct* header,
2118 bool continueParse
2119 )
2120{
2121 // Step through the header fields trying to print them.
2122 // At each step check if we have not overflowed the buffer, if we have
2123 // not then we can print the field, otherwise we print the left over
2124 // bytes assumed to be corrupted rubbish.
2125 int result = CheckHeaderField(header->fNunlikeAnyPt, buffer, bufferSize, continueParse);
2126 if (result != EXIT_SUCCESS) return result;
2127 cout << " Number of unlike all pt triggers: " << header->fNunlikeAnyPt << endl;
2128
2129 result = CheckHeaderField(header->fNunlikeLowPt, buffer, bufferSize, continueParse);
2130 if (result != EXIT_SUCCESS) return result;
2131 cout << " Number of unlike low pt triggers: " << header->fNunlikeLowPt << endl;
2132
2133 result = CheckHeaderField(header->fNunlikeHighPt, buffer, bufferSize, continueParse);
2134 if (result != EXIT_SUCCESS) return result;
2135 cout << " Number of unlike high pt triggers: " << header->fNunlikeHighPt << endl;
2136
2137 result = CheckHeaderField(header->fNlikeAnyPt, buffer, bufferSize, continueParse);
2138 if (result != EXIT_SUCCESS) return result;
2139 cout << " Number of like any pt triggers: " << header->fNlikeAnyPt << endl;
2140
2141 result = CheckHeaderField(header->fNlikeLowPt, buffer, bufferSize, continueParse);
2142 if (result != EXIT_SUCCESS) return result;
2143 cout << " Number of like low pt triggers: " << header->fNlikeLowPt << endl;
2144
2145 result = CheckHeaderField(header->fNlikeHighPt, buffer, bufferSize, continueParse);
2146 if (result != EXIT_SUCCESS) return result;
2147 cout << " Number of like high pt triggers: " << header->fNlikeHighPt << endl;
2148
2149 result = CheckHeaderField(header->fNmassAny, buffer, bufferSize, continueParse);
2150 if (result != EXIT_SUCCESS) return result;
2151 cout << " Number of all invariant mass triggers: " << header->fNmassAny << endl;
2152
2153 result = CheckHeaderField(header->fNmassLow, buffer, bufferSize, continueParse);
2154 if (result != EXIT_SUCCESS) return result;
2155 cout << " Number of low invariant mass triggers: " << header->fNmassLow << endl;
2156
2157 result = CheckHeaderField(header->fNmassHigh, buffer, bufferSize, continueParse);
2158 if (result != EXIT_SUCCESS) return result;
2159 cout << "Number of high invariant mass triggers: " << header->fNmassHigh << endl;
2160
2161 return result;
2162}
2163
2164
2165int DumpPairDecisionStruct(
2166 const char* buffer, unsigned long bufferSize,
2167 const AliHLTMUONPairDecisionStruct* decision,
2168 bool continueParse
2169 )
2170{
2171 // Step through the fields trying to print them.
2172 // At each step check if we have not overflowed the buffer. If we have
2173 // not, then we can print the field, otherwise we print the left over
2174 // bytes assumed to be corrupted rubbish.
2175 int result = CheckField(decision->fTrackAId, buffer, bufferSize, continueParse);
2176 if (result != EXIT_SUCCESS) return result;
710ddf4b 2177 ios::fmtflags oldflags = cout.flags();
3dad52b0 2178 cout << setw(13) << left << decision->fTrackAId << setw(0);
710ddf4b 2179 cout.flags(oldflags);
3dad52b0 2180
2181 result = CheckField(decision->fTrackBId, buffer, bufferSize, continueParse);
2182 if (result != EXIT_SUCCESS) return result;
710ddf4b 2183 oldflags = cout.flags();
3dad52b0 2184 cout << setw(13) << left << decision->fTrackBId << setw(0);
710ddf4b 2185 cout.flags(oldflags);
3dad52b0 2186
2187 result = CheckField(decision->fTriggerBits, buffer, bufferSize, continueParse);
2188 if (result != EXIT_SUCCESS) return result;
710ddf4b 2189 oldflags = cout.flags();
3dad52b0 2190 cout << setw(12) << left << showbase << hex << decision->fTriggerBits
2191 << setw(0) << dec;
710ddf4b 2192 cout.flags(oldflags);
3dad52b0 2193
2194 // Print the individual trigger bits.
2195 bool highMass, lowMass, unlike;
2196 AliHLTUInt8_t highPtCount, lowPtCount;
2197 AliHLTMUONUtils::UnpackPairDecisionBits(
2198 decision->fTriggerBits,
2199 highMass, lowMass, unlike, highPtCount, lowPtCount
2200 );
710ddf4b 2201 oldflags = cout.flags();
3dad52b0 2202 cout << setw(7) << left << (highMass ? "yes" : "no");
2203 cout << setw(7) << left << (lowMass ? "yes" : "no");
2204 cout << setw(7) << left << (unlike ? "yes" : "no");
2205 cout << setw(6) << left << AliHLTUInt16_t(highPtCount);
2206 cout << setw(8) << left << AliHLTUInt16_t(lowPtCount);
2207 cout << setw(0);
710ddf4b 2208 cout.flags(oldflags);
3dad52b0 2209
2210 result = CheckField(decision->fInvMass, buffer, bufferSize, continueParse);
2211 if (result != EXIT_SUCCESS) return result;
2212 cout << decision->fInvMass << endl;
2213
2214 return EXIT_SUCCESS;
2215}
2216
2217
2218int DumpPairsDecisionBlock(
2219 const char* buffer, unsigned long bufferSize,
2220 bool continueParse
2221 )
2222{
2223 int result = EXIT_SUCCESS;
2224 AliHLTMUONPairsDecisionBlockReader block(buffer, bufferSize);
2225
2226 result = CheckCommonHeader(block, buffer, bufferSize, continueParse);
2227 if (result != EXIT_SUCCESS and not continueParse) return result;
2228
2229 // Dump the rest of the block header.
2230 const AliHLTMUONPairsDecisionBlockStruct* header = &block.BlockHeader();
2231 int subResult = DumpPairsDecisionBlockHeader(buffer, bufferSize, header, continueParse);
2232 if (subResult != EXIT_SUCCESS) return subResult;
2233
2234 AliHLTUInt32_t nentries = CalculateNEntries(block, bufferSize);
2235
2236 // Print the data block record entries.
2237 cout << " | | Trigger Bits |" << endl;
2238 cout << "Track A ID | Track B ID | Raw HiMass LoMass Unlike HiPt# LoPt# | Invariant mass" << endl;
2239 cout << "----------------------------------------------------------------------------------------" << endl;
2240 const AliHLTMUONPairDecisionStruct* entry = block.GetArray();
2241 for(AliHLTUInt32_t i = 0; i < nentries; i++)
2242 {
2243 subResult = DumpPairDecisionStruct(buffer, bufferSize, entry++, continueParse);
2244 if (subResult != EXIT_SUCCESS) return subResult;
2245 }
2246
2247 return result;
2248}
2249
2250
2251int DumpCommonHeader(
2252 const char* buffer, unsigned long bufferSize,
2253 const AliHLTMUONDataBlockHeader* header, bool continueParse
2254 )
2255{
2256 // Step through the header fields trying to print them.
2257 // At each step check if we have not overflowed the buffer, if we have
2258 // not then we can print the field, otherwise we print the left over
2259 // bytes assumed to be corrupted rubbish.
2260 int result = CheckHeaderField(header->fType, buffer, bufferSize, continueParse);
2261 if (result != EXIT_SUCCESS) return result;
2262 AliHLTMUONDataBlockType type = AliHLTMUONDataBlockType(header->fType);
2263 cout << " Block type: " << type << endl;
2264
2265 result = CheckHeaderField(header->fRecordWidth, buffer, bufferSize, continueParse);
2266 if (result != EXIT_SUCCESS) return result;
2267 cout << " Record width: " << header->fRecordWidth << endl;
2268
2269 result = CheckHeaderField(header->fNrecords, buffer, bufferSize, continueParse);
2270 if (result != EXIT_SUCCESS) return result;
2271 cout << "Number of entries: " << header->fNrecords << endl;
2272
2273 return result;
2274}
2275
887a669c 2276/**
2277 * Method to look for a certain data word key in the buffer.
2278 * \param buffer The start location in the buffer to search.
2279 * \param end The end of the buffer.
2280 * \returns the pointer position just past the word found or
2281 * NULL if the word was not found in the buffer.
2282 */
2283const char* FindDataWord(const char* buffer, const char* end, UInt_t word)
2284{
2285 for (const char* current = buffer; (current+1) <= end; current += sizeof(UInt_t))
2286 {
2287 const UInt_t* currentWord = reinterpret_cast<const UInt_t*>(current);
2288 if (*currentWord == word) return current + sizeof(UInt_t);
2289 }
2290 return NULL;
2291}
2292
2293/**
2294 * Method to check if the data buffer is really a raw DDL stream.
2295 * \returns kUnknownDataBlock if this does not look like a raw DDL stream.
2296 * kTrackerDDLRawData if this looks like a raw DDL stream from the tracker.
2297 * kTriggerDDLRawData if this looks like a raw DDL stream from the trigger.
2298 */
2299int CheckIfDDLStream(const char* buffer, unsigned long bufferSize)
2300{
2301 if (bufferSize < sizeof(AliRawDataHeader)) return kUnknownDataBlock;
2302
2303 const AliRawDataHeader* cdhHeader =
2304 reinterpret_cast<const AliRawDataHeader*>(buffer);
2305
2306 // Keep scores of indicators / tests that show this is a raw DDL stream
2307 // either from the trigger or tracker. We will decide if the stream is
2308 // indeed a raw DDL stream if the largest of the two scores is above a
2309 // minimum threshold.
2310 int trackerScore = 0;
2311 int triggerScore = 0;
2312
2313 if (cdhHeader->fSize == UInt_t(-1) or cdhHeader->fSize == bufferSize)
2314 {
2315 trackerScore++;
2316 triggerScore++;
2317 }
2318
2319 if (cdhHeader->GetVersion() == 2)
2320 {
2321 trackerScore++;
2322 triggerScore++;
2323 }
2324
2325 const char* payload = buffer + sizeof(AliRawDataHeader);
2326 const char* payloadEnd = buffer + bufferSize;
2327
2328 typedef AliMUONTrackerDDLDecoder<AliMUONTrackerDDLDecoderEventHandler> AliTrkDecoder;
2329 typedef AliMUONTriggerDDLDecoder<AliMUONTriggerDDLDecoderEventHandler> AliTrgDecoder;
2330
2331 // See if the DDL data has a payload with data word keys as expected by
2332 // AliMUONTrackerDDLDecoder.
2333 const char* current = payload;
2334 while ( (current = FindDataWord(current, payloadEnd, AliTrkDecoder::BlockDataKeyWord())) != NULL )
2335 {
2336 trackerScore++;
2337 }
2338 current = payload;
2339 while ( (current = FindDataWord(current, payloadEnd, AliTrkDecoder::DspDataKeyWord())) != NULL )
2340 {
2341 trackerScore++;
2342 }
2343 current = payload;
2344 while ( (current = FindDataWord(current, payloadEnd, AliTrkDecoder::BusPatchDataKeyWord())) != NULL )
2345 {
2346 trackerScore++;
2347 }
2348 current = payload;
2349 while ( (current = FindDataWord(current, payloadEnd, AliTrkDecoder::EndOfDDLWord())) != NULL )
2350 {
2351 trackerScore++;
2352 }
2353
2354 // See if the DDL data has a payload with data word keys as expected by
2355 // AliMUONTriggerDDLDecoder.
2356 current = payload;
2357 while ( (current = FindDataWord(current, payloadEnd, AliTrgDecoder::EndOfDarcWord())) != NULL )
2358 {
2359 triggerScore++;
2360 }
2361 current = payload;
2362 while ( (current = FindDataWord(current, payloadEnd, AliTrgDecoder::EndOfGlobalWord())) != NULL )
2363 {
2364 triggerScore++;
2365 }
2366 current = payload;
2367 while ( (current = FindDataWord(current, payloadEnd, AliTrgDecoder::EndOfRegionalWord())) != NULL )
2368 {
2369 triggerScore++;
2370 }
2371 current = payload;
2372 while ( (current = FindDataWord(current, payloadEnd, AliTrgDecoder::EndOfLocalWord())) != NULL )
2373 {
2374 triggerScore++;
2375 }
2376
2377 if (triggerScore > trackerScore)
2378 {
2379 if (triggerScore >= 6) return kTriggerDDLRawData;
2380 }
2381 else
2382 {
2383 if (trackerScore >= 6) return kTrackerDDLRawData;
2384 }
2385 return kUnknownDataBlock;
2386}
3dad52b0 2387
887a669c 2388/**
2389 * Parses the buffer and prints the contents to screen.
8984a6aa 2390 * \param [in] buffer The pointer to the buffer to parse.
2391 * \param [in] bufferSize The size of the buffer in bytes.
2392 * \param [in] continueParse If specified then the we try to continue parsing the
887a669c 2393 * buffer as much as possible.
8984a6aa 2394 * \param [in] tryrecover Indicates if the DDL decoders should have special
daab896b 2395 * recovery logic enabled.
8984a6aa 2396 * \param [in,out] type Initialy this should indicate the type of the data block
887a669c 2397 * or kUnknownDataBlock if not known. On exit it will be filled with
2398 * the type of the data block as discovered by this routine if type
2399 * initially contained kUnknownDataBlock.
2400 * \returns The error code indicating the problem. EXIT_SUCCESS is returned
2401 * on success.
2402 */
3dad52b0 2403int ParseBuffer(
daab896b 2404 char* buffer, unsigned long bufferSize,
2405 bool continueParse, bool tryrecover, int& type
3dad52b0 2406 )
2407{
2408 assert( buffer != NULL );
2409 int result = EXIT_SUCCESS;
887a669c 2410 int subResult = EXIT_FAILURE;
3dad52b0 2411
887a669c 2412 // If the -type|-t option was not used in the command line then we need to
2413 // figure out what type of data block this is from the data itself.
2414 bool ddlStream = false;
2415 if (type == kUnknownDataBlock)
3dad52b0 2416 {
887a669c 2417 // First check if this is a raw DDL stream, if not then assume it is
2418 // some kind of internal dHLT raw data block.
2419 int streamType = CheckIfDDLStream(buffer, bufferSize);
2420 if (streamType == kTrackerDDLRawData or streamType == kTriggerDDLRawData)
2421 {
2422 type = streamType;
2423 ddlStream = true;
2424 }
3dad52b0 2425 }
887a669c 2426 else if (type == kTrackerDDLRawData or type == kTriggerDDLRawData)
3dad52b0 2427 {
887a669c 2428 ddlStream = true;
3dad52b0 2429 }
887a669c 2430
2431 if (not ddlStream)
3dad52b0 2432 {
887a669c 2433 if (bufferSize < sizeof(AliHLTMUONDataBlockHeader))
2434 {
2435 cerr << "ERROR: The size of the file is too small to contain a"
2436 " valid data block." << endl;
2437 result = PARSE_ERROR;
2438 if (not continueParse) return result;
2439 }
2440 const AliHLTMUONDataBlockHeader* header =
2441 reinterpret_cast<const AliHLTMUONDataBlockHeader*>(buffer);
2442
2443 subResult = DumpCommonHeader(buffer, bufferSize, header, continueParse);
2444 if (subResult != EXIT_SUCCESS) return subResult;
2445
2446
2447 // Check if the block type in the header corresponds to the type given
2448 // by the '-type' command line parameter. If they do not then print an
2449 // error or big fat warning message and force interpretation of the data
2450 // block with the type given by '-type'.
2451 AliHLTMUONDataBlockType headerType = AliHLTMUONDataBlockType(header->fType);
2452
2453 if (type == kUnknownDataBlock)
2454 {
2455 // -type not used in the command line so just use what is given
2456 // by the data block header.
2457 type = headerType;
2458 }
2459 else if (type != headerType)
2460 {
2461 cerr << "WARNING: The data block header indicates a type"
2462 " different from what was specified on the command line."
2463 " The data could be corrupt."
2464 << endl;
2465 cerr << "WARNING: The type value in the file is "
2466 << showbase << hex << header->fType
2467 << " (" << headerType << "), but on the command line it is "
2468 << showbase << hex << int(type) << dec
2469 << " (" << type << ")."
2470 << endl;
2471 cerr << "WARNING: Will force the interpretation of the data block"
2472 " with a type of " << type << "." << endl;
2473 }
3dad52b0 2474 }
2475
887a669c 2476 // Now we know what type the data block is supposed to be, so we can
3dad52b0 2477 // dump it to screen with the appropriate dump routine.
2478 switch (type)
2479 {
887a669c 2480 case kTrackerDDLRawData:
daab896b 2481 subResult = DumpTrackerDDLRawStream(buffer, bufferSize, continueParse, tryrecover);
887a669c 2482 if (subResult != EXIT_SUCCESS) result = subResult;
2483 break;
2484 case kTriggerDDLRawData:
daab896b 2485 subResult = DumpTriggerDDLRawStream(buffer, bufferSize, continueParse, tryrecover);
887a669c 2486 if (subResult != EXIT_SUCCESS) result = subResult;
2487 break;
3dad52b0 2488 case kTriggerRecordsDataBlock:
2489 subResult = DumpTriggerRecordsBlock(buffer, bufferSize, continueParse);
2490 if (subResult != EXIT_SUCCESS) result = subResult;
2491 break;
2492 case kTrigRecsDebugDataBlock:
2493 subResult = DumpTrigRecsDebugBlock(buffer, bufferSize, continueParse);
2494 if (subResult != EXIT_SUCCESS) result = subResult;
2495 break;
3dad52b0 2496 case kRecHitsDataBlock:
2497 subResult = DumpRecHitsBlock(buffer, bufferSize, continueParse);
2498 if (subResult != EXIT_SUCCESS) result = subResult;
2499 break;
2500 case kClustersDataBlock:
2501 subResult = DumpClustersBlock(buffer, bufferSize, continueParse);
2502 if (subResult != EXIT_SUCCESS) result = subResult;
2503 break;
2504 case kChannelsDataBlock:
329db59d 2505 subResult = DumpChannelsBlock(buffer, bufferSize, continueParse);
3dad52b0 2506 if (subResult != EXIT_SUCCESS) result = subResult;
2507 break;
2508 case kMansoTracksDataBlock:
2509 subResult = DumpMansoTracksBlock(buffer, bufferSize, continueParse);
2510 if (subResult != EXIT_SUCCESS) result = subResult;
2511 break;
2512 case kMansoCandidatesDataBlock:
2513 subResult = DumpMansoCandidatesBlock(buffer, bufferSize, continueParse);
2514 if (subResult != EXIT_SUCCESS) result = subResult;
2515 break;
a8982f78 2516 case kTracksDataBlock:
2517 subResult = DumpTracksBlock(buffer, bufferSize, continueParse);
2518 if (subResult != EXIT_SUCCESS) result = subResult;
2519 break;
3dad52b0 2520 case kSinglesDecisionDataBlock:
2521 subResult = DumpSinglesDecisionBlock(buffer, bufferSize, continueParse);
2522 if (subResult != EXIT_SUCCESS) result = subResult;
2523 break;
2524 case kPairsDecisionDataBlock:
efe258e1 2525 subResult = DumpPairsDecisionBlock(buffer, bufferSize, continueParse);
3dad52b0 2526 if (subResult != EXIT_SUCCESS) result = subResult;
2527 break;
2528 default :
2529 cout << "ERROR: Unknown data block type. Found a type number of "
710ddf4b 2530 << showbase << hex << int(type) << noshowbase << dec
3dad52b0 2531 << " (" << int(type) << ")." << endl;
2532 result = PARSE_ERROR;
2533 }
2534
2535 return result;
2536}
2537
887a669c 2538/**
2539 * Convert the type code to a string.
2540 */
2541const char* TypeToString(int type)
2542{
2543 if (type == kTrackerDDLRawData or type == kTriggerDDLRawData)
2544 {
2545 static char str[kAliHLTComponentDataTypefIDsize+1];
2546 AliHLTComponentDataType t = AliHLTMUONConstants::DDLRawDataType();
2547 memcpy(&str, &t.fID, kAliHLTComponentDataTypefIDsize);
2548 // Must insert the NULL character to make this an ANSI C string.
2549 str[kAliHLTComponentDataTypefIDsize] = '\0';
2550 return &str[0];
2551 }
2552 else
2553 {
2554 return AliHLTMUONUtils::DataBlockTypeToString(AliHLTMUONDataBlockType(type));
2555 }
2556}
2557
2558/**
2559 * Find the data specification from the filename and return it in string format.
2560 */
2561const char* TryDecodeDataSpec(const char* filename)
2562{
2563 TString name = filename;
2564
2565 TRegexp re1("MUONTR[GK]_[0123456789]+\\.ddl$");
2566 Ssiz_t length = 0;
2567 Ssiz_t pos = re1.Index(name, &length);
2568 if (pos != kNPOS)
2569 {
2570 TString substr;
2571 for (Ssiz_t i = 0; i < length; i++)
2572 substr += name[pos+i];
2573 TRegexp re("[0123456789]+");
2574 pos = re.Index(substr, &length);
2575 TString num;
2576 for (Ssiz_t i = 0; i < length; i++)
2577 num += substr[pos+i];
2578 AliHLTUInt32_t spec = AliHLTMUONUtils::EquipIdToSpec(num.Atoi());
2579 static char strbuf[32];
2580 sprintf(&strbuf[0], "0x%8.8X", spec);
2581 return &strbuf[0];
2582 }
2583
2584 TRegexp re2("_MUON\\:.+_0x[0123456789abscefABCDEF]+\\.dat$");
2585 pos = re2.Index(name, &length);
2586 if (pos != kNPOS)
2587 {
2588 TString substr;
2589 for (Ssiz_t i = 0; i < length; i++)
2590 substr += name[pos+i];
2591 TRegexp re("0x[0123456789abscefABCDEF]+");
2592 pos = re.Index(substr, &length);
2593 TString num;
2594 for (Ssiz_t i = 0; i < length; i++)
2595 num += substr[pos+i];
2596 static TString result = num;
2597 return result.Data();
2598 }
2599
2600 return NULL;
2601}
2602
ffb64d3e 2603namespace
2604{
2605 // CDB path and run number to use.
162637e4 2606 const char* gCDBPath = "local://$ALICE_ROOT/OCDB";
ffb64d3e 2607 Int_t gRunNumber = 0;
2608}
2609
887a669c 2610/**
2611 * Performs basic data integrity checks of the data block using the
2612 * AliHLTMUONDataCheckerComponent.
8984a6aa 2613 * \param [in] sys The HLT system framework.
2614 * \param [in] filename The name of the file containing the data block to check.
2615 * \param [in] type Must indicate the type of the data block.
2616 * \param [in] dataspec The data specification of the data block. NULL if none.
2617 * \param [in] maxLogging If set to true then full logging is turned on for AliHLTSystem.
887a669c 2618 * \returns The error code indicating the problem. EXIT_SUCCESS is returned
2619 * on success.
2620 */
2621int CheckDataIntegrity(
2622 AliHLTSystem& sys, const char* filename, int type,
2623 const char* dataspec, bool maxLogging
2624 )
2625{
2626 if (maxLogging)
2627 {
2628 AliLog::SetGlobalLogLevel(AliLog::kMaxType);
2629 sys.SetGlobalLoggingLevel(kHLTLogAll);
2630 }
2631 else
2632 {
2633 AliLog::SetGlobalLogLevel(AliLog::kWarning);
2634 int level = kHLTLogWarning | kHLTLogError | kHLTLogFatal;
2635 sys.SetGlobalLoggingLevel(AliHLTComponentLogSeverity(level));
2636 }
2637
ffb64d3e 2638 sys.LoadComponentLibraries("libAliHLTUtil.so");
2639 sys.LoadComponentLibraries("libAliHLTMUON.so");
887a669c 2640
2641 // Setup the component parameter lists and then the components.
ffb64d3e 2642 TString dcparams = "-return_error -warn_on_unexpected_block -no_global_check -cdbpath ";
2643 dcparams += gCDBPath;
2644 dcparams += " -run ";
2645 dcparams += gRunNumber;
887a669c 2646 TString fpparams = "-datatype '";
2647 fpparams += TypeToString(type);
2648 fpparams += "' 'MUON'";
2649 if (dataspec != NULL)
2650 {
2651 fpparams += " -dataspec ";
2652 fpparams += dataspec;
2653 }
2654 else
2655 {
2656 const char* spec = TryDecodeDataSpec(filename);
2657 if (spec != NULL)
2658 {
2659 fpparams += " -dataspec ";
2660 fpparams += spec;
2661 }
2662 else
2663 {
2664 dcparams += " -ignorespec";
2665 }
2666 }
2667 fpparams += " -datafile ";
2668 fpparams += filename;
2669 TString fpname = "filePublisher_";
2670 fpname += filename;
2671 TString dcname = "checker_";
2672 dcname += filename;
2673
2674 if (maxLogging)
2675 {
2676 cout << "DEBUG: Using the following flags for FilePublisher: \""
2677 << fpparams.Data() << "\""<< endl;
2678 cout << "DEBUG: Using the following flags for "
2679 << AliHLTMUONConstants::DataCheckerComponentId()
2680 << ": \"" << dcparams.Data() << "\""<< endl;
2681 }
2682
2683 AliHLTConfiguration(fpname.Data(), "FilePublisher", NULL, fpparams.Data());
2684 AliHLTConfiguration checker(
2685 dcname.Data(), AliHLTMUONConstants::DataCheckerComponentId(),
2686 fpname.Data(), dcparams.Data()
2687 );
2688
2689 // Build and run the HLT tasks.
2690 if (sys.BuildTaskList(dcname.Data()) != 0) return HLTSYSTEM_ERROR;
2691 if (maxLogging) sys.PrintTaskList();
2692 if (sys.Run() != 1) return HLTSYSTEM_ERROR;
2693
2694 // Now clean up.
2695 if (sys.CleanTaskList() != 0) return HLTSYSTEM_ERROR;
2696
2697 return EXIT_SUCCESS;
2698}
2699
3dad52b0 2700
2701/**
2702 * The caller is responsible for freeing memory allocated for buffer with a call
2703 * to delete [] buffer.
2704 */
2705int ReadFile(const char* filename, char*& buffer, unsigned long& bufferSize)
2706{
2707 assert( filename != NULL );
2708
2709 // Open the file and find its size.
2710 fstream file;
2711 file.open(filename, ios::in);
2712 if (not file)
2713 {
2714 cerr << "ERROR: Could not open the file: " << filename << endl;
2715 return SYSTEM_ERROR;
2716 }
2717 file.seekg(0, ios::end);
2718 if (not file)
2719 {
2720 cerr << "ERROR: Could not seek in the file: " << filename << endl;
2721 return SYSTEM_ERROR;
2722 }
2723 bufferSize = file.tellg();
2724 if (not file)
2725 {
2726 cerr << "ERROR: Could not get file size for the file: " <<
2727 filename << endl;
2728 return SYSTEM_ERROR;
2729 }
2730 file.seekg(0, ios::beg);
2731 if (not file)
2732 {
2733 cerr << "ERROR: Could not seek in the file: " << filename << endl;
2734 return SYSTEM_ERROR;
2735 }
2736
2737 // Allocate the memory for the file.
2738 try
2739 {
2740 buffer = new char[bufferSize];
2741 }
2742 catch (const std::bad_alloc&)
2743 {
2744 cerr << "ERROR: Out of memory. Tried to allocate " << bufferSize
2745 << " bytes." << endl;
2746 return SYSTEM_ERROR;
2747 }
2748
2749 file.read(buffer, bufferSize);
2750 if (not file)
2751 {
2752 delete [] buffer;
2753 buffer = NULL;
2754 bufferSize = 0;
2755 cerr << "ERROR: Could not read from file: " << filename << endl;
2756 return SYSTEM_ERROR;
2757 }
2758
2759 file.close();
2760 if (not file)
2761 {
2762 delete [] buffer;
2763 buffer = NULL;
2764 bufferSize = 0;
2765 cerr << "ERROR: Could not close the file: " << filename << endl;
2766 return SYSTEM_ERROR;
2767 }
2768
2769 return EXIT_SUCCESS;
2770}
2771
2772/**
2773 * Prints the command line usage of this program to standard error.
2774 */
c9537879 2775void PrintUsage(bool asError = true)
3dad52b0 2776{
c9537879 2777 std::ostream& os = asError ? cerr : cout;
887a669c 2778 os << "Usage: dHLTdumpraw [-help|-h] [-continue|-c] [-type|-t <typename>] [-check|-k]" << endl;
ffb64d3e 2779 os << " [-debug|-d] [-dataspec|-s <number>] [-cdbpath|-p <url>] [-run|-r <number>]" << endl;
2780 os << " <filename> [<filename> ...]" << endl;
c9537879 2781 os << "Where <filename> is the name of a file containing a raw data block." << endl;
2782 os << "Options:" << endl;
2783 os << " -help | -h" << endl;
2784 os << " Displays this message." << endl;
f68b4e0e 2785 os << " -continue | -c" << endl;
c9537879 2786 os << " If specified, the program will try to continue parsing the data block" << endl;
2787 os << " as much as possible rather than stopping at the first error." << endl;
daab896b 2788 os << " -recover | -e" << endl;
2789 os << " If specified, special DDL decoder recovery logic is enabled to handle" << endl;
2790 os << " corrupt DDL data." << endl;
f68b4e0e 2791 os << " -type | -t <typename>" << endl;
c9537879 2792 os << " Forces the contents of the subsequent files specified on the command" << endl;
2793 os << " line to be interpreted as a specific type of data block." << endl;
2794 os << " Where <typename> can be one of:" << endl;
887a669c 2795 os << " rawtracker - raw DDL stream from tracker chambers." << endl;
2796 os << " rawtrigger - raw DDL stream from trigger chambers." << endl;
c9537879 2797 os << " trigrecs - trigger records data." << endl;
2798 os << " trigrecsdebug - debugging information about trigger records." << endl;
c9537879 2799 os << " rechits - reconstructed hits data." << endl;
2800 os << " channels - channel debugging information from hit reconstruction." << endl;
2801 os << " clusters - cluster debugging information from hit reconstruction." << endl;
2802 os << " mansotracks - partial tracks from Manso algorithm." << endl;
2803 os << " mansocandidates - track candidates considered in the Manso algorithm." << endl;
a8982f78 2804 os << " tracks - tracks from full tracker component." << endl;
c9537879 2805 os << " singlesdecision - trigger decisions for single tracks." << endl;
2806 os << " pairsdecision - trigger decisions for track pairs." << endl;
2807 os << " autodetect - the type of the data block will be automatically" << endl;
2808 os << " detected." << endl;
887a669c 2809 os << " -check | -k" << endl;
2810 os << " If specified then data integrity checks are performed on the raw data." << endl;
2811 os << " Warnings and errors are printed as problems are found with the data, but" << endl;
2812 os << " the data will still be converted into ROOT objects as best as possible." << endl;
2813 os << " -debug | -d" << endl;
2814 os << " If specified, then the all debug messages are printed by the AliHLTSystem." << endl;
2815 os << " This is only useful if experiencing problems with the -check|-k option." << endl;
2816 os << " -dataspec | -s <number>" << endl;
2817 os << " When specified, then <number> is used as the data specification for the" << endl;
2818 os << " data file that follows. This option is only useful with the -check|-k option." << endl;
ffb64d3e 2819 os << " -cdbpath | -p <url>" << endl;
2820 os << " The path to the CDB to use when running with the -check | -k option." << endl;
2821 os << " -run | -r <number>" << endl;
2822 os << " The run number to use when running with the -check | -k option." << endl;
3dad52b0 2823}
2824
3dad52b0 2825/**
2826 * Parses the command line.
2827 * @param argc Number of arguments as given in main().
2828 * @param argv Array of arguments as given in main().
c9537879 2829 * @param filenames Pointer to buffer storing file name strings.
c9537879 2830 * @param filetypes Array that receives the type of the data block expected, i.e.
2831 * the value of the -type flag for the corresponding file.
887a669c 2832 * @param dataspecs Data specifications to use for the data files.
2833 * @param numOfFiles Receives the number of file name strings that were found
2834 * and added to 'filenames'.
2835 * @param continueParse Set to true if the user requested to continue to parse
2836 * after errors.
daab896b 2837 * @param tryrecover Set to true if the user wants DDL decoder recovery logic enabled.
887a669c 2838 * @param checkData Set to true if data integrity checking was requested.
2839 * @param maxLogging Set to true if maximal logging was requested.
3dad52b0 2840 * @return A status flag suitable for returning from main(), containing either
2841 * EXIT_SUCCESS or CMDLINE_ERROR.
2842 */
2843int ParseCommandLine(
c9537879 2844 int argc,
5b882880 2845 char** argv,
c9537879 2846 const char** filenames,
887a669c 2847 int* filetypes,
2848 const char** dataspecs,
c9537879 2849 int& numOfFiles,
2850 bool& continueParse,
daab896b 2851 bool& tryrecover,
887a669c 2852 bool& checkData,
2853 bool& maxLogging
3dad52b0 2854 )
2855{
c9537879 2856 numOfFiles = 0;
3dad52b0 2857 continueParse = false;
daab896b 2858 tryrecover = false;
887a669c 2859 maxLogging = false;
2860 checkData = false;
2861 int currentType = kUnknownDataBlock;
2862 const char* currentDataSpec = NULL;
ffb64d3e 2863 bool pathSet = false;
2864 bool runSet = false;
3dad52b0 2865
2866 // Parse the command line.
2867 for (int i = 1; i < argc; i++)
2868 {
f68b4e0e 2869 if (strcmp(argv[i], "-help") == 0 or strcmp(argv[i], "-h") == 0)
3dad52b0 2870 {
c9537879 2871 PrintUsage(false);
3dad52b0 2872 return EXIT_SUCCESS;
2873 }
f68b4e0e 2874 else if (strcmp(argv[i], "-continue") == 0 or strcmp(argv[i], "-c") == 0)
3dad52b0 2875 {
2876 continueParse = true;
2877 }
daab896b 2878 else if (strcmp(argv[i], "-recover") == 0 or strcmp(argv[i], "-e") == 0)
2879 {
2880 tryrecover = true;
2881 }
f68b4e0e 2882 else if (strcmp(argv[i], "-type") == 0 or strcmp(argv[i], "-t") == 0)
3dad52b0 2883 {
c9537879 2884 if (++i >= argc)
3dad52b0 2885 {
887a669c 2886 cerr << "ERROR: Missing a type specifier." << endl << endl;
3dad52b0 2887 PrintUsage();
2888 return CMDLINE_ERROR;
2889 }
c9537879 2890 // Now we need to parse the typename in the command line.
2891 if (strcmp(argv[i], "autodetect") == 0)
2892 {
2893 currentType = kUnknownDataBlock;
2894 }
887a669c 2895 else if (strcmp(argv[i], "rawtracker") == 0)
2896 {
2897 currentType = kTrackerDDLRawData;
2898 }
2899 else if (strcmp(argv[i], "rawtrigger") == 0)
2900 {
2901 currentType = kTriggerDDLRawData;
2902 }
3dad52b0 2903 else
c9537879 2904 {
f68b4e0e 2905 currentType = AliHLTMUONUtils::ParseCommandLineTypeString(argv[i]);
2906 if (currentType == kUnknownDataBlock)
2907 {
2908 cerr << "ERROR: Invalid type name '" << argv[i]
2909 << "' specified for argument " << argv[i-1]
2910 << "." << endl << endl;
2911 PrintUsage();
2912 return CMDLINE_ERROR;
2913 }
c9537879 2914 }
2915 }
887a669c 2916 else if (strcmp(argv[i], "-debug") == 0 or strcmp(argv[i], "-d") == 0)
2917 {
2918 maxLogging = true;
2919 }
2920 else if (strcmp(argv[i], "-check") == 0 or strcmp(argv[i], "-k") == 0)
2921 {
2922 checkData = true;
2923 }
2924 else if (strcmp(argv[i], "-dataspec") == 0 or strcmp(argv[i], "-s") == 0)
2925 {
2926 if (++i >= argc)
2927 {
2928 cerr << "ERROR: Missing a specification number." << endl << endl;
2929 PrintUsage();
2930 return CMDLINE_ERROR;
2931 }
2932
2933 char* errPos = NULL;
2934 unsigned long num = strtoul(argv[i], &errPos, 0);
2935 if (errPos == NULL or *errPos != '\0')
2936 {
2937 cerr << "ERROR: Cannot convert '%s' to a data specification number."
2938 << argv[i] << endl;
2939 return CMDLINE_ERROR;
2940 }
2941 if (not AliHLTMUONUtils::IsSpecValid(num))
2942 {
2943 cerr << "ERROR: The data specification number is not a valid format."
2944 << endl;
2945 return CMDLINE_ERROR;
2946 }
2947 currentDataSpec = argv[i];
2948 }
ffb64d3e 2949 else if (strcmp(argv[i], "-cdbpath") == 0 or strcmp(argv[i], "-p") == 0)
2950 {
2951 if (pathSet)
2952 {
2953 cerr << "WARNING: Already used -cdbpath|-p with '" << gCDBPath
2954 << "' before. Will override it with the last value specified with -cdbpath|-p."
2955 << endl;
2956 }
2957 if (++i >= argc)
2958 {
2959 cerr << "ERROR: Missing the URL for the CDB path." << endl << endl;
2960 PrintUsage();
2961 return CMDLINE_ERROR;
2962 }
2963 gCDBPath = argv[i];
2964 pathSet = true;
2965 }
2966 else if (strcmp(argv[i], "-run") == 0 or strcmp(argv[i], "-r") == 0)
2967 {
2968 if (runSet)
2969 {
2970 cerr << "WARNING: Already used -run|-r with " << gRunNumber
2971 << " before. Will override it with the last value specified with -run|-r."
2972 << endl;
2973 }
2974 if (++i >= argc)
2975 {
2976 cerr << "ERROR: Missing the run number." << endl << endl;
2977 PrintUsage();
2978 return CMDLINE_ERROR;
2979 }
2980
2981 char* cpErr = NULL;
2982 Int_t run = Int_t( strtol(argv[i], &cpErr, 0) );
2983 if (cpErr == NULL or *cpErr != '\0' or run < 0)
2984 {
2985 cerr << "ERROR: Cannot convert '" << argv[i] << "' to a valid run number."
2986 " Expected a positive integer value." << endl;
2987 return CMDLINE_ERROR;
2988 }
2989
2990 gRunNumber = run;
2991 runSet = true;
2992 }
c9537879 2993 else
2994 {
2995 assert( numOfFiles < argc );
2996 filenames[numOfFiles] = argv[i];
2997 filetypes[numOfFiles] = currentType;
887a669c 2998 dataspecs[numOfFiles] = currentDataSpec;
2999 currentDataSpec = NULL; // Reset because '-dataspec' option is only valid for one file.
c9537879 3000 numOfFiles++;
3dad52b0 3001 }
3002 }
3003
c9537879 3004 // Now check that we have at least one filename and all the flags we need.
3005 if (numOfFiles == 0)
3dad52b0 3006 {
c9537879 3007 cerr << "ERROR: Missing a file name. You must specify at least one file to process."
3dad52b0 3008 << endl << endl;
3009 PrintUsage();
3010 return CMDLINE_ERROR;
3011 }
3012
3013 return EXIT_SUCCESS;
3014}
3015
3016
5b882880 3017int main(int argc, char** argv)
3dad52b0 3018{
1120d5de 3019 // Test endianess of this machine during runtime and print warning if it is not
3020 // little endian.
3021 union
3022 {
3023 int dword;
3024 char byte[4];
3025 } endianTest;
3026 endianTest.dword = 0x1;
3027 if (endianTest.byte[0] != 0x1)
3028 {
3029 cerr << "!!!! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !!!!!" << endl;
3030 cerr << "!!!! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !!!!!" << endl;
3031 cerr << "!! !!" << endl;
3032 cerr << "!! This is not a little endian machine, but dHLT raw data is normally !!" << endl;
3033 cerr << "!! generated in little endian format. Unless you are looking at localy !!" << endl;
3034 cerr << "!! created simulated data, then this program will not show you correct !!" << endl;
3035 cerr << "!! output. !!" << endl;
3036 cerr << "!! !!" << endl;
3037 cerr << "!!!! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !!!!!" << endl;
3038 cerr << "!!!! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !! WARNING !!!!!" << endl;
3039 cerr << endl;
3040 }
3041
3042
c9537879 3043 int numOfFiles = 0;
3dad52b0 3044 bool continueParse = false;
daab896b 3045 bool tryrecover = false;
887a669c 3046 bool checkData = false;
3047 bool maxLogging = false;
3dad52b0 3048 int returnCode = EXIT_SUCCESS;
3dad52b0 3049 char* buffer = NULL;
887a669c 3050 const char** filename = NULL;
3051 int* filetype = NULL;
3052 const char** dataspec = NULL;
3dad52b0 3053
3054 try
3055 {
887a669c 3056 // Reduce logging now to get rid of informationals from AliHLTSystem constructor.
3057 AliLog::SetGlobalLogLevel(AliLog::kWarning);
3058 AliHLTSystem sys;
3059
c9537879 3060 // There will be at least 'argc' number of filenames.
3061 typedef const char* AnsiString;
887a669c 3062 filename = new AnsiString[argc];
3063 filetype = new int[argc];
3064 dataspec = new AnsiString[argc];
c9537879 3065
887a669c 3066 returnCode = ParseCommandLine(
3067 argc, argv, filename, filetype, dataspec, numOfFiles,
daab896b 3068 continueParse, tryrecover, checkData, maxLogging
887a669c 3069 );
3dad52b0 3070
c9537879 3071 if (returnCode == EXIT_SUCCESS)
3dad52b0 3072 {
c9537879 3073 for (int i = 0; i < numOfFiles; i++)
3074 {
3075 unsigned long bufferSize = 0;
3076 returnCode = ReadFile(filename[i], buffer, bufferSize);
3077 if (returnCode != EXIT_SUCCESS) break;
a2be583d 3078 if (numOfFiles > 1)
3079 {
887a669c 3080 cout << "########## Start of dump for file: "
3081 << filename[i] << " ##########" << endl;
a2be583d 3082 }
daab896b 3083 int result = ParseBuffer(buffer, bufferSize, continueParse, tryrecover, filetype[i]);
c9537879 3084 if (buffer != NULL) delete [] buffer;
3085 if (result != EXIT_SUCCESS)
3086 {
3087 returnCode = result;
3088 if (not continueParse) break;
3089 }
887a669c 3090 if (checkData)
3091 {
3092 result = CheckDataIntegrity(
3093 sys, filename[i], filetype[i],
3094 dataspec[i], maxLogging
3095 );
3096 if (result != EXIT_SUCCESS)
3097 {
3098 returnCode = result;
3099 if (not continueParse) break;
3100 }
3101 }
a2be583d 3102 if (numOfFiles > 1)
3103 {
887a669c 3104 cout << "########## End of dump for file: " <<
3105 filename[i] << " ##########" << endl;
a2be583d 3106 }
c9537879 3107 }
3dad52b0 3108 }
3109
c9537879 3110 delete [] filename;
3111 delete [] filetype;
887a669c 3112 delete [] dataspec;
3dad52b0 3113 }
3114 catch (...)
3115 {
3116 cerr << "FATAL ERROR: An unknown exception occurred!" << endl << endl;
3117 returnCode = FATAL_ERROR;
3118 if (buffer != NULL) delete [] buffer;
887a669c 3119 if (filename != NULL) delete [] filename;
3120 if (filetype != NULL) delete [] filetype;
3121 if (dataspec != NULL) delete [] dataspec;
3dad52b0 3122 }
c9537879 3123
3dad52b0 3124 return returnCode;
3125}