TauPacker.cc

Go to the documentation of this file.

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

#include "CaloTokens.h"

#include "L1TStage2Layer2Constants.h"
#include "TauPacker.h"

namespace l1t {
  namespace stage2 {
    Blocks TauPacker::pack(const edm::Event& event, const PackerTokens* toks) {
      edm::Handle<TauBxCollection> taus;
      event.getByToken(static_cast<const CommonTokens*>(toks)->getTauToken(), taus);

      std::vector<uint32_t> load1, load2;

      int nBx = 0;
      for (int i = taus->getFirstBX(); i <= taus->getLastBX(); ++i) {
        for (auto j = taus->begin(i); j != taus->end(i); ++j) {
          uint32_t packed_eta = abs(j->hwEta()) & 0x7F;
          if (j->hwEta() < 0) {
            packed_eta = (128 - packed_eta) | 1 << 7;
          }

          uint32_t word = std::min(j->hwPt(), 0x1FF) | packed_eta << 9 | (j->hwPhi() & 0xFF) << 17 |
                          (j->hwIso() & 0x3) << 25 |  //two bits for tau isolation
                          (j->hwQual() & 0x7) << 27;

          // Depending on the number, push onto appropriate link
          if (load1.size() - nBx * l1t::stage2::layer2::demux::nTauPerLink < l1t::stage2::layer2::demux::nTauPerLink)
            load1.push_back(word);
          else
            load2.push_back(word);
        }

        // push zeroes if not enough taus must do this for each BX and each link
        while (load1.size() - nBx * l1t::stage2::layer2::demux::nOutputFramePerBX <
               l1t::stage2::layer2::demux::nOutputFramePerBX)
          load1.push_back(0);
        while (load2.size() - nBx * l1t::stage2::layer2::demux::nOutputFramePerBX <
               l1t::stage2::layer2::demux::nOutputFramePerBX)
          load2.push_back(0);

        nBx++;
      }

      return {Block(b1_, load1), Block(b2_, load2)};
    }
  }  // namespace stage2
}  // namespace l1t

DEFINE_L1T_PACKER(l1t::stage2::GTTauPacker);
DEFINE_L1T_PACKER(l1t::stage2::CaloTauPacker);