5 Utility to simulate a DATE raw data stream using a given set of raw
6 data files and a configuration file.
10 V01.00 4/05/2004 RD Created
11 V01.01 25/10/2005 RD Support added for timestamp
12 V01.02 4/04/2006 RD Support for CDH
13 V01.03 24/05/2006 RD Added "Direct disk access" option
26 #define DESCRIPTION "DATE raw data stream simulator"
30 char fileHandlerIdent[]= "@(#)""" __FILE__ """: """ DESCRIPTION \
32 """ compiled """ __DATE__ """ """ __TIME__;
34 #define DBG_BASE if ( debug > 0 )
35 #define DBG_DETAILED if ( debug > 1 )
36 #define DBG_VERBOSE if ( debug > 2 )
39 # define TRUE (0 == 0)
42 # define FALSE (0 == 1)
48 typedef enum { unknown, ldc, gdc } workingAsType;
49 typedef enum { collider, fixedTarget } workingModeType;
50 workingAsType workingAs;
51 workingModeType workingMode;
52 struct ldcDescriptorStruct {
54 struct ldcDescriptorStruct *next;
55 } *ldcsHead, *ldcsTail;
56 void *eventsHead, *eventsTail;
57 struct gdcEventDescriptorStruct {
58 struct ldcEventDescriptorStruct *head;
59 struct ldcEventDescriptorStruct *tail;
60 struct gdcEventDescriptorStruct *next;
61 struct eventHeaderStruct header;
64 struct ldcEventDescriptorStruct {
65 struct equipmentEventDescriptorStruct *head;
66 struct equipmentEventDescriptorStruct *tail;
67 struct ldcEventDescriptorStruct *next;
69 struct eventHeaderStruct header;
72 struct equipmentEventDescriptorStruct {
73 struct equipmentEventDescriptorStruct *next;
75 struct payloadDescriptorStruct *payload;
76 struct equipmentHeaderStruct header;
78 struct payloadDescriptorStruct {
79 struct payloadDescriptorStruct *next;
84 } *payloadsHead, *payloadsTail;
86 eventGdcIdType currGdcId;
87 unsigned long32 currDetPattern;
88 eventLdcIdType currLdcId;
89 equipmentIdType currEquipmentId;
95 eventIdType oneEventDelta;
96 eventIdType currEventId;
100 struct commonDataHeaderStruct *cdhRef = NULL;
102 void dumpPayload( const struct payloadDescriptorStruct *p ) {
107 if ( p->data != NULL ) {
108 for ( i = 0, c = (char *)p->data, printable = TRUE;
109 printable && i != p->size;
111 printable = isascii( *c );
114 for ( i = 0, c = (char *)p->data; i != p->size; c++, i++ ) {
120 if ( *c != '\n' ) printf( "\"\n" );
123 for ( i = 0, v = (long32 *)p->data;
126 if ( i % (4*8) == 0 ) {
127 if ( i != 0 ) printf( "\n" );
130 printf( "%08x ", *v );
136 while ( i != p->size ) {
137 printf( "%02x ", *((char *)p->data + p->size - j - 1) & 0xff );
145 } /* End of dumpPayload */
148 assert( workingAs == ldc || workingAs == gdc );
149 if ( eventsHead != NULL ) {
150 printf( "Events:\n" );
151 if ( workingAs == gdc ) {
152 struct gdcEventDescriptorStruct *gdc;
154 for ( gdc = (struct gdcEventDescriptorStruct *)eventsHead;
157 struct ldcEventDescriptorStruct *ldc;
159 printf( " GDC (%p)\n", (void*)gdc );
160 for ( ldc = gdc->head; ldc != NULL; ldc = ldc->next ) {
161 struct equipmentEventDescriptorStruct *eq;
163 printf( " LDC (%p): %d\n", (void*)ldc, ldc->id );
164 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
165 printf( " EQUIPMENT (%p): %d PAYLOAD (%p):",
168 (void*)eq->payload );
170 printf( "\"%s\" (%d bytes)\n",
171 eq->payload->fileName,
173 dumpPayload( eq->payload );
178 if ( workingAs == ldc ) {
179 struct ldcEventDescriptorStruct *ldc;
181 for ( ldc = (struct ldcEventDescriptorStruct *)eventsHead;
184 struct equipmentEventDescriptorStruct *eq;
187 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
188 printf( " EQUIPMENT (%p): %d PAYLOAD (%p):",
191 (void*)eq->payload );
193 printf( "\"%s\" (%d bytes)\n",
194 eq->payload->fileName,
196 dumpPayload( eq->payload );
201 printf( "Events: EMPTY\n" );
203 } /* End of dumpEvents */
205 void getLine( char *line, const int maxSize ) {
209 for ( read = 0; !feof( stdin ) && !ferror( stdin ) && read != maxSize; read++ ) {
210 if ( (line[read] = getchar()) == '\n' ) break;
212 if ( ferror( stdin ) ) {
214 "%s: failed to read configuration input errno:%d ",
219 if ( feof( stdin ) ) read--;
220 if ( read == maxSize && line[read] != '\n' ) {
222 "%s: Input line # %d too long (%d chars max)\n",
223 myName, lineNmb, maxSize-1 );
228 if ( !( read == 0 && feof( stdin ) ) ) {
229 printf( "%d) [%3d] \"%s\"", lineNmb, read, line );
232 for ( c = 0; c != read; c++ ) {
233 if ( line[c] == '#' ) {
240 printf( " => \"%s\"", line );
242 if ( feof( stdin ) ) printf( "<<< EOF >>>" );
243 if ( ferror( stdin ) ) printf( "<<< FERROR >>>" );
246 } /* End of getLine */
248 void handleLdc( eventLdcIdType ldcId ) {
249 struct ldcDescriptorStruct *ldc;
251 if ( ldcsHead != NULL ) {
252 for ( ldc = ldcsHead; ldc != NULL; ldc = ldc->next ) {
253 if ( ldc->id == ldcId ) {
258 if ( (ldc = (struct ldcDescriptorStruct *)malloc( sizeof( *ldc ) )) == NULL ) {
260 "%s: Failed to malloc for %d bytes (struct ldcDescriptorStruct)\n",
261 myName, (int)sizeof( *ldc ) );
266 if ( ldcsHead == NULL ) {
267 ldcsHead = ldcsTail = ldc;
269 ldcsTail->next = ldc;
273 } /* End of handleLdc */
275 void createNewEvent() {
276 assert( workingAs == ldc || workingAs == gdc );
277 if ( workingAs == ldc ) {
278 struct ldcEventDescriptorStruct *p;
280 if ( (p = (struct ldcEventDescriptorStruct *)malloc( sizeof( *p ) ))
283 "%s: failed to malloc for %d bytes (createNewEvent: struct ldcEventDescriptorStruct)",
284 myName, (int)sizeof( *p ) );
289 p->head = p->tail = NULL;
292 if ( eventsHead == NULL ) {
293 eventsHead = eventsTail = p;
295 struct ldcEventDescriptorStruct *q =
296 (struct ldcEventDescriptorStruct *)eventsTail;
302 } else if ( workingAs == gdc ) {
303 struct gdcEventDescriptorStruct *p;
305 if ( (p = (struct gdcEventDescriptorStruct *)malloc( sizeof( *p ) ))
308 "%s: failed to malloc for %d bytes (createNewEvent: struct gdcEventDescriptorStruct)",
309 myName, (int)sizeof( *p ) );
315 p->head = p->tail = NULL;
317 if ( eventsHead == NULL ) {
318 eventsHead = eventsTail = p;
320 struct gdcEventDescriptorStruct *q =
321 (struct gdcEventDescriptorStruct *)eventsTail;
327 } /* End of createNewEvent */
329 void createNewLdcEvent() {
330 struct gdcEventDescriptorStruct *gdcDesc;
331 struct ldcEventDescriptorStruct *p;
333 if ( (p = (struct ldcEventDescriptorStruct *)malloc( sizeof( *p ) ))
336 "%s: failed to malloc for %d bytes (createNewLdcEvent: struct ldcEventDescriptorStruct)",
337 myName, (int)sizeof( *p ) );
342 p->head = p->tail = NULL;
344 gdcDesc = (struct gdcEventDescriptorStruct *)eventsTail;
345 if ( gdcDesc->head == NULL ) {
346 gdcDesc->head = gdcDesc->tail = p;
348 gdcDesc->tail->next = p;
352 } /* End of createNewLdcEvent */
354 void loadBuffer( struct payloadDescriptorStruct * const payload ) {
358 if ( (f = fopen( payload->fileName, "r" )) == NULL ) {
360 "%s: line:%d payload file \"%s\" not found or not readable, errno:%d. ",
365 perror( "System-dependent error " );
368 if ( (payload->data = malloc( payload->size )) == NULL ) {
370 "%s: line:%d Failed to malloc for payload file \"%s\" size:%d errno:%d ",
376 perror( "System-dependent status " );
379 if ( (bytesRead = fread( payload->data, payload->fileSize, 1, f )) != 1 ) {
381 "%s: line:%d Failed to read payload file \"%s\" size:%d requested:1 got:%d feof:%s ferror:%s errno:%d ",
387 feof(f) ? "TRUE" : "false",
388 ferror(f) ? "TRUE" : "false",
390 perror( "System-dependent status " );
394 if ( payload->size != payload->fileSize ) {
395 memset( (char *)payload->data + payload->fileSize,
397 payload->size - payload->fileSize );
399 } /* End of loadBuffer */
401 void unloadBuffer( struct payloadDescriptorStruct * const payload ) {
402 if ( payload->data != NULL ) {
403 free( payload->data );
404 payload->data = NULL;
406 } /* End of unloadBuffer */
408 void unloadAllBuffers() {
409 struct payloadDescriptorStruct *payload;
411 for ( payload = payloadsHead; payload != NULL; payload = payload->next ) {
412 unloadBuffer( payload );
414 } /* End of unloadAllBuffers */
416 void loadPayload( const char *fileName ) {
417 struct payloadDescriptorStruct *payload;
419 for ( payload = payloadsHead; payload != NULL; payload = payload->next ) {
420 if ( strcmp( fileName, payload->fileName ) == 0 )
423 if ( payload == NULL ) {
426 if ( (payload = (struct payloadDescriptorStruct *)malloc( sizeof( *payload ) ))
429 "%s: failed to malloc for %d bytes (loadPayload/payloadDescriptorStruct)\n",
431 (int)sizeof( *payload ) );
434 if ( (payload->fileName = strdup( fileName )) == NULL ) {
436 "%s: failed to duplicate string \"%s\" (loadPaload/fileName)\n",
441 if ( (f = fopen( fileName, "r" )) == NULL ) {
443 "%s: line:%d payload file \"%s\" not found or not readable, errno:%d. ",
448 perror( "System-dependent error " );
451 if ( fseek( f, 0L, SEEK_END ) != 0 ) {
453 "%s: line:%d Failed to seek payload file \"%s\" errno:%d ",
458 perror( "System-dependent error " );
461 if ( (payload->size = ftell( f )) <= 0 ) {
463 "%s: line:%d Failed to get file \"%s\" size size:%d errno:%d ",
469 perror( "System-dependent status " );
472 payload->fileSize = payload->size;
473 while ( (payload->size & 3) != 0 ) payload->size++;
477 loadBuffer( payload );
479 payload->data = NULL;
482 payload->next = NULL;
483 if ( payloadsHead == NULL ) {
484 payloadsHead = payloadsTail = payload;
486 payloadsTail->next = payload;
487 payloadsTail = payload;
492 printf( "%d) Payload \"%s\" loaded at %p\n",
498 strncmp(fileName,"TRG_",4) != 0 ) {
499 struct commonDataHeaderStruct *cdh =
500 (struct commonDataHeaderStruct *)payload->data;
502 printf( " CDH: blockLenght:%d=0x%08x ",
503 cdh->cdhBlockLength, cdh->cdhBlockLength );
504 if ( cdh->cdhBlockLength < sizeof( *cdh ) ) {
505 printf( "TOO SMALL (minimum:%ld=0x%08lx)\n",
506 (unsigned long)sizeof( *cdh ),
507 (unsigned long)sizeof( *cdh ) );
509 printf( "version:%d=0x%x ", cdh->cdhVersion, cdh->cdhVersion );
510 if ( cdh->cdhVersion != CDH_VERSION ) {
511 printf( "EXPECTED:%d=%x (decoding may be inaccurate) ",
512 CDH_VERSION, CDH_VERSION );
515 printf( "L1TriggerMessage:0x%x", cdh->cdhL1TriggerMessage );
516 if ( cdh->cdhL1TriggerMessage != 0 ) {
517 for ( b = 0, n = 0; b != 10; b++ ) {
518 if ( (cdh->cdhL1TriggerMessage & (1<<b)) != 0 ) {
519 if ( n++ != 0 )printf( "+" );
521 case 0: printf( "L1SwC" ); break;
522 case 1: printf( "ESR" ); break;
523 case 2: printf( "RoC1" ); break;
524 case 3: printf( "RoC2" ); break;
525 case 4: printf( "RoC3" ); break;
526 case 5: printf( "RoC4" ); break;
527 case 6: printf( "ClT" ); break;
528 default: printf( "spare %d", b+14 );
535 if ( cdh->cdhMBZ0 != 0 )
536 printf( "MBZ0:0x%x ",
541 printf( "EventId2(orbit):%d=0x%x ",
542 cdh->cdhEventId2, cdh->cdhEventId2 );
543 printf( "EventId1(bunchCrossing):%d=0x%x ",
544 cdh->cdhEventId1, cdh->cdhEventId1 );
548 if ( cdh->cdhMBZ1 != 0 )
549 printf( "MBZ1:0x%x ",
551 printf( "BlockAttributes:0x%x",
552 cdh->cdhBlockAttributes );
553 if ( cdh->cdhBlockAttributes != 0 ) {
555 for ( b = 0, n = 0; b != 8; b++ ) {
556 if ( (cdh->cdhBlockAttributes & (1<<b)) != 0 ) {
565 printf( "ParticipatingSubDetectors:0x%x ",
566 cdh->cdhParticipatingSubDetectors );
570 if ( cdh->cdhMBZ2 != 0 )
571 printf( "MBZ2:0x%x ",
573 printf( "Status/Error:0x%x", cdh->cdhStatusErrorBits );
574 if ( cdh->cdhStatusErrorBits != 0 ) {
576 for ( b = 0,n = 0; b != 16; b++ ) {
577 if ( (cdh->cdhStatusErrorBits & (1<<b)) != 0 ) {
578 if ( n++ != 0 ) printf( "+" );
580 case 0: printf( "TriggerOverLapError" ); break;
581 case 1: printf( "TriggerMissingError" ); break;
582 case 2: printf( "DataParityError" ); break;
583 case 3: printf( "ControlParityError" ); break;
584 case 4: printf( "TriggerInformationUnavailable" ); break;
585 case 5: printf( "FEEError" ); break;
586 case 6: printf( "HLTDecision" ); break;
587 case 7: printf( "HLTPayload" ); break;
588 case 8: printf( "DDGPayload" ); break;
589 default: printf( "spare %d", b );
596 printf( "MiniEventId(bunchCrossing):%d=0x%x ",
597 cdh->cdhMiniEventId, cdh->cdhMiniEventId );
601 printf( "Trigger classes: 0x(%05x-%08x)",
602 cdh->cdhTriggerClassesHigh,
603 cdh->cdhTriggerClassesLow );
604 if ( cdh->cdhTriggerClassesHigh != 0
605 || cdh->cdhTriggerClassesLow != 0 ) {
607 for ( b=0, n=0; b != 32; b++ ) {
608 if ( (cdh->cdhTriggerClassesLow & (1<<b)) != 0 ) {
609 if ( n++ != 0 ) printf( "+" );
613 for ( b=0; b != 18; b++ ) {
614 if ( (cdh->cdhTriggerClassesHigh & (1<<b)) != 0 ) {
615 if ( n++ != 0 ) printf( "+" );
616 printf( "%d", b+32 );
624 if ( cdh->cdhMBZ3 != 0 ) {
625 printf( "MBZ3:0x%x ",
628 printf( "ROI:0x(%08x-%01x)", cdh->cdhRoiHigh, cdh->cdhRoiLow );
629 if ( cdh->cdhRoiHigh != 0
630 || cdh->cdhRoiLow != 0 ) {
632 for ( b=0, n=0; b != 5; b++ ) {
633 if ( (cdh->cdhRoiLow & (1<<b)) != 0 ) {
634 if ( n++ != 0 ) printf( "+" );
638 for ( b=0; b != 32; b++ ) {
639 if ( (cdh->cdhRoiHigh & (1<<b)) != 0 ) {
640 if ( n++ != 0 ) printf( "+" );
652 printf( "%d) Payload \"%s\" already loaded at %p\n",
658 currEvent->payload = payload;
659 } /* End of loadPayload */
661 void parseEquipment( char * const line ) {
662 struct equipmentEventDescriptorStruct *equipment;
663 int payloadFound = FALSE;
668 (struct equipmentEventDescriptorStruct *)malloc( sizeof( *equipment ) )) == NULL ) {
670 "%s: filed to malloc for %d bytes (parseEquipment/equipmentEventDescriptorStruct) errno:%d ",
672 (int)sizeof( *equipment ),
677 currEvent = equipment;
680 while ( (keyword = strtok_r( p, " \t", &p )) != NULL ) {
681 DBG_VERBOSE printf( "%d) Equipment - Keyword:\"%s\"\n",
684 if ( strcasecmp( "id", keyword ) == 0 ) {
687 if ( (idNum = strtok_r( p, " \t", &p )) == NULL ) {
689 "%s: line:%d EQUIPMENT declaration, ID needed",
694 if ( sscanf( idNum, "%d", &currEquipmentId ) != 1 ) {
696 "%s: line:%d EQUIPMENT declaration, numeric ID needed (%s)",
702 DBG_VERBOSE printf( "%d) EQUIPMENT - ID:%d\n",
705 } else if ( strncasecmp( "pay", keyword, 3 ) == 0 ) {
708 if ( (fileName = strtok_r( p, " \t", &p )) == NULL ) {
710 "%s line:%d Payload without filename found\n",
715 DBG_VERBOSE printf( "%d) Equipment - Payload:\"%s\"\n",
718 if ( payloadFound ) {
720 "%s line:%d Payload with multiple filenames found\n",
725 loadPayload( fileName );
729 "%s: line:%d Equipment declaration, unknown keyword \"%s\"\n",
736 if ( !payloadFound ) {
738 "%s: line:%d Equipment without payload found\n",
744 equipment->id = currEquipmentId;
745 equipment->next = NULL;
746 if ( currLdc->head == NULL ) {
747 currLdc->head = currLdc->tail = equipment;
749 currLdc->tail->next = equipment;
750 currLdc->tail = equipment;
752 } /* End of parseEquipment */
754 void parseGdc( char * const line ) {
759 while ( (keyword = strtok_r( p, " \t", &p )) != NULL ) {
760 if ( strcasecmp( "id", keyword ) == 0 ) {
763 if ( (idNum = strtok_r( p, " \t", &p )) == NULL ) {
765 "%s: line:%d GDC declaration, ID needed",
770 if ( sscanf( idNum, "%d", (int*)&currGdcId ) != 1 ) {
772 "%s: line:%d GDC declaration, numeric ID needed (%s)",
778 DBG_VERBOSE printf( "%d) GDC - ID:%d\n",
781 } else if ( strcasecmp( "DetectorPattern", keyword ) == 0 ) {
784 if ( (detPattern = strtok_r( p, " \t", &p )) == NULL ) {
786 "%s: line:%d GDC declaration, DetectorPattern needed",
791 if ( sscanf( detPattern, "%u", &currDetPattern ) != 1 ) {
793 "%s: line:%d GDC declaration, numeric DetectorPattern needed (%s)",
799 DBG_VERBOSE printf( "%d) GDC - DetectorPattern:%u\n",
804 "%s: line:%d GDC declaration, unknown keyword \"%s\"\n",
811 } /* End of parseGdc */
813 void parseLdc( char * const line ) {
818 while ( (keyword = strtok_r( p, " \t", &p )) != NULL ) {
819 if ( strcasecmp( "id", keyword ) == 0 ) {
822 if ( (idNum = strtok_r( p, " \t", &p )) == NULL ) {
824 "%s: line:%d LDC declaration, ID needed",
829 if ( sscanf( idNum, "%d", (int*)&currLdcId ) != 1 ) {
831 "%s: line:%d LDC declaration, numeric ID needed (%s)",
837 DBG_VERBOSE printf( "%d) LDC - ID:%d\n",
842 "%s: line:%d LDC declaration, unknown keyword \"%s\"\n",
849 } /* End of parseLdc */
854 currLdcId = HOST_ID_MIN;
855 currGdcId = HOST_ID_MIN;
858 for ( lineNmb = 1; !feof( stdin ); lineNmb++ ) {
859 getLine( line, sizeof(line) );
860 if ( strlen(line) != 0 ) {
864 if ( (keyword = strtok_r( line, " \t", &p )) != NULL ) {
865 DBG_VERBOSE printf( "%d) Keyword:\"%s\"\n", lineNmb, keyword );
866 if ( strcasecmp( "gdc", keyword ) == 0 ) {
867 if ( workingAs != gdc && workingAs != unknown ) {
869 "%s: line:%d GDC found when working in non-GDC mode (e.g. as a LDC)\n",
876 currLdcId = HOST_ID_MIN;
879 } else if ( strcasecmp( "ldc", keyword ) == 0 ) {
880 if ( workingAs != gdc && workingAs != ldc && workingAs != unknown ) {
882 "%s: line:%d LDC found when working in non-LDC/GDC mode\n",
886 if ( workingAs == unknown ) workingAs = ldc;
888 if ( workingAs == ldc ) {
893 handleLdc( currLdcId );
897 } else if ( strncasecmp( "equ", keyword, 3 ) == 0 ) {
898 if ( workingAs == unknown
899 || (workingAs == ldc && currLdc == NULL )
900 || (workingAs == gdc && currGdc == NULL ) ) {
902 "%s: line:%d Unexpected EQUIPMENT declaration (LDC or GDC needed first)\n",
911 "%s: line:%d Parse error in \"%s\" unknown keyword\n",
919 } while ( !feof( stdin ) ) {}
923 printf( "End of parse: %d line%s found\n",
925 lineNmb != 1 ? "s" : "" );
926 printf( "Working as %s\n",
927 workingAs == gdc ? "GDC" :
928 workingAs == ldc ? "LDC" :
930 if ( workingAs == gdc ) {
931 struct ldcDescriptorStruct *ldc;
933 printf( "LDCs (%d):", numOfLdcs );
934 for ( ldc = ldcsHead; ldc != NULL; ldc = ldc->next ) {
935 printf( " %d", ldc->id );
942 if ( workingAs == ldc ) {
943 assert( ldcsHead == ldcsTail );
944 assert( ldcsTail == NULL );
947 if ( workingAs == gdc ) {
948 struct ldcDescriptorStruct *ldc;
950 assert( ldcsHead != NULL );
951 assert( ldcsTail != NULL );
952 assert( ldcsTail->next == NULL );
953 for ( ldc = ldcsHead; ldc->next != NULL; ldc = ldc->next ) {}
954 assert ( ldc == ldcsTail );
957 if ( workingAs == unknown ) {
958 DBG_VERBOSE printf( "Empty configuration: nothing to do!\n" );
962 assert( (eventsHead == NULL && eventsTail == NULL)
963 || (eventsHead != NULL && eventsTail != NULL) );
964 } /* End of parseRules */
966 void loadTimestamp( struct eventHeaderStruct * const ev ) {
969 if ( time( &t ) == (time_t)-1 ) {
971 "%s: failed to get system time errno:%d (%s)\n",
972 myName, errno, strerror( errno ) );
975 ev->eventTimestamp = (eventTimestampType)t;
976 } /* End of loadTimestamp */
978 void initEvent( struct eventHeaderStruct * const ev ) {
979 memset( ev, 0, sizeof( *ev ) );
981 ev->eventMagic = EVENT_MAGIC_NUMBER;
982 ev->eventHeadSize = EVENT_HEAD_BASE_SIZE;
983 ev->eventVersion = EVENT_CURRENT_VERSION;
984 ev->eventRunNb = currRunNb;
985 ZERO_EVENT_ID( ev->eventId );
986 ZERO_TRIGGER_PATTERN( ev->eventTriggerPattern );
987 ZERO_DETECTOR_PATTERN( ev->eventDetectorPattern );
988 RESET_ATTRIBUTES( ev->eventTypeAttribute );
989 if ( workingMode == collider )
990 SET_SYSTEM_ATTRIBUTE( ev->eventTypeAttribute, ATTR_ORBIT_BC );
991 ev->eventLdcId = VOID_ID;
992 ev->eventGdcId = VOID_ID;
994 } /* End of initEvent */
998 // Swap the endianess of the integer value 'x'
1000 return (((x & 0x000000ffU) << 24) | ((x & 0x0000ff00U) << 8) |
1001 ((x & 0x00ff0000U) >> 8) | ((x & 0xff000000U) >> 24));
1004 void outputEvent( const void * const ev,
1009 const long32 * const v = (long32 *)ev;
1010 printf( "Writing %d bytes @ %p (%d)\n", size, ev, *v );
1013 // .............................Test endianess..............................
1015 char* ptemp = (char*) &temp;
1017 if (ptemp[0]!=1) { // Mac platform: ptemp != 1..............................................................................
1018 int bufSize= size; if (bufSize > (int) sizeof(eventHeaderStruct)) { bufSize = sizeof(eventHeaderStruct); }
1019 char* evTemp = (char*) malloc (bufSize);
1020 memcpy(evTemp, ev, bufSize);
1022 if ((bufSize % sizeof(int)) != 0) {
1023 fprintf( stderr, "%s: size of the input buffer ev is not multiple of 4 (size = %d)\n", myName, bufSize);
1027 // Invert header to evTemp.....................................................
1028 int* buf = (int*) evTemp;
1029 for (int i=0; i < (int) (bufSize / sizeof(int)); i++, buf++) {
1030 int value = Swap(*buf);
1031 memcpy(evTemp + (i * sizeof(int)), &value, sizeof(int));
1034 // Write inverted header to file...............................................
1035 if ((done = fwrite( evTemp, bufSize, 1, outF )) != 1 ) {
1036 fprintf( stderr, "%s: failed to write inverted header. event size:%d bytes, errno:%d (%s)\n", myName, size, errno, strerror( errno ) );
1040 if (size > bufSize) { // Still theraw-data payload to write (but not inverted, since it is inverted eariler).............
1041 if ((done = fwrite( (char*)ev + bufSize, size - bufSize, 1, outF )) != 1 ) {
1042 fprintf( stderr, "%s: failed to write additional event size:%d bytes, errno:%d (%s)\n", myName, size, errno, strerror( errno ) );
1049 else { // Intel platform: ptemp == 1............................................................................
1050 if ((done = fwrite( ev, size, 1, outF )) != 1 ) {
1051 fprintf( stderr, "%s: failed to write event size:%d bytes, errno:%d (%s)\n", myName, size, errno, strerror( errno ) );
1055 } /* End of outputEvent */
1057 void createSorAndEor( const int sor ) {
1058 unsigned char event[ 1000 ];
1059 struct eventHeaderStruct *ev;
1060 struct eventHeaderStruct sev;
1062 assert( workingAs == ldc || workingAs == gdc );
1064 if ( !createSorEor ) return;
1065 ev = (struct eventHeaderStruct *)event;
1067 ev->eventSize = sizeof( event );
1068 ev->eventType = sor ? START_OF_RUN : END_OF_RUN;
1069 if ( workingMode == fixedTarget )
1070 LOAD_RAW_EVENT_ID( ev->eventId, 0, 0, 0 );
1072 LOAD_EVENT_ID( ev->eventId, 0, 0, 0 );
1073 SET_SYSTEM_ATTRIBUTE( ev->eventTypeAttribute, ATTR_P_START );
1075 if ( workingAs == ldc ) {
1076 currLdc = (struct ldcEventDescriptorStruct *)eventsHead;
1078 if ( workingAs == gdc ) {
1080 sev.eventGdcId = currGdcId;
1081 ev->eventGdcId = currGdcId;
1082 currGdc = (struct gdcEventDescriptorStruct *)eventsHead;
1083 currLdc = currGdc->head;
1085 ev->eventLdcId = currLdc->id;
1087 if ( workingAs == ldc ) {
1088 loadTimestamp( ev );
1089 outputEvent( ev, ev->eventSize );
1091 if ( workingAs == gdc ) {
1092 struct ldcDescriptorStruct *ldc;
1094 loadTimestamp( ev );
1096 sev.eventSize = sizeof( sev ) + numOfLdcs * ev->eventSize;
1097 sev.eventType = sor ? START_OF_RUN : END_OF_RUN ;
1098 COPY_EVENT_ID( ev->eventId, sev.eventId );
1099 COPY_SYSTEM_ATTRIBUTES( ev->eventTypeAttribute, sev.eventTypeAttribute );
1100 SET_SYSTEM_ATTRIBUTE( sev.eventTypeAttribute, ATTR_SUPER_EVENT );
1101 loadTimestamp( &sev );
1102 outputEvent( &sev, sizeof( sev ) );
1104 ev->eventGdcId = currGdcId;
1105 for ( ldc = ldcsHead; ldc != NULL; ldc = ldc->next ) {
1106 ev->eventLdcId = ldc->id;
1107 outputEvent( ev, ev->eventSize );
1111 ADD_EVENT_ID( ev->eventId, oneEventDelta );
1112 ev->eventSize = ev->eventSize / 2;
1113 ev->eventType = sor ? START_OF_RUN_FILES : END_OF_RUN_FILES;
1114 CLEAR_SYSTEM_ATTRIBUTE( ev->eventTypeAttribute, ATTR_P_START );
1115 if ( workingAs == ldc ) {
1116 loadTimestamp( ev );
1117 outputEvent( ev, ev->eventSize );
1119 if ( workingAs == gdc ) {
1120 struct ldcDescriptorStruct *ldc;
1122 loadTimestamp( ev );
1124 sev.eventSize = ev->eventSize;
1125 sev.eventType = sor ? START_OF_RUN_FILES : END_OF_RUN_FILES;
1126 COPY_EVENT_ID( ev->eventId, sev.eventId );
1127 COPY_SYSTEM_ATTRIBUTES( ev->eventTypeAttribute, sev.eventTypeAttribute );
1128 CLEAR_SYSTEM_ATTRIBUTE( sev.eventTypeAttribute, ATTR_SUPER_EVENT );
1129 outputEvent( &sev, sizeof( sev ) );
1130 outputEvent( ev, ev->eventSize - sizeof( sev ) );
1132 sev.eventSize = sizeof( sev ) + ev->eventSize;
1133 sev.eventType = sor ? START_OF_RUN_FILES : END_OF_RUN_FILES;
1134 COPY_EVENT_ID( ev->eventId, sev.eventId );
1135 COPY_SYSTEM_ATTRIBUTES( ev->eventTypeAttribute, sev.eventTypeAttribute );
1136 SET_SYSTEM_ATTRIBUTE( sev.eventTypeAttribute, ATTR_SUPER_EVENT );
1138 loadTimestamp( &sev );
1140 ev->eventGdcId = currGdcId;
1141 for ( ldc = ldcsHead; ldc != NULL; ldc = ldc->next ) {
1142 loadTimestamp( &sev );
1143 outputEvent( &sev, sizeof( sev ) );
1144 ev->eventLdcId = ldc->id;
1145 outputEvent( ev, ev->eventSize );
1149 ADD_EVENT_ID( ev->eventId, oneEventDelta );
1150 ev->eventSize = sizeof( *ev );
1151 ev->eventType = sor ? START_OF_RUN : END_OF_RUN;
1152 SET_SYSTEM_ATTRIBUTE( ev->eventTypeAttribute, ATTR_P_END );
1153 if ( workingAs == ldc ) {
1154 loadTimestamp( ev );
1155 outputEvent( ev, ev->eventSize );
1157 if ( workingAs == gdc ) {
1158 struct ldcDescriptorStruct *ldc;
1160 loadTimestamp( ev );
1162 sev.eventSize = sizeof( sev ) + numOfLdcs * ev->eventSize;
1163 sev.eventType = sor ? START_OF_RUN : END_OF_RUN;
1164 COPY_EVENT_ID( ev->eventId, sev.eventId );
1165 COPY_SYSTEM_ATTRIBUTES( ev->eventTypeAttribute, sev.eventTypeAttribute );
1166 SET_SYSTEM_ATTRIBUTE( sev.eventTypeAttribute, ATTR_SUPER_EVENT );
1167 loadTimestamp( &sev );
1169 outputEvent( &sev, sizeof( sev ) );
1171 for ( ldc = ldcsHead; ldc != NULL; ldc = ldc->next ) {
1172 ev->eventLdcId = ldc->id;
1173 outputEvent( ev, ev->eventSize );
1176 } /* End of createSorEor */
1179 createSorAndEor( TRUE );
1180 } /* End of createSor */
1183 createSorAndEor( FALSE );
1184 } /* End of createEor */
1186 void loadCdh( struct commonDataHeaderStruct * const cdh,
1187 eventIdType * const eventId,
1188 equipmentIdType id ) {
1189 if ( !handleCDH ) return;
1191 // CTP raw-data does not contain CDH
1192 if ( id == 4352) return;
1194 if ( gotAliceTrigger ) {
1195 cdh->cdhEventId1 = EVENT_ID_GET_BUNCH_CROSSING( *eventId );
1196 cdh->cdhEventId2 = EVENT_ID_GET_ORBIT( *eventId );
1198 cdh->cdhEventId1 = 0;
1199 cdh->cdhEventId2 = EVENT_ID_GET_NB_IN_RUN( *eventId );
1201 cdh->cdhMiniEventId = cdh->cdhEventId1;
1203 void decodeCDH( struct ldcEventDescriptorStruct * const ldc,
1204 const struct payloadDescriptorStruct * const payloadDesc,
1205 equipmentIdType id );
1207 void createEvent( void ) {
1208 assert( workingAs == ldc || workingAs == gdc );
1210 /* Step 1: load all buffers (if needed) and compose the GDC/LDC headers */
1211 if ( workingAs == gdc ) {
1212 struct ldcEventDescriptorStruct *ldc;
1214 for( ldc = currGdc->head; ldc != NULL; ldc = ldc->next ) {
1215 COPY_EVENT_ID( currEventId, ldc->header.eventId );
1216 loadTimestamp( &ldc->header );
1218 COPY_EVENT_ID( currEventId, currGdc->header.eventId );
1219 loadTimestamp( &currGdc->header );
1221 for( ldc = currGdc->head; ldc != NULL; ldc = ldc->next ) {
1222 struct equipmentEventDescriptorStruct *eq;
1225 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
1226 if ( !bufferData ) {
1227 loadBuffer( eq->payload );
1228 decodeCDH( ldc, eq->payload, eq->id );
1230 loadCdh( (struct commonDataHeaderStruct*)eq->payload->data,
1235 if ( !currGdc->loaded ) {
1236 for ( n = 0; n != EVENT_TRIGGER_PATTERN_WORDS; n++ )
1237 currGdc->header.eventTriggerPattern[n] |= ldc->header.eventTriggerPattern[n];
1238 for ( n = 0; n != EVENT_DETECTOR_PATTERN_WORDS; n++ )
1239 currGdc->header.eventDetectorPattern[n] |= ldc->header.eventDetectorPattern[n];
1240 for ( n = 0; n != ALL_ATTRIBUTE_WORDS; n++ )
1241 currGdc->header.eventTypeAttribute[n] |= ldc->header.eventTypeAttribute[n];
1242 currGdc->loaded = TRUE;
1246 } else if ( workingAs == ldc ) {
1247 struct equipmentEventDescriptorStruct *eq;
1249 COPY_EVENT_ID( currEventId, currLdc->header.eventId );
1250 loadTimestamp( &currLdc->header );
1252 for ( eq = currLdc->head; eq != NULL; eq = eq->next ) {
1253 if ( !bufferData ) {
1254 loadBuffer( eq->payload );
1255 decodeCDH( currLdc, eq->payload, eq->id );
1257 loadCdh( (struct commonDataHeaderStruct*)eq->payload->data,
1260 currLdc->loaded = TRUE;
1264 ADD_EVENT_ID( currEventId, oneEventDelta );
1266 /* Step 2: output the event */
1267 if ( workingAs == gdc ) {
1268 struct ldcEventDescriptorStruct *ldc;
1270 outputEvent( &currGdc->header, sizeof( currGdc->header ) );
1272 for( ldc = currGdc->head; ldc != NULL; ldc = ldc->next ) {
1273 struct equipmentEventDescriptorStruct *eq;
1275 outputEvent( &ldc->header, sizeof( ldc->header ) );
1277 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
1278 outputEvent( &eq->header, sizeof( eq->header ) );
1279 outputEvent( eq->payload->data, eq->payload->size );
1280 if ( !bufferData ) unloadBuffer( eq->payload );
1283 if ( (currGdc = currGdc->next) == NULL )
1284 currGdc = (struct gdcEventDescriptorStruct *)eventsHead;
1285 } else if ( workingAs == ldc ) {
1286 struct equipmentEventDescriptorStruct *eq;
1288 outputEvent( &currLdc->header, sizeof( currLdc->header ) );
1290 for ( eq = currLdc->head; eq != NULL; eq = eq->next ) {
1291 outputEvent( &eq->header, sizeof( eq->header ) );
1292 outputEvent( eq->payload->data, eq->payload->size );
1293 if ( !bufferData ) unloadBuffer( eq->payload );
1295 if ( (currLdc = currLdc->next) == NULL )
1296 currLdc = (struct ldcEventDescriptorStruct *)eventsHead;
1298 } /* End of createEvent */
1300 void createEvents() {
1303 currGdc = (struct gdcEventDescriptorStruct *)eventsHead;
1304 currLdc = (struct ldcEventDescriptorStruct *)eventsHead;
1309 eventNum != numOfEvents && numOfEvents != 0;
1314 } /* End of createEvents */
1318 "Usage: %s [-?][-d][-i definitionFile][-o outputFile][-# numOfEvents][-s][-F|-C]\n\
1320 -v Print version ID and exit\n\
1321 -d Enable debug (repeat for more verbosity)\n\
1322 -i definitionFile File with the description of the events to create (default: stdin)\n\
1323 -o outputFile File used to store events (default: stdout)\n\
1324 -# numOfEvents Number of events to generate (default: 1 event)\n\
1325 -s Do not generate SOR/EOR files (valid only for GDCs)\n\
1326 -F/-C Working in Fixed Target (F) or Collider (C) mode\n\
1328 -D Direct disc access (no buffering)\n",
1331 } /* End of usage */
1333 void parseArgs( int argc, char **argv ) {
1335 int inFileName = -1;
1336 int outFileName = -1;
1339 while ( arg < argc ) {
1340 if ( strcmp( "-?", argv[ arg ] ) == 0 ) {
1344 if ( strcmp( "-i", argv[ arg ] ) == 0 ) {
1345 if ( ++arg == argc ) exit( usage() );
1347 if ( freopen( argv[arg], "r", stdin ) == NULL ){
1349 "%s: failed to open input definition \"%s\" errno:%d ",
1350 myName, argv[arg], errno );
1354 } else if ( strcmp( "-v", argv[ arg ] ) == 0 ) {
1355 printf( "%s\n", fileHandlerIdent );
1357 } else if ( strcmp( "-o", argv[ arg ] ) == 0 ) {
1358 if ( ++arg == argc ) exit( usage() );
1360 } else if ( strcmp( "-#", argv[ arg ] ) == 0 ) {
1363 if ( ++arg == argc ) exit( usage() );
1364 if ( sscanf( argv[ arg ], "%d", &n ) != 1 ) exit( usage() );
1365 if ( n < 0 ) exit( usage() );
1367 } else if ( strcmp( "-s", argv[ arg ] ) == 0 ) {
1368 createSorEor = FALSE;
1369 } else if ( strcmp( "-F", argv[ arg ] ) == 0 ) {
1370 workingMode = fixedTarget;
1371 } else if ( strcmp( "-C", argv[ arg ] ) == 0 ) {
1372 workingMode = collider;
1373 } else if ( strcmp( "-d", argv[ arg ] ) == 0 ) {
1375 } else if ( strcmp( "-c", argv[ arg ] ) == 0 ) {
1377 } else if ( strcmp( "-D", argv[ arg ] ) == 0 ) {
1379 } else if ( strcmp( "-run", argv[ arg ] ) == 0 ) {
1381 if ( ++arg == argc ) exit( usage() );
1382 if ( sscanf( argv[ arg ], "%d", &runnumber ) != 1 ) exit( usage() );
1383 if ( runnumber < 0 ) exit( usage() );
1384 currRunNb = runnumber;
1386 fprintf( stderr, "%s: Unknown switch \"%s\"\n", myName, argv[argc] );
1392 if ( workingMode == fixedTarget )
1393 LOAD_RAW_EVENT_ID( oneEventDelta, 1, 0, 1 );
1395 LOAD_EVENT_ID( oneEventDelta, 0, 0, 1 );
1396 ZERO_EVENT_ID( currEventId );
1399 printf( "Configuration:\n" );
1400 printf( " Debug level: %d\n", debug );
1401 printf( " Configuration: %s\n",
1402 inFileName == -1 ? "stdin" : argv[ inFileName ] );
1403 printf( " Output: %s\n",
1404 outFileName == -1 ? "stdout" : argv[ outFileName ] );
1405 printf( " Working mode: %s\n",
1406 workingMode == fixedTarget ? "fixed target" : "collider" );
1407 printf( " Number of events: %d\n", numOfEvents );
1408 printf( " %s SOR/EOR files\n",
1409 createSorEor ? "Create" : "Do not create" );
1410 printf( " CDH handling: %s\n",
1411 handleCDH ? "enabled" : "disabled" );
1412 printf( " data buffering: %s\n",
1413 bufferData ? "enabled" : "DISABLED" );
1416 if ( outFileName == -1 ) {
1418 printf( "No more trace information from this point...\n" );
1422 if ( ( outF = fopen( argv[ outFileName ], "w" ) ) == NULL ) {
1424 "%s: failed to open output file \"%s\" for writing, errno:%d (%s)\n",
1426 argv[ outFileName ],
1428 strerror( errno ) );
1432 printf( "Output file \"%s\" opened OK for writing\n",
1433 argv[ outFileName ] );
1435 } /* End of parseArgs */
1436 void initEquipment( struct equipmentHeaderStruct * const eq ) {
1437 memset( eq, 0, sizeof( *eq ) );
1438 RESET_ATTRIBUTES( eq->equipmentTypeAttribute );
1439 eq->equipmentBasicElementSize = 4;
1440 } /* End of initEquipment */
1442 void decodeCDH( struct ldcEventDescriptorStruct * const ldc,
1443 const struct payloadDescriptorStruct * const payloadDesc,
1444 equipmentIdType id ) {
1447 struct commonDataHeaderStruct *cdh;
1448 static int softwareTriggerIndicator = FALSE;
1452 if ( payloadDesc->size < CDH_SIZE ) {
1454 "%s: payload too small got:%d CDH:%d\n",
1460 if ( (cdh = (struct commonDataHeaderStruct *)payloadDesc->data) != NULL ) {
1461 if ( cdh->cdhVersion != CDH_VERSION ) {
1463 "%s: CDH version mismatch expected:%d got:%d\n",
1469 if ( cdhRef == NULL ) {
1471 #define CDH_TRIGGER_INFORMATION_UNAVAILABLE_MASK (1<<CDH_TRIGGER_INFORMATION_UNAVAILABLE_BIT)
1472 gotAliceTrigger = (cdhRef->cdhStatusErrorBits & CDH_TRIGGER_INFORMATION_UNAVAILABLE_MASK) == 0;
1473 if ( gotAliceTrigger && workingMode == fixedTarget ) {
1475 "%s: ALICE trigger and fixed target mode are not compatible.\n\
1476 Either work in Collider mode or set the trigger unavailable status bit in the CDH.\n",
1480 if ( gotAliceTrigger ) {
1481 if ( (cdh->cdhL1TriggerMessage & 0x40) != 0 ) {
1483 "%s: CDH is a calibration trigger (unsupported) L1TriggerMessage:0x%x\n",
1484 myName, cdh->cdhL1TriggerMessage );
1487 if ( (cdh->cdhL1TriggerMessage & 0x01) != 0 ) {
1488 softwareTriggerIndicator = TRUE;
1490 if ( softwareTriggerIndicator ) {
1491 switch ((cdh->cdhL1TriggerMessage >> 2) & 0xF) {
1499 /* L1SwC bit = on, Clt bit = off, RoC[4..1] = 0xF --> END_OF_DATA */
1501 /* L1SwC bit = on, Clt bit = off, RoC[4..1] = 0xE0 --> START_OF_DATA */
1503 /* L1SwC bit = on, Clt bit = off, RoC[4] = 1, but not 0xE or 0xF
1504 --> SYSTEM_SOFTWARE_TRIGGER_EVENT */
1506 /* L1SwC bit = on, Clt bit = off, RoC[4] = 0
1507 --> DETECTOR_SOFTWARE_TRIGGER_EVENT */
1509 "%s: CDH trigger SOD/EOD/SST/DST (unsupported) \
1510 L1TriggerMessage:0x%x ALICETrigger:%s\n",
1512 cdh->cdhL1TriggerMessage,
1513 gotAliceTrigger ? "yes" : "no" );
1519 if ( (cdh->cdhStatusErrorBits & CDH_TRIGGER_INFORMATION_UNAVAILABLE_MASK) !=
1520 (cdhRef->cdhStatusErrorBits & CDH_TRIGGER_INFORMATION_UNAVAILABLE_MASK) ) {
1522 "%s: CDH coherency check failed. \
1523 Trigger information reference:%savailable current:%savailable\n",
1525 (cdhRef->cdhStatusErrorBits & CDH_TRIGGER_INFORMATION_UNAVAILABLE_MASK) == 0 ? "UN" : "",
1526 (cdh->cdhStatusErrorBits & CDH_TRIGGER_INFORMATION_UNAVAILABLE_MASK) == 0 ? "UN" : "" );
1529 if ( gotAliceTrigger ) {
1530 if ( cdhRef->cdhL1TriggerMessage != cdh->cdhL1TriggerMessage ) {
1532 "%s: CDH coherency check failed. \
1533 L1 trigger message reference:0x%x current:0x%x\n",
1535 cdhRef->cdhL1TriggerMessage,
1536 cdh->cdhL1TriggerMessage );
1539 if ( cdh->cdhParticipatingSubDetectors != cdhRef->cdhParticipatingSubDetectors ) {
1541 "%s: CDH coherency check failed. \
1542 ParticipatingSubDetectors reference:0x%x current:0x%x\n",
1544 cdhRef->cdhParticipatingSubDetectors,
1545 cdh->cdhParticipatingSubDetectors );
1548 if ( cdh->cdhTriggerClassesLow != cdhRef->cdhTriggerClassesLow
1549 || cdh->cdhTriggerClassesHigh != cdhRef->cdhTriggerClassesHigh ) {
1551 "%s: CDH coherency check failed. \
1552 TriggerClassesHigh/Low reference:0x%x-%x current:0x%x-%x\n",
1554 cdhRef->cdhTriggerClassesHigh, cdhRef->cdhTriggerClassesLow,
1555 cdh ->cdhTriggerClassesHigh, cdh ->cdhTriggerClassesLow );
1558 if ( cdh->cdhBlockLength != 0xffffffff ) {
1559 if ( (unsigned)payloadDesc->size != cdh->cdhBlockLength ) {
1561 "%s: CDH coherency check failed. \
1562 Payload size:%d (0x%08x) CDH block length:%d (0x%08x)\n",
1564 payloadDesc->size, payloadDesc->size,
1565 cdh->cdhBlockLength, cdh->cdhBlockLength );
1569 if ( cdh->cdhRoiLow != cdhRef->cdhRoiLow
1570 || cdh->cdhRoiHigh != cdhRef->cdhRoiHigh ) {
1572 "%s: CDH coherency check failed. \
1573 RoiHigh/Low reference:0x%x-%x current:0x%x-%x\n",
1575 cdhRef->cdhRoiHigh, cdhRef->cdhRoiLow,
1576 cdh ->cdhRoiHigh, cdh ->cdhRoiLow );
1580 if ( cdh->cdhMBZ0 != 0
1581 || cdh->cdhMBZ1 != 0
1582 || cdh->cdhMBZ2 != 0
1583 || cdh->cdhMBZ3 != 0 ) {
1585 "%s: CDH check failed. MBZ0:0x%x MBZ1:0x%x MBZ2:0x%x MBZ3:0x%x\n",
1587 cdh->cdhMBZ0, cdh->cdhMBZ1, cdh->cdhMBZ2, cdh->cdhMBZ3 );
1591 for ( attr = 0; attr != 8; attr++ ) {
1592 if ( (cdh->cdhBlockAttributes & (1<<attr)) != 0 ) {
1593 SET_USER_ATTRIBUTE( ldc->header.eventTypeAttribute, attr );
1596 for ( trig = 0; trig != 32; trig++ ) {
1597 if ( (cdh->cdhTriggerClassesLow & (1<<trig)) != 0 ) {
1598 SET_TRIGGER_IN_PATTERN( ldc->header.eventTriggerPattern,
1602 for ( trig = 0; trig != 18; trig++ ) {
1603 if ( (cdh->cdhTriggerClassesHigh & (1<<trig)) != 0 ) {
1604 SET_TRIGGER_IN_PATTERN( ldc->header.eventTriggerPattern,
1608 if ( gotAliceTrigger )
1609 VALIDATE_TRIGGER_PATTERN( ldc->header.eventTriggerPattern );
1612 } /* End of decodeCDH */
1615 assert( workingAs == ldc || workingAs == gdc );
1617 if ( workingAs == gdc ) {
1618 struct gdcEventDescriptorStruct *gdc;
1620 for ( gdc = (struct gdcEventDescriptorStruct *)eventsHead;
1623 struct ldcEventDescriptorStruct *ldc;
1625 initEvent( &gdc->header );
1626 gdc->header.eventSize = gdc->header.eventHeadSize;
1627 gdc->header.eventType = PHYSICS_EVENT;
1628 SET_SYSTEM_ATTRIBUTE( gdc->header.eventTypeAttribute, ATTR_SUPER_EVENT );
1629 gdc->header.eventGdcId = currGdcId;
1630 COPY_DETECTOR_PATTERN(&currDetPattern, gdc->header.eventDetectorPattern);
1631 for ( ldc = gdc->head; ldc != NULL; ldc = ldc->next ) {
1632 struct equipmentEventDescriptorStruct *eq;
1634 initEvent( &ldc->header );
1635 ldc->header.eventSize = ldc->header.eventHeadSize;
1636 ldc->header.eventType = PHYSICS_EVENT;
1637 ldc->header.eventGdcId = currGdcId;
1638 COPY_DETECTOR_PATTERN(&currDetPattern, ldc->header.eventDetectorPattern);
1639 ldc->header.eventLdcId = ldc->id;
1640 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
1641 initEquipment( &eq->header );
1642 eq->header.equipmentId = eq->id;
1643 if ( workingMode == collider )
1644 SET_SYSTEM_ATTRIBUTE( eq->header.equipmentTypeAttribute,
1646 eq->header.equipmentSize = eq->payload->size + sizeof( eq->header );
1647 ldc->header.eventSize += eq->header.equipmentSize;
1648 decodeCDH( ldc, eq->payload, eq->id );
1649 OR_ALL_ATTRIBUTES( eq->header.equipmentTypeAttribute,
1650 ldc->header.eventTypeAttribute );
1651 OR_ALL_ATTRIBUTES( eq->header.equipmentTypeAttribute,
1652 gdc->header.eventTypeAttribute );
1654 gdc->header.eventSize += ldc->header.eventSize;
1660 printf( "Headers:\n" );
1661 for ( gdc = (struct gdcEventDescriptorStruct *)eventsHead;
1664 struct ldcEventDescriptorStruct *ldc;
1666 printf( " GDC:%d size:%d vers:%08x\n",
1668 gdc->header.eventSize,
1669 gdc->header.eventVersion);
1670 for ( ldc = gdc->head; ldc != NULL; ldc = ldc->next ) {
1671 struct equipmentEventDescriptorStruct *eq;
1673 printf( " LDC:%d size:%d vers:%08x\n",
1675 ldc->header.eventSize,
1676 ldc->header.eventVersion );
1677 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
1678 printf( " EQ:%d size:%d %spayload:%d\n",
1680 eq->header.equipmentSize,
1681 eq->header.equipmentSize - sizeof( struct equipmentHeaderStruct ) == (unsigned)eq->payload->size ? "" : "-ERROR",
1682 eq->payload->size );
1687 } else if ( workingAs == ldc ) {
1688 struct ldcEventDescriptorStruct *ldc;
1690 for ( ldc = (struct ldcEventDescriptorStruct *)eventsHead;
1693 struct equipmentEventDescriptorStruct *eq;
1695 initEvent( &ldc->header );
1696 ldc->header.eventSize = ldc->header.eventHeadSize;
1697 ldc->header.eventType = PHYSICS_EVENT;
1698 ldc->header.eventGdcId = VOID_ID;
1699 ldc->header.eventLdcId = ldc->id;
1700 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
1701 initEquipment( &eq->header );
1702 eq->header.equipmentId = eq->id;
1703 if ( workingMode == collider )
1704 SET_SYSTEM_ATTRIBUTE( eq->header.equipmentTypeAttribute,
1706 eq->header.equipmentSize = eq->payload->size + sizeof( eq->header );
1707 ldc->header.eventSize += eq->header.equipmentSize;
1708 decodeCDH( ldc, eq->payload, eq->id );
1709 OR_ALL_ATTRIBUTES( eq->header.equipmentTypeAttribute,
1710 ldc->header.eventTypeAttribute );
1715 printf( "Headers:\n" );
1716 for ( ldc = (struct ldcEventDescriptorStruct *)eventsHead;
1719 struct equipmentEventDescriptorStruct *eq;
1721 printf( " LDC:%d size:%d vers:%08x\n",
1723 ldc->header.eventSize,
1724 ldc->header.eventVersion );
1725 for ( eq = ldc->head; eq != NULL; eq = eq->next ) {
1726 printf( " EQ:%d size:%d %spayload:%d\n",
1728 eq->header.equipmentSize,
1729 eq->header.equipmentSize - sizeof( struct equipmentHeaderStruct ) == (unsigned)eq->payload->size ? "" : "-ERROR",
1730 eq->payload->size );
1735 } /* End of initEvents */
1739 workingAs = unknown;
1740 workingMode = fixedTarget;
1741 ldcsHead = ldcsTail = NULL;
1742 eventsHead = eventsTail = NULL;
1746 payloadsHead = payloadsTail = NULL;
1750 createSorEor = TRUE;
1752 gotAliceTrigger = TRUE;
1754 } /* End of initVars */
1756 int main( int argc, char **argv ) {
1758 parseArgs( argc, argv );