HGCalConcentratorProcessorSelection.cc

Go to the documentation of this file.

#include "L1Trigger/L1THGCal/interface/concentrator/HGCalConcentratorProcessorSelection.h"
#include <limits>

#include "DataFormats/HGCDigi/interface/HGCDigiCollections.h"

DEFINE_EDM_PLUGIN(HGCalConcentratorFactory, HGCalConcentratorProcessorSelection, "HGCalConcentratorProcessorSelection");

HGCalConcentratorProcessorSelection::HGCalConcentratorProcessorSelection(const edm::ParameterSet& conf)
    : HGCalConcentratorProcessorBase(conf),
      fixedDataSizePerHGCROC_(conf.getParameter<bool>("fixedDataSizePerHGCROC")),
      allTrigCellsInTrigSums_(conf.getParameter<bool>("allTrigCellsInTrigSums")),
      coarsenTriggerCells_(conf.getParameter<std::vector<unsigned>>("coarsenTriggerCells")),
      selectionType_(kNSubDetectors_) {
  std::vector<std::string> selectionType(conf.getParameter<std::vector<std::string>>("Method"));
  if (selectionType.size() != kNSubDetectors_ || coarsenTriggerCells_.size() != kNSubDetectors_) {
    throw cms::Exception("HGCTriggerParameterError")
        << "Inconsistent number of sub-detectors (should be " << kNSubDetectors_ << ")";
  }

  for (int subdet = 0; subdet < kNSubDetectors_; subdet++) {
    if (selectionType[subdet] == "thresholdSelect") {
      selectionType_[subdet] = thresholdSelect;
      if (!thresholdImpl_)
        thresholdImpl_ = std::make_unique<HGCalConcentratorThresholdImpl>(conf);
      if (!trigSumImpl_)
        trigSumImpl_ = std::make_unique<HGCalConcentratorTrigSumImpl>(conf);
    } else if (selectionType[subdet] == "bestChoiceSelect") {
      selectionType_[subdet] = bestChoiceSelect;
      if (!bestChoiceImpl_)
        bestChoiceImpl_ = std::make_unique<HGCalConcentratorBestChoiceImpl>(conf);
      if (!trigSumImpl_)
        trigSumImpl_ = std::make_unique<HGCalConcentratorTrigSumImpl>(conf);
    } else if (selectionType[subdet] == "superTriggerCellSelect") {
      selectionType_[subdet] = superTriggerCellSelect;
      if (!superTriggerCellImpl_)
        superTriggerCellImpl_ = std::make_unique<HGCalConcentratorSuperTriggerCellImpl>(conf);
    } else if (selectionType[subdet] == "autoEncoder") {
      selectionType_[subdet] = autoEncoderSelect;
      if (!autoEncoderImpl_)
        autoEncoderImpl_ = std::make_unique<HGCalConcentratorAutoEncoderImpl>(conf);
    } else if (selectionType[subdet] == "noSelection") {
      selectionType_[subdet] = noSelection;
    } else {
      throw cms::Exception("HGCTriggerParameterError")
          << "Unknown type of concentrator selection '" << selectionType[subdet] << "'";
    }
  }

  if (std::find(coarsenTriggerCells_.begin(), coarsenTriggerCells_.end(), true) != coarsenTriggerCells_.end() ||
      fixedDataSizePerHGCROC_) {
    coarsenerImpl_ = std::make_unique<HGCalConcentratorCoarsenerImpl>(conf);
  }
}

void HGCalConcentratorProcessorSelection::run(const edm::Handle<l1t::HGCalTriggerCellBxCollection>& triggerCellCollInput,
                                              std::tuple<l1t::HGCalTriggerCellBxCollection,
                                                         l1t::HGCalTriggerSumsBxCollection,
                                                         l1t::HGCalConcentratorDataBxCollection>& triggerCollOutput) {
  if (thresholdImpl_)
    thresholdImpl_->setGeometry(geometry());
  if (bestChoiceImpl_)
    bestChoiceImpl_->setGeometry(geometry());
  if (superTriggerCellImpl_)
    superTriggerCellImpl_->setGeometry(geometry());
  if (autoEncoderImpl_)
    autoEncoderImpl_->setGeometry(geometry());
  if (coarsenerImpl_)
    coarsenerImpl_->setGeometry(geometry());
  if (trigSumImpl_)
    trigSumImpl_->setGeometry(geometry());
  triggerTools_.setGeometry(geometry());

  auto& triggerCellCollOutput = std::get<0>(triggerCollOutput);
  auto& triggerSumCollOutput = std::get<1>(triggerCollOutput);
  auto& autoEncoderCollOutput = std::get<2>(triggerCollOutput);

  const l1t::HGCalTriggerCellBxCollection& collInput = *triggerCellCollInput;

  std::unordered_map<uint32_t, std::vector<l1t::HGCalTriggerCell>> tc_modules;
  for (const auto& trigCell : collInput) {
    uint32_t module = geometry()->getModuleFromTriggerCell(trigCell.detId());
    tc_modules[module].push_back(trigCell);
  }

  for (const auto& module_trigcell : tc_modules) {
    std::vector<l1t::HGCalTriggerCell> trigCellVecOutput;
    std::vector<l1t::HGCalTriggerCell> trigCellVecCoarsened;
    std::vector<l1t::HGCalTriggerCell> trigCellVecNotSelected;
    std::vector<l1t::HGCalTriggerSums> trigSumsVecOutput;
    std::vector<l1t::HGCalConcentratorData> ae_EncodedLayerOutput;

    int thickness = triggerTools_.thicknessIndex(module_trigcell.second.at(0).detId());

    HGCalTriggerTools::SubDetectorType subdet = triggerTools_.getSubDetectorType(module_trigcell.second.at(0).detId());

    if (coarsenTriggerCells_[subdet] || (fixedDataSizePerHGCROC_ && thickness > kHighDensityThickness_)) {
      coarsenerImpl_->coarsen(module_trigcell.second, trigCellVecCoarsened);

      switch (selectionType_[subdet]) {
        case thresholdSelect:
          thresholdImpl_->select(trigCellVecCoarsened, trigCellVecOutput, trigCellVecNotSelected);
          break;
        case bestChoiceSelect:
          if (triggerTools_.isEm(module_trigcell.first)) {
            bestChoiceImpl_->select(geometry()->getLinksInModule(module_trigcell.first),
                                    geometry()->getModuleSize(module_trigcell.first),
                                    module_trigcell.second,
                                    trigCellVecOutput,
                                    trigCellVecNotSelected);
          } else {
            bestChoiceImpl_->select(geometry()->getLinksInModule(module_trigcell.first),
                                    geometry()->getModuleSize(module_trigcell.first),
                                    trigCellVecCoarsened,
                                    trigCellVecOutput,
                                    trigCellVecNotSelected);
          }
          break;
        case superTriggerCellSelect:
          superTriggerCellImpl_->select(trigCellVecCoarsened, trigCellVecOutput);
          break;
        case autoEncoderSelect:
          autoEncoderImpl_->select(geometry()->getLinksInModule(module_trigcell.first),
                                   trigCellVecCoarsened,
                                   trigCellVecOutput,
                                   ae_EncodedLayerOutput);
          break;
        case noSelection:
          trigCellVecOutput = trigCellVecCoarsened;
          break;
        default:
          // Should not happen, selection type checked in constructor
          break;
      }

    } else {
      switch (selectionType_[subdet]) {
        case thresholdSelect:
          thresholdImpl_->select(module_trigcell.second, trigCellVecOutput, trigCellVecNotSelected);
          break;
        case bestChoiceSelect:
          bestChoiceImpl_->select(geometry()->getLinksInModule(module_trigcell.first),
                                  geometry()->getModuleSize(module_trigcell.first),
                                  module_trigcell.second,
                                  trigCellVecOutput,
                                  trigCellVecNotSelected);
          break;
        case superTriggerCellSelect:
          superTriggerCellImpl_->select(module_trigcell.second, trigCellVecOutput);
          break;
        case autoEncoderSelect:
          autoEncoderImpl_->select(geometry()->getLinksInModule(module_trigcell.first),
                                   module_trigcell.second,
                                   trigCellVecOutput,
                                   ae_EncodedLayerOutput);
          break;
        case noSelection:
          trigCellVecOutput = module_trigcell.second;
          break;
        default:
          // Should not happen, selection type checked in constructor
          break;
      }
    }

    // trigger sum
    if (trigSumImpl_) {
      if (allTrigCellsInTrigSums_) {  // using all TCs
        trigSumImpl_->doSum(module_trigcell.first, module_trigcell.second, trigSumsVecOutput);
      } else {  // using only unselected TCs
        trigSumImpl_->doSum(module_trigcell.first, trigCellVecNotSelected, trigSumsVecOutput);
      }
    }

    for (const auto& trigCell : trigCellVecOutput) {
      triggerCellCollOutput.push_back(0, trigCell);
    }
    for (const auto& trigSums : trigSumsVecOutput) {
      triggerSumCollOutput.push_back(0, trigSums);
    }
    for (const auto& aeVal : ae_EncodedLayerOutput) {
      autoEncoderCollOutput.push_back(0, aeVal);
    }
  }
}