RCTEmRegionUnpacker.cc

Go to the documentation of this file.

#include "PhysicsToBitConverter.h"
#include "rctDataBase.h"

#include "DataFormats/L1CaloTrigger/interface/L1CaloEmCand.h"
#include "DataFormats/L1TCalorimeter/interface/CaloEmCand.h"
#include "DataFormats/L1CaloTrigger/interface/L1CaloRegion.h"
#include "DataFormats/L1TCalorimeter/interface/CaloRegion.h"
#include "DataFormats/L1CaloTrigger/interface/L1CaloRegionDetId.h"
#include "DataFormats/L1CaloTrigger/interface/L1CaloCollections.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "EventFilter/L1TRawToDigi/plugins/UnpackerFactory.h"

#include <iostream>
#include <fstream>

#include "EventFilter/L1TRawToDigi/interface/Block.h"

#include "CaloCollections.h"
#include "RCTEmRegionUnpacker.h"

namespace l1t {
  namespace stage1 {
    void unpack_em(const Block& block, UnpackerCollections* coll) {
      int nBX, firstBX, lastBX;
      nBX = int(ceil(block.header().getSize() / 6.));
      getBXRange(nBX, firstBX, lastBX);

      auto resRCTEMCands_ = static_cast<CaloCollections*>(coll)->getCaloEmCands();

      // Initialise index
      int unsigned i = 0;

      for (int bx = firstBX; bx <= lastBX; bx++) {
        unsigned int crate;
        bool even = false;

        std::vector<uint32_t> uint;
        uint.reserve(6);

        PhysicsToBitConverter converter;
        rctDataBase database;
        int mp7link = (int)(block.header().getID() / 2);
        database.GetLinkRCT(mp7link, crate, even);

        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);

        LogDebug("L1T") << "--------------- mp7 link =" << mp7link << "RCT crate id=" << crate
                        << ", RCT crate even=" << even << std::endl;

        if (!even) {
          for (int i = 0; i < 6; i++)
            converter.Set32bitWordLinkOdd(i, uint[i]);
          converter.Convert();

          for (int j = 0; j < 4; j++) {
            unsigned int rank = (unsigned int)converter.GetNEEt(j);
            unsigned int reg = (unsigned int)converter.GetNEReg(j);
            unsigned int card = (unsigned int)converter.GetNECard(j);

            LogDebug("L1T") << "UNPACKER, CRATE" << crate << "NON ISO em rank=" << rank << ", region=" << reg
                            << ", card=" << card << std::endl;

            L1CaloEmCand em = L1CaloEmCand(rank, reg, card, crate, false, j, bx);
            resRCTEMCands_->push_back(em);
          }

          for (int j = 0; j < 4; j++) {
            unsigned int rank = converter.GetIEEt(j);
            unsigned int reg = converter.GetIEReg(j);
            unsigned int card = converter.GetIECard(j);

            LogDebug("L1T") << "UNPACKER, CRATE" << crate << "ISO em rank=" << rank << ", region=" << reg
                            << ", card=" << card << std::endl;
            L1CaloEmCand em = L1CaloEmCand(rank, reg, card, crate, true, j, bx);
            resRCTEMCands_->push_back(em);
          }
        }  // end if odd
      }  // end of loop over BX
    }

    void unpack_region(const Block& block, UnpackerCollections* coll) {
      int nBX, firstBX, lastBX;
      nBX = int(ceil(block.header().getSize() / 6.));
      getBXRange(nBX, firstBX, lastBX);

      auto resRCTRegions_ = static_cast<CaloCollections*>(coll)->getCaloRegions();

      // Initialise index
      int unsigned i = 0;
      std::vector<uint32_t> uint;
      uint.reserve(6);

      for (int bx = firstBX; bx <= lastBX; bx++) {
        unsigned int crate;
        bool even = false;

        PhysicsToBitConverter converter;
        rctDataBase database;
        int mp7link = (int)(block.header().getID() / 2);
        database.GetLinkRCT(mp7link, crate, even);

        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);
        uint.push_back(block.payload()[i++]);

        LogDebug("L1T") << "--------------- mp7 link =" << mp7link << "RCT crate id=" << crate
                        << ", RCT crate even=" << even << std::endl;

        if (!even) {
          for (int i = 0; i < 6; i++)
            converter.Set32bitWordLinkOdd(i, uint[i]);
          converter.Convert();

          for (int j = 0; j < 8; j++) {
            unsigned int hfet = (unsigned int)converter.GetHFEt(j);
            bool hfgrain = (bool)converter.GetHFFg(j);

            LogDebug("L1T") << "UNPACKER, CRATE" << crate << "region=" << j << ", rgnEt=" << hfet << std::endl;
            L1CaloRegion rgn = L1CaloRegion(hfet, hfgrain, crate, j);
            rgn.setBx(bx);
            resRCTRegions_->push_back(rgn);
          }
        }  // end if odd

        else {
          for (int i = 0; i < 6; i++)
            converter.Set32bitWordLinkEven(i, uint[i]);
          converter.Convert();

          for (int j = 0; j < 7; j++) {
            for (int k = 0; k < 2; k++) {
              unsigned int RCet = (unsigned int)converter.GetRCEt(j, k);
              bool overflow = (bool)converter.GetRCOf(j, k);
              bool tauveto = (bool)converter.GetRCTau(j, k);
              bool hadveto = (bool)converter.GetRCHad(j, k);
              bool quiet = false;  //to be fixed

              LogDebug("L1T") << "UNPACKER, CRATE=" << crate << ",region=" << k << ", card=" << j << ", rgnEt=" << RCet
                              << ", overflow=" << overflow << ", tauveto=" << tauveto << ", hadveto=" << hadveto
                              << std::endl;
              L1CaloRegion rgn = L1CaloRegion(RCet, overflow, tauveto, hadveto, quiet, crate, j, k);
              rgn.setBx(bx);
              resRCTRegions_->push_back(rgn);
            }
          }
        }  // end if even
      }  // end of loop over BX
    }

    bool RCTEmRegionUnpacker::unpack(const Block& block, UnpackerCollections* coll) {
      if (block.header().getCapID() == 0) {
        unpack_region(block, coll);
      } else if (block.header().getCapID() == 1) {
        unpack_em(block, coll);
      } else {
        return false;
      }
      return true;
    }
  }  // namespace stage1
}  // namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage1::RCTEmRegionUnpacker);