/[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 6.10 by mocchiut, Mon Nov 13 11:45:49 2006 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: mocchiut $
3   * $Date: 2004/12/22 11:38:35 $   * $Date: 2006/09/29 10:19:22 $
4   * $Revision: 2.7 $   * $Revision: 6.9 $
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 43  EventReader::EventReader(int packetsLimi Line 41  EventReader::EventReader(int packetsLimi
41    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth,    new CalibTrkBothReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth,    new CalibTrkBothReader));
42    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));
43    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));  
44    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));
45    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));
46    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));
# Line 51  EventReader::EventReader(int packetsLimi Line 48  EventReader::EventReader(int packetsLimi
48    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac1,      new Calib2_Ac1Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac1,      new Calib2_Ac1Reader));
49    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac2,      new Calib1_Ac2Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac2,      new Calib1_Ac2Reader));
50    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new Calib2_Ac2Reader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new Calib2_Ac2Reader));
51      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new CalibCalReader));
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 58  EventReader::EventReader(int packetsLimi Line 56  EventReader::EventReader(int packetsLimi
56    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader));
57    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader));
58    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));  
59    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader));
60    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader));
61    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader));
# Line 76  EventReader::EventReader(int packetsLimi Line 72  EventReader::EventReader(int packetsLimi
72    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init,          new S4InitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init,          new S4InitReader));
73    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Init,         new Ac2InitReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Init,         new Ac2InitReader));
74    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader));
75    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac1Alarm,        new Ac1AlarmReader));
76    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader));
77    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader));
78    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader));
79    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Alarm,         new S4AlarmReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Alarm,         new S4AlarmReader));
80      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Alarm,        new Ac2AlarmReader));  
81      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbT,            new TsbTReader));
82      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbB,            new TsbBReader));
83  }  }
84    
85  /**  /**
# Line 88  EventReader::EventReader(int packetsLimi Line 87  EventReader::EventReader(int packetsLimi
87   */   */
88  std::string EventReader::GetVersionInfo(void) const {  std::string EventReader::GetVersionInfo(void) const {
89    return    return
90      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 2.7 2004/12/22 11:38:35 kusanagi Exp $\n";      "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/EventReader.cpp,v 6.9 2006/09/29 10:19:22 mocchiut Exp $\n";
91  }  }
92    
93  /**  /**
# Line 97  std::string EventReader::GetVersionInfo( Line 96  std::string EventReader::GetVersionInfo(
96   */   */
97  void EventReader::Init(PamelaRun *run) {  void EventReader::Init(PamelaRun *run) {
98      SetInputStream(run);      SetInputStream(run);
       
99      //Create the structure of directories and create xxx.Header.root files      //Create the structure of directories and create xxx.Header.root files
100      run->WriteHeaders(this, &Header);      run->WriteHeaders(this, &Header);
101    
# Line 118  static void SkipToNextHeader(ifstream *) Line 116  static void SkipToNextHeader(ifstream *)
116   * Read the next event header, call the reader algorithms that   * Read the next event header, call the reader algorithms that
117   * correspond to its packet type, and read the event trailer.   * correspond to its packet type, and read the event trailer.
118   */   */
   
 // 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------------------  
119  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()  
120    stringstream oss;    stringstream oss;
121    int step = 0;    int step = 0;
122    const PacketType* type;    const PacketType* type;
# Line 147  void EventReader::RunEvent(int EventNumb Line 137  void EventReader::RunEvent(int EventNumb
137                  oss << "\n No way to read events of type  " << type->GetName().c_str() << Header->GetPscuHeader()->Print();                  oss << "\n No way to read events of type  " << type->GetName().c_str() << Header->GetPscuHeader()->Print();
138                  throw NotExistingAlgorithmException(oss.str().c_str()); //to exctract to an higher level and delete the logger!                  throw NotExistingAlgorithmException(oss.str().c_str()); //to exctract to an higher level and delete the logger!
139              }              }
140    
141                //
142                // E. Mocchiutti
143                //
144    
145                //
146                // Start reading the next event jumping the CPU padding to 64 bits
147                //      
148                Float_t pd0 = (16.+Header->GetPscuHeader()->GetPacketLenght())/64.;
149                Float_t pd1 =  pd0 - (Float_t)int(pd0);
150                Float_t padfrac = 64. - pd1 * 64.;
151                //
152                UInt_t padbytes = (UInt_t)padfrac;
153                if ( padbytes > 0 ){
154                  oss.str("");
155                  oss << " padbytes is " << padbytes << " \n skipping these bytes\n";
156                  logger->info(oss.str().c_str());
157                  InputFile->seekg(padbytes,std::ios::cur);
158                };
159    
160                //
161                // end E. Mocchiutti
162                //
163    
164        }        }
165      // 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??
166      } catch (NotExistingAlgorithmException exc) {      } catch (NotExistingAlgorithmException exc) {
# Line 161  void EventReader::RunEvent(int EventNumb Line 175  void EventReader::RunEvent(int EventNumb
175        oss.str("");        oss.str("");
176        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
177        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
178          archiveCorruptedPacket(Header->GetPscuHeader()->FileOffset, Header->GetPscuHeader()->PacketLenght);
179          
180      } catch (UnidentifiedPacketException exc) {      } catch (UnidentifiedPacketException exc) {
181        oss.str("");        oss.str("");
182        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
# Line 173  void EventReader::RunEvent(int EventNumb Line 189  void EventReader::RunEvent(int EventNumb
189        oss.str("");        oss.str("");
190        oss << exc.print() << " " << Header->GetPscuHeader()->Print();        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
191        logger->error(oss.str().c_str());        logger->error(oss.str().c_str());
192        } catch (BackwardCounterException exc) {
193          oss.str("");
194          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
195          logger->error(oss.str().c_str());
196      } catch (...) {      } catch (...) {
197        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");
198      }      }
199      if ((step%1000) == 0) std::cout << step/1000 << "K \n";      if ((step%1000) == 0) std::cout << step/1000 << "K \n";
200        oss.str("");
201        oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";
202        logger->info(oss.str().c_str());
203    }    }
204      Header->GetCounter()->PrintCounters();      Header->GetCounter()->PrintCounters();
205        //if (corruptedPacketFile.is_open()) corruptedPacketFile.close();
206  }  }
207    
 static void SkipToNextHeader(ifstream* TechmodelFile) {  
   //  yeah.  
 }  
   
   
   
208  /**  /**
209   * Unpack the PSCU header from a file into the structure.   * Unpack the PSCU header from a file into the structure.
210   */   */
211  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException, BackwardCounterException) {
212    stringstream oss;    stringstream oss;
213    int response = false;    int response = 0;
214    char buff[16];    char buff[16];
215      memset(buff,  0, 16*sizeof(char));
216    InputFile->read(buff, sizeof(buff));    InputFile->read(buff, sizeof(buff));
217        
218        
# Line 205  void EventReader::UnpackPscuHeader(void) Line 224  void EventReader::UnpackPscuHeader(void)
224    unsigned char CRC          = buff[15];    unsigned char CRC          = buff[15];
225    unsigned char FileOffset   = buff[15];    unsigned char FileOffset   = buff[15];
226    
227    
228    if (Counter < prevPckCounter){    if (Counter < prevPckCounter){
229        response = prevPckCounter - Counter;
230        //oss.str("");
231        //oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
232        //throw BackwardCounterException(oss.str().c_str());
233        //logger->error(oss.str().c_str());
234      }
235    
236      if (Counter > prevPckCounter + 1){
237      oss.str("");      oss.str("");
238      oss << "Packet counter is less than before of " << (prevPckCounter - Counter);      oss << "Packet counter is greater than before of " << (Counter - prevPckCounter);
239      logger->warn(oss.str().c_str());      logger->error(oss.str().c_str());
240      }
241      
242      if ((OrbitalTime == prevPckOBT) & (PacketId1 == 0x10)){
243        oss.str("");
244        oss << "Onboard Time of this packet is equal to the previous packet OBT";
245        logger->error(oss.str().c_str());
246        logger->error(Header->GetPscuHeader()->Print());
247    }    }
248    
249    if (OrbitalTime < prevPckOBT){    if (OrbitalTime < prevPckOBT){
250      oss.str("");      oss.str("");
251      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
252      logger->warn(oss.str().c_str());      logger->error(oss.str().c_str());
253    }    }
254    
255    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 265  void EventReader::UnpackPscuHeader(void)
265      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
266      Header->GetPscuHeader()->SetCRC(CRC);      Header->GetPscuHeader()->SetCRC(CRC);
267      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);  
268    } else {    } else {
269       /*Here i should extract the block of Data for later analysis */       /*Here i should extract the block of Data for later analysis */
270      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
# Line 264  void EventReader::UnpackPscuHeader(void) Line 273  void EventReader::UnpackPscuHeader(void)
273      throw WrongCRCHeaderException(oss.str().c_str());      throw WrongCRCHeaderException(oss.str().c_str());
274    }    }
275    
276  /*    char tmpId1[4];    if (response > 0){
277      char tmpId2[4];      oss.str("");
278      char tmpLength[12];      oss << "Packet counter is less than before of " << response;
279      char tmpStart[100];      throw BackwardCounterException(oss.str().c_str());
280      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());  
 */  
281  }  }
282    
283  /**  /**
# Line 308  bool EventReader::FindStart(void) throw Line 303  bool EventReader::FindStart(void) throw
303  //------------------------------------------  //------------------------------------------
304    int subIndex = 0;    int subIndex = 0;
305    char dataByte;    char dataByte;
306    
307    int  buffSize = 100;    int  buffSize = 64;
308    int  index = 0;    int  index = 0;
309    int  loop = -1;    int  loop = -1;
310    char buffer[buffSize];    char buffer[buffSize];
311      bool flagOverPad = false;
312      memset(buffer,  0, buffSize*sizeof(char));
313    while (!InputFile->eof()) {    while (!InputFile->eof()) {
314      InputFile->read(buffer, sizeof(buffer));      InputFile->read(buffer, sizeof(buffer));
315      index = 0;      index = 0;
# Line 323  bool EventReader::FindStart(void) throw Line 319  bool EventReader::FindStart(void) throw
319          if (dataByte == (char)(*(subSign+subIndex))){          if (dataByte == (char)(*(subSign+subIndex))){
320              if (subIndex++ == (subSignDim-1)) {              if (subIndex++ == (subSignDim-1)) {
321                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
322                    if (flagOverPad){
323                        oss.str("");
324                        oss << "\n This packet beginning is farther than 64 byte from the end of the previous."
325                            << "\n Below the is the last already unpacked packet";
326                        logger->error(oss.str().c_str());
327                        logger->error(Header->GetPscuHeader()->Print());
328                    }
329                  return true;                  return true;
330              }              }
331          } else {          } else {
# Line 331  bool EventReader::FindStart(void) throw Line 334  bool EventReader::FindStart(void) throw
334          }          }
335      }      }
336      //Needs to guarantee the overap of the buffer(s) in several loop      //Needs to guarantee the overap of the buffer(s) in several loop
337        flagOverPad = true;
338      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
339    }    }
340    return false;    return false;
341  }  }
342    
343    
344    int EventReader::archiveCorruptedPacket(long int offset, long int length) {
345            if (!corruptedPacketFile.is_open()) {
346                    oss.str("");
347                    oss << Run->GetUnpackPath().c_str() << "Corrupted.dat";
348                    cout << oss.str().c_str() << "\n";
349                    cout << dec << offset << "\n";
350                    cout << dec << length << "\n";
351                    corruptedPacketFile.open(oss.str().c_str(), ios::out);
352            }
353            char *buffer = new char[length];
354            InputFile->seekg(offset, ios_base::beg );
355            InputFile->read(buffer, length);
356            corruptedPacketFile.write(buffer, length);
357            InputFile->seekg(offset + 1, ios_base::beg );
358    }
359    
360  ClassImp(EventReader)  ClassImp(EventReader)
Legend:
Removed from v.2.8  
changed lines
  Added in v.6.10
  ViewVC Help
Powered by ViewVC 1.1.23