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
            
      RegionalMuonCandBxCollection *res;
      if (isKalman)
    res = static_cast<BMTFCollections*>(coll)->getkBMTFMuons();
      else
    res = static_cast<BMTFCollections*>(coll)->getBMTFMuons();

      //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:stage2::BMTFUnpackerOutput::unpack") << "Processor found out of range so 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 (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);

      if (muCand.hwQual() == 0 && !isKalman)
        continue;//though away muons with Zero-Quality (ONLY BMTF)

      muCand.setLink(48 + processor);   //the link corresponds to the uGMT input
      if (isKalman) {
        muCand.setHwPt2((raw_secnd >> 23) & 0xFF);
        muCand.setHwDXY((raw_secnd >> 2) & 0x3);
        LogDebug("L1T") << "Pt = " << muCand.hwPt() << " eta: " << muCand.hwEta() << " phi: " << muCand.hwPhi() << " diplacedPt = " << muCand.hwPt2();
      }
      else
        LogDebug("L1T") << "Pt = " << muCand.hwPt() << " eta: " << muCand.hwEta() << " phi: " << muCand.hwPhi();

      res->push_back(ibx, muCand);

    }//for iw
      }//for ibx

      return true;
    }//unpack

  }//ns stage2
}//ns lt1

DEFINE_L1T_UNPACKER(l1t::stage2::BMTFUnpackerOutput);