l1t::stage1::legacy::EtSumUnpacker Class Reference

#include <LegacyEtSumUnpacker.h>

Inheritance diagram for l1t::stage1::legacy::EtSumUnpacker:

Public Member Functions
bool	unpack (const Block &block, UnpackerCollections *coll) override
Public Member Functions inherited from l1t::Unpacker
unsigned int	getAlgoVersion ()
void	setAlgoVersion (const unsigned int version)
	Unpacker ()
virtual	~Unpacker ()=default

Detailed Description

Definition at line 9 of file LegacyEtSumUnpacker.h.

Member Function Documentation

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

overridevirtual

Implements l1t::Unpacker.

Definition at line 11 of file LegacyEtSumUnpacker.cc.

References coll, DEFINE_L1T_UNPACKER, stringResolutionProvider_cfi::et, l1t::getBXRange(), l1t::BlockHeader::getID(), l1t::BlockHeader::getSize(), l1t::Block::header(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), mps_fire::i, createfilelist::int, l1t::EtSum::kMissingEt, l1t::EtSum::kTotalEt, l1t::EtSum::kTotalHt, LogDebug, RazorAnalyzer::met, l1t::Block::payload(), l1t::L1Candidate::setHwPhi(), l1t::L1Candidate::setHwPt(), l1t::L1Candidate::setHwQual(), l1t::EtSum::setType(), and globals_cff::x1.

         {

            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 res_ = static_cast<CaloCollections*>(coll)->getEtSums();
            res_->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 totet=candbit[0] & 0xFFF;
               int overflowtotet=(candbit[0]>>12) & 0x1;
               int etmiss=candbit[1] & 0xFFF;
               int overflowetmiss=(candbit[1]>>12) & 0x1;
               int totht=candbit[2] & 0xFFF;
               int overflowtotht=(candbit[2]>>12) & 0x1;
               int etmissphi=candbit[3] & 0x7F;

               l1t::EtSum et = l1t::EtSum();
               et.setHwPt(totet);
               et.setType(l1t::EtSum::kTotalEt);      
               int flagtotet=et.hwQual();
               flagtotet|= overflowtotet;
               et.setHwQual(flagtotet);       
               LogDebug("L1T") << "ET: pT " << et.hwPt()<<"is overflow "<<overflowtotet<<std::endl;
               res_->push_back(bx,et);

               l1t::EtSum ht = l1t::EtSum();
               ht.setHwPt(totht);
               ht.setType(l1t::EtSum::kTotalHt);       
               int flagtotht=ht.hwQual();
               flagtotht|= overflowtotht;
               ht.setHwQual(flagtotht);       
               LogDebug("L1T") << "HT: pT " << ht.hwPt()<<"is overflow "<<overflowtotht<<std::endl;
               res_->push_back(bx,ht);

               l1t::EtSum met = l1t::EtSum();
               met.setHwPt(etmiss);
               met.setHwPhi(etmissphi);
               met.setType(l1t::EtSum::kMissingEt);   
               int flagetmiss=met.hwQual();
               flagetmiss|= overflowetmiss;
               met.setHwQual(flagetmiss);       
               LogDebug("L1T") << "MET: pT " << met.hwPt()<<"is overflow "<<overflowetmiss<<std::endl;
               res_->push_back(bx,met);

            }

            return true;

         }