HGCalTowerMap2DImpl.h

Go to the documentation of this file.

#ifndef __L1Trigger_L1THGCal_HGCalTowerMap2DImpl_h__
#define __L1Trigger_L1THGCal_HGCalTowerMap2DImpl_h__

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

#include "DataFormats/L1THGCal/interface/HGCalTriggerSums.h"
#include "DataFormats/L1THGCal/interface/HGCalTriggerCell.h"
#include "DataFormats/L1THGCal/interface/HGCalTowerMap.h"
#include "L1Trigger/L1THGCal/interface/HGCalTriggerTools.h"
#include "L1Trigger/L1THGCal/interface/HGCalTriggerTowerGeometryHelper.h"

class HGCalTowerMap2DImpl {
public:
  HGCalTowerMap2DImpl(const edm::ParameterSet& conf);

  void resetTowerMaps();

  template <class T>
  void buildTowerMap2D(const std::vector<edm::Ptr<T>>& ptrs, l1t::HGCalTowerMapBxCollection& towerMaps) {
    std::unordered_map<int, l1t::HGCalTowerMap> towerMapsTmp = newTowerMaps();

    for (const auto& ptr : ptrs) {
      bool isNose = triggerTools_.isNose(ptr->detId());
      unsigned layer = triggerTools_.layerWithOffset(ptr->detId());

      if (towerMapsTmp.find(layer) == towerMapsTmp.end()) {
        throw cms::Exception("Out of range")
            << "HGCalTowerMap2dImpl: Found trigger sum in layer " << layer << " for which there is no tower map\n";
      }
      // FIXME: should actually sum the energy not the Et...
      double calibPt = ptr->pt();
      if (useLayerWeights_)
        calibPt = layerWeights_[layer] * ptr->mipPt();

      double etEm = layer <= triggerTools_.lastLayerEE(isNose) ? calibPt : 0;
      double etHad = layer > triggerTools_.lastLayerEE(isNose) ? calibPt : 0;

      towerMapsTmp[layer].addEt(towerGeometryHelper_.getTriggerTower(*ptr), etEm, etHad);
    }

    /* store towerMaps in the persistent collection */
    towerMaps.resize(0, towerMapsTmp.size());
    int i = 0;
    for (const auto& towerMap : towerMapsTmp) {
      towerMaps.set(0, i, towerMap.second);
      i++;
    }
  }

  void setGeometry(const HGCalTriggerGeometryBase* const geom) {
    triggerTools_.setGeometry(geom);
    towerGeometryHelper_.setGeometry(geom);
  }

private:
  bool useLayerWeights_;
  std::vector<double> layerWeights_;
  HGCalTriggerTools triggerTools_;
  std::unordered_map<int, l1t::HGCalTowerMap> newTowerMaps();

  HGCalTriggerTowerGeometryHelper towerGeometryHelper_;
};

#endif