yoda/techmodel/CalibS4Reader.cpp

/** @file
 * $Source: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v $
 * $Id: CalibS4Reader.cpp,v 2.5 2005/02/19 17:20:29 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.5 2005/02/19 17:20:29 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[length];    //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[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&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	kusanagi	1.1	/** @file
2	kusanagi	1.2	* $Source: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v $
3	kusanagi	3.0	* $Id: CalibS4Reader.cpp,v 2.5 2005/02/19 17:20:29 kusanagi Exp $
4	kusanagi	1.2	* $Author: kusanagi $
5	kusanagi	1.1	*
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	kusanagi	1.3	#include <log4cxx/logger.h>
15	kusanagi	1.1	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	kusanagi	2.3	#include "event/physics/S4/S4Event.h"
24			//#include "event/CalibS4Event.h"
25	kusanagi	1.1
26
27			using namespace pamela;
28			using namespace pamela::techmodel;
29
30	kusanagi	1.3	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibS4Reader"));
31	kusanagi	1.1
32			/**
33			* Constructor.
34			*/
35			CalibS4Reader::CalibS4Reader(void):
36			TechmodelAlgorithm(PacketType::Log, "TechmodelCalibS4Reader") {
37	kusanagi	1.3	logger->debug(_T("Constructor"));
38	kusanagi	2.3	calibS4 = new CalibS4Event();
39	kusanagi	1.1	}
40
41			/**
42			* Get a string with the version info of the algorithm.
43			*/
44			std::string CalibS4Reader::GetVersionInfo(void) const {
45			return
46	kusanagi	3.0	"$Header: /home/cvsmanager/yoda/techmodel/CalibS4Reader.cpp,v 2.5 2005/02/19 17:20:29 kusanagi Exp $\n";
47	kusanagi	1.1	}
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	kusanagi	1.3	logger->debug(_T("Initialize"));
55	kusanagi	1.1	SetInputStream(run);
56	kusanagi	2.3	run->WriteSubPacket(this, &calibS4, calibS4->Class());
57	kusanagi	1.1	}
58
59			/**
60			* Unpack the CalibS4 event from an input file.
61			*/
62	kusanagi	2.4	void CalibS4Reader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
63	kusanagi	2.3
64			int offset;
65	kusanagi	2.5	char subData[length]; //impulse(11 + 10 + 01 + 00)eventLength + CRC = 12846 + UINT16
66	kusanagi	2.4	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	kusanagi	2.3
71			S4::S4Event* rec;
72			calibS4->Records->Clear();
73			TClonesArray &recs = *(calibS4->Records);
74
75	kusanagi	2.4	InputFile->read(subData, sizeof(unsigned char)*length);
76			subCRC = CM_Compute_CRC16(0, (BYTE*)subData, length - 2);
77	kusanagi	2.5	readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
78	kusanagi	2.4
79			if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for CalibS4 Packet ");
80	kusanagi	2.3
81			for(int i = 0; i < numRecords; i++) {
82			offset = i*6;
83			rec = new(recs[i]) S4::S4Event(); //add a new S4Event
84			rec->S4_REG_STATUS = (((UINT8)subData[offset])&0xF0);
85			rec->S4_DATA = ((((UINT16)subData[offset]<<8)&0x0FFF) + (((UINT16)subData[offset+1])&0x00FF));
86			rec->S4_CMD_NUM = (((UINT8)subData[offset+2])&0xFF);
87			rec->S4_RESP_LENGHT = ((((UINT16)subData[offset+3]<<8)&0xFF00) + (((UINT16)subData[offset+4])&0x00FF));
88			rec->S4_OVERALL_CHKCODE = (((UINT8)subData[offset+5])&0xFF);
89	kusanagi	2.4	}
90	kusanagi	1.1	}
91