yoda/techmodel/McmdReader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v $
 * $Id: McmdReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the McmdReader class.
 */

#define BYTE unsigned char
#include <string>
#include <log4cpp/Category.hh>
extern "C" {
#include <sys/time.h>
#include "CRC.h"
}

#include <fstream>
#include "stdio.h"
#include "ReaderAlgorithms.h"

#include "event/mcmd/McmdRecord.h"

using namespace pamela;
using namespace pamela::techmodel;

static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.McmdReader");  

/**
 * Constructor. 
 */
McmdReader::McmdReader(void): 
  TechmodelAlgorithm(PacketType::Mcmd, "Mcmd") { 
  cat <<  log4cpp::Priority::DEBUG
      <<  "Constructor "
      <<  "\n " << log4cpp::CategoryStream::ENDLINE;
  Mcmd = new McmdEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string McmdReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $\n";
}

/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 */ 
void McmdReader::Init(PamelaRun *run) {
  SetInputStream(run);
  run->WriteSubPacket(this, &Mcmd, Mcmd->Class());
}

/**
 * Unpack the Mcmd event from an input file.
 *  Each subpacket is prafaceded by 4 OBT bytes.
 *  The elementar structure is a kind of 
 * --------CPU - OBT---------------------------------
 *  OBT     - 4 Bytes
 * --------Start Sub-Packet---------------------------------
 *  SeqID  - 2 Bytes
 *  T          - 1 bit       | spare  - 7 bits | Mcmd-ID - 1 byte
 *  Spare  - 4 bits     | Mcmd Lenght 12 bits
 *  Time Tag      - 4 Bytes
 *  DATA            - McmdLength * 2 Bytes
 *  don't care    - 1 Byte
 *  END-ID         - 1 Byte is the sum of subPacket Bytes module 256
 * --------End Sub-Packet---------------------------------
 * --------CPU CRC-16 on OBT + subPacket---------------------------------
 * subCRC        - 2 Bytes
 * this structure is repeated one or more times.
 * The last 2 Bytes are CRC16 computed by CPU on all previous (OBT + subpacket data + subCRC) repetitions

 * @Event Number
 * @length is the size in bytes of the event (or packet)
 */
void McmdReader::RunEvent(int EventNumber, long int length) {
  
  int i = 0;
  char OBT[4];
  char subHeader[10];
  char *subData;
  char subTrailer[2];
  char subPckCRC[2];
  char eventCRC[2];
  long int start;
  unsigned short headerCRC;
  UINT16    subCRC;      //CRC of the subpacket (updated as subPckt is read)
  UINT16    readCRC;     //CRC read from the subpacket
  UINT16    partialCRC;  //partial CRC updated as mcmd packet is read (to compare with CRC read in the total Mcmd header)

  long int lastPosition, dataLength;
  McmdRecord *rec;

  //the 2 bytes subtracted belong to the final event CRC bytes
  lastPosition = (long int)InputFile->tellg() + length - 2;
  Mcmd->Records->Clear();
  TClonesArray &recs = *(Mcmd->Records);
  partialCRC = 0;
  while(InputFile->tellg() < lastPosition) {

//Just to test crc
/*  char tempbuff[length];
  InputFile->read(tempbuff, sizeof(tempbuff));
  subCRC = 0;
  for (int jj = 0; jj < length ; jj++){
    subCRC = CM_Compute_CRC16(subCRC, (BYTE*)(tempbuff+jj), 1);
  }*/

    //Read the OBT preceeding the subpacket then calculate a partial CRC for it
    //and update the partialCRC
    start = InputFile->tellg();
    InputFile->read(OBT, sizeof(OBT));
    subCRC = CM_Compute_CRC16(0, (BYTE*)&OBT, 4);
    partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&OBT, 4);
    
    //Read the header for the subPacket and read mcmdLength 
    //12 is the total lenght of subHeader + subTrailer
    InputFile->read(subHeader, sizeof(subHeader));
    subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subHeader, 10);
    partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subHeader, 10);
    dataLength = (((0x0fff)&((UINT16)subHeader[4]))<<8|((UINT16)subHeader[5])*2) - 12;

    if (dataLength < 0) break; //it should throw an exception ***TBD***
    //read subpacket data according to data length then calculate partial CRC for data
    //and update the partialCRC
    
    // If dataLength > 0 reads the data and compute again a partial CRC
    // else set subData to 0
    //if(dataLength > 0){
        subData = new char[dataLength];
        InputFile->read(subData, sizeof(unsigned char)*dataLength);
        for (int jj = 0; jj < dataLength ; jj++){
            subCRC = CM_Compute_CRC16(subCRC, (BYTE*)(subData+jj), 1);
            partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)(subData+jj), 1);
        }
        //subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subData, dataLength);

    //} else {
    //    subData = "0";
    //}
    
    //Read the CRC inside of MCMD
    //To check sum of all sub packet data, module 256 ****TBD****
    InputFile->read(subTrailer, sizeof(subTrailer));
    subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subTrailer, 2);
    partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subTrailer, 2);

    //Read the CRC of subPacket
    InputFile->read(subPckCRC, sizeof(subPckCRC));
    partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subPckCRC, 2);
    readCRC = (((BYTE)subPckCRC[0])<<8) + ((BYTE)subPckCRC[1]);
    
    //finally check if the RecordCRC is correct
    //and finally update the partialCRC
    //TO DO - if one CRC is wrong also the total one will be corrupted
    if (subCRC == readCRC){
        rec                  = new(recs[i++]) McmdRecord(); //aggiungo un nuovo McmdRecord all'evento
        rec->MCMD_RECORD_OBT = (((UINT32)OBT[0]<<24)&0xFF000000) + (((UINT32)OBT[1]<<16)&0x00FF0000) +  (((UINT32)OBT[2]<<8)&0x0000FF00) + (((UINT32)OBT[3])&0x000000FF);
        rec->SeqID           = (((BYTE)subHeader[0]<<8)&0xFF00) + (((BYTE)subHeader[1])&0x00FF);
        rec->Tbit            = ((BYTE)((subHeader[2]&0x80))>>7);
        rec->ID1             = (BYTE)subHeader[3];
        rec->McmdLength      = (0x0fff)&(((UINT16)(subHeader[4]<<8)) + ((UINT16)subHeader[5]));
        rec->TimeTag         = (((UINT32)OBT[6]<<24)&0xFF000000) + (((UINT32)OBT[7]<<16)&0x00FF0000) +  (((UINT32)OBT[8]<<8)&0x0000FF00) + (((UINT32)OBT[9])&0x000000FF);
        rec->endID           = (BYTE)subTrailer[1];
        rec->setMcmdData((BYTE*)&subData);
        delete [] subData;
    } else {
        cat <<  log4cpp::Priority::ERROR
            <<  "Wrong CRC on Subpacket int MCMD Packet starting at position " << start
            <<  "\n " << log4cpp::CategoryStream::ENDLINE;
    }
  }
    //in the end compare the calculated partial CRC with the MCMD packet CRC
    InputFile->read(eventCRC, sizeof(eventCRC));
    readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);
    
    if(!(partialCRC == readCRC)) {
           cat <<  log4cpp::Priority::ERROR
            <<  "Wrong CRC for CalibS4 Packet"
            <<  "\n " << log4cpp::CategoryStream::ENDLINE;
    }
}


1	kusanagi	1.1	/** @file
2	kusanagi	1.2	* $Source: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v $
3	kusanagi	1.3	* $Id: McmdReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $
4	kusanagi	1.2	* $Author: kusanagi $
5	kusanagi	1.1	*
6			* Implementation of the McmdReader class.
7			*/
8
9			#define BYTE unsigned char
10			#include <string>
11			#include <log4cpp/Category.hh>
12			extern "C" {
13			#include <sys/time.h>
14			#include "CRC.h"
15			}
16
17			#include <fstream>
18			#include "stdio.h"
19			#include "ReaderAlgorithms.h"
20
21			#include "event/mcmd/McmdRecord.h"
22
23			using namespace pamela;
24			using namespace pamela::techmodel;
25
26			static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.McmdReader");
27
28			/**
29			* Constructor.
30			*/
31			McmdReader::McmdReader(void):
32			TechmodelAlgorithm(PacketType::Mcmd, "Mcmd") {
33			cat << log4cpp::Priority::DEBUG
34			<< "Constructor "
35			<< "\n " << log4cpp::CategoryStream::ENDLINE;
36			Mcmd = new McmdEvent();
37			}
38
39			/**
40			* Get a string with the version info of the algorithm.
41			*/
42			std::string McmdReader::GetVersionInfo(void) const {
43			return
44	kusanagi	1.3	"$Header: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $\n";
45	kusanagi	1.1	}
46
47			/**
48			* Initialize the algorithm with a special run. This will initialize the
49			* event reader routines for all packet types.
50			*/
51			void McmdReader::Init(PamelaRun *run) {
52			SetInputStream(run);
53			run->WriteSubPacket(this, &Mcmd, Mcmd->Class());
54			}
55
56			/**
57			* Unpack the Mcmd event from an input file.
58			* Each subpacket is prafaceded by 4 OBT bytes.
59			* The elementar structure is a kind of
60			* --------CPU - OBT---------------------------------
61			* OBT - 4 Bytes
62			* --------Start Sub-Packet---------------------------------
63			* SeqID - 2 Bytes
64			* T - 1 bit \| spare - 7 bits \| Mcmd-ID - 1 byte
65			* Spare - 4 bits \| Mcmd Lenght 12 bits
66			* Time Tag - 4 Bytes
67			* DATA - McmdLength * 2 Bytes
68			* don't care - 1 Byte
69			* END-ID - 1 Byte is the sum of subPacket Bytes module 256
70			* --------End Sub-Packet---------------------------------
71			* --------CPU CRC-16 on OBT + subPacket---------------------------------
72			* subCRC - 2 Bytes
73			* this structure is repeated one or more times.
74			* The last 2 Bytes are CRC16 computed by CPU on all previous (OBT + subpacket data + subCRC) repetitions
75
76			* @Event Number
77			* @length is the size in bytes of the event (or packet)
78			*/
79			void McmdReader::RunEvent(int EventNumber, long int length) {
80
81			int i = 0;
82			char OBT[4];
83			char subHeader[10];
84			char *subData;
85			char subTrailer[2];
86			char subPckCRC[2];
87			char eventCRC[2];
88			long int start;
89			unsigned short headerCRC;
90			UINT16 subCRC; //CRC of the subpacket (updated as subPckt is read)
91			UINT16 readCRC; //CRC read from the subpacket
92			UINT16 partialCRC; //partial CRC updated as mcmd packet is read (to compare with CRC read in the total Mcmd header)
93
94			long int lastPosition, dataLength;
95			McmdRecord *rec;
96
97			//the 2 bytes subtracted belong to the final event CRC bytes
98			lastPosition = (long int)InputFile->tellg() + length - 2;
99			Mcmd->Records->Clear();
100			TClonesArray &recs = *(Mcmd->Records);
101			partialCRC = 0;
102			while(InputFile->tellg() < lastPosition) {
103
104			//Just to test crc
105			/* char tempbuff[length];
106			InputFile->read(tempbuff, sizeof(tempbuff));
107			subCRC = 0;
108			for (int jj = 0; jj < length ; jj++){
109			subCRC = CM_Compute_CRC16(subCRC, (BYTE*)(tempbuff+jj), 1);
110			}*/
111
112			//Read the OBT preceeding the subpacket then calculate a partial CRC for it
113			//and update the partialCRC
114			start = InputFile->tellg();
115			InputFile->read(OBT, sizeof(OBT));
116			subCRC = CM_Compute_CRC16(0, (BYTE*)&OBT, 4);
117			partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&OBT, 4);
118
119			//Read the header for the subPacket and read mcmdLength
120			//12 is the total lenght of subHeader + subTrailer
121			InputFile->read(subHeader, sizeof(subHeader));
122			subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subHeader, 10);
123			partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subHeader, 10);
124			dataLength = (((0x0fff)&((UINT16)subHeader[4]))<<8\|((UINT16)subHeader[5])*2) - 12;
125
126			if (dataLength < 0) break; //it should throw an exception *TBD*
127			//read subpacket data according to data length then calculate partial CRC for data
128			//and update the partialCRC
129
130			// If dataLength > 0 reads the data and compute again a partial CRC
131			// else set subData to 0
132			//if(dataLength > 0){
133			subData = new char[dataLength];
134			InputFile->read(subData, sizeof(unsigned char)*dataLength);
135			for (int jj = 0; jj < dataLength ; jj++){
136			subCRC = CM_Compute_CRC16(subCRC, (BYTE*)(subData+jj), 1);
137			partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)(subData+jj), 1);
138			}
139			//subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subData, dataLength);
140
141			//} else {
142			// subData = "0";
143			//}
144
145			//Read the CRC inside of MCMD
146			//To check sum of all sub packet data, module 256 **TBD**
147			InputFile->read(subTrailer, sizeof(subTrailer));
148			subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subTrailer, 2);
149			partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subTrailer, 2);
150
151			//Read the CRC of subPacket
152			InputFile->read(subPckCRC, sizeof(subPckCRC));
153			partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subPckCRC, 2);
154			readCRC = (((BYTE)subPckCRC[0])<<8) + ((BYTE)subPckCRC[1]);
155
156			//finally check if the RecordCRC is correct
157			//and finally update the partialCRC
158			//TO DO - if one CRC is wrong also the total one will be corrupted
159			if (subCRC == readCRC){
160			rec = new(recs[i++]) McmdRecord(); //aggiungo un nuovo McmdRecord all'evento
161			rec->MCMD_RECORD_OBT = (((UINT32)OBT[0]<<24)&0xFF000000) + (((UINT32)OBT[1]<<16)&0x00FF0000) + (((UINT32)OBT[2]<<8)&0x0000FF00) + (((UINT32)OBT[3])&0x000000FF);
162			rec->SeqID = (((BYTE)subHeader[0]<<8)&0xFF00) + (((BYTE)subHeader[1])&0x00FF);
163			rec->Tbit = ((BYTE)((subHeader[2]&0x80))>>7);
164			rec->ID1 = (BYTE)subHeader[3];
165			rec->McmdLength = (0x0fff)&(((UINT16)(subHeader[4]<<8)) + ((UINT16)subHeader[5]));
166			rec->TimeTag = (((UINT32)OBT[6]<<24)&0xFF000000) + (((UINT32)OBT[7]<<16)&0x00FF0000) + (((UINT32)OBT[8]<<8)&0x0000FF00) + (((UINT32)OBT[9])&0x000000FF);
167			rec->endID = (BYTE)subTrailer[1];
168			rec->setMcmdData((BYTE*)&subData);
169	kusanagi	1.3	delete [] subData;
170	kusanagi	1.1	} else {
171			cat << log4cpp::Priority::ERROR
172			<< "Wrong CRC on Subpacket int MCMD Packet starting at position " << start
173			<< "\n " << log4cpp::CategoryStream::ENDLINE;
174			}
175			}
176			//in the end compare the calculated partial CRC with the MCMD packet CRC
177			InputFile->read(eventCRC, sizeof(eventCRC));
178			readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);
179
180			if(!(partialCRC == readCRC)) {
181			cat << log4cpp::Priority::ERROR
182	kusanagi	1.2	<< "Wrong CRC for CalibS4 Packet"
183	kusanagi	1.1	<< "\n " << log4cpp::CategoryStream::ENDLINE;
184			}
185			}
186
187