BMTFUnpackerOutput.cc

Go to the documentation of this file.

#include "EventFilter/L1TRawToDigi/plugins/UnpackerFactory.h"
#include "L1Trigger/L1TMuon/interface/RegionalMuonRawDigiTranslator.h"
#include "BMTFUnpackerOutput.h"

namespace l1t {
  namespace stage2 {

    bool BMTFUnpackerOutput::unpack(const Block &block, UnpackerCollections *coll) {
      unsigned int blockId = block.header().getID();
      LogDebug("L1T") << "Block ID: " << blockId << " size: " << block.header().getSize();

      //ZeroSupression Handler
      BxBlocks bxBlocks;
      bool ZS_enabled =
          (bool)((block.header().getFlags() >> 1) & 0x01);  //getFlags() returns first 8-bits from the amc header
      if (ZS_enabled)
        bxBlocks =
            block.getBxBlocks((unsigned int)6, true);  //it returnes 7-32bit bxBlocks originated from the amc13 Block
      else
        bxBlocks =
            block.getBxBlocks((unsigned int)6, false);  //it returnes 6-32bit bxBlocks originated from the amc13 Block

      edm::LogInfo("L1T") << "Will use the setup:"
                          << " ZS_enabled->" << ZS_enabled << " isTriggeringAlgo->" << isTriggeringAlgo << " isKalman->"
                          << isKalman;

      RegionalMuonCandBxCollection *res;
      if (isTriggeringAlgo)
        res = static_cast<BMTFCollections *>(coll)->getBMTFMuons();
      else
        res = static_cast<BMTFCollections *>(coll)->getBMTF2Muons();

      //BxBlocks changed the format of the blocks
      int firstBX = 0, lastBX = 0;
      int nBX = 0;
      if (!bxBlocks.empty()) {
        nBX = bxBlocks[0].header().getTotalBx();  //how many BX included in the BxBlock before Suppression
        getBXRange(nBX, firstBX, lastBX);
        res->setBXRange(firstBX, lastBX);
      } else {
        res->setBXRange(-2, 2);
        LogDebug("L1T") << "No BXs included in the given Block. Set the BXRange to be (-2,2).";
        return true;
      }

      LogDebug("L1T") << "nBX = " << nBX << " firstBX = " << firstBX << " lastBX = " << lastBX;

      int processor = block.amc().getBoardID() - 1;
      if (processor < 0 || processor > 11) {
        edm::LogInfo("L1T") << "Processor found out of range, it will be calculated by the old way";
        if (block.amc().getAMCNumber() % 2 != 0)
          processor = block.amc().getAMCNumber() / 2;
        else
          processor = 6 + (block.amc().getAMCNumber() / 2 - 1);
      }

      for (const auto &bxBlock : bxBlocks) {
        int ibx = bxBlock.header().getBx();

        for (auto iw = 0; iw < 6; iw += 2) {
          uint32_t raw_first = bxBlock.payload()[iw];      //payload[ip+(ibx+lastBX)*6];
          uint32_t raw_secnd = bxBlock.payload()[iw + 1];  //payload[ip+(ibx+lastBX)*6];
          if (raw_first == 0) {
            LogDebug("L1T") << "Raw data is zero";
            continue;
          }

          RegionalMuonCand muCand;
          RegionalMuonRawDigiTranslator::fillRegionalMuonCand(
              muCand, raw_first, raw_secnd, processor, tftype::bmtf, isKalman, false, false);

          if (muCand.hwPt() == 0) {
            continue;
          }

          if (isKalman) {
            LogDebug("L1T") << "Pt = " << muCand.hwPt() << " eta: " << muCand.hwEta() << " phi: " << muCand.hwPhi()
                            << " diplacedPt = " << muCand.hwPtUnconstrained();
          } else {
            LogDebug("L1T") << "Pt = " << muCand.hwPt() << " eta: " << muCand.hwEta() << " phi: " << muCand.hwPhi();
          }

          res->push_back(ibx, muCand);

        }  //for iw
      }    //for ibx

      return true;
    }  //unpack

  }  // namespace stage2
}  // namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage2::BMTFUnpackerOutput);