/[PAMELA software]/yoda/techmodel/TabDumpReader.cpp
ViewVC logotype

Diff of /yoda/techmodel/TabDumpReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.3 by kusanagi, Mon Jul 26 23:09:45 2004 UTC revision 4.4 by kusanagi, Sat May 28 10:44:11 2005 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Source: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v $   * $Source: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v $
3   * $Id: TabDumpReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $   * $Id: TabDumpReader.cpp,v 4.0 2005/03/06 04:33:02 kusanagi Exp $
4   * $Author: kusanagi $   * $Author: kusanagi $
5   *   *
6   * Implementation of the TabDumpReader class.   * Implementation of the TabDumpReader class.
7   */   */
8    
9  #include <string>  #include <string>
10  #include <log4cpp/Category.hh>  #include <log4cxx/logger.h>
11  #include <fstream>  #include <fstream>
12  #include "stdio.h"  #include "stdio.h"
13  extern "C" {  extern "C" {
# Line 18  extern "C" { Line 18  extern "C" {
18  using namespace pamela;  using namespace pamela;
19  using namespace pamela::techmodel;  using namespace pamela::techmodel;
20    
21  static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.TabDumpReader");    static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.TabDumpReader"));
22    
23  /**  /**
24   * Constructor.   * Constructor.
# Line 26  static log4cpp::Category& cat = log4cpp: Line 26  static log4cpp::Category& cat = log4cpp:
26  TabDumpReader::TabDumpReader(void):  TabDumpReader::TabDumpReader(void):
27    TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {    TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {
28    TabDump = new TabDumpEvent();    TabDump = new TabDumpEvent();
29      logger->debug(_T("Constructor"));
30  }  }
31    
32  /**  /**
33   * Get a string with the version info of the algorithm.   * Get a string with the version info of the algorithm.
34   */   */
35  std::string TabDumpReader::GetVersionInfo(void) const {  std::string TabDumpReader::GetVersionInfo(void) const {
36    return "$Header: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v 1.2 2004/07/17 20:03:38 kusanagi Exp $\n";    return "$Header: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v 4.0 2005/03/06 04:33:02 kusanagi Exp $\n";
37  }  }
38    
39  /**  /**
# Line 40  std::string TabDumpReader::GetVersionInf Line 41  std::string TabDumpReader::GetVersionInf
41   * event reader routines for all packet types.   * event reader routines for all packet types.
42   */   */
43  void TabDumpReader::Init(PamelaRun *run) {  void TabDumpReader::Init(PamelaRun *run) {
44      logger->debug(_T("Initialize"));
45    SetInputStream(run);    SetInputStream(run);
46    run->WriteSubPacket(this, &TabDump, TabDump->Class());    run->WriteSubPacket(this, &TabDump, TabDump->Class());
47  }  }
# Line 47  void TabDumpReader::Init(PamelaRun *run) Line 49  void TabDumpReader::Init(PamelaRun *run)
49  /**  /**
50   * Unpack the TabDump event from an input file.   * Unpack the TabDump event from an input file.
51   */   */
52  void TabDumpReader::RunEvent(int EventNumber, long int length) {  void TabDumpReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
53        
54    char     *subData;      char      subData[length];
55    char      eventCRC[2];      UINT16    subCRC;      //calculated CRC of the data
56    UINT16    subCRC;      //CRC of the data      UINT16    readCRC;     //CRC read from the end of the subpacket
57    UINT16    readCRC;     //CRC read from the end of the subpacket      long int  dataLength = length - 2; //the block of data
   long int  dataLength;  
58        
59      //the 2 bytes subtracted belong to the final event CRC bytes      InputFile->read(subData, sizeof(subData));
60      dataLength = length - (long int)2;      subCRC  = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
61        readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
62        
63        //if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
64        if (subCRC != readCRC) {
65            logger->error("WRONG CRC FOR TabDump PACKET. Processed anyway for CPU debugging");
66        }
67    
68      subData = new char[dataLength];      TabDumpRecord* rec;
69      InputFile->read(subData, sizeof(unsigned char)*dataLength);      TabDump->PARAMETER_STAMP = (((UINT32)subData[0]<<24)&0xFF000000) + (((UINT32)subData[1]<<16)&0x00FF0000) +  (((UINT32)subData[2]<<8)&0x0000FF00) + (((UINT32)subData[3])&0x000000FF);
70      subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);      long int offset = 4;
71        int i = 0;
72      //took the final CRC to compare it with the previous calculated CRC of the data      TabDump->Records->Clear();
73      InputFile->read(eventCRC, sizeof(eventCRC));      TClonesArray &recs = *(TabDump->Records);
74      readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);      while (offset < dataLength){
75            rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
76      //finally check if the RecordCRC is correct and finally update the partialCRC          rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;
77      if (subCRC == readCRC){          rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
78          TabDumpRecord* rec;          rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
79          long int offset = 0;          rec->Data = new TArrayI((int)((rec->Nrow)*(rec->Ncol)), (int*)(subData+offset+3));
80          int i = 0;          offset = offset + sizeof(UINT32)*(rec->Nrow)*(rec->Ncol) + 3;
         int size;  
         TabDump->Records->Clear();  
         TClonesArray &recs = *(TabDump->Records);  
         while (offset < dataLength){  
             rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump  
             rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;  
             rec->Nrow = ((UINT8)subData[offset+1])&0xFF;  
             rec->Ncol = ((UINT8)subData[offset+2])&0xFF;  
             rec->Data = new UINT32[(rec->Nrow)*(rec->Ncol)];  
             size = sizeof(UINT32)*(rec->Nrow)*(rec->Ncol);  
             memcpy(rec->Data, (UINT32*)(subData+offset+3), size);  
             offset = offset + size + 3;  
         }  
         free(subData);  
     } else {  
         cat <<  log4cpp::Priority::ERROR  
             <<  "Wrong CRC on Subpacket in TabDump Packet "  
             <<  "\n " << log4cpp::CategoryStream::ENDLINE;  
81      }      }
82    }  }
Legend:
Removed from v.1.3  
changed lines
  Added in v.4.4
  ViewVC Help
Powered by ViewVC 1.1.23