/[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 4.0 by kusanagi, Sun Mar 6 04:33:02 2005 UTC revision 6.3 by kusanagi, Tue Apr 25 16:51:03 2006 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: kusanagi $
3   * $Date: 2005/03/06 04:29:20 $   * $Date: 2006/04/05 07:30:54 $
4   * $Revision: 3.1 $   * $Revision: 6.2 $
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.
8   */   */
9    
10  #include <log4cxx/logger.h>  #include <log4cxx/logger.h>
 #include <fstream>  
11  #include "EventReader.h"  #include "EventReader.h"
12  #include "ReaderAlgorithms.h"  #include "ReaderAlgorithms.h"
13    
   
14  extern "C" {  extern "C" {
15      #include "CRC.h"      #include "CRC.h"
16  }  }
# Line 40  EventReader::EventReader(int packetsLimi Line 38  EventReader::EventReader(int packetsLimi
38    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1,  new CalibCalPulse1Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse1,  new CalibCalPulse1Reader));
39    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2,  new CalibCalPulse2Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPulse2,  new CalibCalPulse2Reader));
40    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));  
41    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));
42    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));  
43    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));
44    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));
45    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));
# Line 58  EventReader::EventReader(int packetsLimi Line 54  EventReader::EventReader(int packetsLimi
54    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader));
55    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader));
56    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader));
   TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk,         new TestTrkReader));  
   TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof,         new TestTofReader));  
57    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader));
58    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader));
59    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader));
# Line 90  EventReader::EventReader(int packetsLimi Line 84  EventReader::EventReader(int packetsLimi
84   */   */
85  std::string EventReader::GetVersionInfo(void) const {  std::string EventReader::GetVersionInfo(void) const {
86    return    return
87      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 3.1 2005/03/06 04:29:20 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 6.2 2006/04/05 07:30:54 kusanagi Exp $\n";
88  }  }
89    
90  /**  /**
# Line 99  std::string EventReader::GetVersionInfo( Line 93  std::string EventReader::GetVersionInfo(
93   */   */
94  void EventReader::Init(PamelaRun *run) {  void EventReader::Init(PamelaRun *run) {
95      SetInputStream(run);      SetInputStream(run);
       
96      //Create the structure of directories and create xxx.Header.root files      //Create the structure of directories and create xxx.Header.root files
97      run->WriteHeaders(this, &Header);      run->WriteHeaders(this, &Header);
98    
# Line 120  static void SkipToNextHeader(ifstream *) Line 113  static void SkipToNextHeader(ifstream *)
113   * Read the next event header, call the reader algorithms that   * Read the next event header, call the reader algorithms that
114   * correspond to its packet type, and read the event trailer.   * correspond to its packet type, and read the event trailer.
115   */   */
   
 // 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------------------  
116  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()  
117    stringstream oss;    stringstream oss;
118    int step = 0;    int step = 0;
119    const PacketType* type;    const PacketType* type;
# Line 163  void EventReader::RunEvent(int EventNumb Line 148  void EventReader::RunEvent(int EventNumb
148        oss.str("");        oss.str("");
149        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
150        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
151          //archiveCorruptedPacket(Header->GetPscuHeader()->FileOffset, Header->GetPscuHeader()->PacketLenght);
152          //InputFile->seekg( (-1)*(Header->GetPscuHeader()->GetPacketLenght() + 14) , std::ios::cur);
153          
154      } catch (UnidentifiedPacketException exc) {      } catch (UnidentifiedPacketException exc) {
155        oss.str("");        oss.str("");
156        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
# Line 175  void EventReader::RunEvent(int EventNumb Line 163  void EventReader::RunEvent(int EventNumb
163        oss.str("");        oss.str("");
164        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
165        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
166        } catch (BackwardCounterException exc) {
167          oss.str("");
168          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
169          logger->error(oss.str().c_str());
170      } catch (...) {      } catch (...) {
171        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");
172      }      }
# Line 184  void EventReader::RunEvent(int EventNumb Line 176  void EventReader::RunEvent(int EventNumb
176      logger->info(oss.str().c_str());      logger->info(oss.str().c_str());
177    }    }
178      Header->GetCounter()->PrintCounters();      Header->GetCounter()->PrintCounters();
179        //if (corruptedPacketFile.is_open()) corruptedPacketFile.close();
180  }  }
181    
 static void SkipToNextHeader(ifstream* TechmodelFile) {  
   //  yeah.  
 }  
   
   
   
182  /**  /**
183   * Unpack the PSCU header from a file into the structure.   * Unpack the PSCU header from a file into the structure.
184   */   */
185  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {
186    stringstream oss;    stringstream oss;
187    int response = false;    int response = 0;
188    char buff[16];    char buff[16];
189    InputFile->read(buff, sizeof(buff));    InputFile->read(buff, sizeof(buff));
190        
# Line 210  void EventReader::UnpackPscuHeader(void) Line 197  void EventReader::UnpackPscuHeader(void)
197    unsigned char CRC          = buff[15];    unsigned char CRC          = buff[15];
198    unsigned char FileOffset   = buff[15];    unsigned char FileOffset   = buff[15];
199    
200    
201    if (Counter < prevPckCounter){    if (Counter < prevPckCounter){
202      oss.str("");      response = prevPckCounter - Counter;
203      oss << "Packet counter is less than before of " << (prevPckCounter - Counter);      //oss.str("");
204      logger->error(oss.str().c_str());      //oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
205        //throw BackwardCounterException(oss.str().c_str());
206        //logger->error(oss.str().c_str());
207    }    }
208    
209    if (Counter > prevPckCounter + 1){    if (Counter > prevPckCounter + 1){
# Line 248  void EventReader::UnpackPscuHeader(void) Line 238  void EventReader::UnpackPscuHeader(void)
238      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
239      Header->GetPscuHeader()->SetCRC(CRC);      Header->GetPscuHeader()->SetCRC(CRC);
240      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 << "\n The begin of the next packet is far more than 64 byte from the end of the previous."  
                 << "\n Below the is the now unpacking packet";  
             logger->error(oss.str().c_str());  
             logger->error(Header->GetPscuHeader()->Print());  
             //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);  
241    } else {    } else {
242       /*Here i should extract the block of Data for later analysis */       /*Here i should extract the block of Data for later analysis */
243      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
# Line 283  void EventReader::UnpackPscuHeader(void) Line 246  void EventReader::UnpackPscuHeader(void)
246      throw WrongCRCHeaderException(oss.str().c_str());      throw WrongCRCHeaderException(oss.str().c_str());
247    }    }
248    
249  /*    char tmpId1[4];    if (response > 0){
250      char tmpId2[4];      oss.str("");
251      char tmpLength[12];      oss << "Packet counter is less than before of " << response;
252      char tmpStart[100];      throw BackwardCounterException(oss.str().c_str());
253      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());  
 */  
254  }  }
255    
256  /**  /**
# Line 327  bool EventReader::FindStart(void) throw Line 276  bool EventReader::FindStart(void) throw
276  //------------------------------------------  //------------------------------------------
277    int subIndex = 0;    int subIndex = 0;
278    char dataByte;    char dataByte;
279    
280    int  buffSize = 100;    int  buffSize = 64;
281    int  index = 0;    int  index = 0;
282    int  loop = -1;    int  loop = -1;
283    char buffer[buffSize];    char buffer[buffSize];
284      bool flagOverPad = false;
285    
286    while (!InputFile->eof()) {    while (!InputFile->eof()) {
287      InputFile->read(buffer, sizeof(buffer));      InputFile->read(buffer, sizeof(buffer));
# Line 342  bool EventReader::FindStart(void) throw Line 292  bool EventReader::FindStart(void) throw
292          if (dataByte == (char)(*(subSign+subIndex))){          if (dataByte == (char)(*(subSign+subIndex))){
293              if (subIndex++ == (subSignDim-1)) {              if (subIndex++ == (subSignDim-1)) {
294                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
295                    if (flagOverPad){
296                        oss.str("");
297                        oss << "\n This packet beginning is farther than 64 byte from the end of the previous."
298                            << "\n Below the is the last already unpacked packet";
299                        logger->error(oss.str().c_str());
300                        logger->error(Header->GetPscuHeader()->Print());
301                    }
302                  return true;                  return true;
303              }              }
304          } else {          } else {
# Line 350  bool EventReader::FindStart(void) throw Line 307  bool EventReader::FindStart(void) throw
307          }          }
308      }      }
309      //Needs to guarantee the overap of the buffer(s) in several loop      //Needs to guarantee the overap of the buffer(s) in several loop
310        flagOverPad = true;
311      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
312    }    }
313    return false;    return false;
314  }  }
315    
316    /*
317    int EventReader::archiveCorruptedPacket(long int offset, long int length) {
318            if (!corruptedPacketFile.is_open()) {
319                    oss.str("");
320                    oss << Run->GetUnpackPath().c_str() << "Corrupted.dat";
321                    cout << oss.str().c_str() << "\n";
322                    cout << dec << offset << "\n";
323                    cout << dec << length << "\n";
324                    corruptedPacketFile.open(oss.str().c_str(), ios::out);
325            }
326            char *buffer = new char[length];
327            InputFile->seekg(offset, ios_base::beg );
328            InputFile->read(buffer, length);
329            corruptedPacketFile.write(buffer, length);
330            InputFile->seekg(offset + 1, ios_base::beg );
331    }
332    */
333  ClassImp(EventReader)  ClassImp(EventReader)
Legend:
Removed from v.4.0  
changed lines
  Added in v.6.3
  ViewVC Help
Powered by ViewVC 1.1.23