LegacyEtSumUnpacker.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "EventFilter/L1TRawToDigi/plugins/UnpackerFactory.h"

#include "CaloCollections.h"
#include "LegacyEtSumUnpacker.h"

namespace l1t {
  namespace stage1 {
    namespace legacy {
      bool EtSumUnpacker::unpack(const Block& block, UnpackerCollections* coll) {
        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::kTotalEt};
          et.setHwPt(totet);
          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::kTotalHt};
          ht.setHwPt(totht);
          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::kMissingEt};
          met.setHwPt(etmiss);
          met.setHwPhi(etmissphi);
          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;
      }
    }  // namespace legacy
  }  // namespace stage1
}  // namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage1::legacy::EtSumUnpacker);