1 |
// 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])&0xFF); |
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])&0xFF); |
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 |
// for (Int_t pi=-8; pi<40; pi++){ |
159 |
// printf(" subdata[%u] %X \n",offset+pi,(UINT8)subData[offset+pi]); |
160 |
// }; |
161 |
|
162 |
// UINT32 TB_ALARM_MASK; |
163 |
// UINT32 TB_PMT_MASK_S3S2S12; |
164 |
// UINT32 TB_PMT_MASK_S11CRC; |
165 |
// UINT16 TB_S4_CAL_MASK; |
166 |
// UINT32 TB_BUSY_MASK; |
167 |
// UINT32 TB_TRIG_CONF; |
168 |
|
169 |
physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF); |
170 |
offset += 3; |
171 |
// |
172 |
physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12 = (((UINT32)(subData[offset]<<24))&0xFF000000) + (((UINT32)(subData[offset + 1]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 2]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 3]))&0x000000FF); |
173 |
offset += 4; |
174 |
// |
175 |
physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF); |
176 |
offset += 3; |
177 |
// |
178 |
physEndRun->TB_ENDRUN.TB_S4_CAL_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)(subData[offset + 1]))&0x00FF); |
179 |
offset += 2; |
180 |
// |
181 |
physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)(subData[offset]<<24))&0xFF000000) + (((UINT32)(subData[offset + 1]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 2]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 3]))&0x000000FF); |
182 |
offset += 4; |
183 |
// |
184 |
physEndRun->TB_ENDRUN.TB_TRIG_CONF = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF); |
185 |
offset += 3; |
186 |
// |
187 |
|
188 |
// printf(" physEndRun->TB_ENDRUN.TB_ALARM_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_ALARM_MASK); |
189 |
// printf(" physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12 %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12); |
190 |
// printf(" physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC); |
191 |
// printf(" physEndRun->TB_ENDRUN.TB_S4_CAL_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_S4_CAL_MASK); |
192 |
// printf(" physEndRun->TB_ENDRUN.TB_BUSY_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_BUSY_MASK); |
193 |
// printf(" physEndRun->TB_ENDRUN.TB_TRIG_CONF %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_TRIG_CONF); |
194 |
|
195 |
// I want to use the packet so the exception is thrown here after all .. |
196 |
if (subCRC != readCRC) |
197 |
{ |
198 |
oss.str(""); |
199 |
oss<<"Wrong CRC for PhysEndRun Packet: "<<" CRC COMPUTED= "<< subCRC<<" CRC READ= "<< readCRC; |
200 |
msg=oss.str(); |
201 |
PamOffLineSW::mainLogUtil->logWarning(msg); |
202 |
//the packet is not good but used |
203 |
throw WrongCRCException_PKTUsed(" Wrong CRC for PhysEndRun Packet "); |
204 |
} |
205 |
} |
206 |
|
207 |
/** |
208 |
* Get a string with the version info of the algorithm. |
209 |
*/ |
210 |
std::string PhysEndRunReader::GetVersionInfo(void) const { |
211 |
return |
212 |
"$Header: /afs/ba.infn.it/user/pamela/src/CVS/chewbacca/PamOffLineSW/techmodel/PhysEndRunReader.cpp,v 1.2 2008/09/24 08:47:53 mocchiut Exp $\n"; |
213 |
} |