l1t::stage1::EtSumUnpacker Class Reference

Inheritance diagram for l1t::stage1::EtSumUnpacker:

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

Detailed Description

Definition at line 9 of file EtSumUnpacker.cc.

Member Function Documentation

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

overridevirtual

Implements l1t::Unpacker.

Definition at line 21 of file EtSumUnpacker.cc.

References coll, l1t::getBXRange(), l1t::BlockHeader::getID(), l1t::BlockHeader::getSize(), l1t::Block::header(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), i, l1t::EtSum::kTotalEt, l1t::EtSum::kTotalHt, LogDebug, l1t::Block::payload(), l1t::L1Candidate::setHwPt(), l1t::L1Candidate::setHwQual(), 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 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++){

          res_->resize(bx,4);

          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[2];
          candbit[0] = raw_data0 & 0xFFFF;
          candbit[1] = raw_data1 & 0xFFFF;

          int totet=candbit[0] & 0xFFF;
          int overflowtotet=(candbit[0]>>12) & 0x1;
          int totht=candbit[1] & 0xFFF;
          int overflowtotht=(candbit[1]>>12) & 0x1;

          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_->set(bx, 2,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_->set(bx, 3,ht);

        }

        return true;

      }