/[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 2.1 by kusanagi, Fri Dec 3 22:08:01 2004 UTC revision 3.1 by kusanagi, Sun Mar 6 04:29:20 2005 UTC
# Line 15  extern "C" { Line 15  extern "C" {
15  #include "stdio.h"  #include "stdio.h"
16  #include "ReaderAlgorithms.h"  #include "ReaderAlgorithms.h"
17    
 #include "event/PhysEndRunEvent.h"  
   
18  using namespace pamela;  using namespace pamela;
19  using namespace pamela::techmodel;  using namespace pamela::techmodel;
20    
# Line 36  PhysEndRunReader::PhysEndRunReader(void) Line 34  PhysEndRunReader::PhysEndRunReader(void)
34   */   */
35  std::string PhysEndRunReader::GetVersionInfo(void) const {  std::string PhysEndRunReader::GetVersionInfo(void) const {
36    return    return
37      "$Header: /home/cvsmanager/yoda/techmodel/PhysEndRunReader.cpp,v 2.0 2004/09/21 20:50:54 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/PhysEndRunReader.cpp,v 3.0 2005/03/04 15:54:11 kusanagi Exp $\n";
38  }  }
39    
40  /**  /**
# Line 52  void PhysEndRunReader::Init(PamelaRun *r Line 50  void PhysEndRunReader::Init(PamelaRun *r
50  /**  /**
51   * Unpack the PhysEndRun event from an input file.   * Unpack the PhysEndRun event from an input file.
52   */   */
53    
54  void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){  void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
55        char subData[length];
56        UINT16    subCRC;      //CRC of the data
57        UINT16    readCRC;     //CRC read from the end of the subpacket
58        long int  dataLength;
59        long int  tbInitPos;
60    
61        //int numRecords = ((length - 2)/6);
62        dataLength = length - (long int)2;    
63        tbInitPos = dataLength - 19; // 19 is the length of th TB
64        //physEndRun->CALO_ENDRUN[];
65        //physEndRun->TB_ENDRUN;
66        
67        InputFile->read(subData, sizeof(unsigned char)*length);
68        subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
69        readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
70            
71        //if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for PhysEndRun Packet ");
72        /*if (subCRC != readCRC) {
73            logger->error(" Wrong CRC for PhysEndRun Packet but I process it anyway");
74        }*/
75        //--------------- CALO SECTION------------------------------------
76        long int offset = 0;
77        int j = 0;
78        //physEndRun->CALO_ENDRUN.
79        while(offset < tbInitPos){
80            physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK     = (((UINT8)subData[offset])&0x1F);
81            physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
82            //subData[offset+2--->3] is a fixed word = 0x23
83            //If this bit is ? the subsequent CALO_HK are not valid
84            if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
85                for(int k = 0; k < 11 ; k++){
86                    physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+3+k]<<8))&0xFF00) + (((UINT16)subData[offset+4+k])&0x00FF);
87                    physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+14+k]<<8))&0xFF00) + (((UINT16)subData[offset+15+k])&0x00FF);
88                }
89                offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
90            } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
91            //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
92    
93            for(int k = 0; k < 7 ; k++){
94                physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG     = ((UINT8)(subData[offset])&0x1F);
95                physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
96                //subData[offset+2--->3] is a fixed word = 0x2
97                physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
98                //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
99                offset = offset + 8; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
100            }
101            j++;
102        }
103            
104        //--------------- TB SECTION------------------------------------
105        physEndRun->TB_ENDRUN.TB_ALARM_MASK   = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
106        //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
107        //Here are missing the TB_PMT_MASK;
108        //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
109        physEndRun->TB_ENDRUN.TB_S4_MASK      = ((UINT8)subData[offset+10])>>5;
110        physEndRun->TB_ENDRUN.TB_CALO_MASK    = (((UINT8)subData[offset+10])>>1)&0x0F;
111        //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
112        physEndRun->TB_ENDRUN.TB_BUSY_MASK    = (((UINT32)subData[offset+12]<<16)&0x00FF0000) +  (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
113        //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
114        physEndRun->TB_ENDRUN.TB_CALIB_FLAG   = ((UINT8)subData[offset+16])>>7;
115        physEndRun->TB_ENDRUN.TB_S4_TRIG      = ((UINT8)subData[offset+16])>>6;
116        physEndRun->TB_ENDRUN.TB_CALO_TRIG    = ((UINT8)subData[offset+16])>>5;
117        physEndRun->TB_ENDRUN.TB_TOF_TRIG     = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>5);
118        //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
119  }  }
120    
Legend:
Removed from v.2.1  
changed lines
  Added in v.3.1
  ViewVC Help
Powered by ViewVC 1.1.23