#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_

bool	coarsenTriggerCells_

HGCalCoarseTriggerCellMapping	coarseTCmapping_

EnergyDivisionType	energyDivisionType_

bool	fixedDataSizePerHGCROC_

double	oneBitFractionHighValue_

double	oneBitFractionLowValue_

double	oneBitFractionThreshold_

HGCalCoarseTriggerCellMapping	superTCmapping_

HGCalTriggerTools	triggerTools_

HGCalVFECompressionImpl	vfeCompression_

Static Private Attributes
static int	kHighDensityThickness_ = 0

static int	kOddNumberMask_ = 1

static int	kTriggerCellsForDivision_ = 4

Detailed Description

Definition at line 13 of file HGCalConcentratorSuperTriggerCellImpl.h.

Member Enumeration Documentation

enum HGCalConcentratorSuperTriggerCellImpl::EnergyDivisionType

private

Enumerator
superTriggerCell
oneBitFraction
equalShare

Definition at line 27 of file HGCalConcentratorSuperTriggerCellImpl.h.

                           {
     superTriggerCell,
     oneBitFraction,
     equalShare,
   };

Constructor & Destructor Documentation

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

Definition at line 3 of file HGCalConcentratorSuperTriggerCellImpl.cc.

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

     : fixedDataSizePerHGCROC_(conf.getParameter<bool>("fixedDataSizePerHGCROC")),
       coarsenTriggerCells_(conf.getParameter<bool>("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;
   }
 }

Member Function Documentation

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

private

Definition at line 93 of file HGCalConcentratorSuperTriggerCellImpl.cc.

References HGCalTriggerCellCalibration::calibrateInGeV(), calibration_, coarsenTriggerCells_, coarseTCmapping_, HLTTauDQMOffline_cfi::denominator, l1t::HGCalTriggerCell::detId(), energyDivisionType_, equalShare, PV3DBase< T, PVType, FrameType >::eta(), Exception, fixedDataSizePerHGCROC_, cropTnPTrees::frac, HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getCoarseTriggerCellPosition(), getCompressedSTCEnergy(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getFractionSum(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getMaxId(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getSTCId(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getSumPt(), HGCalTriggerTools::getTCPosition(), HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::getTCpt(), getTriggerCellOneBitFraction(), kHighDensityThickness_, kTriggerCellsForDivision_, oneBitFraction, p4, PV3DBase< T, PVType, FrameType >::phi(), point, reco::LeafCandidate::pt(), l1t::L1Candidate::setHwPt(), reco::LeafCandidate::setP4(), l1t::HGCalTriggerCell::setPosition(), findQualityFiles::size, superTCmapping_, superTriggerCell, Calorimetry_cff::thickness, HGCalTriggerTools::thicknessIndex(), and triggerTools_.

Referenced by createAllTriggerCells(), and HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::size().

                                                                                                                        {
   //Compress and recalibrate STC energy
   uint32_t compressed_value = getCompressedSTCEnergy(stc);
 
   int thickness = triggerTools_.thicknessIndex(c.detId(), true);
 
   GlobalPoint point;
   if ((fixedDataSizePerHGCROC_ && thickness > kHighDensityThickness_) || coarsenTriggerCells_) {
     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_
             ? 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);
   }
 }

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

private

Definition at line 36 of file HGCalConcentratorSuperTriggerCellImpl.cc.

References assignSuperTriggerCellEnergyAndPosition(), coarseTCmapping_, energyDivisionType_, Exception, fixedDataSizePerHGCROC_, HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), HGCalCoarseTriggerCellMapping::getConstituentTriggerCells(), HGCalCoarseTriggerCellMapping::getRepresentativeDetId(), 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(), and HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::size().

                                                                                                                  {
   for (auto& s : STCs) {
     std::vector<uint32_t> output_ids = superTCmapping_.getConstituentTriggerCells(s.second.getSTCId());
 
     int thickness = (!output_ids.empty() ? triggerTools_.thicknessIndex(output_ids.at(0), true) : 0);
 
     for (const auto& id : output_ids) {
       if (fixedDataSizePerHGCROC_ && thickness > kHighDensityThickness_ &&
           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);
   }
 }

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

inline

Definition at line 19 of file HGCalConcentratorSuperTriggerCellImpl.h.

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

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

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

private

Definition at line 29 of file HGCalConcentratorSuperTriggerCellImpl.cc.

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

Referenced by assignSuperTriggerCellEnergyAndPosition(), and HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::size().

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

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

private

Definition at line 149 of file HGCalConcentratorSuperTriggerCellImpl.cc.

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

Referenced by assignSuperTriggerCellEnergyAndPosition(), createAllTriggerCells(), and HGCalConcentratorSuperTriggerCellImpl::SuperTriggerCell::size().

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

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

Definition at line 161 of file HGCalConcentratorSuperTriggerCellImpl.cc.

References createAllTriggerCells(), HGCalCoarseTriggerCellMapping::getCoarseTriggerCellId(), and superTCmapping_.

                                                                                                       {
   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);
 }