yoda/techmodel/TabDumpReader.cpp

/** @file
 * $Source: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/TabDumpReader.cpp,v $
 * $Id: TabDumpReader.cpp,v 6.2 2006/05/30 19:10:03 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the TabDumpReader class.
 */

extern "C" {
    #include "CRC.h"
}
#include "ReaderAlgorithms.h"
using namespace pamela::techmodel;

static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.TabDumpReader"));

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

/**
 * Get a string with the version info of the algorithm.
 */
std::string TabDumpReader::GetVersionInfo(void) const {
  return "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/TabDumpReader.cpp,v 6.2 2006/05/30 19:10:03 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) {
  logger->debug(_T("Initialize"));
  SetInputStream(run);
  run->WriteSubPacket(this, &TabDump, TabDump->Class());
}

/**
 * Unpack the TabDump event from an input file.
 */
void TabDumpReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
  
    char      subData[length];
    UINT16    subCRC;      //calculated CRC of the data
    UINT16    readCRC;     //CRC read from the end of the subpacket
    long int  dataLength = length - 2; //the block of data
    memset(subData,  0, length*sizeof(char));  
    InputFile->read(subData, sizeof(subData));
    subCRC  = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
    readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
    
    if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
    /*if (subCRC != readCRC) {
        logger->error("WRONG CRC FOR TabDump PACKET. Processed anyway for CPU debugging");
    }*/

    TabDumpRecord* rec;
    TabDump->PARAMETER_STAMP = (((UINT32)subData[0]<<24)&0xFF000000) + (((UINT32)subData[1]<<16)&0x00FF0000) +  (((UINT32)subData[2]<<8)&0x0000FF00) + (((UINT32)subData[3])&0x000000FF);
    long int offset = 4;
    int i = 0;
    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 TArrayI((int)((rec->Nrow)*(rec->Ncol)), (int*)(subData+offset+3));
        offset = offset + sizeof(UINT32)*(rec->Nrow)*(rec->Ncol) + 3;
    }
}
1	/** @file
2	* $Source: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/TabDumpReader.cpp,v $
3	* $Id: TabDumpReader.cpp,v 6.2 2006/05/30 19:10:03 kusanagi Exp $
4	* $Author: kusanagi $
5	*
6	* Implementation of the TabDumpReader class.
7	*/
8
9	extern "C" {
10	#include "CRC.h"
11	}
12	#include "ReaderAlgorithms.h"
13	using namespace pamela::techmodel;
14
15	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.TabDumpReader"));
16
17	/**
18	* Constructor.
19	*/
20	TabDumpReader::TabDumpReader(void):
21	TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {
22	TabDump = new TabDumpEvent();
23	logger->debug(_T("Constructor"));
24	}
25
26	/**
27	* Get a string with the version info of the algorithm.
28	*/
29	std::string TabDumpReader::GetVersionInfo(void) const {
30	return "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/TabDumpReader.cpp,v 6.2 2006/05/30 19:10:03 kusanagi Exp $\n";
31	}
32
33	/**
34	* Initialize the algorithm with a special run. This will initialize the
35	* event reader routines for all packet types.
36	*/
37	void TabDumpReader::Init(PamelaRun *run) {
38	logger->debug(_T("Initialize"));
39	SetInputStream(run);
40	run->WriteSubPacket(this, &TabDump, TabDump->Class());
41	}
42
43	/**
44	* Unpack the TabDump event from an input file.
45	*/
46	void TabDumpReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
47
48	char subData[length];
49	UINT16 subCRC; //calculated CRC of the data
50	UINT16 readCRC; //CRC read from the end of the subpacket
51	long int dataLength = length - 2; //the block of data
52	memset(subData, 0, length*sizeof(char));
53	InputFile->read(subData, sizeof(subData));
54	subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
55	readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
56
57	if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
58	/*if (subCRC != readCRC) {
59	logger->error("WRONG CRC FOR TabDump PACKET. Processed anyway for CPU debugging");
60	}*/
61
62	TabDumpRecord* rec;
63	TabDump->PARAMETER_STAMP = (((UINT32)subData[0]<<24)&0xFF000000) + (((UINT32)subData[1]<<16)&0x00FF0000) + (((UINT32)subData[2]<<8)&0x0000FF00) + (((UINT32)subData[3])&0x000000FF);
64	long int offset = 4;
65	int i = 0;
66	TabDump->Records->Clear();
67	TClonesArray &recs = *(TabDump->Records);
68	while (offset < dataLength){
69	rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
70	rec->Tab_ID = ((UINT8)subData[offset])&0xFF;
71	rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
72	rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
73	rec->Data = new TArrayI((int)((rec->Nrow)(rec->Ncol)), (int)(subData+offset+3));
74	offset = offset + sizeof(UINT32)(rec->Nrow)(rec->Ncol) + 3;
75	}
76	}