#include <HcalDDDSimConstants.h>

Public Member Functions
HcalCellType::HcalCell	cell (const int &det, const int &zside, const int &depth, const int &etaR, const int &iphi) const

int	findDepth (const int &det, const int &eta, const int &phi, const int &zside, const int &lay) const

unsigned int	findLayer (const int &layer, const std::vector< HcalParameters::LayerItem > &layerGroup) const

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

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

int	getDepthEta16M (const int &det) const

int	getDepthEta29 (const int &phi, const int &zside, const int &i) const

int	getDepthEta29M (const int &i, const bool &planOne) const

std::pair< int, double >	getDetEta (const double &eta, const int &depth) const

int	getEta (const int &det, const int &lay, const double &hetaR) const

std::pair< int, int >	getEtaDepth (const int &det, int etaR, const int &phi, const int &zside, int depth, const int &lay) const

double	getEtaHO (const double &etaR, const double &x, const double &y, const double &z) const

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

int	getFrontLayer (const int &det, const int &eta) const

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

const std::vector< HcalDetId > &	getIdHF2QIE () const

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

int	getLastLayer (const int &det, const int &eta) 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 &zside, const int &depth) const

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

int	getLayerMax (const int &eta, const int &depth) const

int	getMaxDepth (const int &type) const

int	getMaxDepth (const int &det, const int &eta, const int &phi, const int &zside, const bool &partialOnly) const

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

int	getMinDepth (const int &det, const int &eta, const int &phi, const int &zside, const bool &partialOnly) const

std::pair< int, int >	getModHalfHBHE (const int &type) const

std::pair< double, double >	getPhiCons (const int &det, const int &ieta) const

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

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

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

std::vector< HcalCellType >	HcalCellTypes () const

std::vector< HcalCellType >	HcalCellTypes (const HcalSubdetector &, int ieta=-1, int depth=-1) const

	HcalDDDSimConstants (const HcalParameters *hp)

bool	isBH () const

bool	isHE () const

const HcalLayerDepthMap *	ldMap () const

int	maxHFDepth (const int &ieta, const int &iphi) const

unsigned int	numberOfCells (const HcalSubdetector &) const

const HcalParameters *	parameter () const

int	phiNumber (const int &phi, const int &unit) const

void	printTiles () const

int	unitPhi (const int &det, const int &etaR) const

int	unitPhi (const double &dphi) const

	~HcalDDDSimConstants ()

Private Member Functions
double	deltaEta (const int &det, const int &eta, const int &depth) const

double	getEta (const int &det, const int &etaR, const int &zside, int depth=1) const

double	getEta (const double &r, const double &z) const

double	getGain (const HcalSubdetector &subdet, const int &depth) const

int	getShift (const HcalSubdetector &subdet, const int &depth) const

void	initialize ()

unsigned int	layerGroup (int eta, int i) const

unsigned int	layerGroup (int det, int eta, int phi, int zside, int i) const

unsigned int	layerGroupSize (int eta) const

void	printTileHB (const int &eta, const int &phi, const int &zside, const int &depth) const

void	printTileHE (const int &eta, const int &phi, const int &zside, const int &depth) const

Private Attributes
int	depthEta16 [2]

int	depthEta29 [2]

std::pair< int, int >	depthMaxDf_

std::pair< int, int >	depthMaxSp_

std::vector< int >	depths [nDepthMax]

double	dlShort

double	dzVcal

double	etaHO [4]

const HcalParameters *	hpar

std::vector< HcalDetId >	idHF2QIE

bool	isBH_

int	layBHB [3]

int	layBHE [4]

int	layFHB [2]

int	layFHE [3]

HcalLayerDepthMap	ldmap_

std::vector< int >	maxDepth

int	nDepth

int	nEta

int	nmodHB

int	nmodHE

int	nPhiF

int	nR

int	nzHB

int	nzHE

double	rminHO

double	zVcal

Static Private Attributes
static const int	maxLayer_ = 18

static const int	maxLayerHB_ = 16

static const int	nDepthMax = 9

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 24 of file HcalDDDSimConstants.h.

Constructor & Destructor Documentation

◆ HcalDDDSimConstants()

HcalDDDSimConstants::HcalDDDSimConstants ( const HcalParameters * hp )

Definition at line 10 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), and initialize().

                                                                  : hpar(hp) {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants::HcalDDDSimConstants (const HcalParameters* hp) constructor";
 #endif
 
   initialize();
 #ifdef EDM_ML_DEBUG
   std::vector<HcalCellType> cellTypes = HcalCellTypes();
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size() << " cells of type HCal (All)";
 #endif
 }

◆ ~HcalDDDSimConstants()

HcalDDDSimConstants::~HcalDDDSimConstants ( )

Definition at line 22 of file HcalDDDSimConstants.cc.

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

Member Function Documentation

◆ cell()

HcalCellType::HcalCell HcalDDDSimConstants::cell	(	const int &	det,
		const int &	zside,
		const int &	depth,
		const int &	etaR,
		const int &	iphi
	)		const

Definition at line 28 of file HcalDDDSimConstants.cc.

References d1, deltaEta(), hcalRecHitTable_cff::depth, HcalParameters::drHB, HcalParameters::dzHE, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, getEta(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, hcalRecHitTable_cff::iphi, maxDepth, maxLayerHB_, nEta, HcalParameters::noff, nR, convertSQLiteXML::ok, phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rHB, HcalParameters::rTable, createJobs::tmp, unitPhi(), HcalParameters::zHE, and ecaldqm::zside().

Referenced by HcalCellTypes(), and HcalTestAnalysis::update().

                                                                                                  {
   double etaMn = hpar->etaMin[0];
   double etaMx = hpar->etaMax[0];
   if (idet == static_cast<int>(HcalEndcap)) {
     etaMn = hpar->etaMin[1];
     etaMx = hpar->etaMax[1];
   } else if (idet == static_cast<int>(HcalForward)) {
     etaMn = hpar->etaMin[2];
     etaMx = hpar->etaMax[2];
   }
   double eta = 0, deta = 0, phi = 0, dphi = 0, rz = 0, drz = 0;
   bool ok = false, flagrz = true;
   if ((idet == static_cast<int>(HcalBarrel) || idet == static_cast<int>(HcalEndcap) ||
        idet == static_cast<int>(HcalOuter) || idet == static_cast<int>(HcalForward)) &&
       etaR >= etaMn && etaR <= etaMx && depth > 0)
     ok = true;
   if (idet == static_cast<int>(HcalEndcap) && depth > static_cast<int>(hpar->zHE.size()))
     ok = false;
   else if (idet == static_cast<int>(HcalBarrel) && depth > maxLayerHB_ + 1)
     ok = false;
   else if (idet == static_cast<int>(HcalOuter) && depth != 4)
     ok = false;
   else if (idet == static_cast<int>(HcalForward) && depth > maxDepth[2])
     ok = false;
   if (ok) {
     eta = getEta(idet, etaR, zside, depth);
     deta = deltaEta(idet, etaR, depth);
     double fibin, fioff;
     if (idet == static_cast<int>(HcalBarrel) || idet == static_cast<int>(HcalOuter)) {
       fioff = hpar->phioff[0];
       fibin = hpar->phibin[etaR - 1];
     } else if (idet == static_cast<int>(HcalEndcap)) {
       fioff = hpar->phioff[1];
       fibin = hpar->phibin[etaR - 1];
     } else {
       fioff = hpar->phioff[2];
       fibin = hpar->phitable[etaR - hpar->etaMin[2]];
       if (unitPhi(fibin) > 2)
         fioff = hpar->phioff[4];
     }
     phi = -fioff + (iphi - 0.5) * fibin;
     dphi = 0.5 * fibin;
     if (idet == static_cast<int>(HcalForward)) {
       int ir = nR + hpar->etaMin[2] - etaR - 1;
       if (ir > 0 && ir < nR) {
         rz = 0.5 * (hpar->rTable[ir] + hpar->rTable[ir - 1]);
         drz = 0.5 * (hpar->rTable[ir] - hpar->rTable[ir - 1]);
       } else {
         ok = false;
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: wrong eta " << etaR << " (" << ir << "/" << nR
                                      << ") Detector " << idet;
 #endif
       }
     } else if (etaR <= nEta) {
       int laymin(depth), laymax(depth);
       if (idet == static_cast<int>(HcalOuter)) {
         laymin =
             (etaR > hpar->noff[2]) ? (static_cast<int>(hpar->zHE.size())) : (static_cast<int>(hpar->zHE.size())) - 1;
         laymax = (static_cast<int>(hpar->zHE.size()));
       }
       double d1 = 0, d2 = 0;
       if (idet == static_cast<int>(HcalEndcap)) {
         flagrz = false;
         d1 = hpar->zHE[laymin - 1] - hpar->dzHE[laymin - 1];
         d2 = hpar->zHE[laymax - 1] + hpar->dzHE[laymax - 1];
       } else {
         d1 = hpar->rHB[laymin - 1] - hpar->drHB[laymin - 1];
         d2 = hpar->rHB[laymax - 1] + hpar->drHB[laymax - 1];
       }
       rz = 0.5 * (d2 + d1);
       drz = 0.5 * (d2 - d1);
     } else {
       ok = false;
       edm::LogWarning("HCalGeom") << "HcalDDDSimConstants: wrong depth " << depth << " or etaR " << etaR
                                   << " for detector " << idet;
     }
   } else {
     ok = false;
     edm::LogWarning("HCalGeom") << "HcalDDDSimConstants: wrong depth " << depth << " det " << idet;
   }
   HcalCellType::HcalCell tmp(ok, eta, deta, phi, dphi, rz, drz, flagrz);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: det/side/depth/etaR/"
                                << "phi " << idet << "/" << zside << "/" << depth << "/" << etaR << "/" << iphi
                                << " Cell Flag " << tmp.ok << " " << tmp.eta << " " << tmp.deta << " phi " << tmp.phi
                                << " " << tmp.dphi << " r(z) " << tmp.rz << " " << tmp.drz << " " << tmp.flagrz;
 #endif
   return tmp;
 }

◆ deltaEta()

double HcalDDDSimConstants::deltaEta	(	const int &	det,
		const int &	eta,
		const int &	depth
	)		const

private

Definition at line 979 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, createJobs::tmp, z, and zVcal.

Referenced by cell().

                                                                                             {
   double tmp = 0;
   if (det == static_cast<int>(HcalForward)) {
     int ir = nR + hpar->etaMin[2] - etaR - 1;
     if (ir > 0 && ir < nR) {
       double z = zVcal;
       if (depth % 2 != 1)
         z += dlShort;
       tmp = 0.5 * (getEta(hpar->rTable[ir - 1], z) - getEta(hpar->rTable[ir], z));
     }
   } else {
     if (etaR > 0 && etaR < nEta) {
       if (etaR == hpar->noff[1] - 1 && depth > depthEta29[0]) {
         tmp = 0.5 * (hpar->etaTable[etaR + 1] - hpar->etaTable[etaR - 1]);
       } else if (det == static_cast<int>(HcalOuter)) {
         if (etaR == hpar->noff[2]) {
           tmp = 0.5 * (etaHO[0] - hpar->etaTable[etaR - 1]);
         } else if (etaR == hpar->noff[2] + 1) {
           tmp = 0.5 * (hpar->etaTable[etaR] - etaHO[1]);
         } else if (etaR == hpar->noff[3]) {
           tmp = 0.5 * (etaHO[2] - hpar->etaTable[etaR - 1]);
         } else if (etaR == hpar->noff[3] + 1) {
           tmp = 0.5 * (hpar->etaTable[etaR] - etaHO[3]);
         } else {
           tmp = 0.5 * (hpar->etaTable[etaR] - hpar->etaTable[etaR - 1]);
         }
       } else {
         tmp = 0.5 * (hpar->etaTable[etaR] - hpar->etaTable[etaR - 1]);
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants::deltaEta " << etaR << " " << depth << " ==> " << tmp;
 #endif
   return tmp;
 }

◆ findDepth()

int HcalDDDSimConstants::findDepth	(	const int &	det,
		const int &	eta,
		const int &	phi,
		const int &	zside,
		const int &	lay
	)		const

Definition at line 121 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, HcalLayerDepthMap::getDepth(), HcalLayerDepthMap::isValid(), ldmap_, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::findDepth(), HcalDDDRecConstants::getHCID(), and layerGroup().

                                                                                             {
   int depth = (ldmap_.isValid(det, phi, zside)) ? ldmap_.getDepth(det, eta, phi, zside, lay) : -1;
   return depth;
 }

◆ findLayer()

unsigned int HcalDDDSimConstants::findLayer	(	const int &	layer,
		const std::vector< HcalParameters::LayerItem > &	layerGroup
	)		const

Definition at line 127 of file HcalDDDSimConstants.cc.

References mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, nano_mu_digi_cff::layer, and layerGroup().

                                                                                                         {
   unsigned int id = layerGroup.size();
   for (unsigned int i = 0; i < layerGroup.size(); i++) {
     if (layer == static_cast<int>(layerGroup[i].layer)) {
       id = i;
       break;
     }
   }
   return id;
 }

◆ getConstHBHE()

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

Definition at line 139 of file HcalDDDSimConstants.cc.

References HcalParameters::drHB, HcalParameters::dzHE, hpar, mps_fire::i, HcalParameters::rHB, and HcalParameters::zHE.

                                                                                            {
   std::vector<std::pair<double, double> > gcons;
   if (type == 0) {
     for (unsigned int i = 0; i < hpar->rHB.size(); ++i) {
       gcons.emplace_back(std::pair<double, double>(hpar->rHB[i], hpar->drHB[i]));
     }
   } else {
     for (unsigned int i = 0; i < hpar->zHE.size(); ++i) {
       gcons.emplace_back(std::pair<double, double>(hpar->zHE[i], hpar->dzHE[i]));
     }
   }
   return gcons;
 }

◆ getDepthEta16()

int HcalDDDSimConstants::getDepthEta16	(	const int &	det,
		const int &	phi,
		const int &	zside
	)		const

Definition at line 153 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta16, HcalLayerDepthMap::getDepth16(), ldmap_, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getDepthEta16(), getEtaDepth(), HcalDDDRecConstants::getHCID(), getMaxDepth(), getMinDepth(), HcalDDDRecConstants::getOneEtaBin(), HcalDDDRecConstants::initialize(), and printTiles().

                                                                                              {
   int depth = ldmap_.getDepth16(det, phi, zside);
   if (depth < 0)
     depth = (det == 2) ? depthEta16[1] : depthEta16[0];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "getDepthEta16: " << det << ":" << depth;
 #endif
   return depth;
 }

◆ getDepthEta16M()

int HcalDDDSimConstants::getDepthEta16M ( const int & det ) const

Definition at line 163 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta16, HcalLayerDepthMap::getDepth16(), hcalRecHitTable_cff::iphi, ldmap_, HcalLayerDepthMap::validDet(), and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getMinDepth(), HcalDDDRecConstants::getOneEtaBin(), and HcalDDDRecConstants::initialize().

                                                             {
   int depth = (det == 2) ? depthEta16[1] : depthEta16[0];
   std::vector<int> phis;
   int detsp = ldmap_.validDet(phis);
   if (detsp == det) {
     int zside = (phis[0] > 0) ? 1 : -1;
     int iphi = (phis[0] > 0) ? phis[0] : -phis[0];
     int depthsp = ldmap_.getDepth16(det, iphi, zside);
     if (det == 1 && depthsp > depth)
       depth = depthsp;
     if (det == 2 && depthsp < depth)
       depth = depthsp;
   }
   return depth;
 }

◆ getDepthEta29()

int HcalDDDSimConstants::getDepthEta29	(	const int &	phi,
		const int &	zside,
		const int &	i
	)		const

Definition at line 179 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), mps_fire::i, ldmap_, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getDepthEta29(), getEtaDepth(), HcalDDDRecConstants::getHCID(), HcalDDDRecConstants::getOneEtaBin(), and HcalDDDRecConstants::mergedDepthList29().

                                                                                            {
   int depth = (i == 0) ? ldmap_.getMaxDepthLastHE(2, phi, zside) : -1;
   if (depth < 0)
     depth = (i == 1) ? depthEta29[1] : depthEta29[0];
   return depth;
 }

◆ getDepthEta29M()

int HcalDDDSimConstants::getDepthEta29M	(	const int &	i,
		const bool &	planOne
	)		const

Definition at line 186 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), mps_fire::i, hcalRecHitTable_cff::iphi, ldmap_, HcalLayerDepthMap::validDet(), and ecaldqm::zside().

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

                                                                                {
   int depth = (i == 1) ? depthEta29[1] : depthEta29[0];
   if (i == 0 && planOne) {
     std::vector<int> phis;
     int detsp = ldmap_.validDet(phis);
     if (detsp == 2) {
       int zside = (phis[0] > 0) ? 1 : -1;
       int iphi = (phis[0] > 0) ? phis[0] : -phis[0];
       int depthsp = ldmap_.getMaxDepthLastHE(2, iphi, zside);
       if (depthsp > depth)
         depth = depthsp;
     }
   }
   return depth;
 }

◆ getDetEta()

std::pair< int, double > HcalDDDSimConstants::getDetEta	(	const double &	eta,
		const int &	depth
	)		const

Definition at line 202 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta16, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, generateTowerEtThresholdLUT::etaRange, HcalParameters::etaTable, JetChargeProducer_cfi::exp, HcalBarrel, HcalEndcap, HcalForward, hpar, mps_fire::i, hcalRecHitTable_cff::ieta, nEta, funct::tan(), tauSpinnerTable_cfi::theta, and zVcal.

Referenced by HcalNumberingFromDDD::unitID().

                                                                                              {
   int hsubdet(0), ieta(0);
   double etaR(0);
   double heta = fabs(eta);
   for (int i = 0; i < nEta; i++)
     if (heta > hpar->etaTable[i])
       ieta = i + 1;
   if (heta <= hpar->etaRange[1]) {
     if (((ieta == hpar->etaMin[1] && depth == depthEta16[1]) || (ieta > hpar->etaMax[0])) &&
         (ieta <= hpar->etaMax[1])) {
       hsubdet = static_cast<int>(HcalEndcap);
     } else {
       hsubdet = static_cast<int>(HcalBarrel);
     }
     etaR = eta;
   } else {
     hsubdet = static_cast<int>(HcalForward);
     double theta = 2. * atan(exp(-heta));
     double hR = zVcal * tan(theta);
     etaR = (eta >= 0. ? hR : -hR);
   }
   return std::pair<int, double>(hsubdet, etaR);
 }

◆ getEta() [1/3]

int HcalDDDSimConstants::getEta	(	const int &	det,
		const int &	lay,
		const double &	hetaR
	)		const

Definition at line 226 of file HcalDDDSimConstants.cc.

References etaHO, HcalParameters::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, HcalParameters::etaMin, HcalParameters::etaTable, HcalBarrel, HcalEndcap, HcalForward, hpar, mps_fire::i, hcalRecHitTable_cff::ieta, maxLayer_, nEta, HcalParameters::noff, and nR.

Referenced by cell(), deltaEta(), getEta(), and HcalNumberingFromDDD::unitID().

                                                                                          {
   int ieta(0);
   if (det == static_cast<int>(HcalForward)) {  // Forward HCal
     ieta = hpar->etaMax[2];
     for (int i = nR - 1; i > 0; i--)
       if (hetaR < hpar->rTable[i])
         ieta = hpar->etaMin[2] + nR - i - 1;
   } else {  // Barrel or Endcap
     ieta = 1;
     for (int i = 0; i < nEta - 1; i++)
       if (hetaR > hpar->etaTable[i])
         ieta = i + 1;
     if (det == static_cast<int>(HcalBarrel)) {
       if (ieta > hpar->etaMax[0])
         ieta = hpar->etaMax[0];
       if (lay == maxLayer_) {
         if (hetaR > etaHO[1] && ieta == hpar->noff[2])
           ieta++;
       }
     } else if (det == static_cast<int>(HcalEndcap)) {
       if (ieta <= hpar->etaMin[1])
         ieta = hpar->etaMin[1];
     }
   }
   return ieta;
 }

◆ getEta() [2/3]

double HcalDDDSimConstants::getEta	(	const int &	det,
		const int &	etaR,
		const int &	zside,
		int	depth = `1`
	)		const

private

Definition at line 1016 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, createJobs::tmp, z, ecaldqm::zside(), and zVcal.

                                                                                                      {
   double tmp = 0;
   if (det == static_cast<int>(HcalForward)) {
     int ir = nR + hpar->etaMin[2] - etaR - 1;
     if (ir > 0 && ir < nR) {
       double z = zVcal;
       if (depth % 2 != 1)
         z += dlShort;
       tmp = 0.5 * (getEta(hpar->rTable[ir - 1], z) + getEta(hpar->rTable[ir], z));
     }
   } else {
     if (etaR > 0 && etaR < nEta) {
       if (etaR == hpar->noff[1] - 1 && depth > depthEta29[0]) {
         tmp = 0.5 * (hpar->etaTable[etaR + 1] + hpar->etaTable[etaR - 1]);
       } else if (det == static_cast<int>(HcalOuter)) {
         if (etaR == hpar->noff[2]) {
           tmp = 0.5 * (etaHO[0] + hpar->etaTable[etaR - 1]);
         } else if (etaR == hpar->noff[2] + 1) {
           tmp = 0.5 * (hpar->etaTable[etaR] + etaHO[1]);
         } else if (etaR == hpar->noff[3]) {
           tmp = 0.5 * (etaHO[2] + hpar->etaTable[etaR - 1]);
         } else if (etaR == hpar->noff[3] + 1) {
           tmp = 0.5 * (hpar->etaTable[etaR] + etaHO[3]);
         } else {
           tmp = 0.5 * (hpar->etaTable[etaR] + hpar->etaTable[etaR - 1]);
         }
       } else {
         tmp = 0.5 * (hpar->etaTable[etaR] + hpar->etaTable[etaR - 1]);
       }
     }
   }
   if (zside == 0)
     tmp = -tmp;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants::getEta " << etaR << " " << zside << " " << depth << " ==> "
                                << tmp;
 #endif
   return tmp;
 }

◆ getEta() [3/3]

double HcalDDDSimConstants::getEta	(	const double &	r,
		const double &	z
	)		const

private

Definition at line 1056 of file HcalDDDSimConstants.cc.

References CrabHelper::log, alignCSCRings::r, funct::tan(), createJobs::tmp, and z.

                                                                          {
   double tmp = 0;
   if (z != 0)
     tmp = -log(tan(0.5 * atan(r / z)));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants::getEta " << r << " " << z << " ==> " << tmp;
 #endif
   return tmp;
 }

◆ getEtaDepth()

std::pair< int, int > HcalDDDSimConstants::getEtaDepth	(	const int &	det,
		int	etaR,
		const int &	phi,
		const int &	zside,
		int	depth,
		const int &	lay
	)		const

Definition at line 253 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, HcalParameters::etaMin, getDepthEta16(), getDepthEta29(), getMaxDepth(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, layerGroup(), HcalParameters::noff, phi, HcalParameters::phibin, HcalParameters::phioff, and ecaldqm::zside().

Referenced by HcalTestNumberingTester::analyze(), and HcalNumberingFromDDD::unitID().

                                                                                                  {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDEsimConstants:getEtaDepth: I/P " << det << ":" << etaR << ":" << phi << ":"
                                << zside << ":" << depth << ":" << lay;
 #endif
   //Modify the depth index
   if ((det == static_cast<int>(HcalEndcap)) && (etaR == 17) && (lay == 1))
     etaR = 18;
   if (det == static_cast<int>(HcalForward)) {  // Forward HCal
   } else if (det == static_cast<int>(HcalOuter)) {
     depth = 4;
   } else {
     if (lay >= 0) {
       depth = layerGroup(det, etaR, phi, zside, lay - 1);
       if (((det == 2) && (etaR == 18)) || ((det == 1) && (etaR == 15)))
         if (etaR == hpar->noff[0] && lay > 1) {
           int kphi = phi + static_cast<int>((hpar->phioff[3] + 0.1) / hpar->phibin[etaR - 1]);
           kphi = (kphi - 1) % 4 + 1;
           if (kphi == 2 || kphi == 3)
             depth = layerGroup(det, etaR, phi, zside, lay - 2);
         }
     } else if (det == static_cast<int>(HcalBarrel)) {
       if (depth > getMaxDepth(det, etaR, phi, zside, false))
         depth = getMaxDepth(det, etaR, phi, zside, false);
     }
     if (etaR >= hpar->noff[1] && depth > getDepthEta29(phi, zside, 0)) {
       etaR -= getDepthEta29(phi, zside, 1);
     } else if (etaR == hpar->etaMin[1]) {
       if (det == static_cast<int>(HcalBarrel)) {
         if (depth > getDepthEta16(det, phi, zside))
           depth = getDepthEta16(det, phi, zside);
       } else {
         if (depth < getDepthEta16(det, phi, zside))
           depth = getDepthEta16(det, phi, zside);
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDEsimConstants:getEtaDepth: O/P " << etaR << ":" << depth;
 #endif
   return std::pair<int, int>(etaR, depth);
 }

◆ getEtaHO()

double HcalDDDSimConstants::getEtaHO	(	const double &	etaR,
		const double &	x,
		const double &	y,
		const double &	z
	)		const

Definition at line 297 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, HcalParameters::etaTable, hpar, alignCSCRings::r, rminHO, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, x, y, z, HcalParameters::zHO, and geometryCSVtoXML::zz.

Referenced by HcalNumberingFromDDD::unitID().

                                                                                                                 {
   if (hpar->zHO.size() > 4) {
     double eta = fabs(etaR);
     double r = std::sqrt(x * x + y * y);
     if (r > rminHO) {
       double zz = fabs(z);
       if (zz > hpar->zHO[3]) {
         if (eta <= hpar->etaTable[10])
           eta = hpar->etaTable[10] + 0.001;
       } else if (zz > hpar->zHO[1]) {
         if (eta <= hpar->etaTable[4])
           eta = hpar->etaTable[4] + 0.001;
       }
     }
     eta = (z >= 0. ? eta : -eta);
 #ifdef EDM_ML_DEBUG
     std::string chk = (eta != etaR) ? " **** Check *****" : " ";
     edm::LogVerbatim("HCalGeom") << "R " << r << " Z " << z << " eta " << etaR << ":" << eta << chk;
 #endif
     return eta;
   } else {
     return etaR;
   }
 }

◆ getEtaTableHF()

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

inline

Definition at line 44 of file HcalDDDSimConstants.h.

References HcalParameters::etaTableHF, and hpar.

44 { return hpar->etaTableHF; }

HcalParameters::etaTableHF

std::vector< double > etaTableHF

Definition: HcalParameters.h:45

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

◆ getFrontLayer()

int HcalDDDSimConstants::getFrontLayer	(	const int &	det,
		const int &	eta
	)		const

Definition at line 322 of file HcalDDDSimConstants.cc.

References funct::abs(), PVValHelper::eta, layFHB, and layFHE.

Referenced by HcalDDDRecConstants::getLayerFront().

                                                                            {
   int lay = 0;
   if (det == 1) {
     if (std::abs(eta) == 16)
       lay = layFHB[1];
     else
       lay = layFHB[0];
   } else {
     if (std::abs(eta) == 16)
       lay = layFHE[1];
     else if (std::abs(eta) == 18)
       lay = layFHE[2];
     else
       lay = layFHE[0];
   }
   return lay;
 }

◆ getGain()

double HcalDDDSimConstants::getGain	(	const HcalSubdetector &	subdet,
		const int &	depth
	)		const

private

Definition at line 1085 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, Hcal_Conditions_forGlobalTag_cff::gain, HcalParameters::HBGains, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEGains, HcalParameters::HFGains, and hpar.

Referenced by HcalCellTypes().

                                                                                          {
   double gain;
   switch (subdet) {
     case HcalEndcap:
       gain = hpar->HEGains[0];
       break;
     case HcalForward:
       gain = hpar->HFGains[(depth - 1) % 2];
       break;
     case HcalOuter:
       gain = hpar->HBGains[3];
       break;
     default:
       gain = hpar->HBGains[0];
       break;
   }
   return gain;
 }

◆ getGparHF()

const std::vector<double>& HcalDDDSimConstants::getGparHF ( ) const

inline

Definition at line 45 of file HcalDDDSimConstants.h.

References HcalParameters::gparHF, and hpar.

Referenced by HCalSD::HCalSD(), HFShower::HFShower(), and HFShowerParam::HFShowerParam().

45 { return hpar->gparHF; }

HcalParameters::gparHF

std::vector< double > gparHF

Definition: HcalParameters.h:40

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

◆ getIdHF2QIE()

const std::vector<HcalDetId>& HcalDDDSimConstants::getIdHF2QIE ( ) const

inline

Definition at line 46 of file HcalDDDSimConstants.h.

References idHF2QIE.

Referenced by HcalDDDRecConstants::getHFCellParameters().

46 { return idHF2QIE; }

HcalDDDSimConstants::idHF2QIE

std::vector< HcalDetId > idHF2QIE

Definition: HcalDDDSimConstants.h:114

◆ getiEtaRange()

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

inline

Definition at line 43 of file HcalDDDSimConstants.h.

References HcalParameters::etaMax, HcalParameters::etaMin, hpar, and mps_fire::i.

43 { return std::pair<int, int>(hpar->etaMin[i], hpar->etaMax[i]); }

mps_fire.i

i

Definition: mps_fire.py:429

HcalParameters::etaMax

std::vector< int > etaMax

Definition: HcalParameters.h:55

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

HcalParameters::etaMin

std::vector< int > etaMin

Definition: HcalParameters.h:54

◆ getLastLayer()

int HcalDDDSimConstants::getLastLayer	(	const int &	det,
		const int &	eta
	)		const

Definition at line 340 of file HcalDDDSimConstants.cc.

References funct::abs(), PVValHelper::eta, layBHB, and layBHE.

Referenced by HcalDDDRecConstants::getLayerBack(), and HcalDDDRecConstants::getThickActive().

                                                                           {
   int lay = 0;
   if (det == 1) {
     if (std::abs(eta) == 15)
       lay = layBHB[1];
     else if (std::abs(eta) == 16)
       lay = layBHB[2];
     else
       lay = layBHB[0];
   } else {
     if (std::abs(eta) == 16)
       lay = layBHE[1];
     else if (std::abs(eta) == 17)
       lay = layBHE[2];
     else if (std::abs(eta) == 18)
       lay = layBHE[3];
     else
       lay = layBHE[0];
   }
   return lay;
 }

◆ getLayer0Wt()

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

Definition at line 362 of file HcalDDDSimConstants.cc.

References HcalLayerDepthMap::getLayer0Wt(), hpar, HcalParameters::Layer0Wt, ldmap_, phi, and ecaldqm::zside().

Referenced by HCalSD::getEnergyDeposit(), and HcalDDDRecConstants::getLayer0Wt().

                                                                                               {
   double wt = ldmap_.getLayer0Wt(det, phi, zside);
   if (wt < 0)
     wt = (det == 2) ? hpar->Layer0Wt[1] : hpar->Layer0Wt[0];
   return wt;
 }

◆ getLayerBack()

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

Definition at line 388 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, depths, PVValHelper::eta, HcalLayerDepthMap::getLayerBack(), nano_mu_digi_cff::layer, ldmap_, maxLayer_, phi, and ecaldqm::zside().

Referenced by printTileHB(), and printTileHE().

                                                                                               {
   int layer = ldmap_.getLayerBack(det, eta, phi, zside, depth);
   if (layer < 0) {
     if (det == 1 || det == 2) {
       layer = depths[depth - 1][eta - 1];
     } else {
       layer = maxLayer_;
     }
   }
   return layer;
 }

◆ getLayerFront()

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

Definition at line 369 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, HcalLayerDepthMap::getLayerFront(), getLayerMax(), hpar, MainPageGenerator::l, nano_mu_digi_cff::layer, layerGroup(), ldmap_, maxLayer_, maxLayerHB_, HcalParameters::noff, phi, and ecaldqm::zside().

Referenced by printTileHB(), and printTileHE().

                                                                                               {
   int layer = ldmap_.getLayerFront(det, eta, phi, zside, depth);
   if (layer < 0) {
     if (det == 1 || det == 2) {
       layer = 1;
       for (int l = 0; l < getLayerMax(eta, depth); ++l) {
         if (static_cast<int>(layerGroup(eta - 1, l)) == depth) {
           layer = l + 1;
           break;
         }
       }
     } else {
       layer = (eta > hpar->noff[2]) ? maxLayerHB_ + 1 : maxLayer_;
     }
   }
   return layer;
 }

◆ getLayerMax()

int HcalDDDSimConstants::getLayerMax	(	const int &	eta,
		const int &	depth
	)		const

Definition at line 401 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, HcalParameters::etaMin, hpar, layerGroupSize(), maxDepth, and maxLayerHB_.

Referenced by getLayerFront(), and initialize().

                                                                            {
   int layermx = ((eta < hpar->etaMin[1]) && depth - 1 < maxDepth[0]) ? maxLayerHB_ + 1
                                                                      : static_cast<int>(layerGroupSize(eta - 1));
   return layermx;
 }

◆ getMaxDepth() [1/2]

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

inline

Definition at line 53 of file HcalDDDSimConstants.h.

References maxDepth.

Referenced by getEtaDepth(), HcalDDDRecConstants::getMaxDepth(), HCalSD::HCalSD(), and printTiles().

53 { return maxDepth[type]; }

HcalDDDSimConstants::maxDepth

std::vector< int > maxDepth

Definition: HcalDDDSimConstants.h:95

gainCalibHelper::gainCalibPI::type

type

Definition: SiPixelGainCalibHelper.h:40

◆ getMaxDepth() [2/2]

int HcalDDDSimConstants::getMaxDepth	(	const int &	det,
		const int &	eta,
		const int &	phi,
		const int &	zside,
		const bool &	partialOnly
	)		const

Definition at line 407 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, HcalParameters::etaMax, getDepthEta16(), HcalLayerDepthMap::getDepths(), hpar, HcalLayerDepthMap::isValid(), layerGroup(), ldmap_, maxDepth, maxHFDepth(), maxLayer_, maxLayerHB_, phi, and ecaldqm::zside().

                                                                                                         {
   int dmax(-1);
   if (partialDetOnly) {
     if (ldmap_.isValid(det, phi, zside)) {
       dmax = ldmap_.getDepths(eta).second;
     }
   } else if (det == 1 || det == 2) {
     if (ldmap_.isValid(det, phi, zside))
       dmax = ldmap_.getDepths(eta).second;
     else if (det == 2)
       dmax = (maxDepth[1] > 0) ? layerGroup(eta - 1, maxLayer_) : 0;
     else if (eta == hpar->etaMax[0])
       dmax = getDepthEta16(det, phi, zside);
     else
       dmax = layerGroup(eta - 1, maxLayerHB_);
   } else if (det == 3) {  // HF
     dmax = maxHFDepth(zside * eta, phi);
   } else if (det == 4) {  // HO
     dmax = maxDepth[3];
   } else {
     dmax = -1;
   }
   return dmax;
 }

◆ getMaxDepthDet()

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

inline

Definition at line 55 of file HcalDDDSimConstants.h.

References depthMaxDf_, depthMaxSp_, and mps_fire::i.

Referenced by HcalDDDRecConstants::initialize().

55 { return ((i == 1) ? depthMaxDf_ : depthMaxSp_); }

mps_fire.i

i

Definition: mps_fire.py:429

HcalDDDSimConstants::depthMaxDf_

std::pair< int, int > depthMaxDf_

Definition: HcalDDDSimConstants.h:115

HcalDDDSimConstants::depthMaxSp_

std::pair< int, int > depthMaxSp_

Definition: HcalDDDSimConstants.h:115

◆ getMinDepth()

int HcalDDDSimConstants::getMinDepth	(	const int &	det,
		const int &	eta,
		const int &	phi,
		const int &	zside,
		const bool &	partialOnly
	)		const

Definition at line 433 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, HcalParameters::etaMin, getDepthEta16(), HcalLayerDepthMap::getDepths(), hpar, HcalLayerDepthMap::isValid(), layerGroup(), layerGroupSize(), ldmap_, maxDepth, phi, and ecaldqm::zside().

Referenced by HcalDDDRecConstants::getMinDepth().

                                                                                                         {
   int lmin(-1);
   if (partialDetOnly) {
     if (ldmap_.isValid(det, phi, zside)) {
       lmin = ldmap_.getDepths(eta).first;
     }
   } else if (det == 3) {  // HF
     lmin = 1;
   } else if (det == 4) {  // HO
     lmin = maxDepth[3];
   } else {
     if (ldmap_.isValid(det, phi, zside)) {
       lmin = ldmap_.getDepths(eta).first;
     } else if (layerGroupSize(eta - 1) > 0) {
       lmin = static_cast<int>(layerGroup(eta - 1, 0));
       unsigned int type = (det == 1) ? 0 : 1;
       if (type == 1 && eta == hpar->etaMin[1])
         lmin = getDepthEta16(det, phi, zside);
     } else {
       lmin = 1;
     }
   }
   return lmin;
 }

◆ getModHalfHBHE()

std::pair< int, int > HcalDDDSimConstants::getModHalfHBHE ( const int & type ) const

Definition at line 459 of file HcalDDDSimConstants.cc.

References nmodHB, nmodHE, nzHB, and nzHE.

                                                                            {
   if (type == 0) {
     return std::pair<int, int>(nmodHB, nzHB);
   } else {
     return std::pair<int, int>(nmodHE, nzHE);
   }
 }

◆ getPhiCons()

std::pair< double, double > HcalDDDSimConstants::getPhiCons	(	const int &	det,
		const int &	ieta
	)		const

Definition at line 467 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, HcalEndcap, HcalForward, hpar, hcalRecHitTable_cff::ieta, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, and unitPhi().

Referenced by getPhis(), HcalDDDRecConstants::getPhis(), and HcalNumberingFromDDD::unitID().

                                                                                              {
   double fioff(0), fibin(0);
   if (det == static_cast<int>(HcalForward)) {  // Forward HCal
     fioff = hpar->phioff[2];
     fibin = hpar->phitable[ieta - hpar->etaMin[2]];
     if (unitPhi(fibin) > 2) {  // HF double-phi
       fioff = hpar->phioff[4];
     }
   } else {  // Barrel or Endcap
     if (det == static_cast<int>(HcalEndcap)) {
       fioff = hpar->phioff[1];
     } else {
       fioff = hpar->phioff[0];
     }
     fibin = hpar->phibin[ieta - 1];
   }
   return std::pair<double, double>(fioff, fibin);
 }

◆ getPhis()

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

Definition at line 486 of file HcalDDDSimConstants.cc.

References angle_units::operators::convertRadToDeg(), getPhiCons(), hcalRecHitTable_cff::ieta, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), createfilelist::int, hcalRecHitTable_cff::iphi, isotrackApplyRegressor::k, nphi, phi, phiNumber(), edm::second(), unitPhi(), and units().

Referenced by HcalCellTypes().

                                                                                                       {
   std::vector<std::pair<int, double> > phis;
   int ietaAbs = (ieta > 0) ? ieta : -ieta;
   std::pair<double, double> ficons = getPhiCons(subdet, ietaAbs);
   int nphi = int((2._pi + 0.1 * ficons.second) / ficons.second);
   int units = unitPhi(subdet, ietaAbs);
   for (int ifi = 0; ifi < nphi; ++ifi) {
     double phi = -ficons.first + (ifi + 0.5) * ficons.second;
     int iphi = phiNumber(ifi + 1, units);
     phis.emplace_back(std::pair<int, double>(iphi, phi));
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "getPhis: 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;
 }

◆ getPhiTableHF()

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

inline

Definition at line 60 of file HcalDDDSimConstants.h.

References hpar, and HcalParameters::phitable.

Referenced by HcalDDDRecConstants::getHFCellParameters().

60 { return hpar->phitable; }

HcalParameters::phitable

std::vector< double > phitable

Definition: HcalParameters.h:38

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

◆ getRTableHF()

const std::vector<double>& HcalDDDSimConstants::getRTableHF ( ) const

inline

Definition at line 61 of file HcalDDDSimConstants.h.

References hpar, and HcalParameters::rTable.

Referenced by HcalDDDRecConstants::getHFCellParameters(), HFShowerFibreBundle::HFShowerFibreBundle(), and HFShowerPMT::HFShowerPMT().

61 { return hpar->rTable; }

HcalParameters::rTable

std::vector< double > rTable

Definition: HcalParameters.h:36

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

◆ getShift()

int HcalDDDSimConstants::getShift	(	const HcalSubdetector &	subdet,
		const int &	depth
	)		const

private

Definition at line 1066 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, HcalParameters::HBShift, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEShift, HcalParameters::HFShift, hpar, and edm::shift.

Referenced by HcalCellTypes().

                                                                                        {
   int shift;
   switch (subdet) {
     case HcalEndcap:
       shift = hpar->HEShift[0];
       break;
     case HcalForward:
       shift = hpar->HFShift[(depth - 1) % 2];
       break;
     case HcalOuter:
       shift = hpar->HBShift[3];
       break;
     default:
       shift = hpar->HBShift[0];
       break;
   }
   return shift;
 }

◆ HcalCellTypes() [1/2]

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

Definition at line 507 of file HcalDDDSimConstants.cc.

References HcalBarrel, HcalEndcap, HcalForward, HcalOuter, and mps_fire::i.

Referenced by HcalDDDRecConstants::HcalCellTypes(), HcalDDDSimConstants(), and numberOfCells().

                                                                  {
   std::vector<HcalCellType> cellTypes = HcalCellTypes(HcalBarrel);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size() << " cells of type HCal Barrel";
   for (unsigned int i = 0; i < cellTypes.size(); i++)
     edm::LogVerbatim("HCalGeom") << "Cell " << i << " " << cellTypes[i];
 #endif
 
   std::vector<HcalCellType> hoCells = HcalCellTypes(HcalOuter);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: " << hoCells.size() << " cells of type HCal Outer";
   for (unsigned int i = 0; i < hoCells.size(); i++)
     edm::LogVerbatim("HCalGeom") << "Cell " << i << " " << hoCells[i];
 #endif
   cellTypes.insert(cellTypes.end(), hoCells.begin(), hoCells.end());
 
   std::vector<HcalCellType> heCells = HcalCellTypes(HcalEndcap);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: " << heCells.size() << " cells of type HCal Endcap";
   for (unsigned int i = 0; i < heCells.size(); i++)
     edm::LogVerbatim("HCalGeom") << "Cell " << i << " " << heCells[i];
 #endif
   cellTypes.insert(cellTypes.end(), heCells.begin(), heCells.end());
 
   std::vector<HcalCellType> hfCells = HcalCellTypes(HcalForward);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: " << hfCells.size() << " cells of type HCal Forward";
   for (unsigned int i = 0; i < hfCells.size(); i++)
     edm::LogVerbatim("HCalGeom") << "Cell " << i << " " << hfCells[i];
 #endif
   cellTypes.insert(cellTypes.end(), hfCells.begin(), hfCells.end());
 
   return cellTypes;
 }

◆ HcalCellTypes() [2/2]

std::vector< HcalCellType > HcalDDDSimConstants::HcalCellTypes	(	const HcalSubdetector &	subdet,
		int	ieta = `-1`,
		int	depth = `-1`
	)		const

Definition at line 542 of file HcalDDDSimConstants.cc.

References cell(), hcalRecHitTable_cff::depth, dlShort, dzVcal, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, Hcal_Conditions_forGlobalTag_cff::gain, getGain(), getPhis(), getShift(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, idHF2QIE, hcalRecHitTable_cff::ieta, GetRecoTauVFromDQM_MC_cff::kk, MainPageGenerator::l, maxDepth, maxLayer_, maxLayerHB_, HcalParameters::noff, nzHB, nzHE, HcalCellType::HcalCell::ok, convertSQLiteXML::ok, phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalCellType::setPhi(), edm::shift, unit(), unitPhi(), and ecaldqm::zside().

                                                                                {
   std::vector<HcalCellType> cellTypes;
   if (subdet == HcalForward) {
     if (dzVcal < 0)
       return cellTypes;
   }
 
   int dmin, dmax, indx, nz;
   double hsize = 0;
   switch (subdet) {
     case HcalEndcap:
       dmin = 1;
       dmax = (maxDepth[1] > 0) ? maxLayer_ + 1 : 0;
       indx = 1;
       nz = nzHE;
       break;
     case HcalForward:
       dmin = 1;
       dmax = (!idHF2QIE.empty()) ? 2 : maxDepth[2];
       indx = 2;
       nz = 2;
       break;
     case HcalOuter:
       dmin = 4;
       dmax = 4;
       indx = 0;
       nz = nzHB;
       break;
     default:
       dmin = 1;
       dmax = maxLayerHB_ + 1;
       indx = 0;
       nz = nzHB;
       break;
   }
   if (depthl > 0)
     dmin = dmax = depthl;
   int ietamin = (ieta > 0) ? ieta : hpar->etaMin[indx];
   int ietamax = (ieta > 0) ? ieta : hpar->etaMax[indx];
   int phi = (indx == 2) ? 3 : 1;
 
   // Get the Cells
   int subdet0 = static_cast<int>(subdet);
   for (int depth = dmin; depth <= dmax; depth++) {
     int shift = getShift(subdet, depth);
     double gain = getGain(subdet, depth);
     if (subdet == HcalForward) {
       if (depth % 2 == 1)
         hsize = dzVcal;
       else
         hsize = dzVcal - 0.5 * dlShort;
     }
     for (int eta = ietamin; eta <= ietamax; eta++) {
       for (int iz = 0; iz < nz; ++iz) {
         int zside = -2 * iz + 1;
         HcalCellType::HcalCell temp1 = cell(subdet0, zside, depth, eta, phi);
         if (temp1.ok) {
           std::vector<std::pair<int, double> > phis = getPhis(subdet0, eta);
           HcalCellType temp2(subdet, eta, zside, depth, temp1, shift, gain, hsize);
           double dphi, fioff;
           std::vector<int> phiMiss2;
           if ((subdet == HcalBarrel) || (subdet == HcalOuter)) {
             fioff = hpar->phioff[0];
             dphi = hpar->phibin[eta - 1];
             if (subdet == HcalOuter) {
               if (eta == hpar->noff[4]) {
                 int kk = (iz == 0) ? 7 : (7 + hpar->noff[5]);
                 for (int miss = 0; miss < hpar->noff[5 + iz]; miss++) {
                   phiMiss2.emplace_back(hpar->noff[kk]);
                   kk++;
                 }
               }
             }
           } else if (subdet == HcalEndcap) {
             fioff = hpar->phioff[1];
             dphi = hpar->phibin[eta - 1];
           } else {
             fioff = hpar->phioff[2];
             dphi = hpar->phitable[eta - hpar->etaMin[2]];
             if (unitPhi(dphi) > 2)
               fioff = hpar->phioff[4];
           }
           int unit = unitPhi(dphi);
           temp2.setPhi(phis, phiMiss2, fioff, dphi, unit);
           cellTypes.emplace_back(temp2);
           // For HF look at extra cells
           if ((subdet == HcalForward) && (!idHF2QIE.empty())) {
             HcalCellType temp3(subdet, eta, zside + 2, depth, temp1, shift, gain, hsize);
             std::vector<int> phiMiss3;
             for (auto& phi : phis) {
               bool ok(false);
               for (auto l : idHF2QIE) {
                 if ((eta * zside == l.ieta()) && (phi.first == l.iphi())) {
                   ok = true;
                   break;
                 }
               }
               if (!ok)
                 phiMiss3.emplace_back(phi.first);
             }
             dphi = hpar->phitable[eta - hpar->etaMin[2]];
             unit = unitPhi(dphi);
             fioff = (unit > 2) ? hpar->phioff[4] : hpar->phioff[2];
             temp3.setPhi(phis, phiMiss2, fioff, dphi, unit);
             cellTypes.emplace_back(temp3);
           }
         }
       }
     }
   }
   return cellTypes;
 }

◆ initialize()

void HcalDDDSimConstants::initialize ( void )

private

Definition at line 743 of file HcalDDDSimConstants.cc.

References depthEta16, depthEta29, depthMaxDf_, depthMaxSp_, depths, dlShort, isotrackApplyRegressor::dtype, HcalParameters::dzVcal, dzVcal, etaHO, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaTable, cms::cuda::for(), HcalLayerDepthMap::getDepthMax(), getLayerMax(), HcalParameters::gparHF, HcalForward, hpar, mps_fire::i, idHF2QIE, l1tnanotables_cff::iet, HcalLayerDepthMap::initialize(), hcalRecHitTable_cff::iphi, isBH_, isotrackApplyRegressor::k, GetRecoTauVFromDQM_MC_cff::kk, MainPageGenerator::l, layBHB, layBHE, HcalParameters::Layer0Wt, layerGroup(), layerGroupSize(), layFHB, layFHE, ldmap_, HcalParameters::maxDepth, maxDepth, maxLayerHB_, HcalParameters::modHB, HcalParameters::modHE, nEta, nmodHB, nmodHE, HcalParameters::noff, nphi, nPhiF, nR, nzHB, nzHE, HcalParameters::rHB, HcalParameters::rHO, rminHO, HcalParameters::rTable, HcalParameters::zHO, ecaldqm::zside(), and zVcal.

Referenced by HcalDDDSimConstants().

                                          {
   nEta = hpar->etaTable.size();
   nR = hpar->rTable.size();
   nPhiF = nR - 1;
   isBH_ = false;
 
 #ifdef EDM_ML_DEBUG
   for (int i = 0; i < nEta - 1; ++i) {
     edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants:Read LayerGroup" << i << ":";
     for (unsigned int k = 0; k < layerGroupSize(i); k++)
       edm::LogVerbatim("HCalGeom") << " [" << k << "] = " << layerGroup(i, k);
   }
 #endif
 
   // Geometry parameters for HF
   dlShort = hpar->gparHF[0];
   zVcal = hpar->gparHF[4];
   dzVcal = hpar->dzVcal;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: dlShort " << dlShort << " zVcal " << zVcal << " and dzVcal "
                                << dzVcal;
 #endif
 
   //Transform some of the parameters
   maxDepth = hpar->maxDepth;
   maxDepth[0] = maxDepth[1] = 0;
   for (int i = 0; i < nEta - 1; ++i) {
     unsigned int imx = layerGroupSize(i);
     int laymax = (imx > 0) ? layerGroup(i, imx - 1) : 0;
     if (i < hpar->etaMax[0]) {
       int laymax0 = (imx > 16) ? layerGroup(i, 16) : laymax;
       if (i + 1 == hpar->etaMax[0] && laymax0 > 2)
         laymax0 = 2;
       if (maxDepth[0] < laymax0)
         maxDepth[0] = laymax0;
     }
     if (i >= hpar->etaMin[1] - 1 && i < hpar->etaMax[1]) {
       if (maxDepth[1] < laymax)
         maxDepth[1] = laymax;
     }
   }
 #ifdef EDM_ML_DEBUG
   for (int i = 0; i < 4; ++i)
     edm::LogVerbatim("HCalGeom") << "Detector Type [" << i << "] iEta " << hpar->etaMin[i] << ":" << hpar->etaMax[i]
                                  << " MaxDepth " << maxDepth[i];
 #endif
 
   int maxdepth = (maxDepth[1] > maxDepth[0]) ? maxDepth[1] : maxDepth[0];
   for (int i = 0; i < maxdepth; ++i) {
     for (int k = 0; k < nEta - 1; ++k) {
       int layermx = getLayerMax(k + 1, i + 1);
       int ll = layermx;
       for (int l = layermx - 1; l >= 0; --l) {
         if (static_cast<int>(layerGroup(k, l)) == i + 1) {
           ll = l + 1;
           break;
         }
       }
       depths[i].emplace_back(ll);
     }
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HCalGeom") << "Depth " << i + 1 << " with " << depths[i].size() << " etas:";
     for (int k = 0; k < nEta - 1; ++k)
       edm::LogVerbatim("HCalGeom") << " [" << k << "] " << depths[i][k];
 #endif
   }
 
   nzHB = hpar->modHB[1];
   nmodHB = hpar->modHB[0];
   nzHE = hpar->modHE[1];
   nmodHE = hpar->modHE[0];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants:: " << nzHB << ":" << nmodHB << " barrel and " << nzHE << ":"
                                << nmodHE << " endcap half-sectors";
 #endif
 
   if (hpar->rHB.size() > maxLayerHB_ + 1 && hpar->zHO.size() > 4) {
     rminHO = hpar->rHO[0];
     for (int k = 0; k < 4; ++k)
       etaHO[k] = hpar->rHO[k + 1];
   } else {
     rminHO = -1.0;
     etaHO[0] = hpar->etaTable[4];
     etaHO[1] = hpar->etaTable[4];
     etaHO[2] = hpar->etaTable[10];
     etaHO[3] = hpar->etaTable[10];
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HO Eta boundaries " << etaHO[0] << " " << etaHO[1] << " " << etaHO[2] << " "
                                << etaHO[3];
   edm::LogVerbatim("HCalGeom") << "HO Parameters " << rminHO << " " << hpar->zHO.size();
   for (int i = 0; i < 4; ++i)
     edm::LogVerbatim("HCalGeom") << " eta[" << i << "] = " << etaHO[i];
   for (unsigned int i = 0; i < hpar->zHO.size(); ++i)
     edm::LogVerbatim("HCalGeom") << " zHO[" << i << "] = " << hpar->zHO[i];
 #endif
 
   int noffsize = 7 + hpar->noff[5] + hpar->noff[6];
   int noffl(noffsize + 5);
   if (static_cast<int>(hpar->noff.size()) > (noffsize + 3)) {
     depthEta16[0] = hpar->noff[noffsize];
     depthEta16[1] = hpar->noff[noffsize + 1];
     depthEta29[0] = hpar->noff[noffsize + 2];
     depthEta29[1] = hpar->noff[noffsize + 3];
     if (static_cast<int>(hpar->noff.size()) > (noffsize + 4)) {
       noffl += (2 * hpar->noff[noffsize + 4]);
       if (static_cast<int>(hpar->noff.size()) > noffl)
         isBH_ = (hpar->noff[noffl] > 0);
     }
   } else {
     depthEta16[0] = 2;
     depthEta16[1] = 3;
     depthEta29[0] = 2;
     depthEta29[1] = 1;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "isBH_ " << hpar->noff.size() << ":" << noffsize << ":" << noffl << ":" << isBH_;
   edm::LogVerbatim("HCalGeom") << "Depth index at ieta = 16 for HB (max) " << depthEta16[0] << " HE (min) "
                                << depthEta16[1] << "; max depth for itea = 29 : (" << depthEta29[0] << ":"
                                << depthEta29[1] << ")";
 #endif
 
   if (static_cast<int>(hpar->noff.size()) > (noffsize + 4)) {
     int npair = hpar->noff[noffsize + 4];
     int kk = noffsize + 4;
     for (int k = 0; k < npair; ++k) {
       idHF2QIE.emplace_back(HcalDetId(HcalForward, hpar->noff[kk + 1], hpar->noff[kk + 2], 1));
       kk += 2;
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << idHF2QIE.size() << " detector channels having 2 QIE cards:";
   for (unsigned int k = 0; k < idHF2QIE.size(); ++k)
     edm::LogVerbatim("HCalGeom") << " [" << k << "] " << idHF2QIE[k];
 #endif
 
   layFHB[0] = 0;
   layFHB[1] = 1;
   layBHB[0] = 16;
   layBHB[1] = 15;
   layBHB[2] = 8;
   if (maxDepth[1] == 0) {
     layFHE[0] = layFHE[1] = layFHE[2] = 0;
     layBHE[0] = layBHE[1] = layBHE[2] = layBHE[3] = 0;
   } else {
     layFHE[0] = 1;
     layFHE[1] = 4;
     layFHE[2] = 0;
     layBHE[0] = 18;
     layBHE[1] = 9;
     layBHE[2] = 14;
     layBHE[3] = 16;
   }
   depthMaxSp_ = std::pair<int, int>(0, 0);
   int noffk(noffsize + 5);
   if (static_cast<int>(hpar->noff.size()) > (noffsize + 5)) {
     noffk += (2 * hpar->noff[noffsize + 4]);
     if (static_cast<int>(hpar->noff.size()) >= noffk + 7) {
       int dtype = hpar->noff[noffk + 1];
       int nphi = hpar->noff[noffk + 2];
       int ndeps = hpar->noff[noffk + 3];
       int ndp16 = hpar->noff[noffk + 4];
       int ndp29 = hpar->noff[noffk + 5];
       double wt = 0.1 * (hpar->noff[noffk + 6]);
       if (static_cast<int>(hpar->noff.size()) >= (noffk + 7 + nphi + 3 * ndeps)) {
         if (dtype == 1 || dtype == 2) {
           std::vector<int> ifi, iet, ily, idp;
           for (int i = 0; i < nphi; ++i)
             ifi.emplace_back(hpar->noff[noffk + 7 + i]);
           for (int i = 0; i < ndeps; ++i) {
             iet.emplace_back(hpar->noff[noffk + 7 + nphi + 3 * i]);
             ily.emplace_back(hpar->noff[noffk + 7 + nphi + 3 * i + 1]);
             idp.emplace_back(hpar->noff[noffk + 7 + nphi + 3 * i + 2]);
           }
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("HCalGeom") << "Initialize HcalLayerDepthMap for "
                                        << "Detector " << dtype << " etaMax " << hpar->etaMax[dtype] << " with " << nphi
                                        << " sectors";
           for (int i = 0; i < nphi; ++i)
             edm::LogVerbatim("HCalGeom") << " [" << i << "] " << ifi[i];
           edm::LogVerbatim("HCalGeom") << "And " << ndeps << " depth sections";
           for (int i = 0; i < ndeps; ++i)
             edm::LogVerbatim("HCalGeom") << " [" << i << "] " << iet[i] << "  " << ily[i] << "  " << idp[i];
           edm::LogVerbatim("HCalGeom") << "Maximum depth for last HE Eta tower " << depthEta29[0] << ":" << ndp16 << ":"
                                        << ndp29 << " L0 Wt " << hpar->Layer0Wt[dtype - 1] << ":" << wt;
 #endif
           ldmap_.initialize(dtype, hpar->etaMax[dtype - 1], ndp16, ndp29, wt, ifi, iet, ily, idp);
           int zside = (ifi[0] > 0) ? 1 : -1;
           int iphi = (ifi[0] > 0) ? ifi[0] : -ifi[0];
           depthMaxSp_ = std::pair<int, int>(dtype, ldmap_.getDepthMax(dtype, iphi, zside));
         }
       }
       int noffm = (noffk + 7 + nphi + 3 * ndeps);
       if (static_cast<int>(hpar->noff.size()) > noffm) {
         int ndnext = hpar->noff[noffm];
         if (ndnext > 4 && static_cast<int>(hpar->noff.size()) >= noffm + ndnext) {
           for (int i = 0; i < 2; ++i)
             layFHB[i] = hpar->noff[noffm + i + 1];
           for (int i = 0; i < 3; ++i)
             layFHE[i] = hpar->noff[noffm + i + 3];
         }
         if (ndnext > 11 && static_cast<int>(hpar->noff.size()) >= noffm + ndnext) {
           for (int i = 0; i < 3; ++i)
             layBHB[i] = hpar->noff[noffm + i + 6];
           for (int i = 0; i < 4; ++i)
             layBHE[i] = hpar->noff[noffm + i + 9];
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "Front Layer Definition for HB: " << layFHB[0] << ":" << layFHB[1]
                                << " and for HE: " << layFHE[0] << ":" << layFHE[1] << ":" << layFHE[2];
   edm::LogVerbatim("HCalGeom") << "Last Layer Definition for HB: " << layBHB[0] << ":" << layBHB[1] << ":" << layBHB[2]
                                << " and for HE: " << layBHE[0] << ":" << layBHE[1] << ":" << layBHE[2] << ":"
                                << layBHE[3];
 #endif
   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") << "Detector type and maximum depth for all RBX " << depthMaxDf_.first << ":"
                                << depthMaxDf_.second << " and for special RBX " << depthMaxSp_.first << ":"
                                << depthMaxSp_.second;
 #endif
 }

◆ isBH()

bool HcalDDDSimConstants::isBH ( ) const

inline

Definition at line 64 of file HcalDDDSimConstants.h.

References isBH_.

Referenced by HcalDDDRecConstants::isBH().

64 { return isBH_; }

HcalDDDSimConstants::isBH_

bool isBH_

Definition: HcalDDDSimConstants.h:113

◆ isHE()

bool HcalDDDSimConstants::isHE ( ) const

inline

Definition at line 65 of file HcalDDDSimConstants.h.

References HcalParameters::etaMax, HcalParameters::etaMin, and hpar.

Referenced by HcalDDDRecConstants::isHE(), and HCalSD::setDetUnitId().

65 { return (hpar->etaMax[1] > hpar->etaMin[1]); }

HcalParameters::etaMax

std::vector< int > etaMax

Definition: HcalParameters.h:55

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

HcalParameters::etaMin

std::vector< int > etaMin

Definition: HcalParameters.h:54

◆ layerGroup() [1/2]

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

private

Definition at line 1216 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, hpar, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, isotrackApplyRegressor::k, and HcalParameters::layerGroupEtaSim.

Referenced by findLayer(), getEtaDepth(), getLayerFront(), getMaxDepth(), getMinDepth(), initialize(), and layerGroup().

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

◆ layerGroup() [2/2]

unsigned int HcalDDDSimConstants::layerGroup	(	int	det,
		int	eta,
		int	phi,
		int	zside,
		int	i
	)		const

private

Definition at line 1229 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, findDepth(), layerGroup(), phi, and ecaldqm::zside().

                                                                                                 {
   int depth0 = findDepth(det, eta, phi, zside, lay);
   unsigned int depth = (depth0 > 0) ? static_cast<unsigned int>(depth0) : layerGroup(eta - 1, lay);
   return depth;
 }

◆ layerGroupSize()

unsigned int HcalDDDSimConstants::layerGroupSize ( int eta ) const

private

Definition at line 1203 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, hpar, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, isotrackApplyRegressor::k, and HcalParameters::layerGroupEtaSim.

Referenced by getLayerMax(), getMinDepth(), and initialize().

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

◆ ldMap()

const HcalLayerDepthMap* HcalDDDSimConstants::ldMap ( ) const

inline

Definition at line 66 of file HcalDDDSimConstants.h.

References ldmap_.

Referenced by HcalDDDRecConstants::getDepth(), HcalDDDRecConstants::getEtaBins(), HcalDDDRecConstants::getLayerBack(), HcalDDDRecConstants::getLayerFront(), HcalDDDRecConstants::getPhiZOne(), HcalDDDRecConstants::initialize(), and HcalDDDRecConstants::withSpecialRBXHBHE().

66 { return &ldmap_; }

HcalDDDSimConstants::ldmap_

HcalLayerDepthMap ldmap_

Definition: HcalDDDSimConstants.h:93

◆ maxHFDepth()

int HcalDDDSimConstants::maxHFDepth	(	const int &	ieta,
		const int &	iphi
	)		const

Definition at line 657 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, idHF2QIE, hcalRecHitTable_cff::iphi, isotrackApplyRegressor::k, maxDepth, and convertSQLiteXML::ok.

Referenced by getMaxDepth(), HcalDDDRecConstants::maxHFDepth(), FastHFShowerLibrary::modifyDepth(), and HCalSD::modifyDepth().

                                                                          {
   int mxdepth = maxDepth[2];
   if (!idHF2QIE.empty()) {
     bool ok(false);
     for (auto k : idHF2QIE) {
       if ((eta == k.ieta()) && (iphi == k.iphi())) {
         ok = true;
         break;
       }
     }
     if (!ok)
       mxdepth = 2;
   }
   return mxdepth;
 }

◆ numberOfCells()

unsigned int HcalDDDSimConstants::numberOfCells ( const HcalSubdetector & subdet ) const

Definition at line 673 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), and EgammaValidation_cff::num.

Referenced by HcalDDDRecConstants::nCells(), and HcalDDDRecConstants::numberOfCells().

                                                                                    {
   unsigned int num = 0;
   std::vector<HcalCellType> cellTypes = HcalCellTypes(subdet);
   for (auto& cellType : cellTypes) {
     num += static_cast<unsigned int>(cellType.nPhiBins());
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants:numberOfCells " << cellTypes.size() << " " << num
                                << " for subdetector " << subdet;
 #endif
   return num;
 }

◆ parameter()

const HcalParameters* HcalDDDSimConstants::parameter ( ) const

inline

Definition at line 69 of file HcalDDDSimConstants.h.

References hpar.

69 { return hpar; }

HcalDDDSimConstants::hpar

const HcalParameters * hpar

Definition: HcalDDDSimConstants.h:92

◆ phiNumber()

int HcalDDDSimConstants::phiNumber	(	const int &	phi,
		const int &	unit
	)		const

Definition at line 686 of file HcalDDDSimConstants.cc.

References phi, and units().

Referenced by HcalDDDRecConstants::getHCID(), getPhis(), HcalDDDRecConstants::getPhis(), and HcalNumberingFromDDD::unitID().

                                                                          {
   int iphi_skip = phi;
   if (units == 2)
     iphi_skip = (phi - 1) * 2 + 1;
   else if (units == 4)
     iphi_skip = (phi - 1) * 4 - 1;
   if (iphi_skip < 0)
     iphi_skip += 72;
   return iphi_skip;
 }

◆ printTileHB()

void HcalDDDSimConstants::printTileHB	(	const int &	eta,
		const int &	phi,
		const int &	zside,
		const int &	depth
	)		const

private

Definition at line 1104 of file HcalDDDSimConstants.cc.

References custom_jme_cff::area, angle_units::operators::convertRadToDeg(), funct::cos(), hcalRecHitTable_cff::depth, HcalParameters::dxHB, HcalParameters::dyHB, PVValHelper::dz, PVValHelper::eta, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, isotrackApplyRegressor::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHB, SiStripPI::mean, SiStripPI::min, phi, HcalParameters::rhoxHB, funct::sin(), SiStripMonitorCluster_cfi::zmax, SiStripMonitorCluster_cfi::zmin, and ecaldqm::zside().

Referenced by printTiles().

                                                                                                               {
   edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants::printTileHB for eta " << eta << " and depth " << depth;
 
   double etaL = hpar->etaTable.at(eta - 1);
   double thetaL = 2. * atan(exp(-etaL));
   double etaH = hpar->etaTable.at(eta);
   double thetaH = 2. * atan(exp(-etaH));
   int layL = getLayerFront(1, eta, phi, zside, depth);
   int layH = getLayerBack(1, eta, phi, zside, depth);
   edm::LogVerbatim("HCalGeom") << "\ntileHB:: eta|depth " << zside * eta << "|" << depth << " theta "
                                << convertRadToDeg(thetaH) << ":" << convertRadToDeg(thetaL) << " Layer " << layL - 1
                                << ":" << layH - 1;
   for (int lay = layL - 1; lay < layH; ++lay) {
     std::vector<double> area(2, 0);
     int kk(0);
     double mean(0);
     for (unsigned int k = 0; k < hpar->layHB.size(); ++k) {
       if (lay == hpar->layHB[k]) {
         double zmin = hpar->rhoxHB[k] * std::cos(thetaL) / std::sin(thetaL);
         double zmax = hpar->rhoxHB[k] * std::cos(thetaH) / std::sin(thetaH);
         double dz = (std::min(zmax, hpar->dxHB[k]) - zmin);
         if (dz > 0) {
           area[kk] = dz * hpar->dyHB[k];
           mean += area[kk];
           kk++;
         }
       }
     }
     if (area[0] > 0) {
       mean /= (kk * 100);
       edm::LogVerbatim("HCalGeom") << std::setw(2) << lay << " Area " << std::setw(8) << area[0] << " " << std::setw(8)
                                    << area[1] << " Mean " << mean;
     }
   }
 }

◆ printTileHE()

void HcalDDDSimConstants::printTileHE	(	const int &	eta,
		const int &	phi,
		const int &	zside,
		const int &	depth
	)		const

private

Definition at line 1140 of file HcalDDDSimConstants.cc.

References custom_jme_cff::area, angle_units::operators::convertRadToDeg(), hcalRecHitTable_cff::depth, HcalParameters::dx1HE, HcalParameters::dyHE, PVValHelper::eta, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, isotrackApplyRegressor::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHE, WZElectronSkims53X_cff::max, SiStripPI::mean, SiStripPI::min, nphi, phi, HcalParameters::phibin, HcalParameters::rhoxHE, funct::tan(), ecaldqm::zside(), and HcalParameters::zxHE.

Referenced by printTiles().

                                                                                                               {
   double etaL = hpar->etaTable[eta - 1];
   double thetaL = 2. * atan(exp(-etaL));
   double etaH = hpar->etaTable[eta];
   double thetaH = 2. * atan(exp(-etaH));
   int layL = getLayerFront(2, eta, phi, zside, depth);
   int layH = getLayerBack(2, eta, phi, zside, depth);
   double phib = hpar->phibin[eta - 1];
   int nphi = 2;
   if (phib > 6._deg)
     nphi = 1;
   edm::LogVerbatim("HCalGeom") << "\ntileHE:: Eta/depth " << zside * eta << "|" << depth << " theta "
                                << convertRadToDeg(thetaH) << ":" << convertRadToDeg(thetaL) << " Layer " << layL - 1
                                << ":" << layH - 1 << " phi " << nphi;
   for (int lay = layL - 1; lay < layH; ++lay) {
     std::vector<double> area(4, 0);
     int kk(0);
     double mean(0);
     for (unsigned int k = 0; k < hpar->layHE.size(); ++k) {
       if (lay == hpar->layHE[k]) {
         double rmin = hpar->zxHE[k] * std::tan(thetaH);
         double rmax = hpar->zxHE[k] * std::tan(thetaL);
         if ((lay != 0 || eta == 18) &&
             (lay != 1 || (eta == 18 && hpar->rhoxHE[k] - hpar->dyHE[k] > 1000) ||
              (eta != 18 && hpar->rhoxHE[k] - hpar->dyHE[k] < 1000)) &&
             rmin + 30 < hpar->rhoxHE[k] + hpar->dyHE[k] && rmax > hpar->rhoxHE[k] - hpar->dyHE[k]) {
           rmin = std::max(rmin, hpar->rhoxHE[k] - hpar->dyHE[k]);
           rmax = std::min(rmax, hpar->rhoxHE[k] + hpar->dyHE[k]);
           double dx1 = rmin * std::tan(phib);
           double dx2 = rmax * std::tan(phib);
           double ar1 = 0, ar2 = 0;
           if (nphi == 1) {
             ar1 = 0.5 * (rmax - rmin) * (dx1 + dx2 - 4. * hpar->dx1HE[k]);
             mean += ar1;
           } else {
             ar1 = 0.5 * (rmax - rmin) * (dx1 + dx2 - 2. * hpar->dx1HE[k]);
             ar2 = 0.5 * (rmax - rmin) * ((rmax + rmin) * tan(10._deg) - 4 * hpar->dx1HE[k]) - ar1;
             mean += (ar1 + ar2);
           }
           area[kk] = ar1;
           area[kk + 2] = ar2;
           kk++;
         }
       }
     }
     if (area[0] > 0 && area[1] > 0) {
       int lay0 = lay - 1;
       if (eta == 18)
         lay0++;
       if (nphi == 1) {
         mean /= (kk * 100);
         edm::LogVerbatim("HCalGeom") << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " "
                                      << std::setw(8) << area[1] << " Mean " << mean;
       } else {
         mean /= (kk * 200);
         edm::LogVerbatim("HCalGeom") << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " "
                                      << std::setw(8) << area[1] << ":" << std::setw(8) << area[2] << " " << std::setw(8)
                                      << area[3] << " Mean " << mean;
       }
     }
   }
 }

◆ printTiles()

void HcalDDDSimConstants::printTiles ( ) const

Definition at line 697 of file HcalDDDSimConstants.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, getDepthEta16(), getMaxDepth(), hpar, hcalRecHitTable_cff::iphi, ldmap_, printTileHB(), printTileHE(), HcalLayerDepthMap::validDet(), and ecaldqm::zside().

                                            {
   std::vector<int> phis;
   int detsp = ldmap_.validDet(phis);
   int kphi = (detsp > 0) ? phis[0] : 1;
   int zside = (kphi > 0) ? 1 : -1;
   int iphi = (kphi > 0) ? kphi : -kphi;
   edm::LogVerbatim("HCalGeom") << "Tile Information for HB from " << hpar->etaMin[0] << " to " << hpar->etaMax[0];
   for (int eta = hpar->etaMin[0]; eta <= hpar->etaMax[0]; eta++) {
     int dmax = getMaxDepth(1, eta, iphi, -zside, false);
     for (int depth = 1; depth <= dmax; depth++)
       printTileHB(eta, iphi, -zside, depth);
     if (detsp == 1) {
       int dmax = getMaxDepth(1, eta, iphi, zside, false);
       for (int depth = 1; depth <= dmax; depth++)
         printTileHB(eta, iphi, zside, depth);
     }
   }
 
   edm::LogVerbatim("HCalGeom") << "\nTile Information for HE from " << hpar->etaMin[1] << " to " << hpar->etaMax[1];
   for (int eta = hpar->etaMin[1]; eta <= hpar->etaMax[1]; eta++) {
     int dmin = (eta == hpar->etaMin[1]) ? getDepthEta16(2, iphi, -zside) : 1;
     int dmax = getMaxDepth(2, eta, iphi, -zside, false);
     for (int depth = dmin; depth <= dmax; depth++)
       printTileHE(eta, iphi, -zside, depth);
     if (detsp == 2) {
       int dmax = getMaxDepth(2, eta, iphi, zside, false);
       for (int depth = 1; depth <= dmax; depth++)
         printTileHE(eta, iphi, zside, depth);
     }
   }
 }

◆ unitPhi() [1/2]

int HcalDDDSimConstants::unitPhi	(	const int &	det,
		const int &	etaR
	)		const

Definition at line 729 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, HcalForward, hpar, HcalParameters::phibin, and HcalParameters::phitable.

Referenced by cell(), HcalDDDRecConstants::getEtaPhi(), HcalDDDRecConstants::getHCID(), getPhiCons(), getPhis(), HcalDDDRecConstants::getPhis(), HcalCellTypes(), HcalDDDRecConstants::HcalCellTypes(), HcalDDDRecConstants::initialize(), and HcalNumberingFromDDD::unitID().

                                                                       {
   double dphi =
       (det == static_cast<int>(HcalForward)) ? hpar->phitable[etaR - hpar->etaMin[2]] : hpar->phibin[etaR - 1];
   return unitPhi(dphi);
 }

◆ unitPhi() [2/2]

int HcalDDDSimConstants::unitPhi ( const double & dphi ) const

Definition at line 735 of file HcalDDDSimConstants.cc.

References createfilelist::int, M_PI, and units().

                                                          {
   const double fiveDegInRad = 2 * M_PI / 72;
   int units = int(dphi / fiveDegInRad + 0.5);
   if (units < 1)
     units = 1;
   return units;
 }

Member Data Documentation

◆ depthEta16

int HcalDDDSimConstants::depthEta16[2]

private

Definition at line 107 of file HcalDDDSimConstants.h.

Referenced by getDepthEta16(), getDepthEta16M(), getDetEta(), and initialize().

◆ depthEta29

int HcalDDDSimConstants::depthEta29[2]

private

Definition at line 108 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getDepthEta29(), getDepthEta29M(), getEta(), and initialize().

◆ depthMaxDf_

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

private

Definition at line 115 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

◆ depthMaxSp_

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

private

Definition at line 115 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

◆ depths

std::vector<int> HcalDDDSimConstants::depths[nDepthMax]

private

Definition at line 99 of file HcalDDDSimConstants.h.

Referenced by getLayerBack(), and initialize().

◆ dlShort

double HcalDDDSimConstants::dlShort

private

Definition at line 106 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getEta(), HcalCellTypes(), and initialize().

◆ dzVcal

double HcalDDDSimConstants::dzVcal

private

Definition at line 105 of file HcalDDDSimConstants.h.

Referenced by HcalCellTypes(), and initialize().

◆ etaHO

double HcalDDDSimConstants::etaHO[4]

private

Definition at line 103 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getEta(), and initialize().

◆ hpar

const HcalParameters* HcalDDDSimConstants::hpar

private

Definition at line 92 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getConstHBHE(), getDetEta(), getEta(), getEtaDepth(), getEtaHO(), getEtaTableHF(), getGain(), getGparHF(), getiEtaRange(), getLayer0Wt(), getLayerFront(), getLayerMax(), getMaxDepth(), getMinDepth(), getPhiCons(), getPhiTableHF(), getRTableHF(), getShift(), HcalCellTypes(), initialize(), isHE(), layerGroup(), layerGroupSize(), parameter(), printTileHB(), printTileHE(), printTiles(), and unitPhi().

◆ idHF2QIE

std::vector<HcalDetId> HcalDDDSimConstants::idHF2QIE

private

Definition at line 114 of file HcalDDDSimConstants.h.

Referenced by getIdHF2QIE(), HcalCellTypes(), initialize(), and maxHFDepth().

◆ isBH_

bool HcalDDDSimConstants::isBH_

private

Definition at line 113 of file HcalDDDSimConstants.h.

Referenced by initialize(), and isBH().

◆ layBHB

int HcalDDDSimConstants::layBHB[3]

private

Definition at line 111 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

◆ layBHE

int HcalDDDSimConstants::layBHE[4]

private

Definition at line 112 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

◆ layFHB

int HcalDDDSimConstants::layFHB[2]

private

Definition at line 109 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

◆ layFHE

int HcalDDDSimConstants::layFHE[3]

private

Definition at line 110 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

◆ ldmap_

HcalLayerDepthMap HcalDDDSimConstants::ldmap_

private

Definition at line 93 of file HcalDDDSimConstants.h.

Referenced by findDepth(), getDepthEta16(), getDepthEta16M(), getDepthEta29(), getDepthEta29M(), getLayer0Wt(), getLayerBack(), getLayerFront(), getMaxDepth(), getMinDepth(), initialize(), ldMap(), and printTiles().

◆ maxDepth

std::vector<int> HcalDDDSimConstants::maxDepth

private

Definition at line 95 of file HcalDDDSimConstants.h.

Referenced by cell(), getLayerMax(), getMaxDepth(), getMinDepth(), HcalCellTypes(), initialize(), and maxHFDepth().

◆ maxLayer_

const int HcalDDDSimConstants::maxLayer_ = 18

staticprivate

Definition at line 77 of file HcalDDDSimConstants.h.

Referenced by getEta(), getLayerBack(), getLayerFront(), getMaxDepth(), and HcalCellTypes().

◆ maxLayerHB_

const int HcalDDDSimConstants::maxLayerHB_ = 16

staticprivate

Definition at line 78 of file HcalDDDSimConstants.h.

Referenced by cell(), getLayerFront(), getLayerMax(), getMaxDepth(), HcalCellTypes(), and initialize().

◆ nDepth

int HcalDDDSimConstants::nDepth

private

Definition at line 100 of file HcalDDDSimConstants.h.

◆ nDepthMax

const int HcalDDDSimConstants::nDepthMax = 9

staticprivate

Definition at line 76 of file HcalDDDSimConstants.h.

◆ nEta

int HcalDDDSimConstants::nEta

private

Definition at line 96 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getDetEta(), getEta(), and initialize().

◆ nmodHB

int HcalDDDSimConstants::nmodHB

private

Definition at line 101 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

◆ nmodHE

int HcalDDDSimConstants::nmodHE

private

Definition at line 102 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

◆ nPhiF

int HcalDDDSimConstants::nPhiF

private

Definition at line 98 of file HcalDDDSimConstants.h.

Referenced by initialize().

◆ nR

int HcalDDDSimConstants::nR

private

Definition at line 97 of file HcalDDDSimConstants.h.

Referenced by cell(), deltaEta(), getEta(), and initialize().

◆ nzHB

int HcalDDDSimConstants::nzHB

private

Definition at line 101 of file HcalDDDSimConstants.h.

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

◆ nzHE

int HcalDDDSimConstants::nzHE

private

Definition at line 102 of file HcalDDDSimConstants.h.

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

◆ rminHO

double HcalDDDSimConstants::rminHO

private

Definition at line 103 of file HcalDDDSimConstants.h.

Referenced by getEtaHO(), and initialize().

◆ zVcal

double HcalDDDSimConstants::zVcal

private

Definition at line 104 of file HcalDDDSimConstants.h.

Referenced by deltaEta(), getDetEta(), getEta(), and initialize().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HcalDDDSimConstants()

◆ ~HcalDDDSimConstants()

Member Function Documentation

◆ cell()

◆ deltaEta()

◆ findDepth()

◆ findLayer()

◆ getConstHBHE()

◆ getDepthEta16()

◆ getDepthEta16M()

◆ getDepthEta29()

◆ getDepthEta29M()

◆ getDetEta()

◆ getEta() [1/3]

◆ getEta() [2/3]

◆ getEta() [3/3]

◆ getEtaDepth()

◆ getEtaHO()

◆ getEtaTableHF()

◆ getFrontLayer()

◆ getGain()

◆ getGparHF()

◆ getIdHF2QIE()

◆ getiEtaRange()

◆ getLastLayer()

◆ getLayer0Wt()

◆ getLayerBack()

◆ getLayerFront()

◆ getLayerMax()

◆ getMaxDepth() [1/2]

◆ getMaxDepth() [2/2]

◆ getMaxDepthDet()

◆ getMinDepth()

◆ getModHalfHBHE()

◆ getPhiCons()

◆ getPhis()

◆ getPhiTableHF()

◆ getRTableHF()

◆ getShift()

◆ HcalCellTypes() [1/2]

◆ HcalCellTypes() [2/2]

◆ initialize()

◆ isBH()

◆ isHE()

◆ layerGroup() [1/2]

◆ layerGroup() [2/2]

◆ layerGroupSize()

◆ ldMap()

◆ maxHFDepth()

◆ numberOfCells()

◆ parameter()

◆ phiNumber()

◆ printTileHB()

◆ printTileHE()

◆ printTiles()

◆ unitPhi() [1/2]

◆ unitPhi() [2/2]

Member Data Documentation

◆ depthEta16

◆ depthEta29

◆ depthMaxDf_

◆ depthMaxSp_

◆ depths

◆ dlShort

◆ dzVcal

◆ etaHO

◆ hpar

◆ idHF2QIE

◆ isBH_

◆ layBHB

◆ layBHE

◆ layFHB

◆ layFHE

◆ ldmap_