#include <HcalDDDRecConstants.h>

Classes
struct	HcalActiveLength

struct	HcalEtaBin

struct	HcalID

struct	HFCellParameters

Public Member Functions
const HcalDDDSimConstants *	dddConstants () const

std::vector< std::pair< double, double > >	getConstHBHE (const int &type) const

std::vector< int >	getDepth (const int &det, const int &phi, const int &zside, const unsigned int &eta) const

std::vector< int >	getDepth (const unsigned int &eta, const bool &extra) const

int	getDepthEta16 (const int &det, const int &iphi, const int &zside) const

int	getDepthEta29 (const int &iphi, const int &zside, const int &type) const

std::vector< HcalEtaBin >	getEtaBins (const int &itype) const

std::pair< double, double >	getEtaLimit (const int &i) const

std::pair< double, double >	getEtaPhi (const int &subdet, const int &ieta, const int &iphi) const

std::pair< int, int >	getEtaRange (const int &i) const

const std::vector< double > &	getEtaTable () const

const std::vector< double > &	getEtaTableHF () const

HcalID	getHCID (int subdet, int ieta, int iphi, int lay, int idepth) const

std::vector< HFCellParameters >	getHFCellParameters () const

double	getLayer0Wt (const int &det, const int &phi, const int &zside) const

int	getLayerBack (const int &det, const int &eta, const int &phi, const int &depth) const

void	getLayerDepth (const int &ieta, std::map< int, int > &layers) const

int	getLayerFront (const int &det, const int &eta, const int &phi, const int &depth) const

int	getMaxDepth (const int &itype, const int &ieta, const int &iphi, const int &zside) const

int	getMaxDepth (const int &type) const

int	getMinDepth (const int &itype, const int &ieta, const int &iphi, const int &zside) const

int	getNEta () const

int	getNoff (const int &i) const

int	getNPhi (const int &type) const

double	getPhiBin (const int &i) const

double	getPhiOff (const int &i) const

const std::vector< double > &	getPhiOffs () const

std::vector< std::pair< int, double > >	getPhis (const int &subdet, const int &ieta) const

const std::vector< double > &	getPhiTable () const

const std::vector< double > &	getPhiTableHF () const

int	getPhiZOne (std::vector< std::pair< int, int >> &phiz) const

std::pair< double, double >	getRZ (const HcalDetId &id) const

double	getRZ (const int &subdet, const int &ieta, const int &depth) const

double	getRZ (const int &subdet, const int &ieta, const int &iphi, const int &depth) const

double	getRZ (const int &subdet, const int &layer) const

std::vector< HcalActiveLength >	getThickActive (const int &type) const

int	getTopoMode () const

int	getTriggerMode () const

std::vector< HcalCellType >	HcalCellTypes (HcalSubdetector) const

	HcalDDDRecConstants (const HcalParameters *hp, const HcalDDDSimConstants &hc)

HcalDetId	idBack (const HcalDetId &id) const

HcalDetId	idFront (const HcalDetId &id) const

bool	isBH () const

bool	isPlan1 (const HcalDetId &id) const

bool	isPlan1MergedId (const HcalDetId &id) const

bool	isPlan1ToBeMergedId (const HcalDetId &id) const

int	maxHFDepth (int ieta, int iphi) const

HcalDetId	mergedDepthDetId (const HcalDetId &id) const

std::vector< int >	mergedDepthList29 (int ieta, int iphi) const

bool	mergedDepthList29 (int ieta, int iphi, int depth) const

unsigned int	nCells () const

unsigned int	nCells (HcalSubdetector) const

unsigned int	numberOfCells (HcalSubdetector) const

bool	specialRBXHBHE (bool flag, std::vector< HcalDetId > &) const

void	specialRBXHBHE (const std::vector< HcalDetId > &, std::vector< HcalDetId > &) const

void	unmergeDepthDetId (const HcalDetId &id, std::vector< HcalDetId > &ids) const

bool	withSpecialRBXHBHE () const

	~HcalDDDRecConstants ()

Private Member Functions
void	getOneEtaBin (HcalSubdetector subdet, int ieta, int zside, std::vector< std::pair< int, double >> &phis, std::map< int, int > &layers, bool planOne, std::vector< HcalDDDRecConstants::HcalEtaBin > &bins) const

void	initialize (void)

unsigned int	layerGroup (int eta, int i) const

unsigned int	layerGroupSize (int eta) const

Private Attributes
std::pair< int, int >	depthMaxDf_

std::pair< int, int >	depthMaxSp_

std::map< HcalDetId, HcalDetId >	detIdSp_

std::map< HcalDetId, std::vector< HcalDetId > >	detIdSpR_

std::vector< std::pair< int, int > >	etaSimValu

std::vector< double >	etaTable

std::vector< std::pair< double, double > >	gconsHB

std::vector< std::pair< double, double > >	gconsHE

const HcalDDDSimConstants &	hcons

const HcalParameters *	hpar

std::vector< int >	ietaMap

std::vector< int >	iEtaMax

std::vector< int >	iEtaMin

std::vector< int >	maxDepth

int	nHalves [2]

int	nModule [2]

std::vector< int >	nPhiBins

std::vector< double >	phibin

std::vector< int >	phiUnitS

Static Private Attributes
static const int	maxLayer_ = 18

static const int	maxLayerHB_ = 16

Detailed Description

this class reads the constant section of the hcal-sim-numbering xml-file

Author: Sunanda Banerjee, SINP sunan.nosp@m.da.b.nosp@m.anerj.nosp@m.ee@c.nosp@m.ern.c.nosp@m.h

Definition at line 23 of file HcalDDDRecConstants.h.

Constructor & Destructor Documentation

◆ HcalDDDRecConstants()

HcalDDDRecConstants::HcalDDDRecConstants	(	const HcalParameters *	hp,
		const HcalDDDSimConstants &	hc
	)

Definition at line 15 of file HcalDDDRecConstants.cc.

     : hpar(hp), hcons(hc) {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants::HcalDDDRecConstants (const HcalParameters* hp) constructor";
 #endif
   initialize();
 }

References initialize().

◆ ~HcalDDDRecConstants()

HcalDDDRecConstants::~HcalDDDRecConstants ( )

Definition at line 23 of file HcalDDDRecConstants.cc.

                                           {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants::destructed!!!";
 #endif
 }

Member Function Documentation

◆ dddConstants()

const HcalDDDSimConstants* HcalDDDRecConstants::dddConstants ( ) const

inline

Definition at line 126 of file HcalDDDRecConstants.h.

126 { return &hcons; }

References hcons.

◆ getConstHBHE()

std::vector<std::pair<double, double> > HcalDDDRecConstants::getConstHBHE ( const int & type ) const

inline

Definition at line 54 of file HcalDDDRecConstants.h.

                                                                          {
     if (type == 0)
       return gconsHB;
     else if (type == 1)
       return gconsHE;
     else {
       std::vector<std::pair<double, double>> gcons;
       return gcons;
     }
   }

References gconsHB, and gconsHE.

Referenced by HcalFlexiHardcodeGeometryLoader::makeHBCells(), and HcalFlexiHardcodeGeometryLoader::makeHECells().

◆ getDepth() [1/2]

std::vector< int > HcalDDDRecConstants::getDepth	(	const int &	det,
		const int &	phi,
		const int &	zside,
		const unsigned int &	eta
	)		const

Definition at line 52 of file HcalDDDRecConstants.cc.

                                                                               {
   std::map<int, int> layers;
   hcons.ldMap()->getLayerDepth(det, eta + 1, phi, zside, layers);
   if (layers.empty()) {
     return getDepth(eta, false);
   } else {
     std::vector<int> depths;
     for (unsigned int lay = 0; lay < layers.size(); ++lay)
       depths.emplace_back(layers[lay + 1]);
     return depths;
   }
 }

References HLT_FULL_cff::depths, PVValHelper::eta, HcalLayerDepthMap::getLayerDepth(), hcons, hgcalTopologyTester_cfi::layers, HcalDDDSimConstants::ldMap(), phi, and ecaldqm::zside().

Referenced by HcalTopology::HcalTopology().

◆ getDepth() [2/2]

std::vector< int > HcalDDDRecConstants::getDepth	(	const unsigned int &	eta,
		const bool &	extra
	)		const

Definition at line 29 of file HcalDDDRecConstants.cc.

                                                                                              {
   if (!extra) {
     std::vector<HcalParameters::LayerItem>::const_iterator last = hpar->layerGroupEtaRec.begin();
     for (std::vector<HcalParameters::LayerItem>::const_iterator it = hpar->layerGroupEtaRec.begin();
          it != hpar->layerGroupEtaRec.end();
          ++it) {
       if (it->layer == eta + 1)
         return it->layerGroup;
       if (it->layer > eta + 1)
         return last->layerGroup;
       last = it;
     }
     return last->layerGroup;
   } else {
     std::map<int, int> layers;
     hcons.ldMap()->getLayerDepth(eta + 1, layers);
     std::vector<int> depths;
     for (unsigned int lay = 0; lay < layers.size(); ++lay)
       depths.emplace_back(layers[lay + 1]);
     return depths;
   }
 }

References HLT_FULL_cff::depths, PVValHelper::eta, HcalLayerDepthMap::getLayerDepth(), hcons, hpar, dqmdumpme::last, HcalParameters::layerGroupEtaRec, hgcalTopologyTester_cfi::layers, and HcalDDDSimConstants::ldMap().

◆ getDepthEta16()

int HcalDDDRecConstants::getDepthEta16	(	const int &	det,
		const int &	iphi,
		const int &	zside
	)		const

inline

Definition at line 66 of file HcalDDDRecConstants.h.

                                                                              {
     return hcons.getDepthEta16(det, iphi, zside);
   }

References HcalDDDSimConstants::getDepthEta16(), hcons, LEDCalibrationChannels::iphi, and ecaldqm::zside().

Referenced by HcalTopology::depthBinInformation(), HcalTopology::incrementDepth(), and HcalTopology::validRaw().

◆ getDepthEta29()

int HcalDDDRecConstants::getDepthEta29	(	const int &	iphi,
		const int &	zside,
		const int &	type
	)		const

inline

Definition at line 69 of file HcalDDDRecConstants.h.

                                                                               {
     return hcons.getDepthEta29(iphi, zside, type);
   }

References HcalDDDSimConstants::getDepthEta29(), hcons, LEDCalibrationChannels::iphi, and ecaldqm::zside().

Referenced by HcalTopology::decrementDepth().

◆ getEtaBins()

std::vector< HcalDDDRecConstants::HcalEtaBin > HcalDDDRecConstants::getEtaBins ( const int & itype ) const

Definition at line 68 of file HcalDDDRecConstants.cc.

                                                                                                {
   std::vector<HcalDDDRecConstants::HcalEtaBin> bins;
   unsigned int type = (itype == 0) ? 0 : 1;
   HcalSubdetector subdet = HcalSubdetector(type + 1);
   std::vector<int> phiSp;
   HcalSubdetector subdetSp = HcalSubdetector(hcons.ldMap()->validDet(phiSp));
   std::map<int, int> layers;
   for (int iz = 0; iz < 2; ++iz) {
     int zside = 2 * iz - 1;
     for (int ieta = iEtaMin[type]; ieta <= iEtaMax[type]; ++ieta) {
       std::vector<std::pair<int, double>> phis = getPhis(subdet, ieta);
       std::vector<std::pair<int, double>> phiUse;
       getLayerDepth(ieta, layers);
       if (subdet == subdetSp) {
         for (auto& phi : phis) {
           if (std::find(phiSp.begin(), phiSp.end(), (zside * phi.first)) == phiSp.end()) {
             phiUse.emplace_back(phi);
           }
         }
       } else {
         phiUse.insert(phiUse.end(), phis.begin(), phis.end());
       }
       if (!phiUse.empty())
         getOneEtaBin(subdet, ieta, zside, phiUse, layers, false, bins);
     }
   }
   if (subdetSp == subdet) {
     for (int ieta = iEtaMin[type]; ieta <= iEtaMax[type]; ++ieta) {
       std::vector<std::pair<int, double>> phis = getPhis(subdet, ieta);
       for (int iz = 0; iz < 2; ++iz) {
         int zside = 2 * iz - 1;
         std::vector<std::pair<int, double>> phiUse;
         for (int i : phiSp) {
           for (auto& phi : phis) {
             if (i == zside * phi.first) {
               phiUse.emplace_back(phi);
               break;
             }
           }
         }
         if (!phiUse.empty()) {
           hcons.ldMap()->getLayerDepth(subdet, ieta, phiUse[0].first, zside, layers);
           getOneEtaBin(subdet, ieta, zside, phiUse, layers, true, bins);
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "Prepares " << bins.size() << " eta bins for type " << type;
   for (unsigned int i = 0; i < bins.size(); ++i) {
     edm::LogVerbatim("HcalGeom") << "Bin[" << i << "]: Eta = (" << bins[i].ieta << ":" << bins[i].etaMin << ":"
                                  << bins[i].etaMax << "), Zside = " << bins[i].zside << ", phis = ("
                                  << bins[i].phis.size() << ":" << bins[i].dphi << ") and " << bins[i].layer.size()
                                  << " depths (start) " << bins[i].depthStart;
     for (unsigned int k = 0; k < bins[i].layer.size(); ++k)
       edm::LogVerbatim("HcalGeom") << " [" << k << "] " << bins[i].layer[k].first << ":" << bins[i].layer[k].second;
     edm::LogVerbatim("HcalGeom") << "Phi sets";
     for (unsigned int k = 0; k < bins[i].phis.size(); ++k)
       edm::LogVerbatim("HcalGeom") << "[" << k << "] " << bins[i].phis[k].first << ":" << bins[i].phis[k].second;
   }
 #endif
   return bins;
 }

References trigObjTnPSource_cfi::bins, spr::find(), dqmdumpme::first, HcalLayerDepthMap::getLayerDepth(), getLayerDepth(), getOneEtaBin(), getPhis(), hcons, mps_fire::i, LEDCalibrationChannels::ieta, iEtaMax, iEtaMin, dqmdumpme::k, hgcalTopologyTester_cfi::layers, HcalDDDSimConstants::ldMap(), phi, HcalLayerDepthMap::validDet(), and ecaldqm::zside().

Referenced by getThickActive(), HcalCellTypes(), HcalTopology::HcalTopology(), HcalFlexiHardcodeGeometryLoader::makeHBCells(), HcalFlexiHardcodeGeometryLoader::makeHECells(), and nCells().

◆ getEtaLimit()

std::pair<double, double> HcalDDDRecConstants::getEtaLimit ( const int & i ) const

inline

Definition at line 77 of file HcalDDDRecConstants.h.

                                                           {
     return std::pair<double, double>(etaTable[i], etaTable[i + 1]);
   }

References etaTable, and mps_fire::i.

◆ getEtaPhi()

std::pair< double, double > HcalDDDRecConstants::getEtaPhi	(	const int &	subdet,
		const int &	ieta,
		const int &	iphi
	)		const

Definition at line 132 of file HcalDDDRecConstants.cc.

                                                                                                                 {
   int ietaAbs = (ieta > 0) ? ieta : -ieta;
   double eta(0), phi(0);
   if ((subdet == static_cast<int>(HcalBarrel)) || (subdet == static_cast<int>(HcalEndcap)) ||
       (subdet == static_cast<int>(HcalOuter))) {  // Use Eta Table
     int unit = hcons.unitPhi(phibin[ietaAbs - 1]);
     int kphi = (unit == 2) ? ((iphi - 1) / 2 + 1) : iphi;
     double foff = (ietaAbs <= iEtaMax[0]) ? hpar->phioff[0] : hpar->phioff[1];
     eta = 0.5 * (etaTable[ietaAbs - 1] + etaTable[ietaAbs]);
     phi = foff + (kphi - 0.5) * phibin[ietaAbs - 1];
   } else {
     ietaAbs -= iEtaMin[2];
     int unit = hcons.unitPhi(hpar->phitable[ietaAbs - 1]);
     int kphi = (unit == 4) ? ((iphi - 3) / 4 + 1) : ((iphi - 1) / 2 + 1);
     double foff = (unit > 2) ? hpar->phioff[4] : hpar->phioff[2];
     eta = 0.5 * (hpar->etaTableHF[ietaAbs - 1] + hpar->etaTableHF[ietaAbs]);
     phi = foff + (kphi - 0.5) * hpar->phitable[ietaAbs - 1];
   }
   if (ieta < 0)
     eta = -eta;
   if (phi > M_PI)
     phi -= (2 * M_PI);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getEtaPhi: subdet|ieta|iphi " << subdet << "|" << ieta << "|" << iphi << " eta|phi "
                                << eta << "|" << phi;
 #endif
   return std::pair<double, double>(eta, phi);
 }

References PVValHelper::eta, etaTable, HcalParameters::etaTableHF, HcalBarrel, HcalEndcap, HcalOuter, hcons, hpar, LEDCalibrationChannels::ieta, iEtaMax, iEtaMin, LEDCalibrationChannels::iphi, M_PI, phi, phibin, HcalParameters::phioff, HcalParameters::phitable, unit(), and HcalDDDSimConstants::unitPhi().

Referenced by HGCalSimHitStudy::analyzeHits(), and HGCalSimHitValidation::analyzeHits().

◆ getEtaRange()

std::pair<int, int> HcalDDDRecConstants::getEtaRange ( const int & i ) const

inline

Definition at line 74 of file HcalDDDRecConstants.h.

74 { return std::pair<int, int>(iEtaMin[i], iEtaMax[i]); }

References mps_fire::i, iEtaMax, and iEtaMin.

Referenced by HcalTestNumberingTester::analyze(), HcalSimHitStudy::bookHistograms(), SimHitsValidationHcal::bookHistograms(), HcalSimHitsValidation::bookHistograms(), CaloTowersAnalyzer::dqmBeginRun(), HcalRecHitsAnalyzer::dqmBeginRun(), HGCalTriggerTools::getLayerZ(), SimHitsValidationHcal::getLimits(), HcalTopology::HcalTopology(), and HcalTopology::nPhiBins().

◆ getEtaTable()

const std::vector<double>& HcalDDDRecConstants::getEtaTable ( ) const

inline

Definition at line 75 of file HcalDDDRecConstants.h.

75 { return etaTable; }

References etaTable.

Referenced by HcalTopology::HcalTopology(), and CaloTowerHardcodeGeometryLoader::load().

◆ getEtaTableHF()

const std::vector<double>& HcalDDDRecConstants::getEtaTableHF ( ) const

inline

Definition at line 76 of file HcalDDDRecConstants.h.

76 { return hpar->etaTableHF; }

References HcalParameters::etaTableHF, and hpar.

Referenced by HcalTopology::HcalTopology(), and CaloTowerHardcodeGeometryLoader::load().

◆ getHCID()

HcalDDDRecConstants::HcalID HcalDDDRecConstants::getHCID	(	int	subdet,
		int	ieta,
		int	iphi,
		int	lay,
		int	idepth
	)		const

Definition at line 161 of file HcalDDDRecConstants.cc.

                                                                                                                 {
   int ieta = (keta > 0) ? keta : -keta;
   int zside = (keta > 0) ? 1 : -1;
   int eta(ieta), phi(iphi), depth(idepth);
   if ((subdet == static_cast<int>(HcalOuter)) ||
       ((subdet == static_cast<int>(HcalBarrel)) && (lay > maxLayerHB_ + 1))) {
     subdet = static_cast<int>(HcalOuter);
     depth = 4;
   } else if (subdet == static_cast<int>(HcalBarrel) || subdet == static_cast<int>(HcalEndcap)) {
     eta = ietaMap[ieta - 1];
     int unit = phiUnitS[ieta - 1];
     int phi0 = (iphi - 1) / (hpar->phigroup[eta - 1]);
     if (unit == 2) {
       phi0 = (iphi + 1) / 2;
       phi0 = (phi0 - 1) / (hpar->phigroup[eta - 1]);
     } else if (unit == 4) {
       phi0 = (iphi + 1) / 4;
       phi0 = (phi0 - 1) / (hpar->phigroup[eta - 1]);
     }
     ++phi0;
     unit = hcons.unitPhi(phibin[eta - 1]);
     phi = hcons.phiNumber(phi0, unit);
     depth = hcons.findDepth(subdet, eta, phi, zside, lay - 1);
     if (depth <= 0)
       depth = layerGroup(eta - 1, lay - 1);
     if (eta == iEtaMin[1]) {
       if (subdet == static_cast<int>(HcalBarrel)) {
         if (depth > hcons.getDepthEta16(subdet, phi, zside))
           depth = hcons.getDepthEta16(subdet, phi, zside);
       } else {
         if (depth < hcons.getDepthEta16(subdet, phi, zside))
           depth = hcons.getDepthEta16(subdet, phi, zside);
       }
     } else if (eta == hpar->noff[0] && lay > 1) {
       int kphi = phi + int((hpar->phioff[3] + 0.1) / phibin[eta - 1]);
       kphi = (kphi - 1) % 4 + 1;
       if (kphi == 2 || kphi == 3)
         depth = layerGroup(eta - 1, lay - 2);
     } else if (eta == hpar->noff[1] && depth > hcons.getDepthEta29(phi, zside, 0)) {
       eta -= hcons.getDepthEta29(phi, zside, 1);
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getHCID: input " << subdet << ":" << ieta << ":" << iphi << ":" << idepth << ":"
                                << lay << " output " << eta << ":" << phi << ":" << depth;
 #endif
   return HcalDDDRecConstants::HcalID(subdet, eta, phi, depth);
 }

References LEDCalibrationChannels::depth, PVValHelper::eta, HcalDDDSimConstants::findDepth(), HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta29(), HcalBarrel, HcalEndcap, HcalOuter, hcons, hpar, LEDCalibrationChannels::ieta, ietaMap, iEtaMin, createfilelist::int, LEDCalibrationChannels::iphi, layerGroup(), maxLayerHB_, HcalParameters::noff, phi, phibin, HcalParameters::phigroup, HcalDDDSimConstants::phiNumber(), HcalParameters::phioff, phiUnitS, unit(), HcalDDDSimConstants::unitPhi(), and ecaldqm::zside().

Referenced by HGCalHitValidation::analyze(), HcalGeometry::getBackPosition(), HcalGeometry::getPosition(), and HcalHitRelabeller::relabel().

◆ getHFCellParameters()

std::vector< HcalDDDRecConstants::HFCellParameters > HcalDDDRecConstants::getHFCellParameters ( ) const

Definition at line 210 of file HcalDDDRecConstants.cc.

                                                                                               {
   std::vector<HcalDDDRecConstants::HFCellParameters> cells;
   unsigned int nEta = hcons.getPhiTableHF().size();
   if (maxDepth[2] > 0) {
     for (unsigned int k = 0; k < nEta; ++k) {
       int ieta = iEtaMin[2] + k;
       int dphi = (int)(0.001 + hcons.getPhiTableHF()[k] / (5._deg));
       int iphi = (dphi == 4) ? 3 : 1;
       int nphi = 72 / dphi;
       double rMin = hcons.getRTableHF()[nEta - k - 1] / CLHEP::cm;
       double rMax = hcons.getRTableHF()[nEta - k] / CLHEP::cm;
       HcalDDDRecConstants::HFCellParameters cell1(ieta, 1, iphi, dphi, nphi, rMin, rMax);
       cells.emplace_back(cell1);
       HcalDDDRecConstants::HFCellParameters cell2(-ieta, 1, iphi, dphi, nphi, rMin, rMax);
       cells.emplace_back(cell2);
     }
   }
   if (maxDepth[2] > 2) {
     if (!hcons.getIdHF2QIE().empty()) {
       for (unsigned int k = 0; k < hcons.getIdHF2QIE().size(); ++k) {
         int ieta = hcons.getIdHF2QIE()[k].ieta();
         int ind = std::abs(ieta) - iEtaMin[2];
         int dphi = (int)(0.001 + hcons.getPhiTableHF()[ind] / (5._deg));
         int iphi = hcons.getIdHF2QIE()[k].iphi();
         double rMin = hcons.getRTableHF()[nEta - ind - 1] / CLHEP::cm;
         double rMax = hcons.getRTableHF()[nEta - ind] / CLHEP::cm;
         HcalDDDRecConstants::HFCellParameters cell1(ieta, 3, iphi, dphi, 1, rMin, rMax);
         cells.emplace_back(cell1);
       }
     } else {
       for (unsigned int k = 0; k < nEta; ++k) {
         int ieta = iEtaMin[2] + k;
         int dphi = (int)(0.001 + hcons.getPhiTableHF()[k] / (5._deg));
         int iphi = (dphi == 4) ? 3 : 1;
         int nphi = 72 / dphi;
         double rMin = hcons.getRTableHF()[nEta - k - 1] / CLHEP::cm;
         double rMax = hcons.getRTableHF()[nEta - k] / CLHEP::cm;
         HcalDDDRecConstants::HFCellParameters cell1(ieta, 3, iphi, dphi, nphi, rMin, rMax);
         cells.emplace_back(cell1);
         HcalDDDRecConstants::HFCellParameters cell2(-ieta, 3, iphi, dphi, nphi, rMin, rMax);
         cells.emplace_back(cell2);
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants returns " << cells.size() << " HF cell parameters";
   for (unsigned int k = 0; k < cells.size(); ++k)
     edm::LogVerbatim("HcalGeom") << "Cell[" << k << "] : (" << cells[k].ieta << ", " << cells[k].depth << ", "
                                  << cells[k].firstPhi << ", " << cells[k].stepPhi << ", " << cells[k].nPhi << ", "
                                  << cells[k].rMin << ", " << cells[k].rMax << ")";
 #endif
   return cells;
 }

References funct::abs(), hgcalTopologyTester_cfi::cell1, hgcalTopologyTester_cfi::cell2, postprocess-scan-build::cells, HcalDDDSimConstants::getIdHF2QIE(), HcalDDDSimConstants::getPhiTableHF(), HcalDDDSimConstants::getRTableHF(), hcons, LEDCalibrationChannels::ieta, iEtaMin, createfilelist::int, LEDCalibrationChannels::iphi, dqmdumpme::k, maxDepth, HLT_FULL_cff::nEta, nphi, photonAnalyzer_cfi::rMax, and photonAnalyzer_cfi::rMin.

Referenced by HcalFlexiHardcodeGeometryLoader::makeHFCells().

◆ getLayer0Wt()

double HcalDDDRecConstants::getLayer0Wt	(	const int &	det,
		const int &	phi,
		const int &	zside
	)		const

inline

Definition at line 85 of file HcalDDDRecConstants.h.

                                                                              {
     return hcons.getLayer0Wt(det, phi, zside);
   }

References HcalDDDSimConstants::getLayer0Wt(), hcons, phi, and ecaldqm::zside().

Referenced by HcalHitRelabeller::energyWt().

◆ getLayerBack()

int HcalDDDRecConstants::getLayerBack	(	const int &	det,
		const int &	eta,
		const int &	phi,
		const int &	depth
	)		const

Definition at line 277 of file HcalDDDRecConstants.cc.

                                                                                                                {
   int subdet = (idet == 1) ? 1 : 2;
   int zside = (ieta > 0) ? 1 : -1;
   int eta = zside * ieta;
   int layBack = hcons.ldMap()->getLayerBack(subdet, eta, iphi, zside, depth);
   int laymax = hcons.getLastLayer(subdet, ieta);
   if (layBack < 0 && eta <= hpar->etaMax[1]) {
     for (unsigned int k = 0; k < layerGroupSize(eta - 1); ++k) {
       if (depth + 1 == (int)layerGroup(eta - 1, k)) {
         layBack = k - 1;
         break;
       }
     }
   }
   if (layBack < 0 || layBack > laymax)
     layBack = laymax;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getLayerBack::Input " << idet << ":" << ieta << ":" << iphi << ":" << depth
                                << " Output " << layBack;
 #endif
   return layBack;
 }

References LEDCalibrationChannels::depth, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalDDDSimConstants::getLastLayer(), HcalLayerDepthMap::getLayerBack(), hcons, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, dqmdumpme::k, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), and ecaldqm::zside().

Referenced by getRZ().

◆ getLayerDepth()

void HcalDDDRecConstants::getLayerDepth	(	const int &	ieta,
		std::map< int, int > &	layers
	)		const

Definition at line 264 of file HcalDDDRecConstants.cc.

                                                                                        {
   layers.clear();
   for (unsigned int l = 0; l < layerGroupSize(ieta - 1); ++l) {
     int lay = l + 1;
     layers[lay] = layerGroup(ieta - 1, l);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getLayerDepth::Input " << ieta << " Output " << layers.size() << " entries";
   for (std::map<int, int>::iterator itr = layers.begin(); itr != layers.end(); ++itr)
     edm::LogVerbatim("HcalGeom") << " [" << itr->first << "] " << itr->second;
 #endif
 }

References LEDCalibrationChannels::ieta, cmsLHEtoEOSManager::l, layerGroup(), layerGroupSize(), and hgcalTopologyTester_cfi::layers.

Referenced by getEtaBins().

◆ getLayerFront()

int HcalDDDRecConstants::getLayerFront	(	const int &	det,
		const int &	eta,
		const int &	phi,
		const int &	depth
	)		const

Definition at line 300 of file HcalDDDRecConstants.cc.

                                                                                                                 {
   int subdet = (idet == 1) ? 1 : 2;
   int zside = (ieta > 0) ? 1 : -1;
   int eta = zside * ieta;
   int layFront = hcons.ldMap()->getLayerFront(subdet, eta, iphi, zside, depth);
   int laymin = hcons.getFrontLayer(subdet, ieta);
   if ((layFront < 0) || ((subdet == static_cast<int>(HcalEndcap)) && (eta == 16))) {
     if ((subdet == static_cast<int>(HcalEndcap)) && (eta == 16)) {
       layFront = laymin;
     } else if (eta <= hpar->etaMax[1]) {
       for (unsigned int k = 0; k < layerGroupSize(eta - 1); ++k) {
         if (depth == (int)layerGroup(eta - 1, k)) {
           if ((int)(k) >= laymin) {
             layFront = k;
             break;
           }
         }
       }
     }
   } else {
     if (layFront < laymin)
       layFront = laymin;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getLayerFront::Input " << idet << ":" << ieta << ":" << iphi << ":" << depth
                                << " Output " << layFront;
 #endif
   return layFront;
 }

References LEDCalibrationChannels::depth, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalDDDSimConstants::getFrontLayer(), HcalLayerDepthMap::getLayerFront(), HcalEndcap, hcons, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, dqmdumpme::k, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), and ecaldqm::zside().

Referenced by getRZ(), and getThickActive().

◆ getMaxDepth() [1/2]

int HcalDDDRecConstants::getMaxDepth	(	const int &	itype,
		const int &	ieta,
		const int &	iphi,
		const int &	zside
	)		const

Definition at line 330 of file HcalDDDRecConstants.cc.

                                                                                                                {
   unsigned int type = (itype == 0) ? 0 : 1;
   int lmax = hcons.getMaxDepth(type + 1, ieta, iphi, zside, true);
   if (lmax < 0) {
     unsigned int lymax = (type == 0) ? maxLayerHB_ + 1 : maxLayer_ + 1;
     lmax = 0;
     if (layerGroupSize(ieta - 1) > 0) {
       if (layerGroupSize(ieta - 1) < lymax)
         lymax = layerGroupSize(ieta - 1);
       lmax = (int)(layerGroup(ieta - 1, lymax - 1));
       if (type == 0 && ieta == iEtaMax[type])
         lmax = hcons.getDepthEta16M(1);
       if (type == 1 && ieta >= hpar->noff[1])
         lmax = hcons.getDepthEta29M(0, false);
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getMaxDepth::Input " << itype << ":" << ieta << ":" << iphi << ":" << zside
                                << " Output " << lmax;
 #endif
   return lmax;
 }

References HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29M(), HcalDDDSimConstants::getMaxDepth(), hcons, hpar, LEDCalibrationChannels::ieta, iEtaMax, createfilelist::int, LEDCalibrationChannels::iphi, layerGroup(), layerGroupSize(), maxLayer_, maxLayerHB_, HcalParameters::noff, and ecaldqm::zside().

◆ getMaxDepth() [2/2]

int HcalDDDRecConstants::getMaxDepth ( const int & type ) const

inline

Definition at line 88 of file HcalDDDRecConstants.h.

88 { return maxDepth[type]; }

References maxDepth.

Referenced by HGCalDigiClient::beginRun(), HGCalRecHitsClient::beginRun(), HGCalSimHitsClient::beginRun(), HcalSimHitsClient::beginRun(), HGCalSimHitStudy::beginRun(), HcalRecHitsDQMClient::beginRun(), HcalSimHitStudy::bookHistograms(), SimHitsValidationHcal::bookHistograms(), HcalSimHitsValidation::bookHistograms(), HcalTopology::depthBinInformation(), HGCalRecHitValidation::dqmBeginRun(), HcalDigisValidation::dqmBeginRun(), HcalRecHitsAnalyzer::dqmBeginRun(), HGCalDigiValidation::dqmBeginRun(), HGCalSimHitValidation::dqmBeginRun(), HGCalTriggerTools::eventSetup(), HcalTopology::HcalTopology(), HGCalTriggerGeometryHexLayerBasedImp1::initialize(), mergedDepthList29(), HcalTopologyIdealEP::produce(), and HcalTopology::validRaw().

◆ getMinDepth()

int HcalDDDRecConstants::getMinDepth	(	const int &	itype,
		const int &	ieta,
		const int &	iphi,
		const int &	zside
	)		const

Definition at line 353 of file HcalDDDRecConstants.cc.

                                                                                                                {
   int lmin = hcons.getMinDepth(itype + 1, ieta, iphi, zside, true);
   if (lmin < 0) {
     if (itype == 2) {  // HFn
       lmin = 1;
     } else if (itype == 3) {  //HO
       lmin = maxDepth[3];
     } else {
       unsigned int type = (itype == 0) ? 0 : 1;
       if (layerGroupSize(ieta - 1) > 0) {
         if (type == 1 && ieta == iEtaMin[type])
           lmin = hcons.getDepthEta16M(2);
         else
           lmin = (int)(layerGroup(ieta - 1, 0));
       }
     }
   }
   return lmin;
 }

References HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getMinDepth(), hcons, LEDCalibrationChannels::ieta, iEtaMin, createfilelist::int, LEDCalibrationChannels::iphi, layerGroup(), layerGroupSize(), maxDepth, and ecaldqm::zside().

Referenced by HcalHBHEMuonSimAnalyzer::analyze(), HcalRaddamMuon::analyze(), HcalHBHEMuonAnalyzer::depth16HE(), HcalHBHEMuonHighEtaAnalyzer::depth16HE(), HcalTopology::depthBinInformation(), and HcalGeometry::getClosestCell().

◆ getNEta()

int HcalDDDRecConstants::getNEta ( ) const

inline

Definition at line 91 of file HcalDDDRecConstants.h.

91 { return hpar->etagroup.size(); }

References HcalParameters::etagroup, and hpar.

Referenced by HcalTopology::HcalTopology().

◆ getNoff()

int HcalDDDRecConstants::getNoff ( const int & i ) const

inline

Definition at line 92 of file HcalDDDRecConstants.h.

92 { return hpar->noff[i]; }

References hpar, mps_fire::i, and HcalParameters::noff.

Referenced by HcalTopology::validRaw().

◆ getNPhi()

int HcalDDDRecConstants::getNPhi ( const int & type ) const

inline

Definition at line 93 of file HcalDDDRecConstants.h.

93 { return nPhiBins[type]; }

References nPhiBins.

Referenced by HcalSimHitStudy::bookHistograms(), SimHitsValidationHcal::bookHistograms(), HcalSimHitsValidation::bookHistograms(), CaloTowersAnalyzer::dqmBeginRun(), HcalRecHitsAnalyzer::dqmBeginRun(), and HcalTopology::HcalTopology().

◆ getOneEtaBin()

void HcalDDDRecConstants::getOneEtaBin	(	HcalSubdetector	subdet,
		int	ieta,
		int	zside,
		std::vector< std::pair< int, double >> &	phis,
		std::map< int, int > &	layers,
		bool	planOne,
		std::vector< HcalDDDRecConstants::HcalEtaBin > &	bins
	)		const

private

Definition at line 707 of file HcalDDDRecConstants.cc.

                                                                                              {
   unsigned int lymax = (subdet == HcalBarrel) ? maxLayerHB_ + 1 : maxLayer_ + 1;
   int type = (subdet == HcalBarrel) ? 0 : 1;
   double dphi = phibin[ieta - 1];
   HcalDDDRecConstants::HcalEtaBin etabin =
       HcalDDDRecConstants::HcalEtaBin(ieta, zside, dphi, etaTable[ieta - 1], etaTable[ieta]);
   etabin.phis.insert(etabin.phis.end(), phis.begin(), phis.end());
   int n = (ieta == iEtaMax[type]) ? 0 : 1;
   HcalDDDRecConstants::HcalEtaBin etabin0 =
       HcalDDDRecConstants::HcalEtaBin(ieta, zside, dphi, etaTable[ieta - 1], etaTable[ieta + n]);
   etabin0.depthStart = hcons.getDepthEta29(phis[0].first, zside, 0) + 1;
   int dstart = -1;
   int lmin(0), lmax(0);
  
   std::map<int, int>::iterator itr = layers.begin();
   if (!layers.empty()) {
     int dep = itr->second;
     if (subdet == HcalEndcap && ieta == iEtaMin[type])
       dep = hcons.getDepthEta16(subdet, phis[0].first, zside);
     unsigned lymx0 = (layers.size() > lymax) ? lymax : layers.size();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalGeom") << "Eta " << ieta << ":" << hpar->noff[1] << " zside " << zside << " lymax " << lymx0
                                  << ":" << lymax << " Depth " << dep << ":" << itr->second;
     unsigned int l(0);
     for (itr = layers.begin(); itr != layers.end(); ++itr, ++l)
       edm::LogVerbatim("HcalGeom") << "Layer [" << l << "] " << itr->first << ":" << itr->second;
     edm::LogVerbatim("HcalGeom") << "With " << phis.size() << " phis";
     for (unsigned int l = 0; l < phis.size(); ++l)
       edm::LogVerbatim("HcalGeom") << "[" << l << "] " << phis[l].first << ":" << convertRadToDeg(phis[l].second);
 #endif
     for (itr = layers.begin(); itr != layers.end(); ++itr) {
       if (itr->first <= (int)(lymx0)) {
         if (itr->second == dep) {
           if (lmin == 0)
             lmin = itr->first;
           lmax = itr->first;
         } else if (itr->second > dep) {
           if (dstart < 0)
             dstart = dep;
           int lmax0 = (lmax >= lmin) ? lmax : lmin;
           if (subdet == HcalEndcap && ieta + 1 == hpar->noff[1] && dep > hcons.getDepthEta29(phis[0].first, zside, 0)) {
             etabin0.layer.emplace_back(std::pair<int, int>(lmin, lmax0));
           } else {
             etabin.layer.emplace_back(std::pair<int, int>(lmin, lmax0));
           }
           lmin = itr->first;
           lmax = lmin - 1;
           dep = itr->second;
         }
         if (subdet == HcalBarrel && ieta == iEtaMax[type] && dep > hcons.getDepthEta16M(1))
           break;
         if (subdet == HcalEndcap && ieta == hpar->noff[1] && dep > hcons.getDepthEta29M(0, planOne)) {
           lmax = lymx0;
           break;
         }
         if (itr->first == (int)(lymx0))
           lmax = lymx0;
       }
     }
     if (lmax >= lmin) {
       if (ieta + 1 == hpar->noff[1]) {
         etabin0.layer.emplace_back(std::pair<int, int>(lmin, lmax));
         etabin0.phis.insert(etabin0.phis.end(), phis.begin(), phis.end());
         bins.emplace_back(etabin0);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HcalGeom") << "etabin0: dStatrt " << etabin0.depthStart << " layers " << etabin0.layer.size()
                                      << ":" << lmin << ":" << lmax << " phis " << phis.size();
         for (unsigned int k = 0; k < etabin0.layer.size(); ++k)
           edm::LogVerbatim("HcalGeom") << " [" << k << "] " << etabin0.layer[k].first << ":" << etabin0.layer[k].second;
 #endif
       } else if (ieta == hpar->noff[1]) {
       } else {
         etabin.layer.emplace_back(std::pair<int, int>(lmin, lmax));
         if (dstart < 0)
           dstart = dep;
       }
     }
   }
   etabin.depthStart = dstart;
   bins.emplace_back(etabin);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "etabin: dStatrt " << etabin.depthStart << " layers " << etabin.layer.size() << ":"
                                << lmin << ":" << lmax << " phis " << etabin.phis.size();
   for (unsigned int k = 0; k < etabin.layer.size(); ++k)
     edm::LogVerbatim("HcalGeom") << "[" << k << "] " << etabin.layer[k].first << ":" << etabin.layer[k].second;
 #endif
 }

References trigObjTnPSource_cfi::bins, angle_units::operators::convertRadToDeg(), HcalDDDRecConstants::HcalEtaBin::depthStart, etaTable, dqmdumpme::first, HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29(), HcalDDDSimConstants::getDepthEta29M(), HcalBarrel, HcalEndcap, hcons, hpar, LEDCalibrationChannels::ieta, iEtaMax, iEtaMin, dqmdumpme::k, cmsLHEtoEOSManager::l, HcalDDDRecConstants::HcalEtaBin::layer, hgcalTopologyTester_cfi::layers, maxLayer_, maxLayerHB_, dqmiodumpmetadata::n, HcalParameters::noff, phibin, HcalDDDRecConstants::HcalEtaBin::phis, edm::second(), and ecaldqm::zside().

Referenced by getEtaBins().

◆ getPhiBin()

double HcalDDDRecConstants::getPhiBin ( const int & i ) const

inline

Definition at line 94 of file HcalDDDRecConstants.h.

94 { return phibin[i]; }

References mps_fire::i, and phibin.

◆ getPhiOff()

double HcalDDDRecConstants::getPhiOff ( const int & i ) const

inline

Definition at line 95 of file HcalDDDRecConstants.h.

95 { return hpar->phioff[i]; }

References hpar, mps_fire::i, and HcalParameters::phioff.

◆ getPhiOffs()

const std::vector<double>& HcalDDDRecConstants::getPhiOffs ( ) const

inline

Definition at line 96 of file HcalDDDRecConstants.h.

96 { return hpar->phioff; }

References hpar, and HcalParameters::phioff.

Referenced by HcalTopology::HcalTopology().

◆ getPhis()

std::vector< std::pair< int, double > > HcalDDDRecConstants::getPhis	(	const int &	subdet,
		const int &	ieta
	)		const

Definition at line 373 of file HcalDDDRecConstants.cc.

                                                                                                      {
   std::vector<std::pair<int, double>> phis;
   int ietaAbs = (ieta > 0) ? ieta : -ieta;
   int keta = (subdet != HcalForward) ? etaSimValu[ietaAbs - 1].first : ietaAbs;
   std::pair<double, double> ficons = hcons.getPhiCons(subdet, keta);
   double fioff = ficons.first;
   double dphi = (subdet != HcalForward) ? phibin[ietaAbs - 1] : ficons.second;
   int nphi = int((2._pi + 0.1 * dphi) / dphi);
   int units = hcons.unitPhi(subdet, keta);
   for (int ifi = 0; ifi < nphi; ++ifi) {
     double phi = -fioff + (ifi + 0.5) * dphi;
     int iphi = hcons.phiNumber(ifi + 1, units);
     phis.emplace_back(std::pair<int, double>(iphi, phi));
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getEtaPhi: subdet|ieta|iphi " << subdet << "|" << ieta << " with " << phis.size()
                                << " phi bins";
   for (unsigned int k = 0; k < phis.size(); ++k)
     edm::LogVerbatim("HcalGeom") << "[" << k << "] iphi " << phis[k].first << " phi "
                                  << convertRadToDeg(phis[k].second);
 #endif
   return phis;
 }

References angle_units::operators::convertRadToDeg(), etaSimValu, dqmdumpme::first, HcalDDDSimConstants::getPhiCons(), HcalForward, hcons, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, dqmdumpme::k, nphi, phi, phibin, HcalDDDSimConstants::phiNumber(), edm::second(), HcalDDDSimConstants::unitPhi(), and units().

Referenced by HcalTestNumberingTester::analyze(), and getEtaBins().

◆ getPhiTable()

const std::vector<double>& HcalDDDRecConstants::getPhiTable ( ) const

inline

Definition at line 98 of file HcalDDDRecConstants.h.

98 { return phibin; }

References phibin.

Referenced by HcalTopology::HcalTopology().

◆ getPhiTableHF()

const std::vector<double>& HcalDDDRecConstants::getPhiTableHF ( ) const

inline

Definition at line 99 of file HcalDDDRecConstants.h.

99 { return hpar->phitable; }

References hpar, and HcalParameters::phitable.

Referenced by HcalTopology::HcalTopology().

◆ getPhiZOne()

int HcalDDDRecConstants::getPhiZOne ( std::vector< std::pair< int, int >> & phiz ) const

Definition at line 397 of file HcalDDDRecConstants.cc.

                                                                             {
   phiz.clear();
   int subdet = hcons.ldMap()->getSubdet();
   if (subdet > 0) {
     std::vector<int> phis = hcons.ldMap()->getPhis();
     for (int k : phis) {
       int zside = (k > 0) ? 1 : -1;
       int phi = (k > 0) ? k : -k;
       phiz.emplace_back(std::pair<int, int>(phi, zside));
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "Special RBX for detector " << subdet << " with " << phiz.size() << " phi/z bins";
   for (unsigned int k = 0; k < phiz.size(); ++k)
     edm::LogVerbatim("HcalGeom") << " [" << k << "] " << phiz[k].first << ":" << phiz[k].second;
 #endif
   return subdet;
 }

References HcalLayerDepthMap::getPhis(), HcalLayerDepthMap::getSubdet(), hcons, dqmdumpme::k, HcalDDDSimConstants::ldMap(), phi, and ecaldqm::zside().

Referenced by HcalTopology::getPhiZOne().

◆ getRZ() [1/4]

std::pair< double, double > HcalDDDRecConstants::getRZ ( const HcalDetId & id ) const

Definition at line 441 of file HcalDDDRecConstants.cc.

                                                                             {
   int subdet = id.subdetId();
   int ieta = id.ieta();
   int iphi = id.iphi();
   int depth = id.depth();
   int zside = (subdet == static_cast<int>(HcalBarrel)) ? 1 : id.zside();
   int layf = getLayerFront(subdet, ieta, iphi, depth);
   double rzf = (layf < 0)
                    ? 0.0
                    : ((subdet == static_cast<int>(HcalBarrel)) ? zside * (gconsHB[layf].first - gconsHB[layf].second)
                                                                : zside * (gconsHE[layf].first - gconsHE[layf].second));
   int layb = getLayerBack(subdet, ieta, iphi, depth);
   double rzb = (layb < 0)
                    ? 0.0
                    : ((subdet == static_cast<int>(HcalBarrel)) ? zside * (gconsHB[layb].first + gconsHB[layb].second)
                                                                : zside * (gconsHE[layb].first + gconsHE[layb].second));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getRZ: subdet|ieta|ipho|depth " << subdet << "|" << ieta << "|" << iphi << "|"
                                << depth << " lay|rz (front) " << layf << "|" << rzf << " lay|rz (back) " << layb << "|"
                                << rzb;
 #endif
   return std::pair<double, double>(rzf, rzb);
 }

References LEDCalibrationChannels::depth, dqmdumpme::first, gconsHB, gconsHE, getLayerBack(), getLayerFront(), HcalBarrel, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, edm::second(), and ecaldqm::zside().

◆ getRZ() [2/4]

double HcalDDDRecConstants::getRZ	(	const int &	subdet,
		const int &	ieta,
		const int &	depth
	)		const

Definition at line 416 of file HcalDDDRecConstants.cc.

                                                                                             {
   return getRZ(subdet, ieta, 1, depth);
 }

References LEDCalibrationChannels::depth, and LEDCalibrationChannels::ieta.

Referenced by HGCalSimHitStudy::analyzeHits(), HGCalSimHitValidation::analyzeHits(), HGCalTriggerTools::getLayerZ(), HcalHBHEMuonAnalyzer::goodCell(), and HcalHBHEMuonHighEtaAnalyzer::goodCell().

◆ getRZ() [3/4]

double HcalDDDRecConstants::getRZ	(	const int &	subdet,
		const int &	ieta,
		const int &	iphi,
		const int &	depth
	)		const

Definition at line 420 of file HcalDDDRecConstants.cc.

                                                                                                              {
   int layf = getLayerFront(subdet, ieta, iphi, depth);
   double rz =
       (layf < 0) ? 0.0 : ((subdet == static_cast<int>(HcalBarrel)) ? (gconsHB[layf].first) : (gconsHE[layf].first));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getRZ: subdet|ieta|ipho|depth " << subdet << "|" << ieta << "|" << iphi << "|"
                                << depth << " lay|rz " << layf << "|" << rz;
 #endif
   return rz;
 }

References LEDCalibrationChannels::depth, dqmdumpme::first, gconsHB, gconsHE, getLayerFront(), HcalBarrel, LEDCalibrationChannels::ieta, and LEDCalibrationChannels::iphi.

◆ getRZ() [4/4]

double HcalDDDRecConstants::getRZ	(	const int &	subdet,
		const int &	layer
	)		const

Definition at line 431 of file HcalDDDRecConstants.cc.

                                                                            {
   double rz(0);
   if (layer > 0 && layer <= (int)(layerGroupSize(0)))
     rz = ((subdet == static_cast<int>(HcalBarrel)) ? (gconsHB[layer - 1].first) : (gconsHE[layer - 1].first));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "getRZ: subdet|layer " << subdet << "|" << layer << " rz " << rz;
 #endif
   return rz;
 }

References dqmdumpme::first, gconsHB, gconsHE, HcalBarrel, and layerGroupSize().

◆ getThickActive()

std::vector< HcalDDDRecConstants::HcalActiveLength > HcalDDDRecConstants::getThickActive ( const int & type ) const

Definition at line 465 of file HcalDDDRecConstants.cc.

                                                                                                         {
   std::vector<HcalDDDRecConstants::HcalActiveLength> actives;
   std::vector<HcalDDDRecConstants::HcalEtaBin> bins = getEtaBins(type);
 #ifdef EDM_ML_DEBUG
   unsigned int kount(0);
 #endif
   for (auto& bin : bins) {
     int ieta = bin.ieta;
     int zside = bin.zside;
     int stype = (bin.phis.size() > 4) ? 0 : 1;
     int layf = getLayerFront(type + 1, zside * ieta, bin.phis[0].first, bin.depthStart) + 1;
     int layl = hcons.getLastLayer(type + 1, zside * ieta) + 1;
     double eta = 0.5 * (bin.etaMin + bin.etaMax);
     double theta = 2 * atan(exp(-eta));
     double scale = 1.0 / ((type == 0) ? sin(theta) : cos(theta));
     int depth = bin.depthStart;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalGeom") << "Eta " << ieta << " zside " << zside << " depth " << depth << " Layers " << layf
                                  << ":" << layl << ":" << bin.layer.size();
     for (auto ll : bin.layer)
       edm::LogVerbatim("HcalGeom") << "Layer " << ll.first << ":" << ll.second;
     for (auto phi : bin.phis)
       edm::LogVerbatim("HcalGeom") << "Phi " << phi.first << ":" << convertRadToDeg(phi.second);
 #endif
     for (unsigned int i = 0; i < bin.layer.size(); ++i) {
       double thick(0);
       int lmin = (type == 1 && ieta == iEtaMin[1]) ? layf : std::max(bin.layer[i].first, layf);
       int lmax = std::min(bin.layer[i].second, layl);
       for (int j = lmin; j <= lmax; ++j) {
         double t = ((type == 0) ? gconsHB[j - 1].second : gconsHE[j - 1].second);
         if ((type == 1) && (ieta <= 18))
           t = gconsHE[j].second;
         if (t > 0)
           thick += t;
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalGeom") << "Type " << type << " L " << lmin << ":" << lmax << " T " << thick;
 #endif
       thick *= (2. * scale);
       HcalDDDRecConstants::HcalActiveLength active(ieta, depth, zside, stype, zside * eta, thick);
       for (auto phi : bin.phis)
         active.iphis.emplace_back(phi.first);
       actives.emplace_back(active);
       ++depth;
 #ifdef EDM_ML_DEBUG
       kount++;
       edm::LogVerbatim("HcalGeom") << "getThickActive: [" << kount << "] eta:" << active.ieta << ":" << active.eta
                                    << " zside " << active.zside << " depth " << active.depth << " type " << active.stype
                                    << " thick " << active.thick;
 #endif
     }
   }
   return actives;
 }

References newFWLiteAna::bin, trigObjTnPSource_cfi::bins, angle_units::operators::convertRadToDeg(), funct::cos(), HcalDDDRecConstants::HcalActiveLength::depth, LEDCalibrationChannels::depth, HcalDDDRecConstants::HcalActiveLength::eta, PVValHelper::eta, JetChargeProducer_cfi::exp, gconsHB, gconsHE, getEtaBins(), HcalDDDSimConstants::getLastLayer(), getLayerFront(), hcons, mps_fire::i, HcalDDDRecConstants::HcalActiveLength::ieta, LEDCalibrationChannels::ieta, iEtaMin, HcalDDDRecConstants::HcalActiveLength::iphis, dqmiolumiharvest::j, SiStripPI::max, min(), phi, L1EGammaCrystalsEmulatorProducer_cfi::scale, edm::second(), funct::sin(), HcalDDDRecConstants::HcalActiveLength::stype, submitPVValidationJobs::t, theta(), HcalDDDRecConstants::HcalActiveLength::thick, HcalDDDRecConstants::HcalActiveLength::zside, and ecaldqm::zside().

Referenced by HcalHBHEMuonSimAnalyzer::beginRun(), HcalRaddamMuon::beginRun(), HcalHBHEMuonAnalyzer::beginRun(), and HcalHBHEMuonHighEtaAnalyzer::beginRun().

◆ getTopoMode()

int HcalDDDRecConstants::getTopoMode ( ) const

inline

Definition at line 106 of file HcalDDDRecConstants.h.

106 { return ((hpar->topologyMode) & 0xFF); }

References hpar, and HcalParameters::topologyMode.

Referenced by HcalTopology::HcalTopology(), and HcalTopologyIdealEP::produce().

◆ getTriggerMode()

int HcalDDDRecConstants::getTriggerMode ( ) const

inline

Definition at line 107 of file HcalDDDRecConstants.h.

107 { return (((hpar->topologyMode) >> 8) & 0xFF); }

References hpar, and HcalParameters::topologyMode.

Referenced by HcalTopology::HcalTopology().

◆ HcalCellTypes()

std::vector< HcalCellType > HcalDDDRecConstants::HcalCellTypes ( HcalSubdetector subdet ) const

Definition at line 520 of file HcalDDDRecConstants.cc.

                                                                                        {
   if (subdet == HcalBarrel || subdet == HcalEndcap) {
     std::vector<HcalCellType> cells;
     int isub = (subdet == HcalBarrel) ? 0 : 1;
     std::vector<HcalDDDRecConstants::HcalEtaBin> etabins = getEtaBins(isub);
     std::vector<int> missPhi;
     for (const auto& etabin : etabins) {
       std::vector<HcalCellType> temp;
       std::vector<int> count;
       std::vector<double> dmin, dmax;
       for (unsigned int il = 0; il < etabin.layer.size(); ++il) {
         HcalCellType cell(subdet, etabin.ieta, etabin.zside, 0, HcalCellType::HcalCell());
         temp.emplace_back(cell);
         count.emplace_back(0);
         dmin.emplace_back(0);
         dmax.emplace_back(0);
       }
       int ieta = etabin.ieta;
       for (int keta = etaSimValu[ieta - 1].first; keta <= etaSimValu[ieta - 1].second; ++keta) {
         std::vector<HcalCellType> cellsm = hcons.HcalCellTypes(subdet, keta, -1);
         for (unsigned int il = 0; il < etabin.layer.size(); ++il) {
           for (auto& ic : cellsm) {
             if (ic.depthSegment() >= etabin.layer[il].first && ic.depthSegment() <= etabin.layer[il].second &&
                 ic.etaBin() == temp[il].etaBin() && ic.zside() == temp[il].zside()) {
               if (count[il] == 0) {
                 temp[il] = ic;
                 dmin[il] = ic.depthMin();
                 dmax[il] = ic.depthMax();
               }
               ++count[il];
               if (ic.depthMin() < dmin[il])
                 dmin[il] = ic.depthMin();
               if (ic.depthMax() > dmax[il])
                 dmax[il] = ic.depthMax();
             }
           }
         }
       }
       for (unsigned int il = 0; il < etabin.layer.size(); ++il) {
         int depth = etabin.depthStart + (int)(il);
         temp[il].setEta(ieta, etabin.etaMin, etabin.etaMax);
         temp[il].setDepth(depth, dmin[il], dmax[il]);
         double foff = (etabin.ieta <= iEtaMax[0]) ? hpar->phioff[0] : hpar->phioff[1];
         int unit = hcons.unitPhi(etabin.dphi);
         temp[il].setPhi(etabin.phis, missPhi, foff, etabin.dphi, unit);
         cells.emplace_back(temp[il]);
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants: found " << cells.size() << " cells for sub-detector type "
                                  << isub;
     for (unsigned int ic = 0; ic < cells.size(); ++ic)
       edm::LogVerbatim("HcalGeom") << "Cell[" << ic << "] " << cells[ic];
 #endif
     return cells;
   } else {
     return hcons.HcalCellTypes(subdet, -1, -1);
   }
 }

References postprocess-scan-build::cells, submitPVResolutionJobs::count, LEDCalibrationChannels::depth, defaults_cfi::etabins, etaSimValu, dqmdumpme::first, getEtaBins(), HcalBarrel, HcalDDDSimConstants::HcalCellTypes(), HcalEndcap, hcons, hpar, LEDCalibrationChannels::ieta, iEtaMax, createfilelist::int, HcalParameters::phioff, groupFilesInBlocks::temp, unit(), and HcalDDDSimConstants::unitPhi().

Referenced by HcalDDDGeometryLoader::fill(), and numberOfCells().

◆ idBack()

HcalDetId HcalDDDRecConstants::idBack ( const HcalDetId & id ) const

Definition at line 663 of file HcalDDDRecConstants.cc.

                                                                {
   HcalDetId hid(id);
   std::map<HcalDetId, std::vector<HcalDetId>>::const_iterator itr = detIdSpR_.find(id);
   if (itr != detIdSpR_.end())
     hid = HcalDetId(id.subdet(), id.ieta(), id.iphi(), (itr->second).back().depth());
   return hid;
 }

References detIdSpR_.

Referenced by HcalTopology::idBack().

◆ idFront()

HcalDetId HcalDDDRecConstants::idFront ( const HcalDetId & id ) const

Definition at line 655 of file HcalDDDRecConstants.cc.

                                                                 {
   HcalDetId hid(id);
   std::map<HcalDetId, std::vector<HcalDetId>>::const_iterator itr = detIdSpR_.find(id);
   if (itr != detIdSpR_.end())
     hid = HcalDetId(id.subdet(), id.ieta(), id.iphi(), (itr->second)[0].depth());
   return hid;
 }

References LEDCalibrationChannels::depth, and detIdSpR_.

Referenced by HcalTopology::idFront().

◆ initialize()

void HcalDDDRecConstants::initialize ( void )

private

Definition at line 801 of file HcalDDDRecConstants.cc.

                                          {
   //Eta grouping
   int nEta = (int)(hpar->etagroup.size());
   if (nEta != (int)(hpar->phigroup.size())) {
     edm::LogError("HCalGeom") << "HcalDDDRecConstants: sizes of the vectors "
                               << " etaGroup (" << nEta << ") and phiGroup (" << hpar->phigroup.size()
                               << ") do not match";
     throw cms::Exception("DDException") << "HcalDDDRecConstants: inconsistent array sizes" << nEta << ":"
                                         << hpar->phigroup.size();
   }
  
   // First eta table
   iEtaMin = hpar->etaMin;
   iEtaMax = hpar->etaMax;
   etaTable.clear();
   ietaMap.clear();
   etaSimValu.clear();
   int ieta(0), ietaHB(0), ietaHE(0), ietaHEM(0);
   etaTable.emplace_back(hpar->etaTable[ieta]);
   for (int i = 0; i < nEta; ++i) {
     int ef = ieta + 1;
     ieta += (hpar->etagroup[i]);
     if (ieta >= (int)(hpar->etaTable.size())) {
       edm::LogError("HCalGeom") << "HcalDDDRecConstants: Going beyond the array boundary " << hpar->etaTable.size()
                                 << " at index " << i << " of etaTable from SimConstant";
       throw cms::Exception("DDException")
           << "HcalDDDRecConstants: Going beyond the array boundary " << hpar->etaTable.size() << " at index " << i
           << " of etaTable from SimConstant";
     } else {
       etaTable.emplace_back(hpar->etaTable[ieta]);
       etaSimValu.emplace_back(std::pair<int, int>(ef, ieta));
     }
     for (int k = 0; k < (hpar->etagroup[i]); ++k)
       ietaMap.emplace_back(i + 1);
     if (ieta <= hpar->etaMax[0])
       ietaHB = i + 1;
     if (ieta <= hpar->etaMin[1])
       ietaHE = i + 1;
     if (ieta <= hpar->etaMax[1])
       ietaHEM = i + 1;
   }
   iEtaMin[1] = ietaHE;
   iEtaMax[0] = ietaHB;
   iEtaMax[1] = ietaHEM;
  
   // Then Phi bins
   nPhiBins.clear();
   for (unsigned int k = 0; k < 4; ++k)
     nPhiBins.emplace_back(0);
   ieta = 0;
   phibin.clear();
   phiUnitS.clear();
   for (int i = 0; i < nEta; ++i) {
     double dphi = (hpar->phigroup[i]) * (hpar->phibin[ieta]);
     phibin.emplace_back(dphi);
     int nphi = (int)((2._pi + 0.001) / dphi);
     if (ieta <= iEtaMax[0]) {
       if (nphi > nPhiBins[0])
         nPhiBins[3] = nPhiBins[0] = nphi;
     }
     if (ieta >= iEtaMin[1]) {
       if (nphi > nPhiBins[1])
         nPhiBins[1] = nphi;
     }
     ieta += (hpar->etagroup[i]);
   }
   for (unsigned int i = 1; i < hpar->etaTable.size(); ++i) {
     int unit = hcons.unitPhi(hpar->phibin[i - 1]);
     phiUnitS.emplace_back(unit);
   }
   for (double i : hpar->phitable) {
     int nphi = (int)((2._pi + 0.001) / i);
     if (nphi > nPhiBins[2])
       nPhiBins[2] = nphi;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants: Modified eta/deltaphi table for " << nEta << " bins";
   for (int i = 0; i < nEta; ++i)
     edm::LogVerbatim("HcalGeom") << "Eta[" << i << "] = " << etaTable[i] << ":" << etaTable[i + 1] << ":"
                                  << etaSimValu[i].first << ":" << etaSimValu[i].second << " PhiBin[" << i
                                  << "] = " << convertRadToDeg(phibin[i]);
   for (unsigned int i = 0; i < phiUnitS.size(); ++i)
     edm::LogVerbatim("HcalGeom") << " PhiUnitS[" << i << "] = " << phiUnitS[i];
   for (unsigned int i = 0; i < nPhiBins.size(); ++i)
     edm::LogVerbatim("HcalGeom") << " nPhiBins[" << i << "] = " << nPhiBins[i];
   for (unsigned int i = 0; i < hpar->etaTableHF.size(); ++i)
     edm::LogVerbatim("HcalGeom") << " EtaTableHF[" << i << "] = " << hpar->etaTableHF[i];
   for (unsigned int i = 0; i < hpar->phitable.size(); ++i)
     edm::LogVerbatim("HcalGeom") << " PhiBinHF[" << i << "] = " << hpar->phitable[i];
 #endif
  
   //Now the depths
   maxDepth = hpar->maxDepth;
   maxDepth[0] = maxDepth[1] = 0;
   for (int i = 0; i < nEta; ++i) {
     unsigned int imx = layerGroupSize(i);
     int laymax = (imx > 0) ? layerGroup(i, imx - 1) : 0;
     if (i < iEtaMax[0]) {
       int laymax0 = (imx > 16) ? layerGroup(i, 16) : laymax;
       if (i + 1 == iEtaMax[0])
         laymax0 = hcons.getDepthEta16M(1);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:HB " << i << " " << imx << " " << laymax << " " << laymax0;
 #endif
       if (maxDepth[0] < laymax0)
         maxDepth[0] = laymax0;
     }
     if (i >= iEtaMin[1] - 1 && i < iEtaMax[1]) {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:HE " << i << " " << imx << " " << laymax;
 #endif
       if (maxDepth[1] < laymax)
         maxDepth[1] = laymax;
     }
   }
 #ifdef EDM_ML_DEBUG
   for (int i = 0; i < 4; ++i)
     edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:Detector Type[" << i << "] iEta " << iEtaMin[i] << ":"
                                  << iEtaMax[i] << " MaxDepth " << maxDepth[i];
 #endif
  
   //Now the geometry constants
   nModule[0] = hpar->modHB[0];
   nHalves[0] = hpar->modHB[1];
   for (unsigned int i = 0; i < hpar->rHB.size(); ++i) {
     gconsHB.emplace_back(std::pair<double, double>(hpar->rHB[i] / CLHEP::cm, hpar->drHB[i] / CLHEP::cm));
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:HB with " << nModule[0] << " modules and " << nHalves[0]
                                << " halves and " << gconsHB.size() << " layers";
   for (unsigned int i = 0; i < gconsHB.size(); ++i)
     edm::LogVerbatim("HcalGeom") << "rHB[" << i << "] = " << gconsHB[i].first << " +- " << gconsHB[i].second;
 #endif
   nModule[1] = hpar->modHE[0];
   nHalves[1] = hpar->modHE[1];
   for (unsigned int i = 0; i < hpar->zHE.size(); ++i) {
     gconsHE.emplace_back(std::pair<double, double>(hpar->zHE[i] / CLHEP::cm, hpar->dzHE[i] / CLHEP::cm));
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:HE with " << nModule[1] << " modules and " << nHalves[1]
                                << " halves and " << gconsHE.size() << " layers";
   for (unsigned int i = 0; i < gconsHE.size(); ++i)
     edm::LogVerbatim("HcalGeom") << "zHE[" << i << "] = " << gconsHE[i].first << " +- " << gconsHE[i].second;
 #endif
  
   //Special RBX
   depthMaxSp_ = hcons.getMaxDepthDet(0);
   if (depthMaxSp_.first == 0) {
     depthMaxSp_ = depthMaxDf_ = std::pair<int, int>(2, maxDepth[1]);
   } else if (depthMaxSp_.first == 1) {
     depthMaxDf_ = std::pair<int, int>(1, maxDepth[0]);
     if (depthMaxSp_.second > maxDepth[0])
       maxDepth[0] = depthMaxSp_.second;
   } else {
     depthMaxDf_ = std::pair<int, int>(2, maxDepth[1]);
     if (depthMaxSp_.second > maxDepth[1])
       maxDepth[1] = depthMaxSp_.second;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:Detector type and maximum depth for all RBX "
                                << depthMaxDf_.first << ":" << depthMaxDf_.second << " and for special RBX "
                                << depthMaxSp_.first << ":" << depthMaxSp_.second;
 #endif
  
   //Map of special DetId's
   std::vector<int> phis;
   HcalSubdetector subdet = HcalSubdetector(hcons.ldMap()->validDet(phis));
   detIdSp_.clear();
   detIdSpR_.clear();
   if ((subdet == HcalBarrel) || (subdet == HcalEndcap)) {
     int phi = (phis[0] > 0) ? phis[0] : -phis[0];
     int zside = (phis[0] > 0) ? 1 : -1;
     int lymax = (subdet == HcalBarrel) ? maxLayerHB_ + 1 : maxLayer_ + 1;
     std::pair<int, int> etas = hcons.ldMap()->validEta();
     for (int eta = etas.first; eta <= etas.second; ++eta) {
       std::map<int, std::pair<int, int>> oldDep;
       int depth(0);
       int lmin = layerGroup(eta - 1, 0);
       for (int lay = 0; lay < lymax; ++lay) {
         int depc = layerGroup(eta - 1, lay);
         if (depth != depc) {
           if (depth != 0)
             oldDep[depth] = std::pair<int, int>(lmin, lay - 1);
           depth = depc;
           lmin = lay;
         }
       }
       if (depth != 0)
         oldDep[depth] = std::pair<int, int>(lmin, lymax - 1);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:Eta|Phi|Zside " << eta << ":" << phi << ":" << zside
                                    << " with " << oldDep.size() << " old Depths";
       unsigned int kk(0);
       for (std::map<int, std::pair<int, int>>::const_iterator itr = oldDep.begin(); itr != oldDep.end(); ++itr, ++kk)
         edm::LogVerbatim("HcalGeom") << "[" << kk << "] " << itr->first << " --> " << itr->second.first << ":"
                                      << itr->second.second;
 #endif
       std::pair<int, int> depths = hcons.ldMap()->getDepths(eta);
       for (int ndepth = depths.first; ndepth <= depths.second; ++ndepth) {
         bool flag = ((subdet == HcalBarrel && eta == iEtaMax[0] && ndepth > hcons.getDepthEta16(subdet, phi, zside)) ||
                      (subdet == HcalEndcap && eta == iEtaMin[1] && ndepth < hcons.getDepthEta16(subdet, phi, zside)));
         if (!flag) {
           std::vector<int> count(oldDep.size(), 0);
           int layFront = hcons.ldMap()->getLayerFront(subdet, eta, phi, zside, ndepth);
           int layBack = hcons.ldMap()->getLayerBack(subdet, eta, phi, zside, ndepth);
           for (int lay = layFront; lay <= layBack; ++lay) {
             unsigned int l(0);
             for (std::map<int, std::pair<int, int>>::iterator itr = oldDep.begin(); itr != oldDep.end(); ++itr, ++l) {
               if (lay >= (itr->second).first && lay <= (itr->second).second) {
                 ++count[l];
                 break;
               }
             }
           }
           int odepth(0), maxlay(0);
           unsigned int l(0);
           for (std::map<int, std::pair<int, int>>::iterator itr = oldDep.begin(); itr != oldDep.end(); ++itr, ++l) {
             if (count[l] > maxlay) {
               odepth = itr->first;
               maxlay = count[l];
             }
           }
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:New Depth " << ndepth << " old Depth " << odepth
                                        << " max " << maxlay;
 #endif
           for (int k : phis) {
             zside = (k > 0) ? 1 : -1;
             phi = (k > 0) ? k : -k;
             if (subdet == HcalEndcap && eta == hpar->noff[1] && ndepth > hcons.getDepthEta29M(0, true))
               break;
             HcalDetId newId(subdet, zside * eta, phi, ndepth);
             HcalDetId oldId(subdet, zside * eta, phi, odepth);
             detIdSp_[newId] = oldId;
             std::vector<HcalDetId> ids;
             std::map<HcalDetId, std::vector<HcalDetId>>::iterator itr = detIdSpR_.find(oldId);
             if (itr != detIdSpR_.end())
               ids = itr->second;
             ids.emplace_back(newId);
             detIdSpR_[oldId] = ids;
           }
         }
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:Map for merging new channels to old channel"
                                  << " IDs with " << detIdSp_.size() << " entries";
     int l(0);
     for (auto itr : detIdSp_) {
       edm::LogVerbatim("HcalGeom") << "[" << l << "] Special " << itr.first << " Standard " << itr.second;
       ++l;
     }
     edm::LogVerbatim("HcalGeom") << "HcalDDDRecConstants:Reverse Map for mapping old to new IDs with "
                                  << detIdSpR_.size() << " entries";
     l = 0;
     for (auto itr : detIdSpR_) {
       edm::LogVerbatim("HcalGeom") << "[" << l << "] Standard " << itr.first << " Special";
       for (auto itr1 : itr.second)
         edm::LogVerbatim("HcalGeom") << "ID " << (itr1);
       ++l;
     }
 #endif
   }
 }

References angle_units::operators::convertRadToDeg(), submitPVResolutionJobs::count, LEDCalibrationChannels::depth, depthMaxDf_, depthMaxSp_, HLT_FULL_cff::depths, detIdSp_, detIdSpR_, HcalParameters::drHB, HcalParameters::dzHE, PVValHelper::eta, HcalParameters::etagroup, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, HcalParameters::etaMin, EnergyCorrector::etas, etaSimValu, HcalParameters::etaTable, etaTable, HcalParameters::etaTableHF, Exception, RemoveAddSevLevel::flag, gconsHB, gconsHE, HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29M(), HcalLayerDepthMap::getDepths(), HcalLayerDepthMap::getLayerBack(), HcalLayerDepthMap::getLayerFront(), HcalDDDSimConstants::getMaxDepthDet(), HcalBarrel, HcalEndcap, hcons, hpar, mps_fire::i, LEDCalibrationChannels::ieta, ietaMap, iEtaMax, iEtaMin, createfilelist::int, dqmdumpme::k, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::l, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), genParticles_cff::map, HcalParameters::maxDepth, maxDepth, maxLayer_, maxLayerHB_, HcalParameters::modHB, HcalParameters::modHE, ndepth, HLT_FULL_cff::nEta, nHalves, nModule, HcalParameters::noff, nphi, nPhiBins, phi, HcalParameters::phibin, phibin, HcalParameters::phigroup, HcalParameters::phitable, phiUnitS, HcalParameters::rHB, edm::second(), unit(), HcalDDDSimConstants::unitPhi(), HcalLayerDepthMap::validDet(), HcalLayerDepthMap::validEta(), HcalParameters::zHE, and ecaldqm::zside().

Referenced by HcalDDDRecConstants().

◆ isBH()

bool HcalDDDRecConstants::isBH ( ) const

inline

Definition at line 109 of file HcalDDDRecConstants.h.

109 { return hcons.isBH(); }

References hcons, and HcalDDDSimConstants::isBH().

Referenced by HcalDDDGeometryLoader::HcalDDDGeometryLoader(), HcalTopology::isBH(), and HcalFlexiHardcodeGeometryLoader::load().

◆ isPlan1()

bool HcalDDDRecConstants::isPlan1 ( const HcalDetId & id ) const

inline

Definition at line 110 of file HcalDDDRecConstants.h.

110 { return detIdSp_.find(id) != detIdSp_.end(); };

References detIdSp_.

Referenced by CaloTPGTranscoderULUT::setup(), and HcaluLUTTPGCoder::update().

◆ isPlan1MergedId()

bool HcalDDDRecConstants::isPlan1MergedId ( const HcalDetId & id ) const

inline

Definition at line 125 of file HcalDDDRecConstants.h.

125 { return detIdSpR_.find(id) != detIdSpR_.end(); };

References detIdSpR_.

Referenced by HcalTopology::validHcal().

◆ isPlan1ToBeMergedId()

bool HcalDDDRecConstants::isPlan1ToBeMergedId ( const HcalDetId & id ) const

inline

Definition at line 124 of file HcalDDDRecConstants.h.

124 { return detIdSp_.find(id) != detIdSp_.end(); };

References detIdSp_.

Referenced by HcalTopology::validHcal().

◆ layerGroup()

unsigned int HcalDDDRecConstants::layerGroup	(	int	eta,
		int	i
	)		const

private

Definition at line 1079 of file HcalDDDRecConstants.cc.

                                                                  {
   unsigned int k = 0;
   for (auto const& it : hpar->layerGroupEtaRec) {
     if (it.layer == (unsigned int)(eta + 1)) {
       return it.layerGroup.at(i);
     }
     if (it.layer > (unsigned int)(eta + 1))
       break;
     k = it.layerGroup.at(i);
   }
   return k;
 }

References PVValHelper::eta, hpar, mps_fire::i, dqmdumpme::k, and HcalParameters::layerGroupEtaRec.

Referenced by getHCID(), getLayerBack(), getLayerDepth(), getLayerFront(), getMaxDepth(), getMinDepth(), and initialize().

◆ layerGroupSize()

unsigned int HcalDDDRecConstants::layerGroupSize ( int eta ) const

private

Definition at line 1066 of file HcalDDDRecConstants.cc.

                                                               {
   unsigned int k = 0;
   for (auto const& it : hpar->layerGroupEtaRec) {
     if (it.layer == (unsigned int)(eta + 1)) {
       return it.layerGroup.size();
     }
     if (it.layer > (unsigned int)(eta + 1))
       break;
     k = it.layerGroup.size();
   }
   return k;
 }

References PVValHelper::eta, hpar, dqmdumpme::k, and HcalParameters::layerGroupEtaRec.

Referenced by getLayerBack(), getLayerDepth(), getLayerFront(), getMaxDepth(), getMinDepth(), getRZ(), and initialize().

◆ maxHFDepth()

int HcalDDDRecConstants::maxHFDepth	(	int	ieta,
		int	iphi
	)		const

inline

Definition at line 111 of file HcalDDDRecConstants.h.

111 { return hcons.maxHFDepth(ieta, iphi); }

References hcons, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and HcalDDDSimConstants::maxHFDepth().

Referenced by HcalTopology::validRaw().

◆ mergedDepthDetId()

HcalDetId HcalDDDRecConstants::mergedDepthDetId ( const HcalDetId & id ) const

Definition at line 647 of file HcalDDDRecConstants.cc.

                                                                          {
   std::map<HcalDetId, HcalDetId>::const_iterator itr = detIdSp_.find(id);
   if (itr == detIdSp_.end())
     return id;
   else
     return itr->second;
 }

References detIdSp_.

Referenced by HcalHBHEMuonAnalyzer::analyze(), HcalHBHEMuonHighEtaAnalyzer::analyzeTracks(), HcalHBHEMuonAnalyzer::beginRun(), HcalHBHEMuonHighEtaAnalyzer::beginRun(), HcalTopology::mergedDepthDetId(), and HcalIsoTrkAnalyzer::storeEnergy().

◆ mergedDepthList29() [1/2]

std::vector< int > HcalDDDRecConstants::mergedDepthList29	(	int	ieta,
		int	iphi
	)		const

Definition at line 593 of file HcalDDDRecConstants.cc.

                                                                               {
   std::vector<int> depths;
   int eta = std::abs(ieta);
   int zside = (ieta > 0) ? 1 : -1;
   int etamin = iEtaMax[1] - hcons.getDepthEta29(iphi, zside, 1);
   if ((eta >= etamin) && (eta <= iEtaMax[1])) {
     int depthMax = getMaxDepth(1, etamin, iphi, zside);
     int depthMin = hcons.getDepthEta29(iphi, zside, 0) + 1;
     depths.reserve(depthMax - depthMin + 1);
     for (int depth = depthMin; depth <= depthMax; ++depth)
       depths.emplace_back(depth);
   }
   return depths;
 }

References funct::abs(), LEDCalibrationChannels::depth, HLT_FULL_cff::depths, PVValHelper::eta, muonTiming_cfi::etamin, HcalDDDSimConstants::getDepthEta29(), getMaxDepth(), hcons, LEDCalibrationChannels::ieta, iEtaMax, LEDCalibrationChannels::iphi, and ecaldqm::zside().

◆ mergedDepthList29() [2/2]

bool HcalDDDRecConstants::mergedDepthList29	(	int	ieta,
		int	iphi,
		int	depth
	)		const

Definition at line 580 of file HcalDDDRecConstants.cc.

                                                                                {
   int eta = std::abs(ieta);
   int zside = (ieta > 0) ? 1 : -1;
   int etamin = iEtaMax[1] - hcons.getDepthEta29(iphi, zside, 1);
   if ((eta >= etamin) && (eta <= iEtaMax[1])) {
     int depthMax = getMaxDepth(1, etamin, iphi, zside);
     int depthMin = hcons.getDepthEta29(iphi, zside, 0) + 1;
     if (depth >= depthMin && depth <= depthMax)
       return true;
   }
   return false;
 }

References funct::abs(), LEDCalibrationChannels::depth, PVValHelper::eta, muonTiming_cfi::etamin, HcalDDDSimConstants::getDepthEta29(), getMaxDepth(), hcons, LEDCalibrationChannels::ieta, iEtaMax, LEDCalibrationChannels::iphi, and ecaldqm::zside().

Referenced by HcalTopology::mergedDepth29(), and HcalTopology::mergedDepthList29().

◆ nCells() [1/2]

unsigned int HcalDDDRecConstants::nCells ( ) const

Definition at line 643 of file HcalDDDRecConstants.cc.

                                                {
   return (nCells(HcalBarrel) + nCells(HcalEndcap) + nCells(HcalOuter) + nCells(HcalForward));
 }

References HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

◆ nCells() [2/2]

unsigned int HcalDDDRecConstants::nCells ( HcalSubdetector subdet ) const

Definition at line 625 of file HcalDDDRecConstants.cc.

                                                                      {
   if (subdet == HcalBarrel || subdet == HcalEndcap) {
     int isub = (subdet == HcalBarrel) ? 0 : 1;
     std::vector<HcalDDDRecConstants::HcalEtaBin> etabins = getEtaBins(isub);
     unsigned int ncell(0);
     for (auto& etabin : etabins) {
       ncell += ((etabin.phis.size()) * (etabin.layer.size()));
     }
     return ncell;
   } else if (subdet == HcalOuter) {
     return kHOSizePreLS1;
   } else if (subdet == HcalForward) {
     return (unsigned int)(hcons.numberOfCells(subdet));
   } else {
     return 0;
   }
 }

References defaults_cfi::etabins, getEtaBins(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, kHOSizePreLS1, and HcalDDDSimConstants::numberOfCells().

◆ numberOfCells()

unsigned int HcalDDDRecConstants::numberOfCells ( HcalSubdetector subdet ) const

Definition at line 608 of file HcalDDDRecConstants.cc.

                                                                             {
   if (subdet == HcalBarrel || subdet == HcalEndcap) {
     unsigned int num = 0;
     std::vector<HcalCellType> cellTypes = HcalCellTypes(subdet);
     for (auto& cellType : cellTypes) {
       num += (unsigned int)(cellType.nPhiBins());
     }
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("HCalGeom") << "HcalDDDRecConstants:numberOfCells " << cellTypes.size() << " " << num
                              << " for subdetector " << subdet;
 #endif
     return num;
   } else {
     return hcons.numberOfCells(subdet);
   }
 }

References HcalBarrel, HcalCellTypes(), HcalEndcap, hcons, createfilelist::int, EgammaValidation_cff::num, and HcalDDDSimConstants::numberOfCells().

Referenced by HcalDDDGeometryLoader::load().

◆ specialRBXHBHE() [1/2]

bool HcalDDDRecConstants::specialRBXHBHE	(	bool	flag,
		std::vector< HcalDetId > &	ids
	)		const

Definition at line 694 of file HcalDDDRecConstants.cc.

                                                                                          {
   if (tobemerged) {
     std::map<HcalDetId, HcalDetId>::const_iterator itr;
     for (itr = detIdSp_.begin(); itr != detIdSp_.end(); ++itr)
       ids.emplace_back(itr->first);
   } else {
     std::map<HcalDetId, std::vector<HcalDetId>>::const_iterator itr;
     for (itr = detIdSpR_.begin(); itr != detIdSpR_.end(); ++itr)
       ids.emplace_back(itr->first);
   }
   return (!ids.empty());
 }

References detIdSp_, and detIdSpR_.

◆ specialRBXHBHE() [2/2]

void HcalDDDRecConstants::specialRBXHBHE	(	const std::vector< HcalDetId > &	idsOld,
		std::vector< HcalDetId > &	idsNew
	)		const

Definition at line 684 of file HcalDDDRecConstants.cc.

                                                                                                                {
   for (auto k : idsOld) {
     std::map<HcalDetId, HcalDetId>::const_iterator itr = detIdSp_.find(k);
     if (itr == detIdSp_.end())
       idsNew.emplace_back(k);
     else
       idsNew.emplace_back(itr->second);
   }
 }

References detIdSp_, and dqmdumpme::k.

◆ unmergeDepthDetId()

void HcalDDDRecConstants::unmergeDepthDetId	(	const HcalDetId &	id,
		std::vector< HcalDetId > &	ids
	)		const

Definition at line 671 of file HcalDDDRecConstants.cc.

                                                                                                 {
   ids.clear();
   std::map<HcalDetId, std::vector<HcalDetId>>::const_iterator itr = detIdSpR_.find(id);
   if (itr == detIdSpR_.end()) {
     ids.emplace_back(id);
   } else {
     for (auto k : itr->second) {
       HcalDetId hid(id.subdet(), id.ieta(), id.iphi(), k.depth());
       ids.emplace_back(hid);
     }
   }
 }

References detIdSpR_, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and dqmdumpme::k.

Referenced by HcalHBHEMuonAnalyzer::activeLength(), HcalHBHEMuonHighEtaAnalyzer::activeLength(), and HcalTopology::unmergeDepthDetId().

◆ withSpecialRBXHBHE()

bool HcalDDDRecConstants::withSpecialRBXHBHE ( ) const

inline

Definition at line 123 of file HcalDDDRecConstants.h.

123 { return (hcons.ldMap()->getSubdet() != 0); }

References HcalLayerDepthMap::getSubdet(), hcons, and HcalDDDSimConstants::ldMap().

Referenced by HcalTopology::withSpecialRBXHBHE().

Member Data Documentation

◆ depthMaxDf_

std::pair<int, int> HcalDDDRecConstants::depthMaxDf_

private

Definition at line 155 of file HcalDDDRecConstants.h.

Referenced by initialize().

◆ depthMaxSp_

std::pair<int, int> HcalDDDRecConstants::depthMaxSp_

private

Definition at line 155 of file HcalDDDRecConstants.h.

Referenced by initialize().

◆ detIdSp_

std::map<HcalDetId, HcalDetId> HcalDDDRecConstants::detIdSp_

private

Definition at line 156 of file HcalDDDRecConstants.h.

Referenced by initialize(), isPlan1(), isPlan1ToBeMergedId(), mergedDepthDetId(), and specialRBXHBHE().

◆ detIdSpR_

std::map<HcalDetId, std::vector<HcalDetId> > HcalDDDRecConstants::detIdSpR_

private

Definition at line 157 of file HcalDDDRecConstants.h.

Referenced by idBack(), idFront(), initialize(), isPlan1MergedId(), specialRBXHBHE(), and unmergeDepthDetId().

◆ etaSimValu

std::vector<std::pair<int, int> > HcalDDDRecConstants::etaSimValu

private

Definition at line 144 of file HcalDDDRecConstants.h.

Referenced by getPhis(), HcalCellTypes(), and initialize().

◆ etaTable

std::vector<double> HcalDDDRecConstants::etaTable

private

Definition at line 145 of file HcalDDDRecConstants.h.

Referenced by getEtaLimit(), getEtaPhi(), getEtaTable(), getOneEtaBin(), and initialize().

◆ gconsHB

std::vector<std::pair<double, double> > HcalDDDRecConstants::gconsHB

private

Definition at line 152 of file HcalDDDRecConstants.h.

Referenced by getConstHBHE(), getRZ(), getThickActive(), and initialize().

◆ gconsHE

std::vector<std::pair<double, double> > HcalDDDRecConstants::gconsHE

private

Definition at line 153 of file HcalDDDRecConstants.h.

Referenced by getConstHBHE(), getRZ(), getThickActive(), and initialize().

◆ hcons

const HcalDDDSimConstants& HcalDDDRecConstants::hcons

private

Definition at line 143 of file HcalDDDRecConstants.h.

Referenced by dddConstants(), getDepth(), getDepthEta16(), getDepthEta29(), getEtaBins(), getEtaPhi(), getHCID(), getHFCellParameters(), getLayer0Wt(), getLayerBack(), getLayerFront(), getMaxDepth(), getMinDepth(), getOneEtaBin(), getPhis(), getPhiZOne(), getThickActive(), HcalCellTypes(), initialize(), isBH(), maxHFDepth(), mergedDepthList29(), nCells(), numberOfCells(), and withSpecialRBXHBHE().

◆ hpar

const HcalParameters* HcalDDDRecConstants::hpar

private

Definition at line 142 of file HcalDDDRecConstants.h.

Referenced by getDepth(), getEtaPhi(), getEtaTableHF(), getHCID(), getMaxDepth(), getNEta(), getNoff(), getOneEtaBin(), getPhiOff(), getPhiOffs(), getPhiTableHF(), getTopoMode(), getTriggerMode(), HcalCellTypes(), initialize(), layerGroup(), and layerGroupSize().

◆ ietaMap

std::vector<int> HcalDDDRecConstants::ietaMap

private

Definition at line 146 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().

◆ iEtaMax

std::vector<int> HcalDDDRecConstants::iEtaMax

private

Definition at line 147 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), getEtaPhi(), getEtaRange(), getMaxDepth(), getOneEtaBin(), HcalCellTypes(), initialize(), and mergedDepthList29().

◆ iEtaMin

std::vector<int> HcalDDDRecConstants::iEtaMin

private

Definition at line 147 of file HcalDDDRecConstants.h.

Referenced by getEtaBins(), getEtaPhi(), getEtaRange(), getHCID(), getHFCellParameters(), getMinDepth(), getOneEtaBin(), getThickActive(), and initialize().

◆ maxDepth

std::vector<int> HcalDDDRecConstants::maxDepth

private

Definition at line 148 of file HcalDDDRecConstants.h.

Referenced by getHFCellParameters(), getMaxDepth(), getMinDepth(), and initialize().

◆ maxLayer_

const int HcalDDDRecConstants::maxLayer_ = 18

staticprivate

Definition at line 140 of file HcalDDDRecConstants.h.

Referenced by getMaxDepth(), getOneEtaBin(), and initialize().

◆ maxLayerHB_

const int HcalDDDRecConstants::maxLayerHB_ = 16

staticprivate

Definition at line 141 of file HcalDDDRecConstants.h.

Referenced by getHCID(), getMaxDepth(), getOneEtaBin(), and initialize().

◆ nHalves

int HcalDDDRecConstants::nHalves[2]

private

Definition at line 154 of file HcalDDDRecConstants.h.

Referenced by initialize().

◆ nModule

int HcalDDDRecConstants::nModule[2]

private

Definition at line 154 of file HcalDDDRecConstants.h.

Referenced by initialize().

◆ nPhiBins

std::vector<int> HcalDDDRecConstants::nPhiBins

private

Definition at line 149 of file HcalDDDRecConstants.h.

Referenced by getNPhi(), and initialize().

◆ phibin

std::vector<double> HcalDDDRecConstants::phibin

private

Definition at line 150 of file HcalDDDRecConstants.h.

Referenced by getEtaPhi(), getHCID(), getOneEtaBin(), getPhiBin(), getPhis(), getPhiTable(), and initialize().

◆ phiUnitS

std::vector<int> HcalDDDRecConstants::phiUnitS

private

Definition at line 151 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HcalDDDRecConstants()

◆ ~HcalDDDRecConstants()

Member Function Documentation

◆ dddConstants()

◆ getConstHBHE()

◆ getDepth() [1/2]

◆ getDepth() [2/2]

◆ getDepthEta16()

◆ getDepthEta29()

◆ getEtaBins()

◆ getEtaLimit()

◆ getEtaPhi()

◆ getEtaRange()

◆ getEtaTable()

◆ getEtaTableHF()

◆ getHCID()

◆ getHFCellParameters()

◆ getLayer0Wt()

◆ getLayerBack()

◆ getLayerDepth()

◆ getLayerFront()

◆ getMaxDepth() [1/2]

◆ getMaxDepth() [2/2]

◆ getMinDepth()

◆ getNEta()

◆ getNoff()

◆ getNPhi()

◆ getOneEtaBin()

◆ getPhiBin()

◆ getPhiOff()

◆ getPhiOffs()

◆ getPhis()

◆ getPhiTable()

◆ getPhiTableHF()

◆ getPhiZOne()

◆ getRZ() [1/4]

◆ getRZ() [2/4]

◆ getRZ() [3/4]

◆ getRZ() [4/4]

◆ getThickActive()

◆ getTopoMode()

◆ getTriggerMode()

◆ HcalCellTypes()

◆ idBack()

◆ idFront()

◆ initialize()

◆ isBH()

◆ isPlan1()

◆ isPlan1MergedId()

◆ isPlan1ToBeMergedId()

◆ layerGroup()

◆ layerGroupSize()

◆ maxHFDepth()

◆ mergedDepthDetId()

◆ mergedDepthList29() [1/2]

◆ mergedDepthList29() [2/2]

◆ nCells() [1/2]

◆ nCells() [2/2]

◆ numberOfCells()

◆ specialRBXHBHE() [1/2]

◆ specialRBXHBHE() [2/2]

◆ unmergeDepthDetId()

◆ withSpecialRBXHBHE()

Member Data Documentation

◆ depthMaxDf_

◆ depthMaxSp_

◆ detIdSp_

◆ detIdSpR_

◆ etaSimValu

◆ etaTable

◆ gconsHB

◆ gconsHE

◆ hcons