RCTEmRegionPacker.cc

Go to the documentation of this file.

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

#include "EventFilter/L1TRawToDigi/plugins/PackerFactory.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 "PhysicsToBitConverter.h"
#include "rctDataBase.h"

#include "CaloTokens.h"
#include "RCTEmRegionPacker.h"

namespace l1t {
  namespace stage1 {
    Blocks RCTEmRegionPacker::pack(const edm::Event& event, const PackerTokens* toks) {
      edm::Handle<L1CaloRegionCollection> caloregion;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getCaloRegionToken(), caloregion);

      edm::Handle<L1CaloEmCollection> caloemcand;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getCaloEmCandToken(), caloemcand);

      std::vector<uint32_t> load[36];

      for (int i = 0; i <= 0; ++i) {
        int n = 0;
        PhysicsToBitConverter converter[18];
        for (auto j = caloregion->begin(); j != caloregion->end(); ++j, ++n) {
          int et = (int)j->et();
          int overFlow = (int)j->overFlow();
          int fineGrain = (int)j->fineGrain();
          int mip = (int)j->mip();
          bool tauveto = j->tauVeto();
          //int quiet=(int)j->quiet();

          int rctCrate = (int)j->rctCrate();
          int rctCard = (int)j->rctCard();
          int rctRegion = (int)j->rctRegionIndex();
          bool isHf = (bool)j->isHf();

          if (!isHf) {
            converter[rctCrate].SetRCEt(et, rctCard, rctRegion);
            converter[rctCrate].SetRCOf(overFlow, rctCard, rctRegion);
            converter[rctCrate].SetRCTau(tauveto, rctCard, rctRegion);
            converter[rctCrate].SetRCHad(mip, rctCard, rctRegion);
            LogDebug("L1T") << "CRATE" << rctCrate << "region=" << rctRegion << ", card=" << rctCard << ", rgnEt=" << et
                            << ", overflow=" << overFlow << ", tauveto=" << tauveto << ", hadveto=" << mip << std::endl;

          } else {
            converter[rctCrate].SetHFEt(et, rctRegion);
            converter[rctCrate].SetHFFg(fineGrain, rctRegion);
            LogDebug("L1T") << "CRATE" << rctCrate << "region=" << rctRegion << ", rgnEt=" << et << std::endl;
          }
        }  //end calo region

        int m = 0;
        for (auto j = caloemcand->begin(); j != caloemcand->end(); ++j, ++m) {
          int rank = (int)j->rank();
          int index = (int)j->index();
          int rctCrate = (int)j->rctCrate();
          bool isolated = (bool)j->isolated();
          int rctCard = (int)j->rctCard();
          int rctRegion = (int)j->rctRegion();

          if (isolated) {
            converter[rctCrate].SetIEEt(rank, index);
            converter[rctCrate].SetIEReg(rctRegion, index);
            converter[rctCrate].SetIECard(rctCard, index);
            LogDebug("L1T") << "CRATE" << rctCrate << "ISO em rank=" << rank << ", region=" << rctRegion
                            << ", card=" << rctCard << std::endl;
          } else {
            converter[rctCrate].SetNEEt(rank, index);
            converter[rctCrate].SetNEReg(rctRegion, index);
            converter[rctCrate].SetNECard(rctCard, index);
            LogDebug("L1T") << "CRATE" << rctCrate << "NON ISO em rank=" << rank << ", region=" << rctRegion
                            << ", card=" << rctCard << std::endl;
          }
        }  //end of em cand

        for (int in = 0; in < 18; in++) {
          converter[in].Extract32bitwords();
          for (int d = 0; d < 6; d++) {
            load[2 * in].push_back((uint32_t)converter[in].Get32bitWordLinkEven(d));
            load[2 * in + 1].push_back((uint32_t)converter[in].Get32bitWordLinkOdd(d));
          }
        }
      }  // end of BX

      rctDataBase database;

      Blocks res = {};

      for (int i = 0; i < 36; ++i) {
        unsigned int mycrateRCT = (int)(i / 2);
        bool myRCTeven;
        if (i % 2 == 0)
          myRCTeven = true;
        else
          myRCTeven = false;
        int linkMP7 = -1;
        database.GetLinkMP7(mycrateRCT, myRCTeven, linkMP7);
        res.push_back(Block(2 * linkMP7, load[i], 0));
        res.push_back(Block(2 * linkMP7, load[i], 1));
      }
      return res;
    }
  }  // namespace stage1
}  // namespace l1t
DEFINE_L1T_PACKER(l1t::stage1::RCTEmRegionPacker);