yoda/techmodel/TabDumpReader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v $
 * $Id: TabDumpReader.cpp,v 2.1.2.2 2005/03/04 10:42:14 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the TabDumpReader class.
 */

#include <string>
#include <log4cxx/logger.h>
#include <fstream>
#include "stdio.h"
extern "C" {
    #include "CRC.h"
}
#include "ReaderAlgorithms.h"

using namespace pamela;
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: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v 2.1.2.2 2005/03/04 10:42:14 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());
  logger->debug(_T("Initialize"));
}

/**
 * Unpack the TabDump event from an input file.
 */
void TabDumpReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
    char      subData[length];
    char      eventCRC[2];
    UINT16    subCRC;      //CRC of the data
    UINT16    readCRC;     //CRC read from the end of the subpacket
    long int  dataLength = length - 2;
  
    InputFile->read(subData, sizeof(unsigned char)*length);
    subCRC  = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);
    readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
    
    //if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
    /*Pamela CPU group request to disable for a while CRC control to check values because
      the data seems to be well formatted but perhaps there is just a mistake on CRC evaluation
      into the PAMELA CPU program
    */
    if (subCRC != readCRC) {
        logger->error(" Wrong CRC for TabDump Packet but I process it anyway because of CPU DEBUG");
    }

    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;
    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);
        //rec->Data = new TArrayI((int)((rec->Nrow)*(rec->Ncol)), (int*)(subData+offset+3));
        rec->Data = new TArrayI(size, (int*)(subData+offset+3));
        //memcpy(rec->Data, (UINT32*)(subData+offset+3), size);
        offset = offset + size + 3;
    }
   //delete [] subData;
  }
1	/** @file
2	* $Source: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v $
3	* $Id: TabDumpReader.cpp,v 2.1.2.2 2005/03/04 10:42:14 kusanagi Exp $
4	* $Author: kusanagi $
5	*
6	* Implementation of the TabDumpReader class.
7	*/
8
9	#include <string>
10	#include <log4cxx/logger.h>
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 log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.TabDumpReader"));
22
23	/**
24	* Constructor.
25	*/
26	TabDumpReader::TabDumpReader(void):
27	TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {
28	TabDump = new TabDumpEvent();
29	logger->debug(_T("Constructor"));
30	}
31
32	/**
33	* Get a string with the version info of the algorithm.
34	*/
35	std::string TabDumpReader::GetVersionInfo(void) const {
36	return "$Header: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v 2.1.2.2 2005/03/04 10:42:14 kusanagi Exp $\n";
37	}
38
39	/**
40	* Initialize the algorithm with a special run. This will initialize the
41	* event reader routines for all packet types.
42	*/
43	void TabDumpReader::Init(PamelaRun *run) {
44	SetInputStream(run);
45	run->WriteSubPacket(this, &TabDump, TabDump->Class());
46	logger->debug(_T("Initialize"));
47	}
48
49	/**
50	* Unpack the TabDump event from an input file.
51	*/
52	void TabDumpReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
53	char subData[length];
54	char eventCRC[2];
55	UINT16 subCRC; //CRC of the data
56	UINT16 readCRC; //CRC read from the end of the subpacket
57	long int dataLength = length - 2;
58
59	InputFile->read(subData, sizeof(unsigned char)*length);
60	subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);
61	readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
62
63	//if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
64	/*Pamela CPU group request to disable for a while CRC control to check values because
65	the data seems to be well formatted but perhaps there is just a mistake on CRC evaluation
66	into the PAMELA CPU program
67	*/
68	if (subCRC != readCRC) {
69	logger->error(" Wrong CRC for TabDump Packet but I process it anyway because of CPU DEBUG");
70	}
71
72	TabDumpRecord* rec;
73	TabDump->PARAMETER_STAMP = (((UINT32)subData[0]<<24)&0xFF000000) + (((UINT32)subData[1]<<16)&0x00FF0000) + (((UINT32)subData[2]<<8)&0x0000FF00) + (((UINT32)subData[3])&0x000000FF);
74	long int offset = 4;
75	int i = 0;
76	int size;
77	TabDump->Records->Clear();
78	TClonesArray &recs = *(TabDump->Records);
79	while (offset < dataLength){
80	rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
81	rec->Tab_ID = ((UINT8)subData[offset])&0xFF;
82	rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
83	rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
84	//rec->Data = new UINT32[(rec->Nrow)*(rec->Ncol)];
85	size = sizeof(UINT32)(rec->Nrow)(rec->Ncol);
86	//rec->Data = new TArrayI((int)((rec->Nrow)(rec->Ncol)), (int)(subData+offset+3));
87	rec->Data = new TArrayI(size, (int*)(subData+offset+3));
88	//memcpy(rec->Data, (UINT32*)(subData+offset+3), size);
89	offset = offset + size + 3;
90	}
91	//delete [] subData;
92	}