1 |
mocchiut |
1.2 |
// Implementation of the PhysEndRunReader class. |
2 |
|
|
//new version... |
3 |
|
|
extern "C" { |
4 |
|
|
#include "CRC.h" |
5 |
|
|
} |
6 |
|
|
|
7 |
|
|
#include "ReaderAlgorithms.h" |
8 |
|
|
|
9 |
|
|
using namespace pamela::techmodel; |
10 |
|
|
|
11 |
|
|
|
12 |
|
|
/** |
13 |
|
|
* Constructor. |
14 |
|
|
*/ |
15 |
|
|
PhysEndRunReader::PhysEndRunReader(void): |
16 |
|
|
TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") { |
17 |
|
|
physEndRun = new PhysEndRunEvent(); |
18 |
|
|
} |
19 |
|
|
|
20 |
|
|
/** |
21 |
|
|
* Initialize the algorithm with a special run. This will initialize the |
22 |
|
|
* event reader routines for all packet types. |
23 |
|
|
*/ |
24 |
|
|
void PhysEndRunReader::Init(PamelaRun *run) { |
25 |
|
|
run->WriteSubPacket(this, &physEndRun, physEndRun->Class()); |
26 |
|
|
} |
27 |
|
|
|
28 |
|
|
/** |
29 |
|
|
* Unpack the PhysEndRun event from an input file. |
30 |
|
|
* void PhysEndRunReader::PKT_RunEvent(char* subData, long int length) throw (WrongCRCException_PKTUsed){ |
31 |
|
|
* |
32 |
|
|
*/ |
33 |
|
|
void PhysEndRunReader::PKT_RunEvent(char* subData, long int length) throw (Exception){ |
34 |
|
|
std::stringstream oss; |
35 |
|
|
string msg; |
36 |
|
|
UINT16 subCRC; //CRC of the data |
37 |
|
|
UINT16 readCRC; //CRC read from the end of the subpacket |
38 |
|
|
|
39 |
|
|
|
40 |
|
|
// baco /gpfs/wizard/flight/production/preRawreader/05510005.pam |
41 |
|
|
if (length<21) |
42 |
|
|
{ |
43 |
|
|
oss.str(""); |
44 |
|
|
oss<<"Wrong Lenght for PhysEndRun Packet: "<<" length = "<< length<<" too small to be real"; |
45 |
|
|
msg=oss.str(); |
46 |
|
|
PamOffLineSW::mainLogUtil->logError(msg); |
47 |
|
|
//the packet is not good and will be discarded |
48 |
|
|
throw FatalException(" Wrong Lenght for PhysEndRun Packet "); |
49 |
|
|
} |
50 |
|
|
|
51 |
|
|
long int dataLength = length - 2; |
52 |
|
|
long int tbInitPos = dataLength - 19; // 19 is the length of th TB |
53 |
|
|
|
54 |
|
|
subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength); |
55 |
|
|
readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF); |
56 |
|
|
// std::cout<<"computedCRC "<< (UINT16)subCRC<<" readCRC "<< (UINT16)readCRC<<" readCRC_1 "<<(UINT16)subData[length-2]<<" readCRC_2 "<<(UINT16)subData[length-1]<<endl; |
57 |
|
|
|
58 |
|
|
|
59 |
|
|
//--------------- CALO SECTION------------------------------------ |
60 |
|
|
long int offset = 0; |
61 |
|
|
int j = 0; |
62 |
|
|
//physEndRun->CALO_ENDRUN. |
63 |
|
|
while(offset < tbInitPos){ |
64 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F); |
65 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]); |
66 |
|
|
//subData[offset+2--->3] is a fixed word = 0x23 |
67 |
|
|
//If this bit is ? the subsequent CALO_HK are not valid |
68 |
|
|
printf(" j %i calo_board_id_hk %X \n",j,physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK); |
69 |
|
|
printf(" j %i calo_board_status_hk %X \n",j,physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK); |
70 |
|
|
if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) { |
71 |
|
|
UInt_t kk = 0; |
72 |
|
|
for(int k = 0; k < 11 ; k++){ |
73 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]); |
74 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF); |
75 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF); |
76 |
|
|
printf(" j %i k %i calo_hk0 %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_HK0[k]); |
77 |
|
|
printf(" j %i k %i calo_hk1 %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_HK1[k]); |
78 |
|
|
kk += 2; |
79 |
|
|
} |
80 |
|
|
offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes) |
81 |
|
|
} else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes) |
82 |
|
|
//subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24] |
83 |
|
|
// |
84 |
|
|
for(int k = 0; k < 7 ; k++){ |
85 |
|
|
// |
86 |
|
|
// in CALO_BOARD_ID_REG and CALO_BOARD_STATUS_REG we save only the last id and status register of the 7 and we don't save the CRC for each register. |
87 |
|
|
// |
88 |
|
|
// printf(" INLOPP k %i subdata %X \n",k,(UINT8)(subData[offset])); |
89 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k] = ((UINT8)(subData[offset])&0x1F); |
90 |
|
|
offset++; |
91 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k] = (UINT8)(subData[offset]); |
92 |
|
|
offset++; |
93 |
|
|
printf(" j %i k %i calo_board_id_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k]); |
94 |
|
|
printf(" j %i k %i calo_board_status_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k]); |
95 |
|
|
offset++; |
96 |
|
|
//subData[offset+2--->3] is a fixed word = 0x2 that is the data lenght |
97 |
|
|
offset++; |
98 |
|
|
physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset+1])&0x00FF); |
99 |
|
|
printf(" j %i k %i calo_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_REG[k]); |
100 |
|
|
offset++; |
101 |
|
|
// CRC |
102 |
|
|
offset++; |
103 |
|
|
offset++; |
104 |
|
|
// |
105 |
|
|
offset++; |
106 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF); |
107 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0; |
108 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0; |
109 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0; |
110 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0; |
111 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0; |
112 |
|
|
//subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7] |
113 |
|
|
// offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte) |
114 |
|
|
}; |
115 |
|
|
// |
116 |
|
|
j++; |
117 |
|
|
} |
118 |
|
|
|
119 |
|
|
// //--------------- CALO SECTION------------------------------------ |
120 |
|
|
// long int offset = 0; |
121 |
|
|
// int j = 0; |
122 |
|
|
// //physEndRun->CALO_ENDRUN. |
123 |
|
|
// while(offset < tbInitPos){ |
124 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F); |
125 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]); |
126 |
|
|
// //subData[offset+2--->3] is a fixed word = 0x23 |
127 |
|
|
// //If this bit is ? the subsequent CALO_HK are not valid |
128 |
|
|
// if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) { |
129 |
|
|
// UInt_t kk = 0; |
130 |
|
|
// for(int k = 0; k < 11 ; k++){ |
131 |
|
|
// // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]); |
132 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF); |
133 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF); |
134 |
|
|
// kk += 2; |
135 |
|
|
// } |
136 |
|
|
// offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes) |
137 |
|
|
// } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes) |
138 |
|
|
// //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24] |
139 |
|
|
|
140 |
|
|
// // for(int k = 0; k < 7 ; k++){ |
141 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F); |
142 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]); |
143 |
|
|
// //subData[offset+2--->3] is a fixed word = 0x2 |
144 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[0] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF); |
145 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF); |
146 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0; |
147 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0; |
148 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0; |
149 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0; |
150 |
|
|
// physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0; |
151 |
|
|
// //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7] |
152 |
|
|
// offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte) |
153 |
|
|
// // } |
154 |
|
|
// j++; |
155 |
|
|
// } |
156 |
|
|
|
157 |
|
|
//--------------- TB SECTION------------------------------------ |
158 |
|
|
printf(" subdata[%u] %X \n",offset,subData[offset]); |
159 |
|
|
physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF); |
160 |
|
|
//subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2] |
161 |
|
|
physEndRun->TB_ENDRUN.TB_PMT_MASK_S3 = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F); |
162 |
|
|
physEndRun->TB_ENDRUN.TB_PMT_MASK_S2 = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F); |
163 |
|
|
physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF); |
164 |
|
|
physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF); |
165 |
|
|
//subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9] |
166 |
|
|
physEndRun->TB_ENDRUN.TB_CALO_MASK = ((UINT8)subData[offset+10])>>5; |
167 |
|
|
physEndRun->TB_ENDRUN.TB_S4_MASK = (((UINT8)subData[offset+10])>>1)&0x0F; |
168 |
|
|
//subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11] |
169 |
|
|
physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)subData[offset+12]<<16)&0x00FF0000) + (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF); |
170 |
|
|
//subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15] |
171 |
|
|
physEndRun->TB_ENDRUN.TB_CALIB_FLAG = ((UINT8)subData[offset+16])>>7; |
172 |
|
|
physEndRun->TB_ENDRUN.TB_CALO_TRIG = ((UINT8)subData[offset+16])>>6; |
173 |
|
|
physEndRun->TB_ENDRUN.TB_S4_TRIG = ((UINT8)subData[offset+16])>>5; |
174 |
|
|
physEndRun->TB_ENDRUN.TB_TOF_TRIG = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6); |
175 |
|
|
//subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18] |
176 |
|
|
|
177 |
|
|
// I want to use the packet so the exception is thrown here after all .. |
178 |
|
|
if (subCRC != readCRC) |
179 |
|
|
{ |
180 |
|
|
oss.str(""); |
181 |
|
|
oss<<"Wrong CRC for PhysEndRun Packet: "<<" CRC COMPUTED= "<< subCRC<<" CRC READ= "<< readCRC; |
182 |
|
|
msg=oss.str(); |
183 |
|
|
PamOffLineSW::mainLogUtil->logWarning(msg); |
184 |
|
|
//the packet is not good but used |
185 |
|
|
throw WrongCRCException_PKTUsed(" Wrong CRC for PhysEndRun Packet "); |
186 |
|
|
} |
187 |
|
|
} |
188 |
|
|
|
189 |
|
|
/** |
190 |
|
|
* Get a string with the version info of the algorithm. |
191 |
|
|
*/ |
192 |
|
|
std::string PhysEndRunReader::GetVersionInfo(void) const { |
193 |
|
|
return |
194 |
|
|
"$Header: /afs/ba.infn.it/user/pamela/src/CVS/chewbacca/PamOffLineSW/techmodel/PhysEndRunReader.cpp,v 1.1.1.1 2008/09/23 07:20:24 mocchiut Exp $\n"; |
195 |
|
|
} |