/[PAMELA software]/chewbacca/PamOffLineSW/techmodel/PhysEndRunReader.cpp
ViewVC logotype

Contents of /chewbacca/PamOffLineSW/techmodel/PhysEndRunReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (show annotations) (download)
Wed Sep 24 08:47:53 2008 UTC (16 years, 2 months ago) by mocchiut
Branch: MAIN
Changes since 1.1: +195 -133 lines
Calorimeter PhysEndRun updated + broken compilation fixed

1 // Implementation of the PhysEndRunReader class.
2 //new version...
3 extern "C" {
4 #include "CRC.h"
5 }
6
7 #include "ReaderAlgorithms.h"
8
9 using namespace pamela::techmodel;
10
11
12 /**
13 * Constructor.
14 */
15 PhysEndRunReader::PhysEndRunReader(void):
16 TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") {
17 physEndRun = new PhysEndRunEvent();
18 }
19
20 /**
21 * Initialize the algorithm with a special run. This will initialize the
22 * event reader routines for all packet types.
23 */
24 void PhysEndRunReader::Init(PamelaRun *run) {
25 run->WriteSubPacket(this, &physEndRun, physEndRun->Class());
26 }
27
28 /**
29 * Unpack the PhysEndRun event from an input file.
30 * void PhysEndRunReader::PKT_RunEvent(char* subData, long int length) throw (WrongCRCException_PKTUsed){
31 *
32 */
33 void PhysEndRunReader::PKT_RunEvent(char* subData, long int length) throw (Exception){
34 std::stringstream oss;
35 string msg;
36 UINT16 subCRC; //CRC of the data
37 UINT16 readCRC; //CRC read from the end of the subpacket
38
39
40 // baco /gpfs/wizard/flight/production/preRawreader/05510005.pam
41 if (length<21)
42 {
43 oss.str("");
44 oss<<"Wrong Lenght for PhysEndRun Packet: "<<" length = "<< length<<" too small to be real";
45 msg=oss.str();
46 PamOffLineSW::mainLogUtil->logError(msg);
47 //the packet is not good and will be discarded
48 throw FatalException(" Wrong Lenght for PhysEndRun Packet ");
49 }
50
51 long int dataLength = length - 2;
52 long int tbInitPos = dataLength - 19; // 19 is the length of th TB
53
54 subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
55 readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
56 // std::cout<<"computedCRC "<< (UINT16)subCRC<<" readCRC "<< (UINT16)readCRC<<" readCRC_1 "<<(UINT16)subData[length-2]<<" readCRC_2 "<<(UINT16)subData[length-1]<<endl;
57
58
59 //--------------- CALO SECTION------------------------------------
60 long int offset = 0;
61 int j = 0;
62 //physEndRun->CALO_ENDRUN.
63 while(offset < tbInitPos){
64 physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F);
65 physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
66 //subData[offset+2--->3] is a fixed word = 0x23
67 //If this bit is ? the subsequent CALO_HK are not valid
68 printf(" j %i calo_board_id_hk %X \n",j,physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK);
69 printf(" j %i calo_board_status_hk %X \n",j,physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK);
70 if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
71 UInt_t kk = 0;
72 for(int k = 0; k < 11 ; k++){
73 // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]);
74 physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF);
75 physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF);
76 printf(" j %i k %i calo_hk0 %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_HK0[k]);
77 printf(" j %i k %i calo_hk1 %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_HK1[k]);
78 kk += 2;
79 }
80 offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
81 } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
82 //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
83 //
84 for(int k = 0; k < 7 ; k++){
85 //
86 // 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.
87 //
88 // printf(" INLOPP k %i subdata %X \n",k,(UINT8)(subData[offset]));
89 physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k] = ((UINT8)(subData[offset])&0x1F);
90 offset++;
91 physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k] = (UINT8)(subData[offset]);
92 offset++;
93 printf(" j %i k %i calo_board_id_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k]);
94 printf(" j %i k %i calo_board_status_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k]);
95 offset++;
96 //subData[offset+2--->3] is a fixed word = 0x2 that is the data lenght
97 offset++;
98 physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset+1])&0x00FF);
99 printf(" j %i k %i calo_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_REG[k]);
100 offset++;
101 // CRC
102 offset++;
103 offset++;
104 //
105 offset++;
106 // physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF);
107 // physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0;
108 // physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0;
109 // physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0;
110 // physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0;
111 // physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0;
112 //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
113 // offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
114 };
115 //
116 j++;
117 }
118
119 // //--------------- CALO SECTION------------------------------------
120 // long int offset = 0;
121 // int j = 0;
122 // //physEndRun->CALO_ENDRUN.
123 // while(offset < tbInitPos){
124 // physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F);
125 // physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
126 // //subData[offset+2--->3] is a fixed word = 0x23
127 // //If this bit is ? the subsequent CALO_HK are not valid
128 // if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
129 // UInt_t kk = 0;
130 // for(int k = 0; k < 11 ; k++){
131 // // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]);
132 // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF);
133 // physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF);
134 // kk += 2;
135 // }
136 // offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
137 // } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
138 // //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
139
140 // // for(int k = 0; k < 7 ; k++){
141 // physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F);
142 // physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
143 // //subData[offset+2--->3] is a fixed word = 0x2
144 // physEndRun->CALO_ENDRUN[j].CALO_REG[0] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
145 // physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF);
146 // physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0;
147 // physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0;
148 // physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0;
149 // physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0;
150 // physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0;
151 // //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
152 // offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
153 // // }
154 // j++;
155 // }
156
157 //--------------- TB SECTION------------------------------------
158 printf(" subdata[%u] %X \n",offset,subData[offset]);
159 physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
160 //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
161 physEndRun->TB_ENDRUN.TB_PMT_MASK_S3 = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
162 physEndRun->TB_ENDRUN.TB_PMT_MASK_S2 = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
163 physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
164 physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
165 //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
166 physEndRun->TB_ENDRUN.TB_CALO_MASK = ((UINT8)subData[offset+10])>>5;
167 physEndRun->TB_ENDRUN.TB_S4_MASK = (((UINT8)subData[offset+10])>>1)&0x0F;
168 //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
169 physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)subData[offset+12]<<16)&0x00FF0000) + (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
170 //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
171 physEndRun->TB_ENDRUN.TB_CALIB_FLAG = ((UINT8)subData[offset+16])>>7;
172 physEndRun->TB_ENDRUN.TB_CALO_TRIG = ((UINT8)subData[offset+16])>>6;
173 physEndRun->TB_ENDRUN.TB_S4_TRIG = ((UINT8)subData[offset+16])>>5;
174 physEndRun->TB_ENDRUN.TB_TOF_TRIG = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
175 //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
176
177 // I want to use the packet so the exception is thrown here after all ..
178 if (subCRC != readCRC)
179 {
180 oss.str("");
181 oss<<"Wrong CRC for PhysEndRun Packet: "<<" CRC COMPUTED= "<< subCRC<<" CRC READ= "<< readCRC;
182 msg=oss.str();
183 PamOffLineSW::mainLogUtil->logWarning(msg);
184 //the packet is not good but used
185 throw WrongCRCException_PKTUsed(" Wrong CRC for PhysEndRun Packet ");
186 }
187 }
188
189 /**
190 * Get a string with the version info of the algorithm.
191 */
192 std::string PhysEndRunReader::GetVersionInfo(void) const {
193 return
194 "$Header: /afs/ba.infn.it/user/pamela/src/CVS/chewbacca/PamOffLineSW/techmodel/PhysEndRunReader.cpp,v 1.1.1.1 2008/09/23 07:20:24 mocchiut Exp $\n";
195 }

  ViewVC Help
Powered by ViewVC 1.1.23