l1t::stage1::legacy::HFRingUnpacker Class Reference

Inheritance diagram for l1t::stage1::legacy::HFRingUnpacker:

Public Member Functions
virtual bool	unpack (const Block &block, UnpackerCollections *coll) override

Detailed Description

Definition at line 10 of file LegacyHFRingUnpacker.cc.

Member Function Documentation

bool l1t::stage1::legacy::HFRingUnpacker::unpack ( const Block &  block, UnpackerCollections *  coll  )

overridevirtual

Implements l1t::Unpacker.

Definition at line 24 of file LegacyHFRingUnpacker.cc.

References coll, l1t::getBXRange(), l1t::BlockHeader::getID(), l1t::BlockHeader::getSize(), l1t::Block::header(), l1t::CaloSpare::HFBitCount, l1t::CaloSpare::HFRingSum, l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), i, l1t::EtSum::kMissingHt, LogDebug, l1t::Block::payload(), l1t::L1Candidate::setHwPhi(), l1t::L1Candidate::setHwPt(), l1t::L1Candidate::setHwQual(), l1t::CaloSpare::setType(), and l1t::EtSum::setType().

         {

            LogDebug("L1T") << "Block ID  = " << block.header().getID() << " size = " << block.header().getSize();

            int nBX, firstBX, lastBX;
            nBX = int(ceil(block.header().getSize() / 2.)); 
            getBXRange(nBX, firstBX, lastBX);

            auto resHFBitCounts_ = static_cast<CaloCollections*>(coll)->getCaloSpareHFBitCounts();
            resHFBitCounts_->setBXRange(firstBX, lastBX);

            auto resHFRingSums_ = static_cast<CaloCollections*>(coll)->getCaloSpareHFRingSums();
            resHFRingSums_->setBXRange(firstBX, lastBX);

            auto reset_ = static_cast<CaloCollections*>(coll)->getEtSums();
            reset_->setBXRange(firstBX, lastBX);

            LogDebug("L1T") << "nBX = " << nBX << " first BX = " << firstBX << " lastBX = " << lastBX;

            // Initialise index
            int unsigned i = 0;

            // Loop over multiple BX and then number of jets filling jet collection
            for (int bx=firstBX; bx<=lastBX; bx++){
               uint32_t raw_data0 = block.payload()[i++];
               uint32_t raw_data1 = block.payload()[i++];        

               /* if (raw_data0 == 0 || raw_data1==0) continue; */

               uint16_t candbit[4];
               candbit[0] = raw_data0 & 0xFFFF;
               candbit[1] = (raw_data0 >> 16) & 0xFFFF;
               candbit[2] = raw_data1 & 0xFFFF;
               candbit[3] = (raw_data1 >> 16) & 0xFFFF;

               int hfbitcount=candbit[0] & 0xFFF;
               int hfringsum=((candbit[0]>>12) & 0x7) | ((candbit[2] & 0x1FF) << 3);
               int htmissphi=candbit[1] & 0x1F;
               int htmiss=(candbit[1]>>5) & 0x7F;
               int overflowhtmiss=(candbit[1]>>12) & 0x1;

               l1t::CaloSpare hfbc= l1t::CaloSpare();
               hfbc.setHwPt(hfbitcount);
               hfbc.setType(l1t::CaloSpare::HFBitCount);  
               LogDebug("L1T") << "hfbc pT " << hfbc.hwPt(); 
               resHFBitCounts_->push_back(bx,hfbc);        

               l1t::CaloSpare hfrs= l1t::CaloSpare();
               hfrs.setHwPt(hfringsum);
               hfrs.setType(l1t::CaloSpare::HFRingSum);  
               LogDebug("L1T") << "hfrs pT " << hfrs.hwPt();  
               resHFRingSums_->push_back(bx,hfrs);       

               l1t::EtSum mht = l1t::EtSum();
               mht.setHwPt(htmiss);
               mht.setHwPhi(htmissphi);
               mht.setType(l1t::EtSum::kMissingHt); 
               int flaghtmiss=mht.hwQual();
               flaghtmiss|= overflowhtmiss;
               mht.setHwQual(flaghtmiss);       
               LogDebug("L1T") << "MHT: pT " << mht.hwPt()<<"is overflow "<<overflowhtmiss<<std::endl;
               reset_->push_back(bx,mht);       

            }

            return true;

         }