yoda/techmodel/PhysEndRunReader.cpp


// Implementation of the PhysEndRunReader class.
 

#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"

using namespace pamela;
using namespace pamela::techmodel;

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

/**
 * Constructor. 
 */
PhysEndRunReader::PhysEndRunReader(void): 
  TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") { 
  logger->debug(_T("Constructor"));
  physEndRun = new PhysEndRunEvent();
}

/**
 * Get a string with the version info of the algorithm.
 */
std::string PhysEndRunReader::GetVersionInfo(void) const {
  return 
    "$Header: /home/cvsmanager/yoda/techmodel/PhysEndRunReader.cpp,v 4.5 2005/08/29 09:29:03 kusanagi Exp $\n";
}

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

/**
 * Unpack the PhysEndRun event from an input file.
 */

void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
    char subData[length];
    UINT16    subCRC;      //CRC of the data
    UINT16    readCRC;     //CRC read from the end of the subpacket
    long int  dataLength;
    long int  tbInitPos;

    //int numRecords = ((length - 2)/6);
    dataLength = length - (long int)2;    
    tbInitPos = dataLength - 19; // 19 is the length of th TB
    //physEndRun->CALO_ENDRUN[];
    //physEndRun->TB_ENDRUN;
    
    InputFile->read(subData, sizeof(unsigned char)*length);
    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 PhysEndRun Packet ");
    /*if (subCRC != readCRC) {
        logger->error(" Wrong CRC for PhysEndRun Packet but I process it anyway");
    }*/
    //--------------- CALO SECTION------------------------------------
    long int offset = 0;
    int j = 0;
    //physEndRun->CALO_ENDRUN.
    while(offset < tbInitPos){
        physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK     = (((UINT8)subData[offset])&0x1F);
        physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
        //subData[offset+2--->3] is a fixed word = 0x23
        //If this bit is ? the subsequent CALO_HK are not valid
        if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
            for(int k = 0; k < 11 ; k++){
                physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+3+k]<<8))&0xFF00) + (((UINT16)subData[offset+4+k])&0x00FF);
                physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+14+k]<<8))&0xFF00) + (((UINT16)subData[offset+15+k])&0x00FF);
            }
            offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
        } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
        //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]

        for(int k = 0; k < 7 ; k++){
            physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG     = ((UINT8)(subData[offset])&0x1F);
            physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
            //subData[offset+2--->3] is a fixed word = 0x2
            physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
            //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
            offset = offset + 8; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte) 
        }
        j++;
    }
    
    //--------------- TB SECTION------------------------------------
    physEndRun->TB_ENDRUN.TB_ALARM_MASK   = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
    //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
    physEndRun->TB_ENDRUN.TB_PMT_MASK_S3  = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
    physEndRun->TB_ENDRUN.TB_PMT_MASK_S2  = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
    physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
    physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
    //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
    physEndRun->TB_ENDRUN.TB_CALO_MASK    = ((UINT8)subData[offset+10])>>5;
    physEndRun->TB_ENDRUN.TB_S4_MASK      = (((UINT8)subData[offset+10])>>1)&0x0F;
    //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
    physEndRun->TB_ENDRUN.TB_BUSY_MASK    = (((UINT32)subData[offset+12]<<16)&0x00FF0000) +  (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
    //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
    physEndRun->TB_ENDRUN.TB_CALIB_FLAG   = ((UINT8)subData[offset+16])>>7;
    physEndRun->TB_ENDRUN.TB_CALO_TRIG    = ((UINT8)subData[offset+16])>>6; 
    physEndRun->TB_ENDRUN.TB_S4_TRIG      = ((UINT8)subData[offset+16])>>5; 
    physEndRun->TB_ENDRUN.TB_TOF_TRIG     = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
    //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
}

1	kusanagi	1.1
2			// Implementation of the PhysEndRunReader class.
3
4
5			#define UINT unsigned int
6			#define BYTE unsigned char
7			#include <string>
8	kusanagi	1.2	#include <log4cxx/logger.h>
9	kusanagi	1.1	extern "C" {
10			#include <sys/time.h>
11			#include "CRC.h"
12			}
13
14			#include <fstream>
15			#include "stdio.h"
16			#include "ReaderAlgorithms.h"
17
18			using namespace pamela;
19			using namespace pamela::techmodel;
20
21	kusanagi	1.2	static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.PhysEndRunReader"));
22	kusanagi	1.1
23			/**
24			* Constructor.
25			*/
26			PhysEndRunReader::PhysEndRunReader(void):
27			TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") {
28	kusanagi	1.2	logger->debug(_T("Constructor"));
29	kusanagi	1.1	physEndRun = new PhysEndRunEvent();
30			}
31
32			/**
33			* Get a string with the version info of the algorithm.
34			*/
35			std::string PhysEndRunReader::GetVersionInfo(void) const {
36			return
37	kusanagi	5.0	"$Header: /home/cvsmanager/yoda/techmodel/PhysEndRunReader.cpp,v 4.5 2005/08/29 09:29:03 kusanagi Exp $\n";
38	kusanagi	1.1	}
39
40			/**
41			* Initialize the algorithm with a special run. This will initialize the
42			* event reader routines for all packet types.
43			*/
44			void PhysEndRunReader::Init(PamelaRun *run) {
45			SetInputStream(run);
46			run->WriteSubPacket(this, &physEndRun, physEndRun->Class());
47	kusanagi	1.2	logger->debug(_T("Initialize"));
48	kusanagi	1.1	}
49
50			/**
51			* Unpack the PhysEndRun event from an input file.
52			*/
53	kusanagi	2.3
54	kusanagi	2.1	void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
55	kusanagi	2.3	char subData[length];
56			UINT16 subCRC; //CRC of the data
57			UINT16 readCRC; //CRC read from the end of the subpacket
58			long int dataLength;
59			long int tbInitPos;
60
61			//int numRecords = ((length - 2)/6);
62			dataLength = length - (long int)2;
63			tbInitPos = dataLength - 19; // 19 is the length of th TB
64			//physEndRun->CALO_ENDRUN[];
65			//physEndRun->TB_ENDRUN;
66	kusanagi	1.1
67	kusanagi	2.3	InputFile->read(subData, sizeof(unsigned char)*length);
68			subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
69			readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
70
71	kusanagi	4.5	if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for PhysEndRun Packet ");
72	kusanagi	3.1	/*if (subCRC != readCRC) {
73	kusanagi	2.3	logger->error(" Wrong CRC for PhysEndRun Packet but I process it anyway");
74	kusanagi	3.1	}*/
75	kusanagi	2.3	//--------------- CALO SECTION------------------------------------
76			long int offset = 0;
77			int j = 0;
78			//physEndRun->CALO_ENDRUN.
79			while(offset < tbInitPos){
80			physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F);
81			physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
82			//subData[offset+2--->3] is a fixed word = 0x23
83			//If this bit is ? the subsequent CALO_HK are not valid
84			if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
85			for(int k = 0; k < 11 ; k++){
86			physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+3+k]<<8))&0xFF00) + (((UINT16)subData[offset+4+k])&0x00FF);
87			physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+14+k]<<8))&0xFF00) + (((UINT16)subData[offset+15+k])&0x00FF);
88			}
89			offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
90			} else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
91			//subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
92
93			for(int k = 0; k < 7 ; k++){
94			physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F);
95			physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
96			//subData[offset+2--->3] is a fixed word = 0x2
97			physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
98			//subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
99			offset = offset + 8; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
100			}
101			j++;
102			}
103	kusanagi	1.1
104	kusanagi	2.3	//--------------- TB SECTION------------------------------------
105			physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
106			//subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
107	kusanagi	4.2	physEndRun->TB_ENDRUN.TB_PMT_MASK_S3 = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
108			physEndRun->TB_ENDRUN.TB_PMT_MASK_S2 = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
109			physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
110			physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
111	kusanagi	2.3	//subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
112	kusanagi	4.4	physEndRun->TB_ENDRUN.TB_CALO_MASK = ((UINT8)subData[offset+10])>>5;
113			physEndRun->TB_ENDRUN.TB_S4_MASK = (((UINT8)subData[offset+10])>>1)&0x0F;
114	kusanagi	2.3	//subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
115			physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)subData[offset+12]<<16)&0x00FF0000) + (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
116			//subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
117			physEndRun->TB_ENDRUN.TB_CALIB_FLAG = ((UINT8)subData[offset+16])>>7;
118	kusanagi	4.1	physEndRun->TB_ENDRUN.TB_CALO_TRIG = ((UINT8)subData[offset+16])>>6;
119			physEndRun->TB_ENDRUN.TB_S4_TRIG = ((UINT8)subData[offset+16])>>5;
120			physEndRun->TB_ENDRUN.TB_TOF_TRIG = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
121	kusanagi	2.3	//subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
122	kusanagi	1.1	}
123