yoda/techmodel/CalibS4Reader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v $
 * $Id: CalibS4Reader.cpp,v 2.3 2004/12/16 17:32:57 kusanagi Exp $
 * $Author: kusanagi $
 * 
 * Implementation of the LogReader class.
* ToBeDone:
* Control the CRC for the entire data Packet not just for single records
 */

#define UINT unsigned int
#define BYTE  unsigned char
#include <string>
#include <log4cxx/logger.h>
extern "C" {
#include <sys/time.h>
#include "CRC.h"
}

#include <fstream>
#include "stdio.h"
#include "ReaderAlgorithms.h"
#include "event/physics/S4/S4Event.h"
//#include "event/CalibS4Event.h"


using namespace pamela;
using namespace pamela::techmodel;

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

/**
 * Constructor. 
 */
CalibS4Reader::CalibS4Reader(void): 
  TechmodelAlgorithm(PacketType::Log, "TechmodelCalibS4Reader") { 
  logger->debug(_T("Constructor"));
  calibS4 = new CalibS4Event();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string CalibS4Reader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v 2.3 2004/12/16 17:32:57 kusanagi Exp $\n";
}

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

/**
 * Unpack the CalibS4 event from an input file.
 */
void CalibS4Reader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
    
    int offset;
    char subData[3074];    //impulse*(11 + 10 + 01 + 00)*eventLength + CRC = 128*4*6 + UINT16
    UINT16    subCRC;      //CRC of the data
    UINT16    readCRC;     //CRC read from the end of the subpacket

    int numRecords = ((length - 2)/6);
        
    S4::S4Event* rec;
    calibS4->Records->Clear();
    TClonesArray &recs = *(calibS4->Records);

    InputFile->read(subData, sizeof(unsigned char)*length);

    subCRC = CM_Compute_CRC16(0, (BYTE*)subData, length - 2);
    readCRC = (((UINT16)(subData[3072]<<8))&0xFF00) + (((UINT16)subData[3073])&0x00FF);
    
    if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for CalibS4 Packet ");

    for(int i = 0; i < numRecords; i++) {
        offset = i*6; 
        rec = new(recs[i]) S4::S4Event(); //add a new S4Event
        rec->S4_REG_STATUS        = (((UINT8)subData[offset])&0xF0);
        rec->S4_DATA              = ((((UINT16)subData[offset]<<8)&0x0FFF) + (((UINT16)subData[offset+1])&0x00FF));
        rec->S4_CMD_NUM           = (((UINT8)subData[offset+2])&0xFF);
        rec->S4_RESP_LENGHT       = ((((UINT16)subData[offset+3]<<8)&0xFF00) + (((UINT16)subData[offset+4])&0x00FF));
        rec->S4_OVERALL_CHKCODE   = (((UINT8)subData[offset+5])&0xFF);
    }
}

1	/** @file
2	* $Source: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v $
3	* $Id: CalibS4Reader.cpp,v 2.3 2004/12/16 17:32:57 kusanagi Exp $
4	* $Author: kusanagi $
5	*
6	* Implementation of the LogReader class.
7	* ToBeDone:
8	* Control the CRC for the entire data Packet not just for single records
9	*/
10
11	#define UINT unsigned int
12	#define BYTE unsigned char
13	#include <string>
14	#include <log4cxx/logger.h>
15	extern "C" {
16	#include <sys/time.h>
17	#include "CRC.h"
18	}
19
20	#include <fstream>
21	#include "stdio.h"
22	#include "ReaderAlgorithms.h"
23	#include "event/physics/S4/S4Event.h"
24	//#include "event/CalibS4Event.h"
25
26
27	using namespace pamela;
28	using namespace pamela::techmodel;
29
30	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibS4Reader"));
31
32	/**
33	* Constructor.
34	*/
35	CalibS4Reader::CalibS4Reader(void):
36	TechmodelAlgorithm(PacketType::Log, "TechmodelCalibS4Reader") {
37	logger->debug(_T("Constructor"));
38	calibS4 = new CalibS4Event();
39	}
40
41	/**
42	* Get a string with the version info of the algorithm.
43	*/
44	std::string CalibS4Reader::GetVersionInfo(void) const {
45	return
46	"$Header: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v 2.3 2004/12/16 17:32:57 kusanagi Exp $\n";
47	}
48
49	/**
50	* Initialize the algorithm with a special run. This will initialize the
51	* event reader routines for all packet types.
52	*/
53	void CalibS4Reader::Init(PamelaRun *run) {
54	logger->debug(_T("Initialize"));
55	SetInputStream(run);
56	run->WriteSubPacket(this, &calibS4, calibS4->Class());
57	}
58
59	/**
60	* Unpack the CalibS4 event from an input file.
61	*/
62	void CalibS4Reader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
63
64	int offset;
65	char subData[3074]; //impulse(11 + 10 + 01 + 00)eventLength + CRC = 12846 + UINT16
66	UINT16 subCRC; //CRC of the data
67	UINT16 readCRC; //CRC read from the end of the subpacket
68
69	int numRecords = ((length - 2)/6);
70
71	S4::S4Event* rec;
72	calibS4->Records->Clear();
73	TClonesArray &recs = *(calibS4->Records);
74
75	InputFile->read(subData, sizeof(unsigned char)*length);
76
77	subCRC = CM_Compute_CRC16(0, (BYTE*)subData, length - 2);
78	readCRC = (((UINT16)(subData[3072]<<8))&0xFF00) + (((UINT16)subData[3073])&0x00FF);
79
80	if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for CalibS4 Packet ");
81
82	for(int i = 0; i < numRecords; i++) {
83	offset = i*6;
84	rec = new(recs[i]) S4::S4Event(); //add a new S4Event
85	rec->S4_REG_STATUS = (((UINT8)subData[offset])&0xF0);
86	rec->S4_DATA = ((((UINT16)subData[offset]<<8)&0x0FFF) + (((UINT16)subData[offset+1])&0x00FF));
87	rec->S4_CMD_NUM = (((UINT8)subData[offset+2])&0xFF);
88	rec->S4_RESP_LENGHT = ((((UINT16)subData[offset+3]<<8)&0xFF00) + (((UINT16)subData[offset+4])&0x00FF));
89	rec->S4_OVERALL_CHKCODE = (((UINT8)subData[offset+5])&0xFF);
90	}
91	}
92