/[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.4 by kusanagi, Fri Dec 3 22:08:00 2004 UTC revision 2.9 by kusanagi, Thu Jan 13 08:20:26 2005 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Author: kusanagi $   * $Author: kusanagi $
3   * $Date: 2004/10/17 12:28:38 $   * $Date: 2005/01/11 14:59:31 $
4   * $Revision: 2.3 $   * $Revision: 2.8 $
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 36  EventReader::EventReader(int packetsLimi Line 36  EventReader::EventReader(int packetsLimi
36    logger->debug(_T("Constructor"));    logger->debug(_T("Constructor"));
37    Header = new EventHeader();    Header = new EventHeader();
38    
39    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun,      new PhysEndRunReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::PhysEndRun,      new PhysEndRunReader));
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));
43    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth,    new CalibTrkBothReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrkBoth,    new CalibTrkBothReader));
44    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk1,       new CalibTrk1Reader));
45    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2,       new CalibTrk2Reader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrk2,       new CalibTrk2Reader));
46    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd,        new CalibTrdReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrd,        new CalibTrdReader));
47    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTof,        new CalibTofReader));
48    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibS4,         new CalibS4Reader));
49    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibCalPed,     new CalibCalPedReader));
50    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibAc,         new CalibAcReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac1,      new Calib1_Ac1Reader));
51    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac1,      new Calib2_Ac1Reader));
52    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib1_Ac2,      new Calib1_Ac2Reader));
53    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Calib2_Ac2,      new Calib2_Ac2Reader));
54    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer,    new CalibTrailerReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunHeader,       new RunHeaderReader));
55    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::RunTrailer,      new RunTrailerReader));
56    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibHeader,     new CalibHeaderReader));
57    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalibTrailer,    new CalibTrailerReader));
58    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk,         new TestTrkReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitHeader,      new InitHeaderReader));
59    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof,         new TestTofReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::InitTrailer,     new InitTrailerReader));
60    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::EventTrk,        new EventTrkReader));
61    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTrk,         new TestTrkReader));
62    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TestTof,         new TestTofReader));
63    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump,         new TabDumpReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Log,             new LogReader));
64    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::VarDump,         new VarDumpReader));
65    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ArrDump,         new ArrDumpReader));
66    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TabDump,         new TabDumpReader));
67    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcInit,          new AcInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Tmtc,            new TmtcReader));
68    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Mcmd,            new McmdReader));
69    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::ForcedFECmd,     new ForcedFECmdReader));
70    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac1Init,         new Ac1InitReader));
71    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalInit,         new CalInitReader));
72    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new NdInitReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkInit,         new TrkInitReader));
73    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofInit,         new TofInitReader));
74    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgInit,         new TrgInitReader));
75    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::NdInit,          new NdInitReader));
76    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Init,          new S4InitReader));
77    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader()));    TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::Ac2Init,         new Ac2InitReader));
78      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::CalAlarm,        new CalAlarmReader));
79      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::AcAlarm,         new AcAlarmReader));
80      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrkAlarm,        new TrkAlarmReader));
81      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TrgAlarm,        new TrgAlarmReader));
82      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::TofAlarm,        new TofAlarmReader));
83      TechmodelAlgorithmMap.insert(AlgorithmMap::value_type(PacketType::S4Alarm,         new S4AlarmReader));
84  }  }
85    
86  /**  /**
# Line 82  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.3 2004/10/17 12:28:38 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/EventReader.cpp,v 2.8 2005/01/11 14:59:31 kusanagi Exp $\n";
92  }  }
93    
94  /**  /**
# Line 123  void EventReader::RunEvent(int EventNumb Line 129  void EventReader::RunEvent(int EventNumb
129    //to insert a GetPacketStart()    //to insert a GetPacketStart()
130    stringstream oss;    stringstream oss;
131    int step = 0;    int step = 0;
132      const PacketType* type;
133    while (!InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){    while (!InputFile->eof() && ((step++ < maxPackets) || (maxPackets == 0))){
134      try {      try {
135        if (FindStart()) {        if (FindStart()) {
136              UnpackPscuHeader();              UnpackPscuHeader();
137              const PacketType* type = Header->GetPscuHeader()->GetPacketType();              type = Header->GetPscuHeader()->GetPacketType();
138              AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type);              AlgorithmMap::iterator i = TechmodelAlgorithmMap.find(type);
139              if (i != TechmodelAlgorithmMap.end()) {              if (i != TechmodelAlgorithmMap.end()) {
140                  TechmodelAlgorithm *EventAlgorithm = i->second;                  TechmodelAlgorithm *EventAlgorithm(i->second);
141                  EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());                  EventAlgorithm->RunEvent(EventNumber, Header->GetPscuHeader()->GetPacketLenght());
142                  Run->FillTrees(type);                  Run->FillTrees(type);
143                  Header->GetCounter()->Increment(type);                  Header->GetCounter()->Increment(type);
144                  logger->info(Header->GetPscuHeader()->Print());                  logger->info(Header->GetPscuHeader()->Print());
145              } else {              } else {
146                  throw NotExistingAlgorithmException(type->GetName().c_str()); //to exctract to an higher level and delete the logger!                  oss.str("");
147                    oss << "\n No way to read events of type  " << type->GetName().c_str() << Header->GetPscuHeader()->Print();
148                    throw NotExistingAlgorithmException(oss.str().c_str()); //to exctract to an higher level and delete the logger!
149              }              }
150        }        }
151      // 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??
152      } catch (NotExistingAlgorithmException exc) {      } catch (NotExistingAlgorithmException exc) {
153        char *temp;        oss.str("");
154        sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
155        logger->warn(temp);        logger->error(oss.str().c_str());
156      } catch (WrongCRCHeaderException exc) {      } catch (WrongCRCHeaderException exc) {
157        char *temp;        oss.str("");
158        sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());        oss << exc.print();
159        logger->warn(temp);        logger->error(oss.str().c_str());
160      } catch (WrongCRCException exc) {      } catch (WrongCRCException exc) {
161        //char *temp;        oss.str("");
162        //sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
163        //logger->warn(temp);        logger->error(oss.str().c_str());
         logger->warn(Header->GetPscuHeader()->Print());  
164      } catch (UnidentifiedPacketException exc) {      } catch (UnidentifiedPacketException exc) {
165        char *temp;        oss.str("");
166        sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
167        logger->warn(temp);        logger->error(oss.str().c_str());
168      } catch (NotExistingCounterException exc) {      } catch (NotExistingCounterException exc) {
169        char *temp;        oss.str("");
170        sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
171        logger->warn(temp);        logger->error(oss.str().c_str());
172      } catch (LengthException exc) {      } catch (LengthException exc) {
173        char *temp;        oss.str("");
174        sprintf(temp, "%s %s \n", exc.print(), Header->GetPscuHeader()->Print());        oss << exc.print() << " " << Header->GetPscuHeader()->Print();
175        logger->warn(temp);        logger->error(oss.str().c_str());
176      } catch (...) {      } catch (...) {
177        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");
178      }      }
# Line 182  static void SkipToNextHeader(ifstream* T Line 190  static void SkipToNextHeader(ifstream* T
190  /**  /**
191   * Unpack the PSCU header from a file into the structure.   * Unpack the PSCU header from a file into the structure.
192   */   */
193  void EventReader::UnpackPscuHeader(void) throw (WrongCRCException, LengthException) {  void EventReader::UnpackPscuHeader(void) throw (WrongCRCHeaderException, LengthException) {
194    stringstream oss;    stringstream oss;
195    int response = false;    int response = false;
196    char buff[16];    char buff[16];
# Line 225  void EventReader::UnpackPscuHeader(void) Line 233  void EventReader::UnpackPscuHeader(void)
233            
234      //commented out because of the above test code      //commented out because of the above test code
235      InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);      InputFile->seekg(Header->GetPscuHeader()->GetPacketLenght(), std::ios::cur);
236      FindStart();      if (FindStart()) {
237      finalPos =  (long int)InputFile->tellg() - (1 + initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));          finalPos =  (long int)InputFile->tellg() - (initPos + (long int)(Header->GetPscuHeader()->GetPacketLenght()));
238          if(finalPos == 0){          if(finalPos == 0){
239              logger->debug(_T("Correct packet length"));              logger->debug(_T(" Correct packet length"));
240              }              }
241          if (finalPos > 0 && finalPos < 64) {          if (finalPos > 0 && finalPos < 64) {
242              oss.str("");              oss.str("");
243              oss << " Correct packet length: Padded of " << finalPos << " bytes";              oss << "\n Correct packet length: Padded of " << finalPos << " bytes";
244              logger->debug(oss.str().c_str());              logger->debug(oss.str().c_str());
245          }          }
246          if (finalPos > 64){          if (finalPos > 64){
247              oss.str("");              oss.str("");
248              oss << "The begin of the next packet is far more than 64 byte from the end of the previous."              oss << "\n The begin of the next packet is far more than 64 byte from the end of the previous."
249                  << Header->GetPscuHeader()->Print();                  << Header->GetPscuHeader()->Print();
250              logger->warn(oss.str().c_str());              logger->error(oss.str().c_str());
251              //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");              //throw LengthException("The begin of the next packet is far more than 64 byte from the end of the previous.");
252          }          }
253        }
254        else {
255            logger->debug(_T(" END OF FILE"));        
256        }
257            
258          InputFile->seekg(initPos, std::ios::beg);          InputFile->seekg(initPos, std::ios::beg);
259    } else {    } else {
260         /*Here i should extract the block of Data for later analysis */
261      InputFile->seekg(-(13), std::ios::cur);      InputFile->seekg(-(13), std::ios::cur);
262      throw WrongCRCHeaderException();      oss.str("");
263        oss << "CRC Header Error on packet:" << PscuHeader::Print(buff);
264        throw WrongCRCHeaderException(oss.str().c_str());
265    }    }
266    
267  /*    char tmpId1[4];  /*    char tmpId1[4];
# Line 279  void EventReader::UnpackPscuTrailer(void Line 295  void EventReader::UnpackPscuTrailer(void
295  /**  /**
296   * 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
297   */   */
298  int EventReader::FindStart(void) throw (std::exception) {  bool EventReader::FindStart(void) throw (std::exception) {
299  //search an hexadecimal sequence in a file  //search an hexadecimal sequence in a file
300  //subSign    ------> pointer to the sequence buffer  //subSign    ------> pointer to the sequence buffer
301  //subSignDim ------> dimension of the buffer  //subSignDim ------> dimension of the buffer
# Line 307  int EventReader::FindStart(void) throw ( Line 323  int EventReader::FindStart(void) throw (
323          if (dataByte == (char)(*(subSign+subIndex))){          if (dataByte == (char)(*(subSign+subIndex))){
324              if (subIndex++ == (subSignDim-1)) {              if (subIndex++ == (subSignDim-1)) {
325                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);                  InputFile->seekg( (index - (subIndex + buffSize)), std::ios::cur);
326                  return 1;                  return true;
327              }              }
328          } else {          } else {
329              index = index - (subIndex);              index = index - (subIndex);
# Line 316  int EventReader::FindStart(void) throw ( Line 332  int EventReader::FindStart(void) throw (
332      }      }
333      //Needs to guarantee the overap of the buffer(s) in several loop      //Needs to guarantee the overap of the buffer(s) in several loop
334      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);      InputFile->seekg( (-1)*(subSignDim + 1) , std::ios::cur);
335    }    }
336    return 0;    return false;
337  }  }
338    
339  ClassImp(EventReader)  ClassImp(EventReader)
Legend:
Removed from v.2.4  
changed lines
  Added in v.2.9
  ViewVC Help
Powered by ViewVC 1.1.23