#include <HGCalConcentratorSuperTriggerCellImpl.h>

Classes
class	SuperTriggerCell

Public Member Functions
void	eventSetup (const edm::EventSetup &es)

	HGCalConcentratorSuperTriggerCellImpl (const edm::ParameterSet &conf)

void	select (const std::vector< l1t::HGCalTriggerCell > &trigCellVecInput, std::vector< l1t::HGCalTriggerCell > &trigCellVecOutput)

Private Types
enum	EnergyDivisionType { superTriggerCell, oneBitFraction, equalShare }

Private Member Functions
void	assignSuperTriggerCellEnergyAndPosition (l1t::HGCalTriggerCell &c, const SuperTriggerCell &stc) const

void	createAllTriggerCells (std::unordered_map< unsigned, SuperTriggerCell > &STCs, std::vector< l1t::HGCalTriggerCell > &trigCellVecOutput) const

uint32_t	getCompressedSTCEnergy (const SuperTriggerCell &stc) const

float	getTriggerCellOneBitFraction (float tcPt, float sumPt) const

Private Attributes
HGCalTriggerCellCalibration	calibration_

std::vector< unsigned >	coarsenTriggerCells_

HGCalCoarseTriggerCellMapping	coarseTCmapping_

EnergyDivisionType	energyDivisionType_

bool	fixedDataSizePerHGCROC_

double	oneBitFractionHighValue_

double	oneBitFractionLowValue_

double	oneBitFractionThreshold_

HGCalCoarseTriggerCellMapping	superTCmapping_

HGCalTriggerTools	triggerTools_

HGCalVFECompressionImpl	vfeCompression_

Static Private Attributes
static constexpr int	kHighDensityThickness_ = 0

static constexpr int	kOddNumberMask_ = 1

static constexpr int	kTriggerCellsForDivision_ = 4

Detailed Description

Definition at line 13 of file HGCalConcentratorSuperTriggerCellImpl.h.

Member Enumeration Documentation

◆ EnergyDivisionType

enum HGCalConcentratorSuperTriggerCellImpl::EnergyDivisionType

private

Enumerator
superTriggerCell
oneBitFraction
equalShare

Definition at line 27 of file HGCalConcentratorSuperTriggerCellImpl.h.

                           {
     superTriggerCell,
     oneBitFraction,
     equalShare,
   };

Constructor & Destructor Documentation

◆ HGCalConcentratorSuperTriggerCellImpl()

HGCalConcentratorSuperTriggerCellImpl::HGCalConcentratorSuperTriggerCellImpl ( const edm::ParameterSet & conf )

Definition at line 3 of file HGCalConcentratorSuperTriggerCellImpl.cc.

     : fixedDataSizePerHGCROC_(conf.getParameter<bool>("fixedDataSizePerHGCROC")),
       coarsenTriggerCells_(conf.getParameter<std::vector<unsigned>>("coarsenTriggerCells")),
       coarseTCmapping_(conf.getParameter<std::vector<unsigned>>("ctcSize")),
       superTCmapping_(conf.getParameter<std::vector<unsigned>>("stcSize")),
       calibration_(conf.getParameterSet("superTCCalibration")),
       vfeCompression_(conf.getParameterSet("superTCCompression")) {
   std::string energyType(conf.getParameter<string>("type_energy_division"));
  
   if (energyType == "superTriggerCell") {
     energyDivisionType_ = superTriggerCell;
   } else if (energyType == "oneBitFraction") {
     energyDivisionType_ = oneBitFraction;
  
     oneBitFractionThreshold_ = conf.getParameter<double>("oneBitFractionThreshold");
     oneBitFractionLowValue_ = conf.getParameter<double>("oneBitFractionLowValue");
     oneBitFractionHighValue_ = conf.getParameter<double>("oneBitFractionHighValue");
  
   } else if (energyType == "equalShare") {
     energyDivisionType_ = equalShare;
  
   } else {
     energyDivisionType_ = superTriggerCell;
   }
 }

References energyDivisionType_, equalShare, edm::ParameterSet::getParameter(), oneBitFraction, oneBitFractionHighValue_, oneBitFractionLowValue_, oneBitFractionThreshold_, AlCaHLTBitMon_QueryRunRegistry::string, and superTriggerCell.

Member Function Documentation

◆ assignSuperTriggerCellEnergyAndPosition()

void HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition	(	l1t::HGCalTriggerCell &	c,
		const SuperTriggerCell &	stc
	)		const

private

Definition at line 93 of file HGCalConcentratorSuperTriggerCellImpl.cc.

                                                                                                                        {
   //Compress and recalibrate STC energy
   uint32_t compressed_value = getCompressedSTCEnergy(stc);
  
   HGCalTriggerTools::SubDetectorType subdet = triggerTools_.getSubDetectorType(c.detId());
   int thickness = triggerTools_.thicknessIndex(c.detId(), true);
  
   GlobalPoint point;
   if ((fixedDataSizePerHGCROC_ && thickness > kHighDensityThickness_) || coarsenTriggerCells_[subdet]) {
     point = coarseTCmapping_.getCoarseTriggerCellPosition(coarseTCmapping_.getCoarseTriggerCellId(c.detId()));
   } else {
     point = triggerTools_.getTCPosition(c.detId());
   }
   c.setPosition(point);
  
   math::PtEtaPhiMLorentzVector p4(c.pt(), point.eta(), point.phi(), 0.);
   c.setP4(p4);
  
   if (energyDivisionType_ == superTriggerCell) {
     if (c.detId() == stc.getMaxId()) {
       c.setHwPt(compressed_value);
       calibration_.calibrateInGeV(c);
     } else {
       throw cms::Exception("NonMaxIdSuperTriggerCell")
           << "Trigger Cell with detId not equal to the maximum of the superTriggerCell found";
     }
   } else if (energyDivisionType_ == equalShare) {
     double coarseTriggerCellSize =
         coarsenTriggerCells_[subdet]
             ? double(
                   coarseTCmapping_.getConstituentTriggerCells(coarseTCmapping_.getCoarseTriggerCellId(stc.getMaxId()))
                       .size())
             : 1.;
  
     double denominator =
         fixedDataSizePerHGCROC_
             ? double(kTriggerCellsForDivision_)
             : double(superTCmapping_.getConstituentTriggerCells(stc.getSTCId()).size()) / coarseTriggerCellSize;
  
     c.setHwPt(std::round(compressed_value / denominator));
     calibration_.calibrateInGeV(c);
  
   } else if (energyDivisionType_ == oneBitFraction) {
     double frac = 0;
  
     if (c.detId() != stc.getMaxId()) {
       frac = getTriggerCellOneBitFraction(stc.getTCpt(c.detId()), stc.getSumPt());
     } else {
       frac = 1 - stc.getFractionSum();
     }
  
     c.setHwPt(std::round(compressed_value * frac));
     calibration_.calibrateInGeV(c);
   }
 }

References HltBtagPostValidation_cff::c, HGCalTriggerCellCalibration::calibrateInGeV(), calibration_, coarsenTriggerCells_, coarseTCmapping_, HLTTauDQMOffline_cfi::denominator, energyDivisionType_, equalShare, Exception, fixedDataSizePerHGCROC_, cropTnPTrees::frac, HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getCoarseTriggerCellPosition(), getCompressedSTCEnergy(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getFractionSum(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getMaxId(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getSTCId(), HGCalTriggerTools::getSubDetectorType(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getSumPt(), HGCalTriggerTools::getTCPosition(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getTCpt(), getTriggerCellOneBitFraction(), kHighDensityThickness_, kTriggerCellsForDivision_, oneBitFraction, p4, point, findQualityFiles::size, superTCmapping_, superTriggerCell, Calorimetry_cff::thickness, HGCalTriggerTools::thicknessIndex(), and triggerTools_.

Referenced by createAllTriggerCells().

◆ createAllTriggerCells()

void HGCalConcentratorSuperTriggerCellImpl::createAllTriggerCells	(	std::unordered_map< unsigned, SuperTriggerCell > &	STCs,
		std::vector< l1t::HGCalTriggerCell > &	trigCellVecOutput
	)		const

private

Definition at line 36 of file HGCalConcentratorSuperTriggerCellImpl.cc.

                                                                                                                  {
   for (auto& s : STCs) {
     std::vector<uint32_t> output_ids = superTCmapping_.getConstituentTriggerCells(s.second.getSTCId());
     if (output_ids.empty())
       continue;
  
     HGCalTriggerTools::SubDetectorType subdet = triggerTools_.getSubDetectorType(output_ids.at(0));
     int thickness = (!output_ids.empty() ? triggerTools_.thicknessIndex(output_ids.at(0), true) : 0);
  
     for (const auto& id : output_ids) {
       if (((fixedDataSizePerHGCROC_ && thickness > kHighDensityThickness_) || coarsenTriggerCells_[subdet]) &&
           (id != coarseTCmapping_.getRepresentativeDetId(id))) {
         continue;
       }
  
       if (!triggerTools_.getTriggerGeometry()->validTriggerCell(id)) {
         continue;
       }
  
       l1t::HGCalTriggerCell triggerCell;
       triggerCell.setDetId(id);
       if (energyDivisionType_ == superTriggerCell && id != s.second.getMaxId()) {
         continue;
       }
  
       DetId tc_Id(id);
  
       //To guard against the case in v8 geometry where
       //there might be different thicknesses within a module
       //This is a small effect, but in principle energy might
       //be lost.
       if (superTCmapping_.getCoarseTriggerCellId(id) != s.second.getSTCId()) {
         if (triggerTools_.getTriggerGeometry()->isV9Geometry()) {
           throw cms::Exception("NonExistingCoarseTC")
               << "The coarse trigger cell correponsing to the nominal trigger cell does not exist";
         } else {
           continue;
         }
       }
       trigCellVecOutput.push_back(triggerCell);
       if (energyDivisionType_ == oneBitFraction) {  //Get the 1 bit fractions
  
         if (id != s.second.getMaxId()) {
           float tc_fraction = getTriggerCellOneBitFraction(s.second.getTCpt(id), s.second.getSumPt());
           s.second.addToFractionSum(tc_fraction);
         }
       }
     }
   }
   // assign energy
   for (l1t::HGCalTriggerCell& tc : trigCellVecOutput) {
     const auto& stc = STCs[superTCmapping_.getCoarseTriggerCellId(tc.detId())];
     assignSuperTriggerCellEnergyAndPosition(tc, stc);
   }
 }

References assignSuperTriggerCellEnergyAndPosition(), coarsenTriggerCells_, coarseTCmapping_, energyDivisionType_, Exception, fixedDataSizePerHGCROC_, HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), HGCalCoarseTriggerCellMapping::getRepresentativeDetId(), HGCalTriggerTools::getSubDetectorType(), getTriggerCellOneBitFraction(), HGCalTriggerTools::getTriggerGeometry(), HGCalTriggerGeometryBase::isV9Geometry(), kHighDensityThickness_, oneBitFraction, alignCSCRings::s, l1t::HGCalTriggerCell::setDetId(), superTCmapping_, superTriggerCell, Calorimetry_cff::thickness, HGCalTriggerTools::thicknessIndex(), triggerTools_, and HGCalTriggerGeometryBase::validTriggerCell().

Referenced by select().

◆ eventSetup()

void HGCalConcentratorSuperTriggerCellImpl::eventSetup ( const edm::EventSetup & es )

inline

Definition at line 19 of file HGCalConcentratorSuperTriggerCellImpl.h.

                                            {
     triggerTools_.eventSetup(es);
     coarseTCmapping_.eventSetup(es);
     superTCmapping_.eventSetup(es);
     calibration_.eventSetup(es);
   }

References calibration_, coarseTCmapping_, HGCalCoarseTriggerCellMapping::eventSetup(), HGCalTriggerCellCalibration::eventSetup(), HGCalTriggerTools::eventSetup(), superTCmapping_, and triggerTools_.

◆ getCompressedSTCEnergy()

uint32_t HGCalConcentratorSuperTriggerCellImpl::getCompressedSTCEnergy ( const SuperTriggerCell & stc ) const

private

Definition at line 29 of file HGCalConcentratorSuperTriggerCellImpl.cc.

                                                                                                         {
   uint32_t code(0);
   uint32_t compressed_value(0);
   vfeCompression_.compressSingle(stc.getSumHwPt(), code, compressed_value);
   return compressed_value;
 }

References HGCalVFECompressionImpl::compressSingle(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getSumHwPt(), and vfeCompression_.

Referenced by assignSuperTriggerCellEnergyAndPosition().

◆ getTriggerCellOneBitFraction()

float HGCalConcentratorSuperTriggerCellImpl::getTriggerCellOneBitFraction	(	float	tcPt,
		float	sumPt
	)		const

private

Definition at line 150 of file HGCalConcentratorSuperTriggerCellImpl.cc.

                                                                                                        {
   double f = tcPt / sumPt;
   double frac = 0;
   if (f < oneBitFractionThreshold_) {
     frac = oneBitFractionLowValue_;
   } else {
     frac = oneBitFractionHighValue_;
   }
  
   return frac;
 }

References f, cropTnPTrees::frac, oneBitFractionHighValue_, oneBitFractionLowValue_, oneBitFractionThreshold_, and TtFullHadEvtBuilder_cfi::sumPt.

Referenced by assignSuperTriggerCellEnergyAndPosition(), and createAllTriggerCells().

◆ select()

void HGCalConcentratorSuperTriggerCellImpl::select	(	const std::vector< l1t::HGCalTriggerCell > &	trigCellVecInput,
		std::vector< l1t::HGCalTriggerCell > &	trigCellVecOutput
	)

Definition at line 162 of file HGCalConcentratorSuperTriggerCellImpl.cc.

                                                                                                       {
   std::unordered_map<unsigned, SuperTriggerCell> STCs;
   // first pass, fill the "coarse" trigger cells
   for (const l1t::HGCalTriggerCell& tc : trigCellVecInput) {
     uint32_t stcid = superTCmapping_.getCoarseTriggerCellId(tc.detId());
     STCs[stcid].add(tc, stcid);
   }
  
   createAllTriggerCells(STCs, trigCellVecOutput);
 }