yoda/techmodel/TsbTReader.cpp

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

extern "C" {
    #include "CRC.h"
}
#include "ReaderAlgorithms.h"

using namespace pamela::techmodel;

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

/**
 * Constructor. 
 */
TsbTReader::TsbTReader(void): 
  TechmodelAlgorithm(PacketType::TsbT, "TechmodelTsbT") { 
  TsbT = new TsbTEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string TsbTReader::GetVersionInfo(void) const {
  return     "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/TsbTReader.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 TsbTReader::Init(PamelaRun *run) {
  logger->debug(_T("Initialize"));
  SetInputStream(run);
  run->WriteSubPacket(this, &TsbT, TsbT->Class());
}

/**
 * Unpack the TsbT event from an input file.
 */
void TsbTReader::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 TsbT Packet ");
    /*if (subCRC != readCRC) {
        throw WrongCRCException("WRONG CRC FOR TsbT PACKET."); 
    }*/
    
    TsbTRecord* rec;
    TsbT->Records->Clear();
    TClonesArray &recs = *(TsbT->Records);
    
    int i = 0;
    long int offset = 0; //dataLength/31;
    while (offset < dataLength){
        rec               = new(recs[i++]) TsbTRecord(); //add a new TsbT
        rec->RECORD_OBT  = (((UINT32)subData[offset]<<24)&0xFF000000) + (((UINT32)subData[offset+1]<<16)&0x00FF0000) +  (((UINT32)subData[offset+2]<<8)&0x0000FF00) + (((UINT32)subData[offset+3])&0x000000FF);
        rec->STATUS_CODE = subData[offset+4];
        for (int j=0; j < 13; j++){
           rec->TEMPERATURES[j]  = ((((UINT16)subData[offset+5+(2*j)]<<8)&0xFF00) + (((UINT16)subData[offset+6+(2*j)])&0x00FF));
        }
        offset = offset + 31;
    }
}

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