| 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.4 2006/08/23 14:39:44 pam-rm2 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 |
memset(subData, 0, length*sizeof(char)); |
| 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 |
|