yoda/techmodel/TabDumpReader.cpp

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

#include <string>
#include <log4cpp/Category.hh>
#include <fstream>
#include "stdio.h"
extern "C" {
    #include "CRC.h"
}
#include "ReaderAlgorithms.h"

using namespace pamela;
using namespace pamela::techmodel;

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

/**
 * Constructor. 
 */
TabDumpReader::TabDumpReader(void): 
  TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") { 
  TabDump = new TabDumpEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string TabDumpReader::GetVersionInfo(void) const {
  return "$Header: /home/cvsmanager/yoda/techmodel/TabDumpReader.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 TabDumpReader::Init(PamelaRun *run) {
  SetInputStream(run);
  run->WriteSubPacket(this, &TabDump, TabDump->Class());
}

/**
 * Unpack the TabDump event from an input file.
 */
void TabDumpReader::RunEvent(int EventNumber, long int length) {
  
  char     *subData;
  char      eventCRC[2];
  UINT16    subCRC;      //CRC of the data
  UINT16    readCRC;     //CRC read from the end of the subpacket
  long int  dataLength;
  
    //the 2 bytes subtracted belong to the final event CRC bytes
    dataLength = length - (long int)2;

    subData = new char[dataLength];
    InputFile->read(subData, sizeof(unsigned char)*dataLength);
    subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);

    //took the final CRC to compare it with the previous calculated CRC of the data
    InputFile->read(eventCRC, sizeof(eventCRC));
    readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);

    //finally check if the RecordCRC is correct and finally update the partialCRC
    if (subCRC == readCRC){
        TabDumpRecord* rec;
        long int offset = 0;
        int i = 0;
        int size;
        TabDump->Records->Clear();
        TClonesArray &recs = *(TabDump->Records);
        while (offset < dataLength){
            rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
            rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;
            rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
            rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
            rec->Data = new UINT32[(rec->Nrow)*(rec->Ncol)];
            size = sizeof(UINT32)*(rec->Nrow)*(rec->Ncol);
            memcpy(rec->Data, (UINT32*)(subData+offset+3), size);
            offset = offset + size + 3;
        }
        free(subData);
    } else {
        cat <<  log4cpp::Priority::ERROR
            <<  "Wrong CRC on Subpacket in TabDump Packet "
            <<  "\n " << log4cpp::CategoryStream::ENDLINE;
    }
  }
1	kusanagi	1.1	/** @file
2	kusanagi	1.2	* $Source: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v $
3	kusanagi	1.3	* $Id: TabDumpReader.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 TabDumpReader class.
7			*/
8
9			#include <string>
10			#include <log4cpp/Category.hh>
11			#include <fstream>
12			#include "stdio.h"
13			extern "C" {
14			#include "CRC.h"
15			}
16			#include "ReaderAlgorithms.h"
17
18			using namespace pamela;
19			using namespace pamela::techmodel;
20
21			static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.TabDumpReader");
22
23			/**
24			* Constructor.
25			*/
26			TabDumpReader::TabDumpReader(void):
27			TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {
28			TabDump = new TabDumpEvent();
29			}
30
31			/**
32			* Get a string with the version info of the algorithm.
33			*/
34			std::string TabDumpReader::GetVersionInfo(void) const {
35	kusanagi	1.3	return "$Header: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $\n";
36	kusanagi	1.1	}
37
38			/**
39			* Initialize the algorithm with a special run. This will initialize the
40			* event reader routines for all packet types.
41			*/
42			void TabDumpReader::Init(PamelaRun *run) {
43			SetInputStream(run);
44			run->WriteSubPacket(this, &TabDump, TabDump->Class());
45			}
46
47			/**
48			* Unpack the TabDump event from an input file.
49			*/
50			void TabDumpReader::RunEvent(int EventNumber, long int length) {
51
52			char *subData;
53			char eventCRC[2];
54			UINT16 subCRC; //CRC of the data
55			UINT16 readCRC; //CRC read from the end of the subpacket
56			long int dataLength;
57
58			//the 2 bytes subtracted belong to the final event CRC bytes
59			dataLength = length - (long int)2;
60
61			subData = new char[dataLength];
62			InputFile->read(subData, sizeof(unsigned char)*dataLength);
63			subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);
64
65			//took the final CRC to compare it with the previous calculated CRC of the data
66			InputFile->read(eventCRC, sizeof(eventCRC));
67			readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);
68	kusanagi	1.2
69	kusanagi	1.1	//finally check if the RecordCRC is correct and finally update the partialCRC
70			if (subCRC == readCRC){
71	kusanagi	1.2	TabDumpRecord* rec;
72			long int offset = 0;
73			int i = 0;
74	kusanagi	1.3	int size;
75	kusanagi	1.2	TabDump->Records->Clear();
76			TClonesArray &recs = *(TabDump->Records);
77			while (offset < dataLength){
78			rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
79			rec->Tab_ID = ((UINT8)subData[offset])&0xFF;
80	kusanagi	1.3	rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
81			rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
82			rec->Data = new UINT32[(rec->Nrow)*(rec->Ncol)];
83			size = sizeof(UINT32)(rec->Nrow)(rec->Ncol);
84			memcpy(rec->Data, (UINT32*)(subData+offset+3), size);
85			offset = offset + size + 3;
86	kusanagi	1.2	}
87	kusanagi	1.1	free(subData);
88			} else {
89			cat << log4cpp::Priority::ERROR
90			<< "Wrong CRC on Subpacket in TabDump Packet "
91			<< "\n " << log4cpp::CategoryStream::ENDLINE;
92			}
93			}