HGCalVFESummationImpl Class Reference

#include <HGCalVFESummationImpl.h>

Public Member Functions
	HGCalVFESummationImpl (const edm::ParameterSet &conf)
void	setGeometry (const HGCalTriggerGeometryBase *const geom)
void	triggerCellSums (const std::vector< std::pair< DetId, uint32_t > > &, std::unordered_map< uint32_t, uint32_t > &)

Private Attributes
double	lsb_scintillator_MIP_
double	lsb_silicon_fC_
double	threshold_scintillator_
std::vector< double >	thresholds_silicon_
HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 12 of file HGCalVFESummationImpl.h.

Constructor & Destructor Documentation

HGCalVFESummationImpl()

HGCalVFESummationImpl::HGCalVFESummationImpl

(

const edm::ParameterSet &

conf

)

Definition at line 3 of file HGCalVFESummationImpl.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), Exception, edm::ParameterSet::getParameter(), hgchebackDigitizer_cfi::noise, DiMuonV_cfg::threshold, threshold_scintillator_, thresholds_silicon_, and HcalDetIdTransform::transform().

    : lsb_silicon_fC_(conf.getParameter<double>("siliconCellLSB_fC")),
      lsb_scintillator_MIP_(conf.getParameter<double>("scintillatorCellLSB_MIP")) {
  constexpr unsigned nThickness = 3;
  thresholds_silicon_ =
      conf.getParameter<edm::ParameterSet>("noiseSilicon").getParameter<std::vector<double>>("values");
  if (thresholds_silicon_.size() != nThickness) {
    throw cms::Exception("Configuration") << thresholds_silicon_.size() << " silicon thresholds are given instead of "
                                          << nThickness << " (the number of sensor thicknesses)";
  }
  threshold_scintillator_ = conf.getParameter<edm::ParameterSet>("noiseScintillator").getParameter<double>("noise_MIP");
  const auto threshold = conf.getParameter<double>("noiseThreshold");
  std::transform(
      thresholds_silicon_.begin(), thresholds_silicon_.end(), thresholds_silicon_.begin(), [threshold](auto noise) {
        return noise * threshold;
      });
  threshold_scintillator_ *= threshold;
}

Member Function Documentation

setGeometry()

void HGCalVFESummationImpl::setGeometry ( const HGCalTriggerGeometryBase *const  geom )

inline

Definition at line 16 of file HGCalVFESummationImpl.h.

References relativeConstraints::geom, HGCalTriggerTools::setGeometry(), and triggerTools_.

{ triggerTools_.setGeometry(geom); }

triggerCellSums()

void HGCalVFESummationImpl::triggerCellSums	(	const std::vector< std::pair< DetId, uint32_t > > &	,
		std::unordered_map< uint32_t, uint32_t > &
	)

Definition at line 22 of file HGCalVFESummationImpl.cc.

References HGCalTriggerGeometryBase::getTriggerCellFromCell(), HGCalTriggerTools::getTriggerGeometry(), HGCalTriggerTools::isScintillator(), HGCalTriggerTools::isSilicon(), lsb_scintillator_MIP_, lsb_silicon_fC_, Calorimetry_cff::thickness, HGCalTriggerTools::thicknessIndex(), DiMuonV_cfg::threshold, threshold_scintillator_, thresholds_silicon_, triggerTools_, and relativeConstraints::value.

                                                                                                {
  if (input_dataframes.empty())
    return;
  // sum energies in trigger cells
  for (const auto& [cellid, value] : input_dataframes) {
    // Apply noise threshold before summing into trigger cells
    uint32_t value_zero_suppressed = value;
    if (triggerTools_.isSilicon(cellid)) {
      int thickness = triggerTools_.thicknessIndex(cellid);
      double threshold = thresholds_silicon_.at(thickness);
      value_zero_suppressed = (value * lsb_silicon_fC_ > threshold ? value : 0);
    } else if (triggerTools_.isScintillator(cellid)) {
      value_zero_suppressed = (value * lsb_scintillator_MIP_ > threshold_scintillator_ ? value : 0);
    }
    if (value_zero_suppressed == 0)
      continue;

    // find trigger cell associated to cell
    uint32_t tcid = triggerTools_.getTriggerGeometry()->getTriggerCellFromCell(cellid);
    triggercells.emplace(tcid, 0);  // do nothing if key exists already

    // sums energy for the same trigger cell id
    triggercells[tcid] += value_zero_suppressed;  // 32 bits integer should be largely enough
  }
}

Member Data Documentation

lsb_scintillator_MIP_

double HGCalVFESummationImpl::lsb_scintillator_MIP_

private

Definition at line 21 of file HGCalVFESummationImpl.h.

Referenced by triggerCellSums().

lsb_silicon_fC_

double HGCalVFESummationImpl::lsb_silicon_fC_

private

Definition at line 20 of file HGCalVFESummationImpl.h.

Referenced by triggerCellSums().

threshold_scintillator_

double HGCalVFESummationImpl::threshold_scintillator_

private

Definition at line 23 of file HGCalVFESummationImpl.h.

Referenced by HGCalVFESummationImpl(), and triggerCellSums().

thresholds_silicon_

std::vector<double> HGCalVFESummationImpl::thresholds_silicon_

private

Definition at line 22 of file HGCalVFESummationImpl.h.

Referenced by HGCalVFESummationImpl(), and triggerCellSums().

triggerTools_

HGCalTriggerTools HGCalVFESummationImpl::triggerTools_

private

Definition at line 25 of file HGCalVFESummationImpl.h.

Referenced by setGeometry(), and triggerCellSums().