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

Annotation of /yoda/techmodel/PhysEndRunReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 6.7 - (hide annotations) (download)
Mon Feb 11 11:04:27 2008 UTC (16 years, 9 months ago) by mocchiut
Branch: MAIN
CVS Tags: YODA6_3/20, HEAD
Changes since 6.6: +31 -16 lines
Calorimeter registry bug in PhysEndRunReader.cpp fixed

1 kusanagi 1.1 // Implementation of the PhysEndRunReader class.
2    
3     extern "C" {
4     #include "CRC.h"
5     }
6    
7     #include "ReaderAlgorithms.h"
8    
9     using namespace pamela::techmodel;
10    
11 kusanagi 1.2 static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.PhysEndRunReader"));
12 kusanagi 1.1
13     /**
14     * Constructor.
15     */
16     PhysEndRunReader::PhysEndRunReader(void):
17     TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") {
18 kusanagi 1.2 logger->debug(_T("Constructor"));
19 kusanagi 1.1 physEndRun = new PhysEndRunEvent();
20     }
21    
22     /**
23     * Get a string with the version info of the algorithm.
24     */
25     std::string PhysEndRunReader::GetVersionInfo(void) const {
26     return
27 mocchiut 6.7 "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/PhysEndRunReader.cpp,v 6.6 2008/01/10 16:21:38 mocchiut Exp $\n";
28 kusanagi 1.1 }
29    
30     /**
31     * Initialize the algorithm with a special run. This will initialize the
32     * event reader routines for all packet types.
33     */
34     void PhysEndRunReader::Init(PamelaRun *run) {
35     SetInputStream(run);
36     run->WriteSubPacket(this, &physEndRun, physEndRun->Class());
37 kusanagi 1.2 logger->debug(_T("Initialize"));
38 kusanagi 1.1 }
39    
40     /**
41     * Unpack the PhysEndRun event from an input file.
42     */
43 kusanagi 2.3
44 kusanagi 2.1 void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
45 kusanagi 2.3 char subData[length];
46     UINT16 subCRC; //CRC of the data
47     UINT16 readCRC; //CRC read from the end of the subpacket
48 kusanagi 5.1 long int dataLength = length - 2;
49     long int tbInitPos = dataLength - 19; // 19 is the length of th TB
50 kusanagi 2.3
51     //physEndRun->CALO_ENDRUN[];
52     //physEndRun->TB_ENDRUN;
53 mocchiut 6.5 memset(subData, 0, length*sizeof(char));
54 kusanagi 5.1 InputFile->read(subData, sizeof(subData));
55 pam-rm2 6.4 subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
56 kusanagi 2.3 readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
57    
58 kusanagi 5.1 //if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for PhysEndRun Packet ");
59     if (subCRC != readCRC) {
60 kusanagi 2.3 logger->error(" Wrong CRC for PhysEndRun Packet but I process it anyway");
61 pam-rm2 6.4 }
62 kusanagi 2.3 //--------------- 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 mocchiut 6.6 UInt_t kk = 0;
73 kusanagi 2.3 for(int k = 0; k < 11 ; k++){
74 mocchiut 6.6 // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]);
75     physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF);
76     physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF);
77     kk += 2;
78     }
79 kusanagi 2.3 offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
80     } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
81     //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
82 mocchiut 6.7 //
83     for(int k = 0; k < 7 ; k++){
84     //
85     // in CALO_BOARD_ID_REG and CALO_BOARD_STATUS_REG we save only the last id and status register of the 7 and we don't save the CRC for each register.
86     //
87     // printf(" INLOPP k %i subdata %X \n",k,(UINT8)(subData[offset]));
88     physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F);
89     offset++;
90     physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset]);
91     offset++;
92     offset++;
93     //subData[offset+2--->3] is a fixed word = 0x2 that is the data lenght
94     offset++;
95     physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset+1])&0x00FF);
96     offset++;
97     // CRC
98     offset++;
99     offset++;
100     //
101     offset++;
102     // physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF);
103     // physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0;
104     // physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0;
105     // physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0;
106     // physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0;
107     // physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0;
108     //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
109     // offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
110     };
111     //
112 kusanagi 2.3 j++;
113     }
114 kusanagi 1.1
115 kusanagi 2.3 //--------------- TB SECTION------------------------------------
116     physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
117     //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
118 kusanagi 4.2 physEndRun->TB_ENDRUN.TB_PMT_MASK_S3 = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
119     physEndRun->TB_ENDRUN.TB_PMT_MASK_S2 = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
120     physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
121     physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
122 kusanagi 2.3 //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
123 kusanagi 4.4 physEndRun->TB_ENDRUN.TB_CALO_MASK = ((UINT8)subData[offset+10])>>5;
124     physEndRun->TB_ENDRUN.TB_S4_MASK = (((UINT8)subData[offset+10])>>1)&0x0F;
125 kusanagi 2.3 //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
126     physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)subData[offset+12]<<16)&0x00FF0000) + (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
127     //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
128     physEndRun->TB_ENDRUN.TB_CALIB_FLAG = ((UINT8)subData[offset+16])>>7;
129 kusanagi 4.1 physEndRun->TB_ENDRUN.TB_CALO_TRIG = ((UINT8)subData[offset+16])>>6;
130     physEndRun->TB_ENDRUN.TB_S4_TRIG = ((UINT8)subData[offset+16])>>5;
131     physEndRun->TB_ENDRUN.TB_TOF_TRIG = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
132 kusanagi 2.3 //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
133 kusanagi 1.1 }
134    

  ViewVC Help
Powered by ViewVC 1.1.23