/[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.1 - (hide annotations) (download)
Tue Sep 23 07:20:24 2008 UTC (16 years, 2 months ago) by mocchiut
Branch: MAIN
Branch point for: v0r00
Initial revision

1 mocchiut 1.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     //--------------- CALO SECTION------------------------------------
59     long int offset = 0;
60     int j = 0;
61     //physEndRun->CALO_ENDRUN.
62     while(offset < tbInitPos){
63     physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F);
64     physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]);
65     //subData[offset+2--->3] is a fixed word = 0x23
66     //If this bit is ? the subsequent CALO_HK are not valid
67     if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) {
68     UInt_t kk = 0;
69     for(int k = 0; k < 11 ; k++){
70     // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]);
71     physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF);
72     physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF);
73     kk += 2;
74     }
75     offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
76     } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
77     //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
78    
79     // for(int k = 0; k < 7 ; k++){
80     physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F);
81     physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
82     //subData[offset+2--->3] is a fixed word = 0x2
83     physEndRun->CALO_ENDRUN[j].CALO_REG[0] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
84     physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF);
85     physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0;
86     physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0;
87     physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0;
88     physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0;
89     physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0;
90     //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
91     offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
92     // }
93     j++;
94     }
95    
96     //--------------- TB SECTION------------------------------------
97     physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
98     //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
99     physEndRun->TB_ENDRUN.TB_PMT_MASK_S3 = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
100     physEndRun->TB_ENDRUN.TB_PMT_MASK_S2 = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
101     physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
102     physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
103     //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
104     physEndRun->TB_ENDRUN.TB_CALO_MASK = ((UINT8)subData[offset+10])>>5;
105     physEndRun->TB_ENDRUN.TB_S4_MASK = (((UINT8)subData[offset+10])>>1)&0x0F;
106     //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
107     physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)subData[offset+12]<<16)&0x00FF0000) + (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
108     //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
109     physEndRun->TB_ENDRUN.TB_CALIB_FLAG = ((UINT8)subData[offset+16])>>7;
110     physEndRun->TB_ENDRUN.TB_CALO_TRIG = ((UINT8)subData[offset+16])>>6;
111     physEndRun->TB_ENDRUN.TB_S4_TRIG = ((UINT8)subData[offset+16])>>5;
112     physEndRun->TB_ENDRUN.TB_TOF_TRIG = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
113     //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
114    
115     // I want to use the packet so the exception is thrown here after all ..
116     if (subCRC != readCRC)
117     {
118     oss.str("");
119     oss<<"Wrong CRC for PhysEndRun Packet: "<<" CRC COMPUTED= "<< subCRC<<" CRC READ= "<< readCRC;
120     msg=oss.str();
121     PamOffLineSW::mainLogUtil->logWarning(msg);
122     //the packet is not good but used
123     throw WrongCRCException_PKTUsed(" Wrong CRC for PhysEndRun Packet ");
124     }
125     }
126    
127     /**
128     * Get a string with the version info of the algorithm.
129     */
130     std::string PhysEndRunReader::GetVersionInfo(void) const {
131     return
132     "$Header: /repository/PamOffLineSW/techmodel/PhysEndRunReader.cpp,v 1.12 2008-03-05 16:52:18 messineo Exp $\n";
133     }

  ViewVC Help
Powered by ViewVC 1.1.23