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