/[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.3 by mocchiut, Fri Sep 29 10:19:27 2006 UTC
# Line 1  Line 1 
1  /** @file  /** @file
2   * $Source: /home/cvsmanager/yoda/techmodel/TabDumpReader.cpp,v $   * $Source: /afs/ba.infn.it/user/pamela/src/CVS/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.2 2006/05/30 19:10:03 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: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/TabDumpReader.cpp,v 6.2 2006/05/30 19:10:03 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
52    long int  dataLength;      memset(subData,  0, length*sizeof(char));  
53          InputFile->read(subData, sizeof(subData));
54      //the 2 bytes subtracted belong to the final event CRC bytes      subCRC  = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
55      dataLength = length - (long int)2;      readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
56        
57      subData = new char[dataLength];      if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
58      InputFile->read(subData, sizeof(unsigned char)*dataLength);      /*if (subCRC != readCRC) {
59      subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);          logger->error("WRONG CRC FOR TabDump PACKET. Processed anyway for CPU debugging");
60        }*/
61      //took the final CRC to compare it with the previous calculated CRC of the data  
62      InputFile->read(eventCRC, sizeof(eventCRC));      TabDumpRecord* rec;
63      readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);      TabDump->PARAMETER_STAMP = (((UINT32)subData[0]<<24)&0xFF000000) + (((UINT32)subData[1]<<16)&0x00FF0000) +  (((UINT32)subData[2]<<8)&0x0000FF00) + (((UINT32)subData[3])&0x000000FF);
64        long int offset = 4;
65      //finally check if the RecordCRC is correct and finally update the partialCRC      int i = 0;
66      if (subCRC == readCRC){      TabDump->Records->Clear();
67          TabDumpRecord* rec;      TClonesArray &recs = *(TabDump->Records);
68          long int offset = 0;      while (offset < dataLength){
69          int i = 0;          rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
70          int size;          rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;
71          TabDump->Records->Clear();          rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
72          TClonesArray &recs = *(TabDump->Records);          rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
73          while (offset < dataLength){          rec->Data = new TArrayI((int)((rec->Nrow)*(rec->Ncol)), (int*)(subData+offset+3));
74              rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump          offset = offset + sizeof(UINT32)*(rec->Nrow)*(rec->Ncol) + 3;
             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);  
             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.3
  ViewVC Help
Powered by ViewVC 1.1.23