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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.3 - (hide annotations) (download)
Wed Sep 24 15:17:23 2008 UTC (16 years, 2 months ago) by mocchiut
Branch: MAIN
Changes since 1.2: +47 -29 lines
PhysEndRun trigger data unpack upgrade

1 mocchiut 1.2 // 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 mocchiut 1.3 physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0xFF);
65 mocchiut 1.2 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 mocchiut 1.3 // 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 mocchiut 1.2 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 mocchiut 1.3 // 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 mocchiut 1.2 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 mocchiut 1.3 physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k] = ((UINT8)(subData[offset])&0xFF);
90 mocchiut 1.2 offset++;
91     physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k] = (UINT8)(subData[offset]);
92     offset++;
93 mocchiut 1.3 // 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 mocchiut 1.2 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 mocchiut 1.3 // printf(" j %i k %i calo_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_REG[k]);
100 mocchiut 1.2 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 mocchiut 1.3 // for (Int_t pi=-8; pi<40; pi++){
159     // printf(" subdata[%u] %X \n",offset+pi,(UINT8)subData[offset+pi]);
160     // };
161    
162     // UINT32 TB_ALARM_MASK;
163     // UINT32 TB_PMT_MASK_S3S2S12;
164     // UINT32 TB_PMT_MASK_S11CRC;
165     // UINT16 TB_S4_CAL_MASK;
166     // UINT32 TB_BUSY_MASK;
167     // UINT32 TB_TRIG_CONF;
168    
169     physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF);
170     offset += 3;
171     //
172     physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12 = (((UINT32)(subData[offset]<<24))&0xFF000000) + (((UINT32)(subData[offset + 1]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 2]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 3]))&0x000000FF);
173     offset += 4;
174     //
175     physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF);
176     offset += 3;
177     //
178     physEndRun->TB_ENDRUN.TB_S4_CAL_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)(subData[offset + 1]))&0x00FF);
179     offset += 2;
180     //
181     physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)(subData[offset]<<24))&0xFF000000) + (((UINT32)(subData[offset + 1]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 2]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 3]))&0x000000FF);
182     offset += 4;
183     //
184     physEndRun->TB_ENDRUN.TB_TRIG_CONF = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF);
185     offset += 3;
186     //
187    
188     // printf(" physEndRun->TB_ENDRUN.TB_ALARM_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_ALARM_MASK);
189     // printf(" physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12 %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12);
190     // printf(" physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC);
191     // printf(" physEndRun->TB_ENDRUN.TB_S4_CAL_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_S4_CAL_MASK);
192     // printf(" physEndRun->TB_ENDRUN.TB_BUSY_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_BUSY_MASK);
193     // printf(" physEndRun->TB_ENDRUN.TB_TRIG_CONF %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_TRIG_CONF);
194    
195 mocchiut 1.2 // I want to use the packet so the exception is thrown here after all ..
196     if (subCRC != readCRC)
197     {
198     oss.str("");
199     oss<<"Wrong CRC for PhysEndRun Packet: "<<" CRC COMPUTED= "<< subCRC<<" CRC READ= "<< readCRC;
200     msg=oss.str();
201     PamOffLineSW::mainLogUtil->logWarning(msg);
202     //the packet is not good but used
203     throw WrongCRCException_PKTUsed(" Wrong CRC for PhysEndRun Packet ");
204     }
205     }
206    
207     /**
208     * Get a string with the version info of the algorithm.
209     */
210     std::string PhysEndRunReader::GetVersionInfo(void) const {
211     return
212 mocchiut 1.3 "$Header: /afs/ba.infn.it/user/pamela/src/CVS/chewbacca/PamOffLineSW/techmodel/PhysEndRunReader.cpp,v 1.2 2008/09/24 08:47:53 mocchiut Exp $\n";
213 mocchiut 1.2 }

  ViewVC Help
Powered by ViewVC 1.1.23