HcalDDDSimConstants Class Reference

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

(

const HcalParameters *

hp

)

Definition at line 10 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), and initialize().

                                                                 : hpar(hp) {
#ifdef EDM_ML_DEBUG
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::HcalDDDSimConstants (const HcalParameters* hp) constructor\n";
#endif

  initialize();
#ifdef EDM_ML_DEBUG
  std::vector<HcalCellType> cellTypes = HcalCellTypes();
  edm::LogVerbatim("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size() << " cells of type HCal (All)\n";
#endif
}

HcalDDDSimConstants::~HcalDDDSimConstants

(

)

Definition at line 22 of file HcalDDDSimConstants.cc.

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

Member Function Documentation

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(), HcalCellType::HcalCell::deta, HcalCellType::HcalCell::dphi, HcalParameters::drHB, HcalCellType::HcalCell::drz, HcalParameters::dzHE, HcalCellType::HcalCell::eta, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, HcalCellType::HcalCell::flagrz, getEta(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, maxDepth, maxLayerHB_, nEta, HcalParameters::noff, nR, HcalCellType::HcalCell::ok, convertSQLiteXML::ok, phi, HcalCellType::HcalCell::phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rHB, HcalParameters::rTable, HcalCellType::HcalCell::rz, createJobs::tmp, unitPhi(), and HcalParameters::zHE.

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 > (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::LogInfo("HCalGeom") << "HcalDDDSimConstants: wrong eta " << etaR << " (" << ir << "/" << nR
                                 << ") Detector " << idet << std::endl;
#endif
      }
    } else if (etaR <= nEta) {
      int laymin(depth), laymax(depth);
      if (idet == static_cast<int>(HcalOuter)) {
        laymin = (etaR > hpar->noff[2]) ? ((int)(hpar->zHE.size())) : ((int)(hpar->zHE.size())) - 1;
        laymax = ((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;
}

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

private

Definition at line 979 of file HcalDDDSimConstants.cc.

References 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;
}

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

Definition at line 120 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, HcalLayerDepthMap::getDepth(), HcalLayerDepthMap::isValid(), and ldmap_.

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

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

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

Definition at line 126 of file HcalDDDSimConstants.cc.

References mps_fire::i, and gpuClustering::id.

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

std::vector< std::pair< double, double > > HcalDDDSimConstants::getConstHBHE

(

const int &

type

)

const

Definition at line 138 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;
}

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

Definition at line 152 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, depthEta16, HcalLayerDepthMap::getDepth16(), and ldmap_.

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

int HcalDDDSimConstants::getDepthEta16M

(

const int &

det

)

const

Definition at line 162 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, depthEta16, HcalLayerDepthMap::getDepth16(), 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;
}

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

Definition at line 178 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), and ldmap_.

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

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

Definition at line 185 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), 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;
}

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

Definition at line 201 of file HcalDDDSimConstants.cc.

References depthEta16, PVValHelper::eta, Puppi_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, generateTowerEtThresholdLUT::etaRange, HcalParameters::etaTable, funct::exp(), HcalBarrel, HcalEndcap, HcalForward, hpar, mps_fire::i, nEta, funct::tan(), 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);
}

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

Definition at line 225 of file HcalDDDSimConstants.cc.

References etaHO, HcalParameters::etaMax, Puppi_cff::etaMin, HcalParameters::etaMin, HcalParameters::etaTable, HcalBarrel, HcalEndcap, HcalForward, hpar, mps_fire::i, 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;
}

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 depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, createJobs::tmp, z, 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;
}

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

private

Definition at line 1056 of file HcalDDDSimConstants.cc.

References log, funct::tan(), and createJobs::tmp.

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

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 252 of file HcalDDDSimConstants.cc.

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

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 (etaR == hpar->noff[0] && lay > 1) {
        int kphi = phi + 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);
}

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

Definition at line 295 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, HcalParameters::etaTable, hpar, alignCSCRings::r, rminHO, mathSSE::sqrt(), and HcalParameters::zHO.

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
    edm::LogVerbatim("HCalGeom") << "R " << r << " Z " << z << " eta " << etaR << ":" << eta;
    if (eta != etaR)
      edm::LogInfo("HCalGeom") << "**** Check *****";
#endif
    return eta;
  } else {
    return etaR;
  }
}

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

inline

Definition at line 44 of file HcalDDDSimConstants.h.

References HcalParameters::etaTableHF, and hpar.

{ return hpar->etaTableHF; }

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

Definition at line 321 of file HcalDDDSimConstants.cc.

References funct::abs(), 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;
}

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

private

Definition at line 1085 of file HcalDDDSimConstants.cc.

References 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;
}

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(), HFFibre::HFFibre(), and HFShowerParam::HFShowerParam().

{ return hpar->gparHF; }

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

inline

Definition at line 46 of file HcalDDDSimConstants.h.

References idHF2QIE.

Referenced by HcalDDDRecConstants::getHFCellParameters().

{ return idHF2QIE; }

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.

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

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

Definition at line 339 of file HcalDDDSimConstants.cc.

References funct::abs(), 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;
}

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

Definition at line 361 of file HcalDDDSimConstants.cc.

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

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

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

Definition at line 387 of file HcalDDDSimConstants.cc.

References depths, HcalLayerDepthMap::getLayerBack(), phase1PixelTopology::layer, ldmap_, and maxLayer_.

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

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

Definition at line 368 of file HcalDDDSimConstants.cc.

References HcalLayerDepthMap::getLayerFront(), getLayerMax(), hpar, cmsLHEtoEOSManager::l, phase1PixelTopology::layer, layerGroup(), ldmap_, maxLayer_, maxLayerHB_, and HcalParameters::noff.

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 ((int)(layerGroup(eta - 1, l)) == depth) {
          layer = l + 1;
          break;
        }
      }
    } else {
      layer = (eta > hpar->noff[2]) ? maxLayerHB_ + 1 : maxLayer_;
    }
  }
  return layer;
}

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

Definition at line 400 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, hpar, layerGroupSize(), maxDepth, and maxLayerHB_.

Referenced by getLayerFront(), and initialize().

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

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().

{ return maxDepth[type]; }

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

Definition at line 405 of file HcalDDDSimConstants.cc.

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

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

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

inline

Definition at line 55 of file HcalDDDSimConstants.h.

References depthMaxDf_, and depthMaxSp_.

Referenced by HcalDDDRecConstants::initialize().

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

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

Definition at line 431 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, getDepthEta16(), HcalLayerDepthMap::getDepths(), hpar, HcalLayerDepthMap::isValid(), layerGroup(), layerGroupSize(), ldmap_, and maxDepth.

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 = (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;
}

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

(

const int &

type

)

const

Definition at line 457 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);
  }
}

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

Definition at line 465 of file HcalDDDSimConstants.cc.

References HcalParameters::etaMin, HcalEndcap, HcalForward, hpar, 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);
}

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

Definition at line 484 of file HcalDDDSimConstants.cc.

References angle_units::operators::convertRadToDeg(), getPhiCons(), 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;
}

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().

{ return hpar->phitable; }

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(), HFFibre::HFFibre(), HFShowerFibreBundle::HFShowerFibreBundle(), and HFShowerPMT::HFShowerPMT().

{ return hpar->rTable; }

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

private

Definition at line 1066 of file HcalDDDSimConstants.cc.

References 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;
}

std::vector< HcalCellType > HcalDDDSimConstants::HcalCellTypes

(

)

const

Definition at line 505 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::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << cellTypes.size() << " cells of type HCal Barrel\n";
  for (unsigned int i = 0; i < cellTypes.size(); i++)
    edm::LogInfo("HCalGeom") << "Cell " << i << " " << cellTypes[i] << "\n";
#endif

  std::vector<HcalCellType> hoCells = HcalCellTypes(HcalOuter);
#ifdef EDM_ML_DEBUG
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << hoCells.size() << " cells of type HCal Outer\n";
  for (unsigned int i = 0; i < hoCells.size(); i++)
    edm::LogInfo("HCalGeom") << "Cell " << i << " " << hoCells[i] << "\n";
#endif
  cellTypes.insert(cellTypes.end(), hoCells.begin(), hoCells.end());

  std::vector<HcalCellType> heCells = HcalCellTypes(HcalEndcap);
#ifdef EDM_ML_DEBUG
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << heCells.size() << " cells of type HCal Endcap\n";
  for (unsigned int i = 0; i < heCells.size(); i++)
    edm::LogInfo("HCalGeom") << "Cell " << i << " " << heCells[i] << "\n";
#endif
  cellTypes.insert(cellTypes.end(), heCells.begin(), heCells.end());

  std::vector<HcalCellType> hfCells = HcalCellTypes(HcalForward);
#ifdef EDM_ML_DEBUG
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants: " << hfCells.size() << " cells of type HCal Forward\n";
  for (unsigned int i = 0; i < hfCells.size(); i++)
    edm::LogInfo("HCalGeom") << "Cell " << i << " " << hfCells[i] << "\n";
#endif
  cellTypes.insert(cellTypes.end(), hfCells.begin(), hfCells.end());

  return cellTypes;
}

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

Definition at line 540 of file HcalDDDSimConstants.cc.

References cell(), HLT_FULL_cff::depth, dlShort, dzVcal, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, getGain(), getPhis(), getShift(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, idHF2QIE, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::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;
}

void HcalDDDSimConstants::initialize ( void  )

private

Definition at line 743 of file HcalDDDSimConstants.cc.

References depthEta16, depthEta29, depthMaxDf_, depthMaxSp_, depths, dlShort, HcalParameters::dzVcal, dzVcal, etaHO, Puppi_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaTable, for(), HcalLayerDepthMap::getDepthMax(), getLayerMax(), HcalParameters::gparHF, HcalForward, hpar, mps_fire::i, idHF2QIE, HcalLayerDepthMap::initialize(), isBH_, isotrackApplyRegressor::k, GetRecoTauVFromDQM_MC_cff::kk, cmsLHEtoEOSManager::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 ((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 ((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 ((int)(hpar->noff.size()) > (noffsize + 4)) {
      noffl += (2 * hpar->noff[noffsize + 4]);
      if ((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 ((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 ((int)(hpar->noff.size()) > (noffsize + 5)) {
    noffk += (2 * hpar->noff[noffsize + 4]);
    if ((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 ((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 ((int)(hpar->noff.size()) > noffm) {
        int ndnext = hpar->noff[noffm];
        if (ndnext > 4 && (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 && (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
}

bool HcalDDDSimConstants::isBH ( ) const

inline

Definition at line 64 of file HcalDDDSimConstants.h.

References isBH_.

Referenced by HcalDDDRecConstants::isBH().

{ return isBH_; }

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

private

Definition at line 1216 of file HcalDDDSimConstants.cc.

References hpar, isotrackApplyRegressor::k, and HcalParameters::layerGroupEtaSim.

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

                                                                 {
  unsigned int k = 0;
  for (auto const& it : hpar->layerGroupEtaSim) {
    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;
}

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 HLT_FULL_cff::depth, findDepth(), and layerGroup().

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

unsigned int HcalDDDSimConstants::layerGroupSize ( int  eta ) const

private

Definition at line 1203 of file HcalDDDSimConstants.cc.

References hpar, isotrackApplyRegressor::k, and HcalParameters::layerGroupEtaSim.

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

                                                              {
  unsigned int k = 0;
  for (auto const& it : hpar->layerGroupEtaSim) {
    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;
}

const HcalLayerDepthMap* HcalDDDSimConstants::ldMap ( ) const

inline

Definition at line 65 of file HcalDDDSimConstants.h.

References ldmap_.

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

{ return &ldmap_; }

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

Definition at line 655 of file HcalDDDSimConstants.cc.

References idHF2QIE, 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;
}

unsigned int HcalDDDSimConstants::numberOfCells

(

const HcalSubdetector &

subdet

)

const

Definition at line 671 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), and pileupDistInMC::num.

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

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

const HcalParameters* HcalDDDSimConstants::parameter ( ) const

inline

Definition at line 68 of file HcalDDDSimConstants.h.

References hpar.

{ return hpar; }

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

Definition at line 684 of file HcalDDDSimConstants.cc.

References phi.

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

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 angle_units::operators::convertRadToDeg(), funct::cos(), HLT_FULL_cff::depth, HcalParameters::dxHB, HcalParameters::dyHB, PVValHelper::dz, HcalParameters::etaTable, funct::exp(), getLayerBack(), getLayerFront(), hpar, isotrackApplyRegressor::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHB, SiStripPI::mean, SiStripPI::min, HcalParameters::rhoxHB, and funct::sin().

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

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 angle_units::operators::convertRadToDeg(), HcalParameters::dx1HE, HcalParameters::dyHE, PVValHelper::eta, HcalParameters::etaTable, funct::exp(), getLayerBack(), getLayerFront(), hpar, isotrackApplyRegressor::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHE, SiStripPI::max, SiStripPI::mean, SiStripPI::min, nphi, HcalParameters::phibin, HcalParameters::rhoxHE, funct::tan(), 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;
      }
    }
  }
}

void HcalDDDSimConstants::printTiles

(

)

const

Definition at line 695 of file HcalDDDSimConstants.cc.

References HLT_FULL_cff::depth, PVValHelper::eta, Puppi_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, getDepthEta16(), getMaxDepth(), hpar, 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]
                               << "\n";
  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]
                               << "\n";
  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);
    }
  }
}

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

int HcalDDDSimConstants::unitPhi

(

const double &

dphi

)

const

Definition at line 735 of file HcalDDDSimConstants.cc.

References 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

int HcalDDDSimConstants::depthEta16[2]

private

Definition at line 106 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::depthEta29[2]

private

Definition at line 107 of file HcalDDDSimConstants.h.

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

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

private

Definition at line 114 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

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

private

Definition at line 114 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

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

private

Definition at line 98 of file HcalDDDSimConstants.h.

Referenced by getLayerBack(), and initialize().

double HcalDDDSimConstants::dlShort

private

Definition at line 105 of file HcalDDDSimConstants.h.

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

double HcalDDDSimConstants::dzVcal

private

Definition at line 104 of file HcalDDDSimConstants.h.

Referenced by HcalCellTypes(), and initialize().

double HcalDDDSimConstants::etaHO[4]

private

Definition at line 102 of file HcalDDDSimConstants.h.

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

const HcalParameters* HcalDDDSimConstants::hpar

private

Definition at line 91 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(), layerGroup(), layerGroupSize(), parameter(), printTileHB(), printTileHE(), printTiles(), and unitPhi().

std::vector<HcalDetId> HcalDDDSimConstants::idHF2QIE

private

Definition at line 113 of file HcalDDDSimConstants.h.

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

bool HcalDDDSimConstants::isBH_

private

Definition at line 112 of file HcalDDDSimConstants.h.

Referenced by initialize(), and isBH().

int HcalDDDSimConstants::layBHB[3]

private

Definition at line 110 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

int HcalDDDSimConstants::layBHE[4]

private

Definition at line 111 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

int HcalDDDSimConstants::layFHB[2]

private

Definition at line 108 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

int HcalDDDSimConstants::layFHE[3]

private

Definition at line 109 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

HcalLayerDepthMap HcalDDDSimConstants::ldmap_

private

Definition at line 92 of file HcalDDDSimConstants.h.

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

std::vector<int> HcalDDDSimConstants::maxDepth

private

Definition at line 94 of file HcalDDDSimConstants.h.

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

const int HcalDDDSimConstants::maxLayer_ = 18

staticprivate

Definition at line 76 of file HcalDDDSimConstants.h.

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

const int HcalDDDSimConstants::maxLayerHB_ = 16

staticprivate

Definition at line 77 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nDepth

private

Definition at line 99 of file HcalDDDSimConstants.h.

const int HcalDDDSimConstants::nDepthMax = 9

staticprivate

Definition at line 75 of file HcalDDDSimConstants.h.

int HcalDDDSimConstants::nEta

private

Definition at line 95 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nmodHB

private

Definition at line 100 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

int HcalDDDSimConstants::nmodHE

private

Definition at line 101 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

int HcalDDDSimConstants::nPhiF

private

Definition at line 97 of file HcalDDDSimConstants.h.

Referenced by initialize().

int HcalDDDSimConstants::nR

private

Definition at line 96 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nzHB

private

Definition at line 100 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nzHE

private

Definition at line 101 of file HcalDDDSimConstants.h.

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

double HcalDDDSimConstants::rminHO

private

Definition at line 102 of file HcalDDDSimConstants.h.

Referenced by getEtaHO(), and initialize().

double HcalDDDSimConstants::zVcal

private

Definition at line 103 of file HcalDDDSimConstants.h.

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