1 |
kusanagi |
1.1 |
/** @file |
2 |
kusanagi |
1.2 |
* $Author: kusanagi $ |
3 |
kusanagi |
1.4 |
* $Date: 2004/07/08 12:31:42 $ |
4 |
|
|
* $Revision: 1.3 $ |
5 |
kusanagi |
1.1 |
* |
6 |
|
|
* Implementation of the functions of a sample Algorithm class. |
7 |
|
|
* This file can be used as a templace to develop your own algorithm. |
8 |
|
|
*/ |
9 |
|
|
|
10 |
|
|
#include <log4cpp/Category.hh> |
11 |
|
|
#include <fstream> |
12 |
|
|
#include "TechmodelAlgorithm.h" |
13 |
|
|
#include "EventReader.h" |
14 |
|
|
#include "ReaderAlgorithms.h" |
15 |
|
|
#include "Exception.h" |
16 |
|
|
#include "stdio.h" |
17 |
|
|
extern "C" { |
18 |
|
|
#include "CRC.h" |
19 |
|
|
} |
20 |
|
|
|
21 |
|
|
#define BYTE unsigned char |
22 |
|
|
|
23 |
|
|
using namespace pamela; |
24 |
|
|
using namespace pamela::techmodel; |
25 |
|
|
|
26 |
|
|
static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.EventReader"); |
27 |
|
|
static unsigned int prevPckCounter = 0; |
28 |
|
|
static unsigned int prevPckOBT = 0; |
29 |
|
|
/** |
30 |
|
|
* Constructor. |
31 |
|
|
*/ |
32 |
|
|
EventReader::EventReader(void): |
33 |
|
|
TechmodelAlgorithm(0, "TechmodelEventReader") { |
34 |
|
|
cat << log4cpp::Priority::DEBUG |
35 |
|
|
<< "Constructor " |
36 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
37 |
|
|
Header = new EventHeader(); |
38 |
|
|
/* |
39 |
|
|
|
40 |
|
|
|
41 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Alarm, |
42 |
|
|
new AlarmReader())); |
43 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Khb, |
44 |
|
|
new KhbReader())); |
45 |
|
|
*/ |
46 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Physics, new PhysicsReader())); |
47 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedPkt, new ForcedPktReader())); |
48 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader, new RunHeaderReader())); |
49 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer, new RunTrailerReader())); |
50 |
kusanagi |
1.3 |
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log, new LogReader())); |
51 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump, new VarDumpReader())); |
52 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump, new ArrDumpReader())); |
53 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump, new TabDumpReader())); |
54 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc, new TmtcReader())); |
55 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd, new McmdReader())); |
56 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCal, new CalibCalReader())); |
57 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,new CalibCalPedReader())); |
58 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCal, new CalibAcReader())); |
59 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1, new CalibTrk1Reader())); |
60 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2, new CalibTrk2Reader())); |
61 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd, new CalibTrdReader())); |
62 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof, new CalibTofReader())); |
63 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4, new CalibS4Reader())); |
64 |
kusanagi |
1.1 |
/* |
65 |
|
|
TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::HA_Header_E5, |
66 |
|
|
new E5Reader())); |
67 |
|
|
*/ |
68 |
|
|
} |
69 |
|
|
|
70 |
|
|
/** |
71 |
|
|
* Get a string with the version info of the algorithm. |
72 |
|
|
*/ |
73 |
|
|
std::string EventReader::GetVersionInfo(void) const { |
74 |
|
|
return |
75 |
kusanagi |
1.4 |
"$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 1.3 2004/07/08 12:31:42 kusanagi Exp $\n"; |
76 |
kusanagi |
1.1 |
} |
77 |
|
|
|
78 |
|
|
/** |
79 |
|
|
* Initialize the algorithm with a special run. This will initialize the |
80 |
|
|
* event reader routines for all packet types. |
81 |
|
|
*/ |
82 |
|
|
void EventReader::Init(PamelaRun *run) { |
83 |
|
|
SetInputStream(run); |
84 |
|
|
|
85 |
|
|
//Create the structure of directories and create xxx.Header.root files |
86 |
|
|
run->WriteHeaders(this, &Header); |
87 |
|
|
|
88 |
|
|
//Create the xxx.root in it's specific directory |
89 |
|
|
for (AlgorithmMap::iterator i = TechmodelAlgorithmMap.begin(); |
90 |
|
|
i != TechmodelAlgorithmMap.end(); i++) { |
91 |
|
|
cat << log4cpp::Priority::DEBUG |
92 |
|
|
<< " Initializing algo " << i->second->GetAlgorithmName() |
93 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
94 |
|
|
i->second->Init(run); |
95 |
|
|
} |
96 |
|
|
Run = dynamic_cast<TechmodelPamelaRun*>(run); |
97 |
|
|
} |
98 |
|
|
|
99 |
|
|
static void SkipToNextHeader(ifstream *); |
100 |
|
|
|
101 |
|
|
/** |
102 |
|
|
* Read the next event header, call the reader algorithms that |
103 |
|
|
* correspond to its packet type, and read the event trailer. |
104 |
|
|
*/ |
105 |
|
|
|
106 |
|
|
// 15 June 2004 ---note------------------ |
107 |
|
|
// A nice refacoring would require to not read again the packet data (which is |
108 |
|
|
// read in the UnpackPscuHeader() and put directly in the |
109 |
|
|
// EventAlgorithm->RunEvent(....) directly the data inside the packet.... will see later.... |
110 |
|
|
// 15 June 2004 ---note------------------ |
111 |
|
|
void EventReader::RunEvent(int EventNumber) { |
112 |
|
|
|
113 |
|
|
cat << log4cpp::Priority::ERROR << "Start Unpacking........" |
114 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
115 |
|
|
//for now i'll suppose that the raw file starts immediately with the right bytes |
116 |
|
|
//to insert a GetPacketStart() |
117 |
|
|
while (!InputFile->eof()){ |
118 |
|
|
try { |
119 |
|
|
if (FindStart()) { |
120 |
|
|
if(UnpackPscuHeader()){ |
121 |
|
|
const PacketType* type = Header->GetPscuHeader()->GetPacketType(); |
122 |
|
|
AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type); |
123 |
|
|
if (i != TechmodelAlgorithmMap.end()) { |
124 |
|
|
TechmodelAlgorithm *EventAlgorithm = i->second; |
125 |
|
|
EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght()); |
126 |
|
|
Run->FillTrees(type); |
127 |
|
|
Header->GetCounter()->Increment(type); |
128 |
|
|
} else { |
129 |
|
|
cat << log4cpp::Priority::ERROR |
130 |
|
|
<< "No way to read events of type" << type->GetName() |
131 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
132 |
|
|
throw Exception("No way to read events of type " + type->GetName()); |
133 |
|
|
} |
134 |
|
|
} //else { |
135 |
|
|
//**TO BE DONE** <<WRONG CRC l --- LOST PACKET>>// |
136 |
|
|
//Have to be defined a "WrongCRC" packet::Type? |
137 |
|
|
//} |
138 |
|
|
} |
139 |
|
|
} catch (...) { |
140 |
|
|
//This have to be more detailed. More exceptions type are needed. |
141 |
|
|
cat << log4cpp::Priority::ERROR << "Couldn't read the event. Skipping to the next header." |
142 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
143 |
|
|
} |
144 |
|
|
} |
145 |
|
|
Header->GetCounter()->PrintCounters(); |
146 |
|
|
} |
147 |
|
|
|
148 |
|
|
static void SkipToNextHeader(ifstream* TechmodelFile) { |
149 |
|
|
// yeah. |
150 |
|
|
} |
151 |
|
|
|
152 |
|
|
|
153 |
|
|
|
154 |
|
|
/** |
155 |
|
|
* Unpack the PSCU header from a file into the structure. |
156 |
|
|
*/ |
157 |
|
|
int EventReader::UnpackPscuHeader(void) throw (std::exception) { |
158 |
|
|
int response; |
159 |
|
|
char buff[16]; |
160 |
|
|
InputFile->read(buff, sizeof(buff)); |
161 |
|
|
|
162 |
|
|
prevPckCounter = 0; |
163 |
|
|
prevPckOBT = 0; |
164 |
|
|
|
165 |
|
|
unsigned char PacketId1 = buff[3]; |
166 |
|
|
unsigned char PacketId2 = buff[4]; |
167 |
|
|
unsigned int Counter = (((UINT32)buff[5]<<16)&0x00FF0000) + (((UINT32)buff[6]<<8)&0x0000FF00) + (((UINT32)buff[7])&0x000000FF); |
168 |
|
|
unsigned int OrbitalTime = (((UINT32)buff[8]<<24)&0xFF000000) + (((UINT32)buff[9]<<16)&0x00FF0000) + (((UINT32)buff[10]<<8)&0x0000FF00) + (((UINT32)buff[11])&0x000000FF); |
169 |
|
|
unsigned int PacketLenght = (((UINT32)buff[12]<<16)&0x00FF0000) + (((UINT32)buff[13]<<8)&0x0000FF00) + (((UINT32)buff[14])&0x000000FF); |
170 |
|
|
unsigned char CRC = buff[15]; |
171 |
|
|
|
172 |
|
|
if (Counter < prevPckCounter){ |
173 |
|
|
cat << log4cpp::Priority::ERROR |
174 |
|
|
<< " Packet counter is less than before of " << (prevPckCounter - Counter) |
175 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
176 |
|
|
} |
177 |
|
|
|
178 |
|
|
if (OrbitalTime < prevPckOBT){ |
179 |
|
|
cat << log4cpp::Priority::WARN |
180 |
|
|
<< " Orbital Time is less than before of " << (prevPckOBT - OrbitalTime) |
181 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
182 |
|
|
} |
183 |
|
|
|
184 |
|
|
|
185 |
|
|
if ((BYTE)CM_Compute_CRC16(0, (BYTE*)&buff, 15) == (BYTE)buff[15]){ |
186 |
|
|
long int initPos = InputFile->tellg(); |
187 |
|
|
long int finalPos; |
188 |
|
|
Header->GetPscuHeader()->SetPacketId(PacketId1, PacketId2); |
189 |
|
|
Header->GetPscuHeader()->SetCounter(Counter); |
190 |
|
|
Header->GetPscuHeader()->SetOrbitalTime(OrbitalTime); |
191 |
|
|
//PacketLength is the length of the whole DATApacket starting from the first byte after the header |
192 |
|
|
//plus the CRC legth (which varies for each type of packet) |
193 |
|
|
Header->GetPscuHeader()->SetPacketLenght(PacketLenght); |
194 |
|
|
Header->GetPscuHeader()->SetCRC(CRC); |
195 |
|
|
|
196 |
|
|
//commented out because of the above test code |
197 |
|
|
InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur); |
198 |
|
|
FindStart(); |
199 |
|
|
finalPos = (long int)InputFile->tellg() - (1 + initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght())); |
200 |
|
|
if(finalPos == 0){ |
201 |
|
|
cat << log4cpp::Priority::INFO << " Correct packet length \n" |
202 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
203 |
|
|
} |
204 |
|
|
if (finalPos > 0 && finalPos < 64) { |
205 |
|
|
cat << log4cpp::Priority::WARN |
206 |
|
|
<< " Correct packet length: Padded of " << finalPos << " bytes \n" |
207 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
208 |
|
|
} |
209 |
|
|
if (finalPos > 64){ |
210 |
|
|
cat << log4cpp::Priority::ERROR << " Wrong packet length? (because of wrong padding?) \n" |
211 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
212 |
|
|
} |
213 |
|
|
InputFile->seekg(initPos, std::ios::beg); |
214 |
|
|
response = true; |
215 |
|
|
} else { |
216 |
|
|
cat << log4cpp::Priority::ERROR << " WRONG CRC FOR HEADER " |
217 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
218 |
|
|
response = false; |
219 |
|
|
} |
220 |
|
|
|
221 |
|
|
char tmpId1[4]; |
222 |
|
|
char tmpId2[4]; |
223 |
|
|
char tmpLength[12]; |
224 |
|
|
char tmpStart[100]; |
225 |
|
|
char tmpCRC[4]; |
226 |
|
|
sprintf(tmpId1, "%02X", PacketId1); |
227 |
|
|
sprintf(tmpId2, "%02X", PacketId2); |
228 |
kusanagi |
1.4 |
sprintf(tmpLength, "%06X", Header->GetPscuHeader()->GetPacketLenght()); |
229 |
kusanagi |
1.1 |
sprintf(tmpStart, "%X", ((long int)(InputFile->tellg()) - 16)); |
230 |
|
|
sprintf(tmpCRC, "%02X", CRC); |
231 |
|
|
cat << log4cpp::Priority::INFO |
232 |
|
|
<< "\n Packet Counter (decimal) : " << Counter |
233 |
|
|
<< "\n Id1 - Id2 : " << tmpId1 << " - " << tmpId2 |
234 |
|
|
<< "\n Orbital Time (decimal) : " << OrbitalTime |
235 |
|
|
<< "\n Lenght : " << tmpLength |
236 |
|
|
<< "\n CRC : " << tmpCRC |
237 |
|
|
<< "\n Header Start Position : " << tmpStart |
238 |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
239 |
|
|
return response; |
240 |
|
|
} |
241 |
|
|
|
242 |
|
|
/** |
243 |
|
|
* Unpack the trailer of a PSCU event into the structure. |
244 |
|
|
*/ |
245 |
|
|
void EventReader::UnpackPscuTrailer(void) throw (std::exception) { |
246 |
|
|
|
247 |
|
|
} |
248 |
|
|
|
249 |
|
|
/** |
250 |
|
|
* Find the next starting poin for the PSCU event looking for a {0xFA, 0xFE, 0xDE} sequence |
251 |
|
|
*/ |
252 |
|
|
int EventReader::FindStart(void) throw (std::exception) { |
253 |
|
|
//search an hexadecimal sequence in a file |
254 |
|
|
//subSign ------> pointer to the sequence buffer |
255 |
|
|
//subSignDim ------> dimension of the buffer |
256 |
|
|
// at exit |
257 |
|
|
// return true if founds a match (else false) |
258 |
|
|
// subF point rigth after the match, if found. Else EOF. |
259 |
|
|
//Maurizio 15/11/2002----------------------- |
260 |
|
|
const unsigned char subSign[3]={0xFA, 0xFE, 0xDE}; |
261 |
|
|
int subSignDim = 3; |
262 |
|
|
//------------------------------------------ |
263 |
|
|
int subIndex = 0; |
264 |
|
|
char dataByte; |
265 |
|
|
/*InputFile->read(fileByte, 1); |
266 |
|
|
if ( *fileByte == *(subSign+subIndex)){ |
267 |
|
|
if (subIndex++ == (subSignDim-1)) { |
268 |
|
|
InputFile->seekg(-(subIndex), std::ios::cur); |
269 |
|
|
return 1; |
270 |
|
|
} |
271 |
|
|
}*/ |
272 |
|
|
while (!InputFile->eof()) { |
273 |
|
|
InputFile->get(dataByte); |
274 |
|
|
if (dataByte == (char)(*(subSign+subIndex))){ |
275 |
|
|
if (subIndex++ == (subSignDim-1)) { |
276 |
|
|
InputFile->seekg(-(subIndex), std::ios::cur); |
277 |
|
|
return 1; |
278 |
|
|
} |
279 |
|
|
} else { |
280 |
|
|
if (InputFile->eof()) { |
281 |
|
|
throw Exception("Extra bytes over the last packets"); |
282 |
|
|
} else { |
283 |
|
|
return 0; |
284 |
|
|
} |
285 |
|
|
//InputFile->seekg(-(subIndex+1), std::ios::cur); |
286 |
|
|
//subIndex = 0; |
287 |
|
|
} |
288 |
|
|
} |
289 |
|
|
return 0; |
290 |
|
|
} |
291 |
|
|
|
292 |
|
|
ClassImp(EventReader) |