HFRingUnpacker.cc

Go to the documentation of this file.

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

#include "CaloCollections.h"
#include "HFRingUnpacker.h"

namespace l1t {
  namespace stage1 {
    bool HFRingUnpacker::unpack(const Block& block, UnpackerCollections* coll) {
      LogDebug("L1T") << "Block ID  = " << block.header().getID() << " size = " << block.header().getSize();

      int nBX = int(ceil(block.header().getSize() / 2.));

      // Find the first and last BXs
      int firstBX = -(ceil((double)nBX / 2.) - 1);
      int lastBX;
      if (nBX % 2 == 0) {
        lastBX = ceil((double)nBX / 2.) + 1;
      } else {
        lastBX = ceil((double)nBX / 2.);
      }

      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;
    }
  }  // namespace stage1
}  // namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage1::HFRingUnpacker);