HFRingPacker.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Event.h"
#include "EventFilter/L1TRawToDigi/plugins/PackerFactory.h"

#include "CaloTokens.h"
#include "HFRingPacker.h"

namespace l1t {
  namespace stage1 {
    Blocks HFRingPacker::pack(const edm::Event& event, const PackerTokens* toks) {
      edm::Handle<EtSumBxCollection> etSums;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getEtSumToken(), etSums);

      edm::Handle<CaloSpareBxCollection> calosparesHFBitCounts;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getCaloSpareHFBitCountsToken(), calosparesHFBitCounts);

      edm::Handle<CaloSpareBxCollection> calosparesHFRingSums;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getCaloSpareHFRingSumsToken(), calosparesHFRingSums);

      std::vector<uint32_t> load;

      for (int i = etSums->getFirstBX(); i <= etSums->getLastBX(); ++i) {
        int n = 0;

        int hfbitcount = 0;
        int hfringsum = 0;
        int htmissphi = 0;
        int htmiss = 0;

        int flaghtmiss = 0;

        for (auto j = etSums->begin(i); j != etSums->end(i) && n < 4; ++j, ++n) {
          if (j->getType() == l1t::EtSum::kMissingHt) {
            flaghtmiss = j->hwQual() & 0x1;
            htmiss = std::min(j->hwPt(), 0x7F);
            htmissphi = std::min(j->hwPhi(), 0x1F);
          }
        }

        n = 0;

        for (auto j = calosparesHFBitCounts->begin(i); j != calosparesHFBitCounts->end(i) && n < 2; ++j, ++n) {
          hfbitcount = std::min(j->hwPt(), 0xFFF);
        }

        n = 0;

        for (auto j = calosparesHFRingSums->begin(i); j != calosparesHFRingSums->end(i) && n < 2; ++j, ++n) {
          hfringsum = std::min(j->hwPt(), 0xFFF);
        }

        uint16_t object[4] = {0, 0, 0, 0};

        object[0] = hfbitcount | ((hfringsum & 0x7) << 12);
        object[1] = htmissphi | ((htmiss & 0x7F) << 5) | (flaghtmiss << 12) | (0x1 << 14);
        object[2] = ((hfringsum >> 3) & 0x1FF) | (0x1) << 10 | (0x1) << 12 | (0x1) << 14;
        object[3] = 0x1 | (0x1 << 2) | (0x1 << 4) | (0x1 << 6) | (0x1 << 8) | (0x1 << 10) | (0x1 << 12) | (0x1 << 14);

        uint32_t word0 = (object[0] & 0xFFFF) | ((object[1] & 0xFFFF) << 16);
        uint32_t word1 = (object[2] & 0xFFFF) | ((object[3] & 0xFFFF) << 16);

        word0 |= (1 << 31) | (1 << 15);
        word1 |= ((i == 0) << 31) | ((i == 0) << 15);

        load.push_back(word0);
        load.push_back(word1);
      }

      return {Block(7, load)};
    }
  }  // namespace stage1
}  // namespace l1t

DEFINE_L1T_PACKER(l1t::stage1::HFRingPacker);