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

Annotation of /yoda/techmodel/EventReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 4.2 - (hide annotations) (download)
Mon May 23 12:22:00 2005 UTC (19 years, 6 months ago) by kusanagi
Branch: MAIN
CVS Tags: YODA4_3/02, YODA4_3/01
Changes since 4.1: +4 -4 lines
Minor change on the description in case he packets are more distant than 64 bytes.

1 kusanagi 1.1 /** @file
2 kusanagi 1.2 * $Author: kusanagi $
3 kusanagi 4.2 * $Date: 2005/03/22 19:54:22 $
4     * $Revision: 4.1 $
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 kusanagi 1.10 #include <log4cxx/logger.h>
11 kusanagi 1.1 #include <fstream>
12     #include "EventReader.h"
13     #include "ReaderAlgorithms.h"
14 kusanagi 2.4
15    
16 kusanagi 1.1 extern "C" {
17 kusanagi 2.4 #include "CRC.h"
18 kusanagi 1.1 }
19    
20     using namespace pamela;
21     using namespace pamela::techmodel;
22    
23 kusanagi 1.10 static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.EventReader"));
24    
25     unsigned int EventReader::maxPackets = 0;
26 kusanagi 1.7 unsigned int EventReader::prevPckCounter = 0;
27     unsigned int EventReader::prevPckOBT = 0;
28 kusanagi 1.10
29    
30 kusanagi 1.1 /**
31     * Constructor.
32     */
33 kusanagi 1.10 EventReader::EventReader(int packetsLimit = -1):
34     TechmodelAlgorithm(0, "TechmodelEventReader"){
35     EventReader::maxPackets = packetsLimit;
36     logger->debug(_T("Constructor"));
37 kusanagi 1.1 Header = new EventHeader();
38 kusanagi 1.5
39 kusanagi 2.5 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun, new PhysEndRunReader));
40     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1, new CalibCalPulse1Reader));
41     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2, new CalibCalPulse2Reader));
42     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Physics, new PhysicsReader));
43     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth, new CalibTrkBothReader));
44     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1, new CalibTrk1Reader));
45     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2, new CalibTrk2Reader));
46     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd, new CalibTrdReader));
47     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof, new CalibTofReader));
48     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4, new CalibS4Reader));
49     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed, new CalibCalPedReader));
50 kusanagi 2.6 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac1, new Calib1_Ac1Reader));
51     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac1, new Calib2_Ac1Reader));
52     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac2, new Calib1_Ac2Reader));
53     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2, new Calib2_Ac2Reader));
54 kusanagi 2.5 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader, new RunHeaderReader));
55     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer, new RunTrailerReader));
56     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader, new CalibHeaderReader));
57     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer, new CalibTrailerReader));
58     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader, new InitHeaderReader));
59     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer, new InitTrailerReader));
60     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk, new EventTrkReader));
61     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk, new TestTrkReader));
62     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof, new TestTofReader));
63     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log, new LogReader));
64     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump, new VarDumpReader));
65     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump, new ArrDumpReader));
66     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump, new TabDumpReader));
67     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc, new TmtcReader));
68     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd, new McmdReader));
69     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd, new ForcedFECmdReader));
70 kusanagi 2.6 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac1Init, new Ac1InitReader));
71 kusanagi 2.5 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit, new CalInitReader));
72     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit, new TrkInitReader));
73     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit, new TofInitReader));
74     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit, new TrgInitReader));
75     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::NdInit, new NdInitReader));
76     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init, new S4InitReader));
77 kusanagi 2.6 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Init, new Ac2InitReader));
78 kusanagi 2.5 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm, new CalAlarmReader));
79     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm, new AcAlarmReader));
80     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm, new TrkAlarmReader));
81     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm, new TrgAlarmReader));
82     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm, new TofAlarmReader));
83     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Alarm, new S4AlarmReader));
84 kusanagi 2.11 TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbT, new TsbTReader));
85     TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbB, new TsbBReader));
86 kusanagi 1.1 }
87    
88     /**
89     * Get a string with the version info of the algorithm.
90     */
91     std::string EventReader::GetVersionInfo(void) const {
92     return
93 kusanagi 4.2 "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 4.1 2005/03/22 19:54:22 kusanagi Exp $\n";
94 kusanagi 1.1 }
95    
96     /**
97     * Initialize the algorithm with a special run. This will initialize the
98     * event reader routines for all packet types.
99     */
100     void EventReader::Init(PamelaRun *run) {
101     SetInputStream(run);
102    
103     //Create the structure of directories and create xxx.Header.root files
104     run->WriteHeaders(this, &Header);
105    
106     //Create the xxx.root in it's specific directory
107     for (AlgorithmMap::iterator i = TechmodelAlgorithmMap.begin();
108     i != TechmodelAlgorithmMap.end(); i++) {
109 kusanagi 2.4 oss.str("");
110 kusanagi 1.10 oss << "Initializing algo " << i->second->GetAlgorithmName();
111     logger->debug(oss.str().c_str());
112 kusanagi 1.1 i->second->Init(run);
113     }
114     Run = dynamic_cast<TechmodelPamelaRun*>(run);
115     }
116    
117     static void SkipToNextHeader(ifstream *);
118    
119     /**
120     * Read the next event header, call the reader algorithms that
121     * correspond to its packet type, and read the event trailer.
122     */
123    
124     // 15 June 2004 ---note------------------
125     // A nice refacoring would require to not read again the packet data (which is
126 kusanagi 2.4 // read in the UnpackPscuHeader() and put directly in the
127 kusanagi 1.1 // EventAlgorithm->RunEvent(....) directly the data inside the packet.... will see later....
128     // 15 June 2004 ---note------------------
129     void EventReader::RunEvent(int EventNumber) {
130     //for now i'll suppose that the raw file starts immediately with the right bytes
131     //to insert a GetPacketStart()
132 kusanagi 1.10 stringstream oss;
133     int step = 0;
134 kusanagi 2.5 const PacketType* type;
135 kusanagi 1.10 while (!InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){
136 kusanagi 1.1 try {
137     if (FindStart()) {
138 kusanagi 2.4 UnpackPscuHeader();
139 kusanagi 2.5 type = Header->GetPscuHeader()->GetPacketType();
140 kusanagi 1.1 AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type);
141     if (i != TechmodelAlgorithmMap.end()) {
142 kusanagi 2.5 TechmodelAlgorithm *EventAlgorithm(i->second);
143 kusanagi 1.1 EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());
144     Run->FillTrees(type);
145     Header->GetCounter()->Increment(type);
146 kusanagi 2.4 logger->info(Header->GetPscuHeader()->Print());
147 kusanagi 1.1 } else {
148 kusanagi 2.8 oss.str("");
149     oss << "\n No way to read events of type " << type->GetName().c_str() << Header->GetPscuHeader()->Print();
150     throw NotExistingAlgorithmException(oss.str().c_str()); //to exctract to an higher level and delete the logger!
151 kusanagi 2.3 }
152     }
153 kusanagi 2.4 // In case of exception have to save the packet in a specific root file??
154     } catch (NotExistingAlgorithmException exc) {
155 kusanagi 2.5 oss.str("");
156     oss << exc.print() << " " << Header->GetPscuHeader()->Print();
157     logger->error(oss.str().c_str());
158 kusanagi 2.4 } catch (WrongCRCHeaderException exc) {
159 kusanagi 2.5 oss.str("");
160 kusanagi 2.8 oss << exc.print();
161 kusanagi 2.5 logger->error(oss.str().c_str());
162 kusanagi 2.4 } catch (WrongCRCException exc) {
163 kusanagi 2.5 oss.str("");
164     oss << exc.print() << " " << Header->GetPscuHeader()->Print();
165     logger->error(oss.str().c_str());
166 kusanagi 2.4 } catch (UnidentifiedPacketException exc) {
167 kusanagi 2.5 oss.str("");
168     oss << exc.print() << " " << Header->GetPscuHeader()->Print();
169     logger->error(oss.str().c_str());
170 kusanagi 2.4 } catch (NotExistingCounterException exc) {
171 kusanagi 2.5 oss.str("");
172     oss << exc.print() << " " << Header->GetPscuHeader()->Print();
173     logger->error(oss.str().c_str());
174 kusanagi 2.4 } catch (LengthException exc) {
175 kusanagi 2.5 oss.str("");
176     oss << exc.print() << " " << Header->GetPscuHeader()->Print();
177     logger->error(oss.str().c_str());
178 kusanagi 1.1 } catch (...) {
179 kusanagi 2.4 logger->error("Couldn't read the event. Skipping to the next header. \n");
180 kusanagi 1.1 }
181 kusanagi 2.1 if ((step%1000) == 0) std::cout << step/1000 << "K \n";
182 kusanagi 2.10 oss.str("");
183     oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";
184     logger->info(oss.str().c_str());
185 kusanagi 1.1 }
186 kusanagi 1.10 Header->GetCounter()->PrintCounters();
187 kusanagi 1.1 }
188    
189     static void SkipToNextHeader(ifstream* TechmodelFile) {
190     // yeah.
191     }
192    
193    
194    
195     /**
196     * Unpack the PSCU header from a file into the structure.
197     */
198 kusanagi 2.8 void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {
199 kusanagi 1.10 stringstream oss;
200 kusanagi 2.4 int response = false;
201 kusanagi 1.1 char buff[16];
202     InputFile->read(buff, sizeof(buff));
203    
204    
205     unsigned char PacketId1 = buff[3];
206     unsigned char PacketId2 = buff[4];
207     unsigned int Counter = (((UINT32)buff[5]<<16)&0x00FF0000) + (((UINT32)buff[6]<<8)&0x0000FF00) + (((UINT32)buff[7])&0x000000FF);
208     unsigned int OrbitalTime = (((UINT32)buff[8]<<24)&0xFF000000) + (((UINT32)buff[9]<<16)&0x00FF0000) + (((UINT32)buff[10]<<8)&0x0000FF00) + (((UINT32)buff[11])&0x000000FF);
209     unsigned int PacketLenght = (((UINT32)buff[12]<<16)&0x00FF0000) + (((UINT32)buff[13]<<8)&0x0000FF00) + (((UINT32)buff[14])&0x000000FF);
210     unsigned char CRC = buff[15];
211 kusanagi 1.9 unsigned char FileOffset = buff[15];
212 kusanagi 1.1
213     if (Counter < prevPckCounter){
214 kusanagi 2.4 oss.str("");
215 kusanagi 1.10 oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
216 kusanagi 2.10 logger->error(oss.str().c_str());
217     }
218    
219     if (Counter > prevPckCounter + 1){
220     oss.str("");
221     oss << "Packet counter is greater than before of " << (Counter - prevPckCounter);
222     logger->error(oss.str().c_str());
223 kusanagi 1.1 }
224 kusanagi 3.1
225     if ((OrbitalTime == prevPckOBT) & (PacketId1 == 0x10)){
226     oss.str("");
227     oss << "Onboard Time of this packet is equal to the previous packet OBT";
228     logger->error(oss.str().c_str());
229     logger->error(Header->GetPscuHeader()->Print());
230     }
231 kusanagi 1.1
232     if (OrbitalTime < prevPckOBT){
233 kusanagi 2.4 oss.str("");
234     oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
235 kusanagi 2.10 logger->error(oss.str().c_str());
236 kusanagi 1.1 }
237    
238 kusanagi 1.7 if (((BYTE)CM_Compute_CRC16(0, (BYTE*)&buff, 15) == (BYTE)buff[15]) && (PacketId1 == PacketId2)){
239     prevPckCounter = Counter;
240     prevPckOBT = OrbitalTime;
241 kusanagi 1.1 long int initPos = InputFile->tellg();
242     long int finalPos;
243     Header->GetPscuHeader()->SetPacketId(PacketId1, PacketId2);
244     Header->GetPscuHeader()->SetCounter(Counter);
245     Header->GetPscuHeader()->SetOrbitalTime(OrbitalTime);
246     //PacketLength is the length of the whole DATApacket starting from the first byte after the header
247     //plus the CRC legth (which varies for each type of packet)
248     Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
249     Header->GetPscuHeader()->SetCRC(CRC);
250 kusanagi 1.9 Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));
251 kusanagi 1.1
252     //commented out because of the above test code
253 kusanagi 4.1 /*InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);
254 kusanagi 2.7 if (FindStart()) {
255     finalPos = (long int)InputFile->tellg() - (initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));
256 kusanagi 1.1 if(finalPos == 0){
257 kusanagi 2.7 logger->debug(_T(" Correct packet length"));
258 kusanagi 1.1 }
259     if (finalPos > 0 && finalPos < 64) {
260 kusanagi 2.4 oss.str("");
261 kusanagi 2.10 oss << "Correct packet length: Padded of " << finalPos << " bytes";
262 kusanagi 2.4 logger->debug(oss.str().c_str());
263 kusanagi 1.1 }
264     if (finalPos > 64){
265 kusanagi 2.4 oss.str("");
266 kusanagi 2.10 oss << "\n The begin of the next packet is far more than 64 byte from the end of the previous."
267     << "\n Below the is the now unpacking packet";
268 kusanagi 2.5 logger->error(oss.str().c_str());
269 kusanagi 2.10 logger->error(Header->GetPscuHeader()->Print());
270 kusanagi 2.4 //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");
271 kusanagi 1.1 }
272 kusanagi 2.7 }
273     else {
274     logger->debug(_T(" END OF FILE"));
275 kusanagi 4.1 }*/
276 kusanagi 2.7
277 kusanagi 4.1 //InputFile->seekg(initPos, std::ios::beg);
278 kusanagi 1.1 } else {
279 kusanagi 2.8 /*Here i should extract the block of Data for later analysis */
280 kusanagi 1.8 InputFile->seekg(-(13), std::ios::cur);
281 kusanagi 2.8 oss.str("");
282     oss << "CRC Header Error on packet:" << PscuHeader::Print(buff);
283     throw WrongCRCHeaderException(oss.str().c_str());
284 kusanagi 1.1 }
285 kusanagi 2.4
286     /* char tmpId1[4];
287 kusanagi 1.1 char tmpId2[4];
288     char tmpLength[12];
289     char tmpStart[100];
290     char tmpCRC[4];
291     sprintf(tmpId1, "%02X", PacketId1);
292     sprintf(tmpId2, "%02X", PacketId2);
293 kusanagi 1.4 sprintf(tmpLength, "%06X", Header->GetPscuHeader()->GetPacketLenght());
294 kusanagi 1.9 sprintf(tmpStart, "%X", Header->GetPscuHeader()->GetFileOffset());
295 kusanagi 1.1 sprintf(tmpCRC, "%02X", CRC);
296 kusanagi 1.10 oss.flush();
297     oss << "\n Packet Counter (decimal) : " << Counter
298 kusanagi 1.1 << "\n Id1 - Id2 : " << tmpId1 << " - " << tmpId2
299     << "\n Orbital Time (decimal) : " << OrbitalTime
300     << "\n Lenght : " << tmpLength
301     << "\n CRC : " << tmpCRC
302 kusanagi 1.10 << "\n Header Start Position : " << tmpStart;
303     logger->info(oss.str().c_str());
304 kusanagi 2.4 */
305 kusanagi 1.1 }
306    
307     /**
308     * Unpack the trailer of a PSCU event into the structure.
309     */
310     void EventReader::UnpackPscuTrailer(void) throw (std::exception) {
311    
312     }
313    
314     /**
315     * Find the next starting poin for the PSCU event looking for a {0xFA, 0xFE, 0xDE} sequence
316     */
317 kusanagi 2.7 bool EventReader::FindStart(void) throw (std::exception) {
318 kusanagi 1.1 //search an hexadecimal sequence in a file
319     //subSign ------> pointer to the sequence buffer
320     //subSignDim ------> dimension of the buffer
321     // at exit
322     // return true if founds a match (else false)
323     // subF point rigth after the match, if found. Else EOF.
324     //Maurizio 15/11/2002-----------------------
325 kusanagi 2.4 const int subSignDim = 3;
326     const unsigned char subSign[subSignDim]={0xFA, 0xFE, 0xDE};
327 kusanagi 1.1 //------------------------------------------
328     int subIndex = 0;
329     char dataByte;
330 kusanagi 4.1
331     int buffSize = 64;
332 kusanagi 1.10 int index = 0;
333 kusanagi 2.4 int loop = -1;
334 kusanagi 1.10 char buffer[buffSize];
335 kusanagi 4.1 bool flagOverPad = false;
336 kusanagi 1.10
337     while (!InputFile->eof()) {
338     InputFile->read(buffer, sizeof(buffer));
339     index = 0;
340 kusanagi 2.4 loop++;
341 kusanagi 1.10 while (index < buffSize){
342     dataByte = buffer[index++];
343     if (dataByte == (char)(*(subSign+subIndex))){
344     if (subIndex++ == (subSignDim-1)) {
345 kusanagi 2.4 InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
346 kusanagi 4.1 if (flagOverPad){
347     oss.str("");
348     oss << "\n This packet beginning is farther than 64 byte from the end of the previous."
349 kusanagi 4.2 << "\n Below the is the last already unpacked packet";
350 kusanagi 4.1 logger->error(oss.str().c_str());
351     logger->error(Header->GetPscuHeader()->Print());
352     }
353 kusanagi 2.7 return true;
354 kusanagi 2.4 }
355 kusanagi 1.10 } else {
356 kusanagi 2.4 index = index - (subIndex);
357 kusanagi 1.10 subIndex = 0;
358     }
359     }
360 kusanagi 2.4 //Needs to guarantee the overap of the buffer(s) in several loop
361 kusanagi 4.1 flagOverPad = true;
362 kusanagi 2.4 InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
363 kusanagi 2.7 }
364     return false;
365 kusanagi 1.1 }
366    
367     ClassImp(EventReader)

  ViewVC Help
Powered by ViewVC 1.1.23