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

Diff of /yoda/techmodel/PhysEndRunReader.cpp

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

revision 1.1 by kusanagi, Thu Aug 19 15:24:46 2004 UTC revision 6.1 by kusanagi, Tue May 30 19:10:02 2006 UTC
# Line 1  Line 1 
   
1  // Implementation of the PhysEndRunReader class.  // Implementation of the PhysEndRunReader class.
   
2    
 #define UINT unsigned int  
 #define BYTE  unsigned char  
 #include <string>  
 #include <log4cpp/Category.hh>  
3  extern "C" {  extern "C" {
 #include <sys/time.h>  
4  #include "CRC.h"  #include "CRC.h"
5  }  }
6    
 #include <fstream>  
 #include "stdio.h"  
7  #include "ReaderAlgorithms.h"  #include "ReaderAlgorithms.h"
8    
 #include "event/PhysEndRunEvent.h"  
   
 using namespace pamela;  
9  using namespace pamela::techmodel;  using namespace pamela::techmodel;
10    
11  static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.PhysEndRunReader");    static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.PhysEndRunReader"));
12    
13  /**  /**
14   * Constructor.   * Constructor.
15   */   */
16  PhysEndRunReader::PhysEndRunReader(void):  PhysEndRunReader::PhysEndRunReader(void):
17    TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") {    TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") {
18    cat <<  log4cpp::Priority::DEBUG    logger->debug(_T("Constructor"));
       <<  "Constructor "  
       <<  "\n " << log4cpp::CategoryStream::ENDLINE;  
19    physEndRun = new PhysEndRunEvent();    physEndRun = new PhysEndRunEvent();
20  }  }
21    
# Line 38  PhysEndRunReader::PhysEndRunReader(void) Line 24  PhysEndRunReader::PhysEndRunReader(void)
24   */   */
25  std::string PhysEndRunReader::GetVersionInfo(void) const {  std::string PhysEndRunReader::GetVersionInfo(void) const {
26    return    return
27      "$Header: /home/cvsmanager/yoda/techmodel/PhysEndRunReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/PhysEndRunReader.cpp,v 6.0 2006/02/07 17:11:09 kusanagi Exp $\n";
28  }  }
29    
30  /**  /**
# Line 48  std::string PhysEndRunReader::GetVersion Line 34  std::string PhysEndRunReader::GetVersion
34  void PhysEndRunReader::Init(PamelaRun *run) {  void PhysEndRunReader::Init(PamelaRun *run) {
35    SetInputStream(run);    SetInputStream(run);
36    run->WriteSubPacket(this, &physEndRun, physEndRun->Class());    run->WriteSubPacket(this, &physEndRun, physEndRun->Class());
37      logger->debug(_T("Initialize"));
38  }  }
39    
40  /**  /**
41   * Unpack the PhysEndRun event from an input file.   * Unpack the PhysEndRun event from an input file.
42   */   */
43  void PhysEndRunReader::RunEvent(int EventNumber, long int length) {  
44    void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
45        char subData[length];
46        UINT16    subCRC;      //CRC of the data
47        UINT16    readCRC;     //CRC read from the end of the subpacket
48        long int  dataLength = length - 2;    
49        long int  tbInitPos  = dataLength - 19; // 19 is the length of th TB
50    
51        //physEndRun->CALO_ENDRUN[];
52        //physEndRun->TB_ENDRUN;
53        
54        InputFile->read(subData, sizeof(subData));
55        subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
56        readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
57            
58        //if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for PhysEndRun Packet ");
59        if (subCRC != readCRC) {
60            logger->error(" Wrong CRC for PhysEndRun Packet but I process it anyway");
61        }
62        //--------------- CALO SECTION------------------------------------
63        long int offset = 0;
64        int j = 0;
65        //physEndRun->CALO_ENDRUN.
66        while(offset < tbInitPos){
67            physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK     = (((UINT8)subData[offset])&0x1F);
68            physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
69            //subData[offset+2--->3] is a fixed word = 0x23
70            //If this bit is ? the subsequent CALO_HK are not valid
71            if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
72                for(int k = 0; k < 11 ; k++){
73                    physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+3+k]<<8))&0xFF00) + (((UINT16)subData[offset+4+k])&0x00FF);
74                    physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+14+k]<<8))&0xFF00) + (((UINT16)subData[offset+15+k])&0x00FF);
75                }
76                offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
77            } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
78            //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
79    
80            for(int k = 0; k < 7 ; k++){
81                physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG     = ((UINT8)(subData[offset])&0x1F);
82                physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
83                //subData[offset+2--->3] is a fixed word = 0x2
84                physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
85                //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
86                offset = offset + 8; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
87            }
88            j++;
89        }
90            
91        //--------------- TB SECTION------------------------------------
92        physEndRun->TB_ENDRUN.TB_ALARM_MASK   = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
93        //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
94        physEndRun->TB_ENDRUN.TB_PMT_MASK_S3  = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
95        physEndRun->TB_ENDRUN.TB_PMT_MASK_S2  = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
96        physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
97        physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
98        //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
99        physEndRun->TB_ENDRUN.TB_CALO_MASK    = ((UINT8)subData[offset+10])>>5;
100        physEndRun->TB_ENDRUN.TB_S4_MASK      = (((UINT8)subData[offset+10])>>1)&0x0F;
101        //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
102        physEndRun->TB_ENDRUN.TB_BUSY_MASK    = (((UINT32)subData[offset+12]<<16)&0x00FF0000) +  (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
103        //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
104        physEndRun->TB_ENDRUN.TB_CALIB_FLAG   = ((UINT8)subData[offset+16])>>7;
105        physEndRun->TB_ENDRUN.TB_CALO_TRIG    = ((UINT8)subData[offset+16])>>6;
106        physEndRun->TB_ENDRUN.TB_S4_TRIG      = ((UINT8)subData[offset+16])>>5;
107        physEndRun->TB_ENDRUN.TB_TOF_TRIG     = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
108        //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
109  }  }
110    
Legend:
Removed from v.1.1  
changed lines
  Added in v.6.1
  ViewVC Help
Powered by ViewVC 1.1.23