/[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 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.2 2004/07/17 20:03:38 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 <log4cpp/Category.hh>  
 #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 log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.TabDumpReader");    static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.TabDumpReader"));
16    
17  /**  /**
18   * Constructor.   * Constructor.
# Line 26  static log4cpp::Category& cat = log4cpp: Line 20  static log4cpp::Category& cat = log4cpp:
20  TabDumpReader::TabDumpReader(void):  TabDumpReader::TabDumpReader(void):
21    TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {    TechmodelAlgorithm(PacketType::TabDump, "TechmodelTabDump") {
22    TabDump = new TabDumpEvent();    TabDump = new TabDumpEvent();
23      logger->debug(_T("Constructor"));
24  }  }
25    
26  /**  /**
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.2 2004/07/17 20:03:38 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 40  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());
41  }  }
# Line 47  void TabDumpReader::Init(PamelaRun *run) Line 43  void TabDumpReader::Init(PamelaRun *run)
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){
     
   char     *subData;  
   char      eventCRC[2];  
   UINT16    subCRC;      //CRC of the data  
   UINT16    readCRC;     //CRC read from the end of the subpacket  
   long int  dataLength;  
47        
48      //the 2 bytes subtracted belong to the final event CRC bytes      char      subData[length];
49      dataLength = length - (long int)2;      UINT16    subCRC;      //calculated CRC of the data
50        UINT16    readCRC;     //CRC read from the end of the subpacket
51      subData = new char[dataLength];      long int  dataLength = length - 2; //the block of data
52      InputFile->read(subData, sizeof(unsigned char)*dataLength);      memset(subData,  0, length*sizeof(char));  
53      subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);      InputFile->read(subData, sizeof(subData));
54        subCRC  = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
55      //took the final CRC to compare it with the previous calculated CRC of the data      readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
56      InputFile->read(eventCRC, sizeof(eventCRC));      
57      readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);      if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for TabDump Packet ");
58        /*if (subCRC != readCRC) {
59      //finally check if the RecordCRC is correct and finally update the partialCRC          logger->error("WRONG CRC FOR TabDump PACKET. Processed anyway for CPU debugging");
60      if (subCRC == readCRC){      }*/
61          TabDumpRecord* rec;  
62          long int offset = 0;      TabDumpRecord* rec;
63          int i = 0;      TabDump->PARAMETER_STAMP = (((UINT32)subData[0]<<24)&0xFF000000) + (((UINT32)subData[1]<<16)&0x00FF0000) +  (((UINT32)subData[2]<<8)&0x0000FF00) + (((UINT32)subData[3])&0x000000FF);
64          int size;      long int offset = 4;
65          TabDump->Records->Clear();      int i = 0;
66          TClonesArray &recs = *(TabDump->Records);      TabDump->Records->Clear();
67          while (offset < dataLength){      TClonesArray &recs = *(TabDump->Records);
68              rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump      while (offset < dataLength){
69              rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;          rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
70              rec->Nrow = ((UINT8)subData[offset+1])&0xFF;          rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;
71              rec->Ncol = ((UINT8)subData[offset+2])&0xFF;          rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
72              rec->Data = new UINT32[(rec->Nrow)*(rec->Ncol)];          rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
73              size = sizeof(UINT32)*(rec->Nrow)*(rec->Ncol);          rec->Data = new TArrayI((int)((rec->Nrow)*(rec->Ncol)), (int*)(subData+offset+3));
74              memcpy(rec->Data, (UINT32*)(subData+offset+3), size);          offset = offset + sizeof(UINT32)*(rec->Nrow)*(rec->Ncol) + 3;
             offset = offset + size + 3;  
         }  
         free(subData);  
     } else {  
         cat <<  log4cpp::Priority::ERROR  
             <<  "Wrong CRC on Subpacket in TabDump Packet "  
             <<  "\n " << log4cpp::CategoryStream::ENDLINE;  
75      }      }
76    }  }

Legend:
Removed from v.1.3  
changed lines
  Added in v.6.3

  ViewVC Help
Powered by ViewVC 1.1.23