/[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 5.1 by kusanagi, Fri Sep 16 12:05:49 2005 UTC revision 6.14 by mocchiut, Thu Nov 16 12:28:17 2006 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: mocchiut $
3   * $Date: 2005/08/29 09:46:13 $   * $Date: 2006/11/16 11:01:34 $
4   * $Revision: 5.0 $   * $Revision: 6.13 $
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));
41      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::CalibTof,        new CalibTofReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));
# Line 49  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 56  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 74  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));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbT,            new TsbTReader));
82    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbB,            new TsbBReader));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TsbB,            new TsbBReader));
83  }  }
# 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 5.0 2005/08/29 09:46:13 kusanagi Exp $\n";      "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/EventReader.cpp,v 6.13 2006/11/16 11:01:34 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 122  void EventReader::RunEvent(int EventNumb Line 120  void EventReader::RunEvent(int EventNumb
120    stringstream oss;    stringstream oss;
121    int step = 0;    int step = 0;
122    const PacketType* type;    const PacketType* type;
123    while (!InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){    while ( InputFile->good() && !InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){
124      try {      try {
125        if (FindStart()) {        if (FindStart()) {
126              UnpackPscuHeader();              UnpackPscuHeader();
# Line 139  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 && !InputFile->eof() ){
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 153  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        InputFile->seekg( (-1)*(Header->GetPscuHeader()->GetPacketLenght() + 15) , std::ios::cur);        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 178  void EventReader::RunEvent(int EventNumb Line 201  void EventReader::RunEvent(int EventNumb
201      oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";      oss << "----endPck " << Header->GetPscuHeader()->GetCounter() << "\n";
202      logger->info(oss.str().c_str());      logger->info(oss.str().c_str());
203    }    }
204      if ( !InputFile->eof() && ( InputFile->bad() || InputFile->fail() ) ) printf("\n\n ERROR READING RAW FILE! \n\n"); // E. Mocchiutti
205    
206      Header->GetCounter()->PrintCounters();      Header->GetCounter()->PrintCounters();
207        //if (corruptedPacketFile.is_open()) corruptedPacketFile.close();
208  }  }
209    
210  /**  /**
# Line 188  void EventReader::UnpackPscuHeader(void) Line 214  void EventReader::UnpackPscuHeader(void)
214    stringstream oss;    stringstream oss;
215    int response = 0;    int response = 0;
216    char buff[16];    char buff[16];
217      memset(buff,  0, 16*sizeof(char));
218    InputFile->read(buff, sizeof(buff));    InputFile->read(buff, sizeof(buff));
219        
220        
# Line 197  void EventReader::UnpackPscuHeader(void) Line 224  void EventReader::UnpackPscuHeader(void)
224    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);
225    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);
226    unsigned char CRC          = buff[15];    unsigned char CRC          = buff[15];
227    unsigned char FileOffset   = buff[15];    unsigned char FileOffset   = 0;
228    
229    
230    if (Counter < prevPckCounter){    if (Counter < prevPckCounter){
# Line 239  void EventReader::UnpackPscuHeader(void) Line 266  void EventReader::UnpackPscuHeader(void)
266      //plus the CRC legth (which varies for each type of packet)      //plus the CRC legth (which varies for each type of packet)
267      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);      Header->GetPscuHeader()->SetPacketLenght(PacketLenght);
268      Header->GetPscuHeader()->SetCRC(CRC);      Header->GetPscuHeader()->SetCRC(CRC);
269      Header->GetPscuHeader()->SetFileOffset(((long int)(InputFile->tellg()) - 16));      Header->GetPscuHeader()->SetFileOffset(((unsigned long int)((unsigned long int)(InputFile->tellg()) - 16)));
270        //(unsigned long int)(InputFile->tellg() - 16));
271    } else {    } else {
272       /*Here i should extract the block of Data for later analysis */       /*Here i should extract the block of Data for later analysis */
273      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
# Line 284  bool EventReader::FindStart(void) throw Line 312  bool EventReader::FindStart(void) throw
312    int  loop = -1;    int  loop = -1;
313    char buffer[buffSize];    char buffer[buffSize];
314    bool flagOverPad = false;    bool flagOverPad = false;
315      memset(buffer,  0, buffSize*sizeof(char));
316    while (!InputFile->eof()) {    while (!InputFile->eof()) {
317      InputFile->read(buffer, sizeof(buffer));      InputFile->read(buffer, sizeof(buffer));
318      index = 0;      index = 0;
# Line 315  bool EventReader::FindStart(void) throw Line 343  bool EventReader::FindStart(void) throw
343    return false;    return false;
344  }  }
345    
346    
347    int EventReader::archiveCorruptedPacket(long int offset, long int length) {
348            if (!corruptedPacketFile.is_open()) {
349                    oss.str("");
350                    oss << Run->GetUnpackPath().c_str() << "Corrupted.dat";
351                    cout << oss.str().c_str() << "\n";
352                    cout << dec << offset << "\n";
353                    cout << dec << length << "\n";
354                    corruptedPacketFile.open(oss.str().c_str(), ios::out);
355            }
356            char *buffer = new char[length];
357            InputFile->seekg(offset, ios_base::beg );
358            InputFile->read(buffer, length);
359            corruptedPacketFile.write(buffer, length);
360            InputFile->seekg(offset + 1, ios_base::beg );
361    }
362    
363  ClassImp(EventReader)  ClassImp(EventReader)
Legend:
Removed from v.5.1  
changed lines
  Added in v.6.14
  ViewVC Help
Powered by ViewVC 1.1.23