EtSumUnpacker_0x10010057.cc

Go to the documentation of this file.

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

#include "L1Trigger/L1TCalorimeter/interface/CaloTools.h"

#include "L1TObjectCollections.h"

#include "L1TStage2Layer2Constants.h"
#include "EtSumUnpacker_0x10010057.h"

namespace l1t {
  namespace stage2 {
    EtSumUnpacker_0x10010057::EtSumUnpacker_0x10010057() : EtSumCopy_(0) {}

    bool EtSumUnpacker_0x10010057::unpack(const Block& block, UnpackerCollections* coll) {
      using namespace l1t::stage2::layer2;

      LogDebug("L1T") << "Block ID  = " << block.header().getID() << " size = " << block.header().getSize();

      int nBX = int(
          ceil(block.header().getSize() / demux::nOutputFramePerBX));  // Since there 6 frames per demux output event
      // expect the first four frames to be the first 4 EtSum objects reported per event (see CMS IN-2013/005)

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

      auto res_ = static_cast<L1TObjectCollections*>(coll)->getEtSums(EtSumCopy_);
      res_->setBXRange(firstBX, lastBX);

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

      // Loop over multiple BX and fill EtSums collection
      for (int bx = firstBX; bx <= lastBX; bx++) {
        // ET
        int iFrame = (bx - firstBX) * demux::nOutputFramePerBX;

        uint32_t raw_data = block.payload().at(iFrame);

        l1t::EtSum et = l1t::EtSum();

        et.setHwPt(raw_data & 0xFFF);
        et.setType(l1t::EtSum::kTotalEt);
        et.setP4(l1t::CaloTools::p4Demux(&et));

        LogDebug("L1T") << "ET: pT " << et.hwPt() << " bx " << bx;

        res_->push_back(bx, et);

        // ET EM

        l1t::EtSum etem = l1t::EtSum();

        etem.setHwPt((raw_data >> 12) & 0xFFF);
        etem.setType(l1t::EtSum::kTotalEtEm);
        etem.setP4(l1t::CaloTools::p4Demux(&etem));

        LogDebug("L1T") << "ETEM: pT " << etem.hwPt() << " bx " << bx;

        res_->push_back(bx, etem);

        // MBHFPT0

        l1t::EtSum mbp0 = l1t::EtSum();
        mbp0.setHwPt((raw_data >> 28) & 0xf);
        mbp0.setType(l1t::EtSum::kMinBiasHFP0);

        res_->push_back(bx, mbp0);

        // HT

        raw_data = block.payload()[iFrame + 1];

        l1t::EtSum ht = l1t::EtSum();

        ht.setHwPt(raw_data & 0xFFF);
        ht.setType(l1t::EtSum::kTotalHt);
        ht.setP4(l1t::CaloTools::p4Demux(&ht));

        LogDebug("L1T") << "HT: pT " << ht.hwPt();

        res_->push_back(bx, ht);

        //MBHFMT0

        l1t::EtSum mbm0 = l1t::EtSum();
        mbm0.setHwPt((raw_data >> 28) & 0xf);
        mbm0.setType(l1t::EtSum::kMinBiasHFM0);

        res_->push_back(bx, mbm0);

        //  MET (no HF)

        raw_data = block.payload()[iFrame + 2];

        l1t::EtSum met = l1t::EtSum();

        met.setHwPt(raw_data & 0xFFF);
        met.setHwPhi((raw_data >> 12) & 0xFF);
        met.setType(l1t::EtSum::kMissingEt);
        met.setP4(l1t::CaloTools::p4Demux(&met));

        LogDebug("L1T") << "MET: phi " << met.hwPhi() << " pT " << met.hwPt() << " bx " << bx;

        res_->push_back(bx, met);

        // MBHFPT1

        l1t::EtSum mbp1 = l1t::EtSum();
        mbp1.setHwPt((raw_data >> 28) & 0xf);
        mbp1.setType(l1t::EtSum::kMinBiasHFP1);

        res_->push_back(bx, mbp1);

        // ET ASYMM

        l1t::EtSum etAsym = l1t::EtSum();
        etAsym.setHwPt((raw_data >> 20) & 0xFF);
        etAsym.setType(l1t::EtSum::kAsymEt);

        res_->push_back(bx, etAsym);

        // MHT

        raw_data = block.payload()[iFrame + 3];

        l1t::EtSum mht = l1t::EtSum();

        mht.setHwPt(raw_data & 0xFFF);
        mht.setHwPhi((raw_data >> 12) & 0xFF);
        mht.setType(l1t::EtSum::kMissingHt);
        mht.setP4(l1t::CaloTools::p4Demux(&mht));

        LogDebug("L1T") << "MHT: phi " << mht.hwPhi() << " pT " << mht.hwPt() << " bx " << bx;

        res_->push_back(bx, mht);

        // MBHFMT1

        l1t::EtSum mbm1 = l1t::EtSum();
        mbm1.setHwPt((raw_data >> 28) & 0xf);
        mbm1.setType(l1t::EtSum::kMinBiasHFM1);

        res_->push_back(bx, mbm1);

        // HT ASYMM

        l1t::EtSum htAsym = l1t::EtSum();
        htAsym.setHwPt((raw_data >> 20) & 0xFF);
        htAsym.setType(l1t::EtSum::kAsymHt);

        res_->push_back(bx, htAsym);

        //  MET (with HF)

        raw_data = block.payload()[iFrame + 4];

        l1t::EtSum methf = l1t::EtSum();

        methf.setHwPt(raw_data & 0xFFF);
        methf.setHwPhi((raw_data >> 12) & 0xFF);
        methf.setType(l1t::EtSum::kMissingEtHF);
        methf.setP4(l1t::CaloTools::p4Demux(&methf));

        LogDebug("L1T") << "METHF: phi " << methf.hwPhi() << " pT " << methf.hwPt() << " bx " << bx;

        res_->push_back(bx, methf);

        // ETHF ASYMM

        l1t::EtSum etHFAsym = l1t::EtSum();
        etHFAsym.setHwPt((raw_data >> 20) & 0xFF);
        etHFAsym.setType(l1t::EtSum::kAsymEtHF);

        res_->push_back(bx, etHFAsym);

        // CENTRALITY LOWER NIB

        uint32_t centLN = ((raw_data >> 28) & 0xF);

        // MHT with HF

        raw_data = block.payload()[iFrame + 5];

        l1t::EtSum mhthf = l1t::EtSum();

        mhthf.setHwPt(raw_data & 0xFFF);
        mhthf.setHwPhi((raw_data >> 12) & 0xFF);
        mhthf.setType(l1t::EtSum::kMissingHtHF);
        mhthf.setP4(l1t::CaloTools::p4Demux(&mhthf));

        LogDebug("L1T") << "MHThf: phi " << mhthf.hwPhi() << " pT " << mhthf.hwPt() << " bx " << bx;

        res_->push_back(bx, mhthf);

        // HTHF ASYMM

        l1t::EtSum htHFAsym = l1t::EtSum();
        htHFAsym.setHwPt((raw_data >> 20) & 0xFF);
        htHFAsym.setType(l1t::EtSum::kAsymHtHF);

        res_->push_back(bx, htHFAsym);

        // CENTRALITY UPPER NIB

        uint32_t centUN = ((raw_data >> 28) & 0xF);

        l1t::EtSum cent = l1t::EtSum();
        cent.setHwPt((centUN << 4) | centLN);
        cent.setType(l1t::EtSum::kCentrality);

        res_->push_back(bx, cent);

        //HI-SUM

        raw_data = block.payload()[iFrame + 1];

        l1t::EtSum towCount = l1t::EtSum();
        towCount.setHwPt((raw_data >> 12) & 0x1FFF);
        towCount.setType((l1t::EtSum::kTowerCount));
        towCount.setP4(l1t::CaloTools::p4Demux(&towCount));

        res_->push_back(bx, towCount);
      }

      return true;
    }
  }  // namespace stage2
}  // namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage2::EtSumUnpacker_0x10010057);