/[PAMELA software]/yoda/techmodel/PhysEndRunReader.cpp
ViewVC logotype

Contents of /yoda/techmodel/PhysEndRunReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 6.1 - (show annotations) (download)
Tue May 30 19:10:02 2006 UTC (18 years, 6 months ago) by kusanagi
Branch: MAIN
Changes since 6.0: +1 -11 lines
Major update.
All the packet officially produced by PAMELA are implemented and unpacked.
The RegistryEvent Packet has been removed and put into another library.
New version, releasd by D.Campana, of tofunpack.

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

  ViewVC Help
Powered by ViewVC 1.1.23