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

Diff of /yoda/techmodel/EventReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2.5 by kusanagi, Thu Dec 9 08:48:41 2004 UTC revision 4.4 by kusanagi, Sat May 28 10:44:11 2005 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: kusanagi $
3   * $Date: 2004/12/03 22:08:00 $   * $Date: 2005/05/28 08:30:11 $
4   * $Revision: 2.4 $   * $Revision: 4.3 $
5   *   *
6   * Implementation of the functions of a sample Algorithm class.   * Implementation of the functions of a sample Algorithm class.
7   * This file can be used as a templace to develop your own algorithm.   * This file can be used as a templace to develop your own algorithm.
# Line 40  EventReader::EventReader(int packetsLimi Line 40  EventReader::EventReader(int packetsLimi
40    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1,  new CalibCalPulse1Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1,  new CalibCalPulse1Reader));
41    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2,  new CalibCalPulse2Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2,  new CalibCalPulse2Reader));
42    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Physics,         new PhysicsReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Physics,         new PhysicsReader));
   TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth,    new CalibTrkBothReader));  
43    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));
44    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2,       new CalibTrk2Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2,       new CalibTrk2Reader));
   TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd,        new CalibTrdReader));  
45    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));
46    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));
47    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));
48    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibAc,         new CalibAcReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac1,      new Calib1_Ac1Reader));
49      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac1,      new Calib2_Ac1Reader));
50      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac2,      new Calib1_Ac2Reader));
51      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new Calib2_Ac2Reader));
52    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader));
53    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader));
54    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader));
# Line 64  EventReader::EventReader(int packetsLimi Line 65  EventReader::EventReader(int packetsLimi
65    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader));
66    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader));
67    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader));
68    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcInit,          new AcInitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac1Init,         new Ac1InitReader));
69    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader));
70    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader));
71    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader));
72    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader));
73    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::NdInit,          new NdInitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::NdInit,          new NdInitReader));
74    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init,          new S4InitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init,          new S4InitReader));
75      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Init,         new Ac2InitReader));
76    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader));
77    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader));
78    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader));
79    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader));
80    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader));
81    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Alarm,         new S4AlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Alarm,         new S4AlarmReader));
82      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbT,            new TsbTReader));
83      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbB,            new TsbBReader));
84  }  }
85    
86  /**  /**
# Line 84  EventReader::EventReader(int packetsLimi Line 88  EventReader::EventReader(int packetsLimi
88   */   */
89  std::string EventReader::GetVersionInfo(void) const {  std::string EventReader::GetVersionInfo(void) const {
90    return    return
91      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 2.4 2004/12/03 22:08:00 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 4.3 2005/05/28 08:30:11 kusanagi Exp $\n";
92  }  }
93    
94  /**  /**
# Line 139  void EventReader::RunEvent(int EventNumb Line 143  void EventReader::RunEvent(int EventNumb
143                  Header->GetCounter()->Increment(type);                  Header->GetCounter()->Increment(type);
144                  logger->info(Header->GetPscuHeader()->Print());                  logger->info(Header->GetPscuHeader()->Print());
145              } else {              } else {
146                  throw NotExistingAlgorithmException(type->GetName().c_str()); //to exctract to an higher level and delete the logger!                  oss.str("");
147                    oss << "\n No way to read events of type  " << type->GetName().c_str() << Header->GetPscuHeader()->Print();
148                    throw NotExistingAlgorithmException(oss.str().c_str()); //to exctract to an higher level and delete the logger!
149              }              }
150        }        }
151      // In case of exception have to save the packet in a specific root file??      // In case of exception have to save the packet in a specific root file??
# Line 149  void EventReader::RunEvent(int EventNumb Line 155  void EventReader::RunEvent(int EventNumb
155        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
156      } catch (WrongCRCHeaderException exc) {      } catch (WrongCRCHeaderException exc) {
157        oss.str("");        oss.str("");
158        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print();
159        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
160      } catch (WrongCRCException exc) {      } catch (WrongCRCException exc) {
161        oss.str("");        oss.str("");
# Line 171  void EventReader::RunEvent(int EventNumb Line 177  void EventReader::RunEvent(int EventNumb
177        logger->error("Couldn't read the event. Skipping to the next header. \n");        logger->error("Couldn't read the event. Skipping to the next header. \n");
178      }      }
179      if ((step%1000) == 0) std::cout << step/1000 << "K \n";      if ((step%1000) == 0) std::cout << step/1000 << "K \n";
180        oss.str("");
181        oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";
182        logger->info(oss.str().c_str());
183    }    }
184      Header->GetCounter()->PrintCounters();      Header->GetCounter()->PrintCounters();
185  }  }
# Line 184  static void SkipToNextHeader(ifstream* T Line 193  static void SkipToNextHeader(ifstream* T
193  /**  /**
194   * Unpack the PSCU header from a file into the structure.   * Unpack the PSCU header from a file into the structure.
195   */   */
196  void EventReader::UnpackPscuHeader(void) throw (WrongCRCException, LengthException) {  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {
197    stringstream oss;    stringstream oss;
198    int response = false;    int response = false;
199    char buff[16];    char buff[16];
# Line 202  void EventReader::UnpackPscuHeader(void) Line 211  void EventReader::UnpackPscuHeader(void)
211    if (Counter < prevPckCounter){    if (Counter < prevPckCounter){
212      oss.str("");      oss.str("");
213      oss << "Packet counter is less than before of " << (prevPckCounter - Counter);      oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
214      logger->warn(oss.str().c_str());      logger->error(oss.str().c_str());
215      }
216    
217      if (Counter > prevPckCounter + 1){
218        oss.str("");
219        oss << "Packet counter is greater than before of " << (Counter - prevPckCounter);
220        logger->error(oss.str().c_str());
221      }
222      
223      if ((OrbitalTime == prevPckOBT) & (PacketId1 == 0x10)){
224        oss.str("");
225        oss << "Onboard Time of this packet is equal to the previous packet OBT";
226        logger->error(oss.str().c_str());
227        logger->error(Header->GetPscuHeader()->Print());
228    }    }
229    
230    if (OrbitalTime < prevPckOBT){    if (OrbitalTime < prevPckOBT){
231      oss.str("");      oss.str("");
232      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
233      logger->warn(oss.str().c_str());      logger->error(oss.str().c_str());
234    }    }
235    
236    if (((BYTE)CM_Compute_CRC16(0, (BYTE*)&buff, 15) == (BYTE)buff[15]) && (PacketId1 == PacketId2)){    if (((BYTE)CM_Compute_CRC16(0, (BYTE*)&buff, 15) == (BYTE)buff[15]) && (PacketId1 == PacketId2)){
# Line 226  void EventReader::UnpackPscuHeader(void) Line 248  void EventReader::UnpackPscuHeader(void)
248      Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));      Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));
249            
250      //commented out because of the above test code      //commented out because of the above test code
251      InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);      /*InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);
252      FindStart();      if (FindStart()) {
253      finalPos =  (long int)InputFile->tellg() - (1 + initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));          finalPos =  (long int)InputFile->tellg() - (initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));
254          if(finalPos == 0){          if(finalPos == 0){
255              logger->debug(_T("Correct packet length"));              logger->debug(_T(" Correct packet length"));
256              }              }
257          if (finalPos > 0 && finalPos < 64) {          if (finalPos > 0 && finalPos < 64) {
258              oss.str("");              oss.str("");
259              oss << " Correct packet length: Padded of " << finalPos << " bytes";              oss << "Correct packet length: Padded of " << finalPos << " bytes";
260              logger->debug(oss.str().c_str());              logger->debug(oss.str().c_str());
261          }          }
262          if (finalPos > 64){          if (finalPos > 64){
263              oss.str("");              oss.str("");
264              oss << "The begin of the next packet is far more than 64 byte from the end of the previous."              oss << "\n The begin of the next packet is far more than 64 byte from the end of the previous."
265                  << Header->GetPscuHeader()->Print();                  << "\n Below the is the now unpacking packet";
266              logger->error(oss.str().c_str());              logger->error(oss.str().c_str());
267                logger->error(Header->GetPscuHeader()->Print());
268              //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");              //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");
269          }          }
270          InputFile->seekg(initPos, std::ios::beg);      }
271        else {
272            logger->debug(_T(" END OF FILE"));        
273        }*/
274            
275            //InputFile->seekg(initPos, std::ios::beg);
276    } else {    } else {
277         /*Here i should extract the block of Data for later analysis */
278      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
279      throw WrongCRCHeaderException();      oss.str("");
280        oss << "CRC Header Error on packet:" << PscuHeader::Print(buff);
281        throw WrongCRCHeaderException(oss.str().c_str());
282    }    }
283    
284  /*    char tmpId1[4];  /*    char tmpId1[4];
# Line 281  void EventReader::UnpackPscuTrailer(void Line 312  void EventReader::UnpackPscuTrailer(void
312  /**  /**
313   * Find the next starting poin for the PSCU event looking for a {0xFA, 0xFE, 0xDE} sequence   * Find the next starting poin for the PSCU event looking for a {0xFA, 0xFE, 0xDE} sequence
314   */   */
315  int EventReader::FindStart(void) throw (std::exception) {  bool EventReader::FindStart(void) throw (std::exception) {
316  //search an hexadecimal sequence in a file  //search an hexadecimal sequence in a file
317  //subSign    ------> pointer to the sequence buffer  //subSign    ------> pointer to the sequence buffer
318  //subSignDim ------> dimension of the buffer  //subSignDim ------> dimension of the buffer
# Line 294  int EventReader::FindStart(void) throw ( Line 325  int EventReader::FindStart(void) throw (
325  //------------------------------------------  //------------------------------------------
326    int subIndex = 0;    int subIndex = 0;
327    char dataByte;    char dataByte;
328    
329    int  buffSize = 100;    int  buffSize = 64;
330    int  index = 0;    int  index = 0;
331    int  loop = -1;    int  loop = -1;
332    char buffer[buffSize];    char buffer[buffSize];
333      bool flagOverPad = false;
334    
335    while (!InputFile->eof()) {    while (!InputFile->eof()) {
336      InputFile->read(buffer, sizeof(buffer));      InputFile->read(buffer, sizeof(buffer));
# Line 309  int EventReader::FindStart(void) throw ( Line 341  int EventReader::FindStart(void) throw (
341          if (dataByte == (char)(*(subSign+subIndex))){          if (dataByte == (char)(*(subSign+subIndex))){
342              if (subIndex++ == (subSignDim-1)) {              if (subIndex++ == (subSignDim-1)) {
343                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
344                  return 1;                  if (flagOverPad){
345                        oss.str("");
346                        oss << "\n This packet beginning is farther than 64 byte from the end of the previous."
347                            << "\n Below the is the last already unpacked packet";
348                        logger->error(oss.str().c_str());
349                        logger->error(Header->GetPscuHeader()->Print());
350                    }
351                    return true;
352              }              }
353          } else {          } else {
354              index = index - (subIndex);              index = index - (subIndex);
# Line 317  int EventReader::FindStart(void) throw ( Line 356  int EventReader::FindStart(void) throw (
356          }          }
357      }      }
358      //Needs to guarantee the overap of the buffer(s) in several loop      //Needs to guarantee the overap of the buffer(s) in several loop
359        flagOverPad = true;
360      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
361    }    }
362    return 0;    return false;
363  }  }
364    
365  ClassImp(EventReader)  ClassImp(EventReader)

Legend:
Removed from v.2.5  
changed lines
  Added in v.4.4

  ViewVC Help
Powered by ViewVC 1.1.23