/[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.8 by kusanagi, Tue Jan 11 14:59:31 2005 UTC revision 5.1 by kusanagi, Fri Sep 16 12:05:49 2005 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: kusanagi $
3   * $Date: 2004/12/22 11:38:35 $   * $Date: 2005/08/29 09:46:13 $
4   * $Revision: 2.7 $   * $Revision: 5.0 $
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));
# Line 81  EventReader::EventReader(int packetsLimi Line 79  EventReader::EventReader(int packetsLimi
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 88  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.7 2004/12/22 11:38:35 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 5.0 2005/08/29 09:46:13 kusanagi Exp $\n";
92  }  }
93    
94  /**  /**
# Line 118  static void SkipToNextHeader(ifstream *) Line 118  static void SkipToNextHeader(ifstream *)
118   * Read the next event header, call the reader algorithms that   * Read the next event header, call the reader algorithms that
119   * correspond to its packet type, and read the event trailer.   * correspond to its packet type, and read the event trailer.
120   */   */
   
 // 15 June 2004 ---note------------------  
 // A nice refacoring would require to not read again the packet data (which is  
 // read in the UnpackPscuHeader() and put directly in  the  
 // EventAlgorithm->RunEvent(....) directly the data inside the packet.... will see later....  
 // 15 June 2004 ---note------------------  
121  void EventReader::RunEvent(int EventNumber) {  void EventReader::RunEvent(int EventNumber) {
   //for now i'll suppose that the raw file starts immediately with the right bytes  
   //to insert a GetPacketStart()  
122    stringstream oss;    stringstream oss;
123    int step = 0;    int step = 0;
124    const PacketType* type;    const PacketType* type;
# Line 161  void EventReader::RunEvent(int EventNumb Line 153  void EventReader::RunEvent(int EventNumb
153        oss.str("");        oss.str("");
154        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
155        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
156          InputFile->seekg( (-1)*(Header->GetPscuHeader()->GetPacketLenght() + 15) , std::ios::cur);
157      } catch (UnidentifiedPacketException exc) {      } catch (UnidentifiedPacketException exc) {
158        oss.str("");        oss.str("");
159        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
# Line 173  void EventReader::RunEvent(int EventNumb Line 166  void EventReader::RunEvent(int EventNumb
166        oss.str("");        oss.str("");
167        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
168        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
169        } catch (BackwardCounterException exc) {
170          oss.str("");
171          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
172          logger->error(oss.str().c_str());
173      } catch (...) {      } catch (...) {
174        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");
175      }      }
176      if ((step%1000) == 0) std::cout << step/1000 << "K \n";      if ((step%1000) == 0) std::cout << step/1000 << "K \n";
177        oss.str("");
178        oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";
179        logger->info(oss.str().c_str());
180    }    }
181      Header->GetCounter()->PrintCounters();      Header->GetCounter()->PrintCounters();
182  }  }
183    
 static void SkipToNextHeader(ifstream* TechmodelFile) {  
   //  yeah.  
 }  
   
   
   
184  /**  /**
185   * Unpack the PSCU header from a file into the structure.   * Unpack the PSCU header from a file into the structure.
186   */   */
187  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException, BackwardCounterException) {
188    stringstream oss;    stringstream oss;
189    int response = false;    int response = 0;
190    char buff[16];    char buff[16];
191    InputFile->read(buff, sizeof(buff));    InputFile->read(buff, sizeof(buff));
192        
# Line 205  void EventReader::UnpackPscuHeader(void) Line 199  void EventReader::UnpackPscuHeader(void)
199    unsigned char CRC          = buff[15];    unsigned char CRC          = buff[15];
200    unsigned char FileOffset   = buff[15];    unsigned char FileOffset   = buff[15];
201    
202    
203    if (Counter < prevPckCounter){    if (Counter < prevPckCounter){
204        response = prevPckCounter - Counter;
205        //oss.str("");
206        //oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
207        //throw BackwardCounterException(oss.str().c_str());
208        //logger->error(oss.str().c_str());
209      }
210    
211      if (Counter > prevPckCounter + 1){
212      oss.str("");      oss.str("");
213      oss << "Packet counter is less than before of " << (prevPckCounter - Counter);      oss << "Packet counter is greater than before of " << (Counter - prevPckCounter);
214      logger->warn(oss.str().c_str());      logger->error(oss.str().c_str());
215      }
216      
217      if ((OrbitalTime == prevPckOBT) & (PacketId1 == 0x10)){
218        oss.str("");
219        oss << "Onboard Time of this packet is equal to the previous packet OBT";
220        logger->error(oss.str().c_str());
221        logger->error(Header->GetPscuHeader()->Print());
222    }    }
223    
224    if (OrbitalTime < prevPckOBT){    if (OrbitalTime < prevPckOBT){
225      oss.str("");      oss.str("");
226      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
227      logger->warn(oss.str().c_str());      logger->error(oss.str().c_str());
228    }    }
229    
230    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 230  void EventReader::UnpackPscuHeader(void) Line 240  void EventReader::UnpackPscuHeader(void)
240      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
241      Header->GetPscuHeader()->SetCRC(CRC);      Header->GetPscuHeader()->SetCRC(CRC);
242      Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));      Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));
       
     //commented out because of the above test code  
     InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);  
     if (FindStart()) {  
         finalPos =  (long int)InputFile->tellg() - (initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));  
         if(finalPos == 0){  
             logger->debug(_T(" Correct packet length"));  
             }  
         if (finalPos > 0 && finalPos < 64) {  
             oss.str("");  
             oss << " Correct packet length: Padded of " << finalPos << " bytes";  
             logger->debug(oss.str().c_str());  
         }  
         if (finalPos > 64){  
             oss.str("");  
             oss << " The begin of the next packet is far more than 64 byte from the end of the previous."  
                 << Header->GetPscuHeader()->Print();  
             logger->error(oss.str().c_str());  
             //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");  
         }  
     }  
     else {  
         logger->debug(_T(" END OF FILE"));          
     }  
           
         InputFile->seekg(initPos, std::ios::beg);  
243    } else {    } else {
244       /*Here i should extract the block of Data for later analysis */       /*Here i should extract the block of Data for later analysis */
245      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
# Line 264  void EventReader::UnpackPscuHeader(void) Line 248  void EventReader::UnpackPscuHeader(void)
248      throw WrongCRCHeaderException(oss.str().c_str());      throw WrongCRCHeaderException(oss.str().c_str());
249    }    }
250    
251  /*    char tmpId1[4];    if (response > 0){
252      char tmpId2[4];      oss.str("");
253      char tmpLength[12];      oss << "Packet counter is less than before of " << response;
254      char tmpStart[100];      throw BackwardCounterException(oss.str().c_str());
255      char tmpCRC[4];    }
     sprintf(tmpId1, "%02X", PacketId1);  
     sprintf(tmpId2, "%02X", PacketId2);  
     sprintf(tmpLength, "%06X", Header->GetPscuHeader()->GetPacketLenght());  
     sprintf(tmpStart, "%X", Header->GetPscuHeader()->GetFileOffset());  
     sprintf(tmpCRC, "%02X", CRC);  
     oss.flush();  
     oss << "\n Packet Counter (decimal) : "  <<  Counter  
         <<  "\n Id1 - Id2                : "  <<  tmpId1 <<  " - " <<  tmpId2  
         <<  "\n Orbital Time (decimal)   : "  <<  OrbitalTime  
         <<  "\n Lenght                   : "  <<  tmpLength  
         <<  "\n CRC                      : "  <<  tmpCRC  
         <<  "\n Header Start Position    : "  <<  tmpStart;  
     logger->info(oss.str().c_str());  
 */  
256  }  }
257    
258  /**  /**
# Line 308  bool EventReader::FindStart(void) throw Line 278  bool EventReader::FindStart(void) throw
278  //------------------------------------------  //------------------------------------------
279    int subIndex = 0;    int subIndex = 0;
280    char dataByte;    char dataByte;
281    
282    int  buffSize = 100;    int  buffSize = 64;
283    int  index = 0;    int  index = 0;
284    int  loop = -1;    int  loop = -1;
285    char buffer[buffSize];    char buffer[buffSize];
286      bool flagOverPad = false;
287    
288    while (!InputFile->eof()) {    while (!InputFile->eof()) {
289      InputFile->read(buffer, sizeof(buffer));      InputFile->read(buffer, sizeof(buffer));
# Line 323  bool EventReader::FindStart(void) throw Line 294  bool EventReader::FindStart(void) throw
294          if (dataByte == (char)(*(subSign+subIndex))){          if (dataByte == (char)(*(subSign+subIndex))){
295              if (subIndex++ == (subSignDim-1)) {              if (subIndex++ == (subSignDim-1)) {
296                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
297                    if (flagOverPad){
298                        oss.str("");
299                        oss << "\n This packet beginning is farther than 64 byte from the end of the previous."
300                            << "\n Below the is the last already unpacked packet";
301                        logger->error(oss.str().c_str());
302                        logger->error(Header->GetPscuHeader()->Print());
303                    }
304                  return true;                  return true;
305              }              }
306          } else {          } else {
# Line 331  bool EventReader::FindStart(void) throw Line 309  bool EventReader::FindStart(void) throw
309          }          }
310      }      }
311      //Needs to guarantee the overap of the buffer(s) in several loop      //Needs to guarantee the overap of the buffer(s) in several loop
312        flagOverPad = true;
313      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
314    }    }
315    return false;    return false;
Legend:
Removed from v.2.8  
changed lines
  Added in v.5.1
  ViewVC Help
Powered by ViewVC 1.1.23