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

Contents of /yoda/techmodel/PhysEndRunReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 6.4 - (show annotations) (download)
Wed Aug 23 14:39:44 2006 UTC (18 years, 3 months ago) by pam-rm2
Branch: MAIN
CVS Tags: YODA6_3/10, YODA6_3/07, YODA6_3/08, YODA6_3/09
Changes since 6.3: +3 -4 lines
YODA6_3/005 + S4/ND upgrade

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 static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.PhysEndRunReader"));
12
13 /**
14 * Constructor.
15 */
16 PhysEndRunReader::PhysEndRunReader(void):
17 TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") {
18 logger->debug(_T("Constructor"));
19 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 "$Header: /afs/ba.infn.it/user/pamela/src/CVS/yoda/techmodel/PhysEndRunReader.cpp,v 6.2 2006/05/30 19:10:03 kusanagi Exp $\n";
28 }
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 logger->debug(_T("Initialize"));
38 }
39
40 /**
41 * Unpack the PhysEndRun event from an input file.
42 */
43
44 void PhysEndRunReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
45 char subData[length];
46 UINT16 subCRC; //CRC of the data
47 UINT16 readCRC; //CRC read from the end of the subpacket
48 long int dataLength = length - 2;
49 long int tbInitPos = dataLength - 19; // 19 is the length of th TB
50
51 //physEndRun->CALO_ENDRUN[];
52 //physEndRun->TB_ENDRUN;
53
54 InputFile->read(subData, sizeof(subData));
55 subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength);
56 readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF);
57
58 //if (subCRC != readCRC) throw WrongCRCException(" Wrong CRC for PhysEndRun Packet ");
59 if (subCRC != readCRC) {
60 logger->error(" Wrong CRC for PhysEndRun Packet but I process it anyway");
61 }
62 //--------------- 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 for(int k = 0; k < 11 ; k++){
73 physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+3+k]<<8))&0xFF00) + (((UINT16)subData[offset+4+k])&0x00FF);
74 physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+14+k]<<8))&0xFF00) + (((UINT16)subData[offset+15+k])&0x00FF);
75 }
76 offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes)
77 } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes)
78 //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24]
79
80 for(int k = 0; k < 7 ; k++){
81 physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F);
82 physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]);
83 //subData[offset+2--->3] is a fixed word = 0x2
84 physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF);
85 //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7]
86 offset = offset + 8; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte)
87 }
88 j++;
89 }
90
91 //--------------- TB SECTION------------------------------------
92 physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset + 1])&0x00FF);
93 //subData[offset+ 0--->1] have a CRC(what?) in subData[offset+2]
94 physEndRun->TB_ENDRUN.TB_PMT_MASK_S3 = ((((UINT16)(subData[offset+2]))<<4)&0x0FF0) + (((UINT16)subData[offset + 3]>>4)&0x000F);
95 physEndRun->TB_ENDRUN.TB_PMT_MASK_S2 = ((((UINT8)(subData[offset+3]))<<4)&0xF0) + (((UINT8)subData[offset + 4]>>4)&0x0F);
96 physEndRun->TB_ENDRUN.TB_PMT_MASK_S12 = ((((UINT16)(subData[offset+4]))<<8)&0x0F00) + (((UINT16)subData[offset + 5])&0x00FF);
97 physEndRun->TB_ENDRUN.TB_PMT_MASK_S11 = ((((UINT16)(subData[offset+6]))<<8)&0xFF00) + (((UINT16)subData[offset + 7])&0x00FF);
98 //subData[offset+ 3--->8] have a CRC(what?) in subData[offset+9]
99 physEndRun->TB_ENDRUN.TB_CALO_MASK = ((UINT8)subData[offset+10])>>5;
100 physEndRun->TB_ENDRUN.TB_S4_MASK = (((UINT8)subData[offset+10])>>1)&0x0F;
101 //subData[tbInitPos+9] have a CRC(what?) in subData[tbInitPos+11]
102 physEndRun->TB_ENDRUN.TB_BUSY_MASK = (((UINT32)subData[offset+12]<<16)&0x00FF0000) + (((UINT32)subData[offset+13]<<8)&0x0000FF00) + (((UINT32)subData[offset+14])&0x000000FF);
103 //subData[tbInitPos+ 12--->14] have a CRC(what?) in subData[tbInitPos+15]
104 physEndRun->TB_ENDRUN.TB_CALIB_FLAG = ((UINT8)subData[offset+16])>>7;
105 physEndRun->TB_ENDRUN.TB_CALO_TRIG = ((UINT8)subData[offset+16])>>6;
106 physEndRun->TB_ENDRUN.TB_S4_TRIG = ((UINT8)subData[offset+16])>>5;
107 physEndRun->TB_ENDRUN.TB_TOF_TRIG = (UINT8)(subData[offset+16]<<3) + (UINT8)(subData[offset+17]>>6);
108 //subData[tbInitPos+ 16--->17] have a CRC(what?) in subData[tbInitPos+18]
109 }
110

  ViewVC Help
Powered by ViewVC 1.1.23