// Implementation of the PhysEndRunReader class. //new version... extern "C" { #include "CRC.h" } #include "ReaderAlgorithms.h" using namespace pamela::techmodel; /** * Constructor. */ PhysEndRunReader::PhysEndRunReader(void): TechmodelAlgorithm(PacketType::PhysEndRun, "TechmodelPhysEndRunReader") { physEndRun = new PhysEndRunEvent(); } /** * Initialize the algorithm with a special run. This will initialize the * event reader routines for all packet types. */ void PhysEndRunReader::Init(PamelaRun *run) { run->WriteSubPacket(this, &physEndRun, physEndRun->Class()); } /** * Unpack the PhysEndRun event from an input file. * void PhysEndRunReader::PKT_RunEvent(char* subData, long int length) throw (WrongCRCException_PKTUsed){ * */ void PhysEndRunReader::PKT_RunEvent(char* subData, long int length) throw (Exception){ std::stringstream oss; string msg; UINT16 subCRC; //CRC of the data UINT16 readCRC; //CRC read from the end of the subpacket // baco /gpfs/wizard/flight/production/preRawreader/05510005.pam if (length<21) { oss.str(""); oss<<"Wrong Lenght for PhysEndRun Packet: "<<" length = "<< length<<" too small to be real"; msg=oss.str(); PamOffLineSW::mainLogUtil->logError(msg); //the packet is not good and will be discarded throw FatalException(" Wrong Lenght for PhysEndRun Packet "); } long int dataLength = length - 2; long int tbInitPos = dataLength - 19; // 19 is the length of th TB subCRC = CM_Compute_CRC16(0, (UINT8*)subData, dataLength); readCRC = (((UINT16)(subData[length - 2]<<8))&0xFF00) + (((UINT16)subData[length - 1])&0x00FF); // std::cout<<"computedCRC "<< (UINT16)subCRC<<" readCRC "<< (UINT16)readCRC<<" readCRC_1 "<<(UINT16)subData[length-2]<<" readCRC_2 "<<(UINT16)subData[length-1]<CALO_ENDRUN. while(offset < tbInitPos){ physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0xFF); physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]); //subData[offset+2--->3] is a fixed word = 0x23 //If this bit is ? the subsequent CALO_HK are not valid // printf(" j %i calo_board_id_hk %X \n",j,physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK); // printf(" j %i calo_board_status_hk %X \n",j,physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK); if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) { UInt_t kk = 0; for(int k = 0; k < 11 ; k++){ // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]); physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF); physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF); // printf(" j %i k %i calo_hk0 %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_HK0[k]); // printf(" j %i k %i calo_hk1 %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_HK1[k]); kk += 2; } offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes) } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes) //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24] // for(int k = 0; k < 7 ; k++){ // // 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. // // printf(" INLOPP k %i subdata %X \n",k,(UINT8)(subData[offset])); physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k] = ((UINT8)(subData[offset])&0xFF); offset++; physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k] = (UINT8)(subData[offset]); offset++; // printf(" j %i k %i calo_board_id_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG[k]); // printf(" j %i k %i calo_board_status_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG[k]); offset++; //subData[offset+2--->3] is a fixed word = 0x2 that is the data lenght offset++; physEndRun->CALO_ENDRUN[j].CALO_REG[k] = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)subData[offset+1])&0x00FF); // printf(" j %i k %i calo_reg %X \n",j,k,physEndRun->CALO_ENDRUN[j].CALO_REG[k]); offset++; // CRC offset++; offset++; // offset++; // physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF); // physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0; //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7] // offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte) }; // j++; } // //--------------- CALO SECTION------------------------------------ // long int offset = 0; // int j = 0; // //physEndRun->CALO_ENDRUN. // while(offset < tbInitPos){ // physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_HK = (((UINT8)subData[offset])&0x1F); // physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK = ((UINT8)subData[offset+1]); // //subData[offset+2--->3] is a fixed word = 0x23 // //If this bit is ? the subsequent CALO_HK are not valid // if (!(physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_HK & 0x02)) { // UInt_t kk = 0; // for(int k = 0; k < 11 ; k++){ // // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (UINT8)(subData[offset+4+kk])+ 255 + (UINT8)(subData[offset+5+kk]); // physEndRun->CALO_ENDRUN[j].CALO_HK0[k] = (((UINT16)(subData[offset+4+kk]<<8))&0xFF00) + (((UINT16)subData[offset+5+kk])&0x00FF); // physEndRun->CALO_ENDRUN[j].CALO_HK1[k] = (((UINT16)(subData[offset+26+kk]<<8))&0xFF00) + (((UINT16)subData[offset+27+kk])&0x00FF); // kk += 2; // } // offset = offset + 50; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + 2*CALO_HKx + CRC (2Bytes) // } else offset = offset + 4; //CALO_BOARD_ID_HK + CALO_BOARD_STATUS_HK + CRC (2Bytes) // //subData[offset+(3+k)--->(15+k)] have a CRC(what?) in subData[offset+ 23--->24] // // for(int k = 0; k < 7 ; k++){ // physEndRun->CALO_ENDRUN[j].CALO_BOARD_ID_REG = ((UINT8)(subData[offset])&0x1F); // physEndRun->CALO_ENDRUN[j].CALO_BOARD_STATUS_REG = (UINT8)(subData[offset+1]); // //subData[offset+2--->3] is a fixed word = 0x2 // physEndRun->CALO_ENDRUN[j].CALO_REG[0] = (((UINT16)(subData[offset+4]<<8))&0xFF00) + (((UINT16)subData[offset+5])&0x00FF); // physEndRun->CALO_ENDRUN[j].CALO_REG[1] = (((UINT16)(subData[offset+6]<<8))&0xFF00) + (((UINT16)subData[offset+7])&0x00FF); // physEndRun->CALO_ENDRUN[j].CALO_REG[2] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[3] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[4] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[5] = 0; // physEndRun->CALO_ENDRUN[j].CALO_REG[6] = 0; // //subData[offset+ 0--->5] have a CRC(what?) in subData[offset+ 6--->7] // offset = offset + 56; //CALO_BOARD_ID_REG + CALO_BOARD_STATUS_REG + CALO_REG + CRC(2Byte) // // } // j++; // } //--------------- TB SECTION------------------------------------ // for (Int_t pi=-8; pi<40; pi++){ // printf(" subdata[%u] %X \n",offset+pi,(UINT8)subData[offset+pi]); // }; // UINT32 TB_ALARM_MASK; // UINT32 TB_PMT_MASK_S3S2S12; // UINT32 TB_PMT_MASK_S11CRC; // UINT16 TB_S4_CAL_MASK; // UINT32 TB_BUSY_MASK; // UINT32 TB_TRIG_CONF; physEndRun->TB_ENDRUN.TB_ALARM_MASK = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF); offset += 3; // 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); offset += 4; // physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF); offset += 3; // physEndRun->TB_ENDRUN.TB_S4_CAL_MASK = (((UINT16)(subData[offset]<<8))&0xFF00) + (((UINT16)(subData[offset + 1]))&0x00FF); offset += 2; // 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); offset += 4; // physEndRun->TB_ENDRUN.TB_TRIG_CONF = (((UINT32)(subData[offset]<<16))&0x00FF0000) + (((UINT32)(subData[offset + 1]<<8))&0x0000FF00) + (((UINT32)(subData[offset + 2]))&0x000000FF); offset += 3; // // printf(" physEndRun->TB_ENDRUN.TB_ALARM_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_ALARM_MASK); // printf(" physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12 %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_PMT_MASK_S3S2S12); // printf(" physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_PMT_MASK_S11CRC); // printf(" physEndRun->TB_ENDRUN.TB_S4_CAL_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_S4_CAL_MASK); // printf(" physEndRun->TB_ENDRUN.TB_BUSY_MASK %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_BUSY_MASK); // printf(" physEndRun->TB_ENDRUN.TB_TRIG_CONF %X \n",(UInt_t)physEndRun->TB_ENDRUN.TB_TRIG_CONF); // I want to use the packet so the exception is thrown here after all .. if (subCRC != readCRC) { oss.str(""); oss<<"Wrong CRC for PhysEndRun Packet: "<<" CRC COMPUTED= "<< subCRC<<" CRC READ= "<< readCRC; msg=oss.str(); PamOffLineSW::mainLogUtil->logWarning(msg); //the packet is not good but used throw WrongCRCException_PKTUsed(" Wrong CRC for PhysEndRun Packet "); } } /** * Get a string with the version info of the algorithm. */ std::string PhysEndRunReader::GetVersionInfo(void) const { return "$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"; }