/[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.0 by kusanagi, Tue Sep 21 20:50:54 2004 UTC revision 6.15 by mocchiut, Mon Nov 20 08:58:04 2006 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: mocchiut $
3   * $Date: 2004/09/21 20:24:33 $   * $Date: 2006/11/16 12:28:17 $
4   * $Revision: 1.10 $   * $Revision: 6.14 $
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>  
 #include "TechmodelAlgorithm.h"  
11  #include "EventReader.h"  #include "EventReader.h"
12  #include "ReaderAlgorithms.h"  #include "ReaderAlgorithms.h"
13  #include "Exception.h"  
 #include "stdio.h"  
14  extern "C" {  extern "C" {
15  #include "CRC.h"      #include "CRC.h"
16  }  }
17    
 #define BYTE unsigned char  
   
18  using namespace pamela;  using namespace pamela;
19  using namespace pamela::techmodel;  using namespace pamela::techmodel;
20    
# Line 35  unsigned int EventReader::prevPckOBT     Line 30  unsigned int EventReader::prevPckOBT    
30   */   */
31  EventReader::EventReader(int packetsLimit = -1):  EventReader::EventReader(int packetsLimit = -1):
32    TechmodelAlgorithm(0, "TechmodelEventReader"){    TechmodelAlgorithm(0, "TechmodelEventReader"){
 //EventReader::EventReader(void):  
 //  TechmodelAlgorithm(0, "TechmodelEventReader") {  
33    EventReader::maxPackets = packetsLimit;    EventReader::maxPackets = packetsLimit;
34    logger->debug(_T("Constructor"));    logger->debug(_T("Constructor"));
35    Header = new EventHeader();    Header = new EventHeader();
36    
37    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun,      new PhysEndRunReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun,      new PhysEndRunReader));
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));
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));
44    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd,        new CalibTrdReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));
45    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));
46    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));
47    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac1,      new Calib1_Ac1Reader));
48    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibAc,         new CalibAcReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac1,      new Calib2_Ac1Reader));
49    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac2,      new Calib1_Ac2Reader));
50    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new Calib2_Ac2Reader));
51    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new CalibCalReader));
52    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer,    new CalibTrailerReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader));
53    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader));
54    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader));
55    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer,    new CalibTrailerReader));
56    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk,         new TestTrkReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader));
57    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof,         new TestTofReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader));
58    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader));
59    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader));
60    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader));
61    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump,         new TabDumpReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader));
62    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump,         new TabDumpReader));
63    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader));
64    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader));
65    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcInit,          new AcInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader));
66    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac1Init,         new Ac1InitReader));
67    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader));
68    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader));
69    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader));
70    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader));
71    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::NdInit,          new NdInitReader));
72    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init,          new S4InitReader));
73    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Init,         new Ac2InitReader));
74    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader));
75      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac1Alarm,        new Ac1AlarmReader));
76      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader));
77      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader));
78      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader));
79      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 86  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 1.10 2004/09/21 20:24:33 kusanagi Exp $\n";      "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/EventReader.cpp,v 6.14 2006/11/16 12:28:17 mocchiut Exp $\n";
91  }  }
92    
93  /**  /**
# Line 94  std::string EventReader::GetVersionInfo( Line 95  std::string EventReader::GetVersionInfo(
95   * event reader routines for all packet types.   * event reader routines for all packet types.
96   */   */
97  void EventReader::Init(PamelaRun *run) {  void EventReader::Init(PamelaRun *run) {
     std::stringstream oss;  
     logger->debug(_T("Initialize"));  
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        fPL = 0;
102        notfirst = false;
103        jgood = true;
104      //Create the xxx.root in it's specific directory      //Create the xxx.root in it's specific directory
105      for (AlgorithmMap::iterator i = TechmodelAlgorithmMap.begin();      for (AlgorithmMap::iterator i = TechmodelAlgorithmMap.begin();
106         i != TechmodelAlgorithmMap.end(); i++) {         i != TechmodelAlgorithmMap.end(); i++) {
107         oss.flush();         oss.str("");
108         oss << "Initializing algo " << i->second->GetAlgorithmName();         oss << "Initializing algo " << i->second->GetAlgorithmName();
        //cout << "Initializing algo " << i->second->GetAlgorithmName();  
109         logger->debug(oss.str().c_str());         logger->debug(oss.str().c_str());
110         i->second->Init(run);         i->second->Init(run);
111      }      }
# Line 119  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    while (!InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){    const PacketType* type;
125      while ( InputFile->good() && !InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){
126      try {      try {
127        if (FindStart()) {        if (FindStart()) {
128          if(UnpackPscuHeader()){              UnpackPscuHeader();
129              const PacketType* type = Header->GetPscuHeader()->GetPacketType();              type = Header->GetPscuHeader()->GetPacketType();
130              AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type);              AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type);
131              if (i != TechmodelAlgorithmMap.end()) {              if (i != TechmodelAlgorithmMap.end()) {
132                  TechmodelAlgorithm *EventAlgorithm = i->second;                  TechmodelAlgorithm *EventAlgorithm(i->second);
133                  EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());                  EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());
134                  Run->FillTrees(type);                  Run->FillTrees(type);
135                  Header->GetCounter()->Increment(type);                  Header->GetCounter()->Increment(type);
136                    logger->info(Header->GetPscuHeader()->Print());
137              } else {              } else {
138                  oss.flush();                  oss.str("");
139                  oss << "No way to read events of type" << type->GetName();                  oss << "\n No way to read events of type  " << type->GetName().c_str() << Header->GetPscuHeader()->Print();
140                  logger->info(oss.str().c_str());                  throw NotExistingAlgorithmException(oss.str().c_str()); //to exctract to an higher level and delete the logger!
141              throw Exception("No way to read events of type " + type->GetName());              }
142          }  
143        } //else {              //
144          //**TO BE DONE**   <<WRONG CRC l --- LOST PACKET>>//              // E. Mocchiutti
145          //Have to be defined a "WrongCRC" packet::Type?              //
146          //}              //
147       }              // Start reading the next event jumping the CPU padding to 64 bits
148                //
149                if ( Header->GetPscuHeader()->GetPacketId1() == 0x10 ) notfirst = true;        
150                //
151                if ( !InputFile->eof() && notfirst ){
152                        const unsigned char fafe[3]={0xFA, 0xFE, 0xDE};
153                        char isfa[3];
154                        memset(isfa, 0, 3*sizeof(char));
155                        InputFile->read(isfa, sizeof(isfa));
156                        InputFile->seekg(-3,std::ios::cur);
157                        //
158                        fPL += 16 + Header->GetPscuHeader()->GetPacketLenght();
159                        //
160                        if ( (0x000000FF & isfa[0]) == fafe[0] && (0x000000FF & isfa[2]) == fafe[2] && (0x000000FF & isfa[2]) == fafe[2] ){
161                               //
162                               fPL += 0;
163                        } else {
164                                if ( !jgood ){
165                                        jgood = true;
166                                } else {
167                                        Float_t pd0 = fPL/64.;
168                                        fPL = 0;
169                                        Float_t pd1 = pd0 - (Float_t)int(pd0);
170                                        Float_t padfrac = 64. - pd1 * 64.;
171                                        //
172                                        UInt_t padbytes = (UInt_t)padfrac;
173                                        if ( padbytes > 0 && padbytes < 64 ){
174                                                  oss.str("");
175                                                  oss << " padbytes is " << padbytes << " \n skipping these bytes\n";
176                                                  logger->info(oss.str().c_str());
177                                                  InputFile->seekg(padbytes,std::ios::cur);
178                                        };
179                                };
180                        };
181                };
182                //
183                // end E. Mocchiutti
184                //
185    
186          }
187        // In case of exception have to save the packet in a specific root file??
188        } catch (NotExistingAlgorithmException exc) {
189          oss.str("");
190          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
191          logger->error(oss.str().c_str());
192          jgood = false;
193        } catch (WrongCRCHeaderException exc) {
194          oss.str("");
195          oss << exc.print();
196          logger->error(oss.str().c_str());
197          jgood = false;
198        } catch (WrongCRCException exc) {
199          oss.str("");
200          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
201          logger->error(oss.str().c_str());
202          archiveCorruptedPacket(Header->GetPscuHeader()->FileOffset, Header->GetPscuHeader()->PacketLenght);
203          jgood = false;      
204        } catch (UnidentifiedPacketException exc) {
205          oss.str("");
206          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
207          logger->error(oss.str().c_str());
208          jgood = false;
209        } catch (NotExistingCounterException exc) {
210          oss.str("");
211          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
212          logger->error(oss.str().c_str());
213          jgood = false;
214        } catch (LengthException exc) {
215          oss.str("");
216          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
217          logger->error(oss.str().c_str());
218          jgood = false;
219        } catch (BackwardCounterException exc) {
220          oss.str("");
221          oss << exc.print() << " " << Header->GetPscuHeader()->Print();
222          logger->error(oss.str().c_str());
223          jgood = false;
224      } catch (...) {      } catch (...) {
225        //This have to be more detailed. More exceptions type are needed.        logger->error("Couldn't read the event. Skipping to the next header. \n");
226        logger->error("Couldn't read the event. Skipping to the next header.");        jgood = false;
227      }      }
228      if ((step%1000) == 0) std::cout << step << "K \n";      if ((step%1000) == 0) std::cout << step/1000 << "K \n";
229        oss.str("");
230        oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";
231        logger->info(oss.str().c_str());
232    }    }
233      Header->GetCounter()->PrintCounters();    if ( !InputFile->eof() && ( InputFile->bad() || InputFile->fail() ) ) printf("\n\n ERROR READING RAW FILE! \n\n"); // E. Mocchiutti
 }  
234    
235  static void SkipToNextHeader(ifstream* TechmodelFile) {      Header->GetCounter()->PrintCounters();
236    //  yeah.      //if (corruptedPacketFile.is_open()) corruptedPacketFile.close();
237  }  }
238    
   
   
239  /**  /**
240   * Unpack the PSCU header from a file into the structure.   * Unpack the PSCU header from a file into the structure.
241   */   */
242  int EventReader::UnpackPscuHeader(void) throw (std::exception) {  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException, BackwardCounterException) {
243    stringstream oss;    stringstream oss;
244    int response;    int response = 0;
245    char buff[16];    char buff[16];
246      memset(buff,  0, 16*sizeof(char));
247    InputFile->read(buff, sizeof(buff));    InputFile->read(buff, sizeof(buff));
248        
249        
# Line 183  int EventReader::UnpackPscuHeader(void) Line 253  int EventReader::UnpackPscuHeader(void)
253    unsigned int  OrbitalTime  = (((UINT32)buff[8]<<24)&0xFF000000) + (((UINT32)buff[9]<<16)&0x00FF0000) +  (((UINT32)buff[10]<<8)&0x0000FF00) + (((UINT32)buff[11])&0x000000FF);    unsigned int  OrbitalTime  = (((UINT32)buff[8]<<24)&0xFF000000) + (((UINT32)buff[9]<<16)&0x00FF0000) +  (((UINT32)buff[10]<<8)&0x0000FF00) + (((UINT32)buff[11])&0x000000FF);
254    unsigned int  PacketLenght = (((UINT32)buff[12]<<16)&0x00FF0000) +  (((UINT32)buff[13]<<8)&0x0000FF00) + (((UINT32)buff[14])&0x000000FF);    unsigned int  PacketLenght = (((UINT32)buff[12]<<16)&0x00FF0000) +  (((UINT32)buff[13]<<8)&0x0000FF00) + (((UINT32)buff[14])&0x000000FF);
255    unsigned char CRC          = buff[15];    unsigned char CRC          = buff[15];
256    unsigned char FileOffset   = buff[15];    unsigned char FileOffset   = 0;
257    
258    
259    if (Counter < prevPckCounter){    if (Counter < prevPckCounter){
260      oss.flush();      response = prevPckCounter - Counter;
261      oss << "Packet counter is less than before of " << (prevPckCounter - Counter);      //oss.str("");
262      logger->info(oss.str().c_str());      //oss << "Packet counter is less than before of " << (prevPckCounter - Counter);
263        //throw BackwardCounterException(oss.str().c_str());
264        //logger->error(oss.str().c_str());
265      }
266    
267      if (Counter > prevPckCounter + 1){
268        oss.str("");
269        oss << "Packet counter is greater than before of " << (Counter - prevPckCounter);
270        logger->error(oss.str().c_str());
271      }
272      
273      if ((OrbitalTime == prevPckOBT) & (PacketId1 == 0x10)){
274        oss.str("");
275        oss << "Onboard Time of this packet is equal to the previous packet OBT";
276        logger->error(oss.str().c_str());
277        logger->error(Header->GetPscuHeader()->Print());
278    }    }
279    
280    if (OrbitalTime < prevPckOBT){    if (OrbitalTime < prevPckOBT){
281      oss.flush();      oss.str("");
282      oss << " Orbital Time is less than before of " << (prevPckOBT - OrbitalTime);      oss << " Onboard Time is less than before of " << (prevPckOBT - OrbitalTime);
283      logger->info(oss.str().c_str());      notfirst = false;
284        logger->error(oss.str().c_str());
285    }    }
286    
   
287    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)){
288      prevPckCounter = Counter;      prevPckCounter = Counter;
289      prevPckOBT     = OrbitalTime;      prevPckOBT     = OrbitalTime;
# Line 210  int EventReader::UnpackPscuHeader(void) Line 296  int EventReader::UnpackPscuHeader(void)
296      //plus the CRC legth (which varies for each type of packet)      //plus the CRC legth (which varies for each type of packet)
297      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
298      Header->GetPscuHeader()->SetCRC(CRC);      Header->GetPscuHeader()->SetCRC(CRC);
299      Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));      Header->GetPscuHeader()->SetFileOffset(((unsigned long int)((unsigned long int)(InputFile->tellg()) - 16)));
300            //(unsigned long int)(InputFile->tellg() - 16));
     //commented out because of the above test code  
     InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);  
     FindStart();  
     finalPos =  (long int)InputFile->tellg() - (1 + initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));  
         if(finalPos == 0){  
             logger->debug(_T("Correct packet length"));  
             }  
         if (finalPos > 0 && finalPos < 64) {  
             oss.flush();  
             oss << " Correct packet length: Padded of " << finalPos << " bytes";  
             logger->warn(oss.str().c_str());  
         }  
         if (finalPos > 64){  
             logger->warn(_T(" Wrong packet length? (because of wrong padding?)"));  
         }  
         InputFile->seekg(initPos, std::ios::beg);  
         response =  true;  
301    } else {    } else {
302      logger->warn(_T(" WRONG CRC FOR HEADER "));       /*Here i should extract the block of Data for later analysis */
303      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
304      response =  false;      oss.str("");
305        oss << "CRC Header Error on packet:" << PscuHeader::Print(buff);
306        throw WrongCRCHeaderException(oss.str().c_str());
307      }
308    
309      if (response > 0){
310        oss.str("");
311        oss << "Packet counter is less than before of " << response;
312        throw BackwardCounterException(oss.str().c_str());
313    }    }
     char tmpId1[4];  
     char tmpId2[4];  
     char tmpLength[12];  
     char tmpStart[100];  
     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());  
   return response;  
314  }  }
315    
316  /**  /**
# Line 265  void EventReader::UnpackPscuTrailer(void Line 323  void EventReader::UnpackPscuTrailer(void
323  /**  /**
324   * 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
325   */   */
326  int EventReader::FindStart(void) throw (std::exception) {  bool EventReader::FindStart(void) throw (std::exception) {
327  //search an hexadecimal sequence in a file  //search an hexadecimal sequence in a file
328  //subSign    ------> pointer to the sequence buffer  //subSign    ------> pointer to the sequence buffer
329  //subSignDim ------> dimension of the buffer  //subSignDim ------> dimension of the buffer
# Line 273  int EventReader::FindStart(void) throw ( Line 331  int EventReader::FindStart(void) throw (
331  // return true if founds a match (else false)  // return true if founds a match (else false)
332  // subF point rigth after the match, if found. Else EOF.  // subF point rigth after the match, if found. Else EOF.
333  //Maurizio 15/11/2002-----------------------  //Maurizio 15/11/2002-----------------------
334    const unsigned char subSign[3]={0xFA, 0xFE, 0xDE};    const int subSignDim = 3;
335    int subSignDim = 3;    const unsigned char subSign[subSignDim]={0xFA, 0xFE, 0xDE};
336  //------------------------------------------  //------------------------------------------
337    int subIndex = 0;    int subIndex = 0;
338    char dataByte;    char dataByte;
339    
340    int  buffSize = 100;    int  buffSize = 64;
341    int  index = 0;    int  index = 0;
342    int  loop = 0;    int  loop = -1;
343    char buffer[buffSize];    char buffer[buffSize];
344      bool flagOverPad = false;
345    /* original version    memset(buffer,  0, buffSize*sizeof(char));
   while (!InputFile->eof()) {  
     InputFile->get(dataByte);  
     if (dataByte == (char)(*(subSign+subIndex))){  
       if (subIndex++ == (subSignDim-1)) {  
          InputFile->seekg(-(subIndex), std::ios::cur);  
         return 1;  
        }  
     } else {  
       if (InputFile->eof()) {  
         throw Exception("Extra bytes over the last packets");  
       } else {  
         subIndex = 0;  
       }  
     }  
   } */  
346    while (!InputFile->eof()) {    while (!InputFile->eof()) {
347      InputFile->read(buffer, sizeof(buffer));      InputFile->read(buffer, sizeof(buffer));
348      index = 0;      index = 0;
349      ++loop;      loop++;
350      while (index < buffSize){      while (index < buffSize){
351          dataByte = buffer[index++];          dataByte = buffer[index++];
352          if (dataByte == (char)(*(subSign+subIndex))){          if (dataByte == (char)(*(subSign+subIndex))){
353              if (subIndex++ == (subSignDim-1)) {              if (subIndex++ == (subSignDim-1)) {
354                  InputFile->seekg( (index - subIndex - (buffSize*loop)), std::ios::cur);                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
355              return 1;                  if (flagOverPad){
356              }                      oss.str("");
357                        oss << "\n This packet beginning is farther than 64 byte from the end of the previous."
358                            << "\n Below the is the last already unpacked packet";
359                        logger->error(oss.str().c_str());
360                        logger->error(Header->GetPscuHeader()->Print());
361                    }
362                    return true;
363                }
364          } else {          } else {
365                index = index - (subIndex);
366              subIndex = 0;              subIndex = 0;
367          }          }
368      }      }
369    }      //Needs to guarantee the overap of the buffer(s) in several loop
370    return 0;      flagOverPad = true;
371        InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
372      }
373      return false;
374  }  }
375    
376    
377    int EventReader::archiveCorruptedPacket(long int offset, long int length) {
378            if (!corruptedPacketFile.is_open()) {
379                    oss.str("");
380                    oss << Run->GetUnpackPath().c_str() << "Corrupted.dat";
381                    cout << oss.str().c_str() << "\n";
382                    cout << dec << offset << "\n";
383                    cout << dec << length << "\n";
384                    corruptedPacketFile.open(oss.str().c_str(), ios::out);
385            }
386            char *buffer = new char[length];
387            InputFile->seekg(offset, ios_base::beg );
388            InputFile->read(buffer, length);
389            corruptedPacketFile.write(buffer, length);
390            InputFile->seekg(offset + 1, ios_base::beg );
391    }
392    
393  ClassImp(EventReader)  ClassImp(EventReader)
Legend:
Removed from v.2.0  
changed lines
  Added in v.6.15
  ViewVC Help
Powered by ViewVC 1.1.23