/[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 2.1.2.2 by kusanagi, Fri Mar 4 10:42:14 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 2.1.2.1 2005/03/04 10:26:41 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 2.1.2.1 2005/03/04 10:26:41 kusanagi Exp $\n";
37  }  }
38    
39  /**  /**
# Line 42  std::string TabDumpReader::GetVersionInf Line 43  std::string TabDumpReader::GetVersionInf
43  void TabDumpReader::Init(PamelaRun *run) {  void TabDumpReader::Init(PamelaRun *run) {
44    SetInputStream(run);    SetInputStream(run);
45    run->WriteSubPacket(this, &TabDump, TabDump->Class());    run->WriteSubPacket(this, &TabDump, TabDump->Class());
46      logger->debug(_T("Initialize"));
47  }  }
48    
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        char      subData[length];
54        char      eventCRC[2];
55        UINT16    subCRC;      //CRC of the data
56        UINT16    readCRC;     //CRC read from the end of the subpacket
57        long int  dataLength = length - 2;
58        
59    char     *subData;      InputFile->read(subData, sizeof(unsigned char)*length);
60    char      eventCRC[2];      subCRC  = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);
61    UINT16    subCRC;      //CRC of the data      readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
62    UINT16    readCRC;     //CRC read from the end of the subpacket      
63    long int  dataLength;      /*Pamela CPU group request to disable for a while CRC control to check values because
64            the data seems to be well formatted but perhaps there is just a mistake on CRC evaluation
65      //the 2 bytes subtracted belong to the final event CRC bytes        into the PAMELA CPU program
66      dataLength = length - (long int)2;      */
67        if (subCRC != readCRC) {
68            logger->error(" Wrong CRC for TabDump Packet but I process it anyway because of CPU DEBUG");
69        }
70    
71      subData = new char[dataLength];      TabDumpRecord* rec;
72      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);
73      subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength);      long int offset = 4;
74        int i = 0;
75      //took the final CRC to compare it with the previous calculated CRC of the data      int size;
76      InputFile->read(eventCRC, sizeof(eventCRC));      TabDump->Records->Clear();
77      readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);      TClonesArray &recs = *(TabDump->Records);
78        while (offset < dataLength){
79      //finally check if the RecordCRC is correct and finally update the partialCRC          rec = new(recs[i++]) TabDumpRecord(); //add a new TabDump
80      if (subCRC == readCRC){          rec->Tab_ID  = ((UINT8)subData[offset])&0xFF;
81          TabDumpRecord* rec;          rec->Nrow = ((UINT8)subData[offset+1])&0xFF;
82          long int offset = 0;          rec->Ncol = ((UINT8)subData[offset+2])&0xFF;
83          int i = 0;          //rec->Data = new UINT32[(rec->Nrow)*(rec->Ncol)];
84          int size;          size = sizeof(UINT32)*(rec->Nrow)*(rec->Ncol);
85          TabDump->Records->Clear();          rec->Data = new TArrayI((int)((rec->Nrow)*(rec->Ncol)), (int*)(subData+offset+3));
86          TClonesArray &recs = *(TabDump->Records);          //memcpy(rec->Data, (UINT32*)(subData+offset+3), size);
87          while (offset < dataLength){          offset = offset + size + 3;
             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;  
88      }      }
89       delete [] subData;
90    }    }
Legend:
Removed from v.1.3  
changed lines
  Added in v.2.1.2.2
  ViewVC Help
Powered by ViewVC 1.1.23