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)
int	getEta (const int &det, const int &lay, const double &hetaR)
std::pair< int, int >	getEtaDepth (const int &det, int etaR, const int &phi, const int &zside, int depth, const int &lay)
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
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 12 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), and initialize().

                                                                 : hpar(hp) {

#ifdef EDM_ML_DEBUG
  edm::LogInfo("HCalGeom") << "HcalDDDSimConstants::HcalDDDSimConstants (const HcalParameter* 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 26 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 32 of file HcalDDDSimConstants.cc.

References 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, createfilelist::int, 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, tmp, unitPhi(), and HcalParameters::zHE.

Referenced by HGCGeometryValidation::analyze(), HGCalHitValidation::analyze(), 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 980 of file HcalDDDSimConstants.cc.

References depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, 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 125 of file HcalDDDSimConstants.cc.

References egammaForCoreTracking_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 133 of file HcalDDDSimConstants.cc.

References mps_fire::i, and triggerObjects_cff::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 145 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 160 of file HcalDDDSimConstants.cc.

References egammaForCoreTracking_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 171 of file HcalDDDSimConstants.cc.

References egammaForCoreTracking_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 185 of file HcalDDDSimConstants.cc.

References egammaForCoreTracking_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 193 of file HcalDDDSimConstants.cc.

References egammaForCoreTracking_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
	)

Definition at line 209 of file HcalDDDSimConstants.cc.

References depthEta16, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, MuonErrorMatrixValues_cff::etaRange, HcalParameters::etaTable, JetChargeProducer_cfi::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])) {
      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
	)

Definition at line 234 of file HcalDDDSimConstants.cc.

References etaHO, HcalParameters::etaMax, ALCARECOTkAlBeamHalo_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 1019 of file HcalDDDSimConstants.cc.

References depthEta29, dlShort, etaHO, HcalParameters::etaMin, HcalParameters::etaTable, getEta(), HcalForward, HcalOuter, hpar, nEta, HcalParameters::noff, nR, HcalParameters::rTable, 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 1059 of file HcalDDDSimConstants.cc.

References cmsBatch::log, funct::tan(), and 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
	)

Definition at line 258 of file HcalDDDSimConstants.cc.

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

Referenced by 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 305 of file HcalDDDSimConstants.cc.

References PVValHelper::eta, HcalParameters::etaTable, hpar, alignCSCRings::r, rminHO, mathSSE::sqrt(), 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
    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 49 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 331 of file HcalDDDSimConstants.cc.

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

Referenced by getIdHF2QIE(), and 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 1091 of file HcalDDDSimConstants.cc.

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

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

inline

Definition at line 50 of file HcalDDDSimConstants.h.

References HcalParameters::gparHF, and hpar.

Referenced by HFFibre::initRun(), HFShower::initRun(), HFShowerLibrary::initRun(), and HCalSD::update().

{return hpar->gparHF;}

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

inline

Definition at line 51 of file HcalDDDSimConstants.h.

References getFrontLayer(), getLastLayer(), getLayer0Wt(), getLayerBack(), getLayerFront(), getLayerMax(), and idHF2QIE.

Referenced by HcalDDDRecConstants::getHFCellParameters().

{return idHF2QIE;}

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

inline

Definition at line 48 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 345 of file HcalDDDSimConstants.cc.

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

Referenced by getIdHF2QIE(), 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(), getIdHF2QIE(), 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 389 of file HcalDDDSimConstants.cc.

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

Referenced by getIdHF2QIE(), 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 369 of file HcalDDDSimConstants.cc.

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

Referenced by getIdHF2QIE(), 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+1) {
      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 404 of file HcalDDDSimConstants.cc.

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

Referenced by getIdHF2QIE(), 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 63 of file HcalDDDSimConstants.h.

References maxDepth.

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

{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 410 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 = layerGroup(eta-1,maxLayer_);
    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 67 of file HcalDDDSimConstants.h.

References depthMaxDf_, depthMaxSp_, getMinDepth(), getModHalfHBHE(), getPhiCons(), and getPhis().

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

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

Referenced by getMaxDepthDet(), and 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 466 of file HcalDDDSimConstants.cc.

References nmodHB, nmodHE, nzHB, and nzHE.

Referenced by getMaxDepthDet().

                                                                          {

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

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

Referenced by getMaxDepthDet(), 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 497 of file HcalDDDSimConstants.cc.

References getPhiCons(), createfilelist::int, gen::k, phi, phiNumber(), unitPhi(), and units().

Referenced by getMaxDepthDet(), and 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((CLHEP::twopi+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 " << phis[k].second/CLHEP::deg;
#endif
  return phis;
}

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

inline

Definition at line 74 of file HcalDDDSimConstants.h.

References hpar, and HcalParameters::phitable.

Referenced by HcalDDDRecConstants::getHFCellParameters(), and HFShowerLibrary::initRun().

{return hpar->phitable;}

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

inline

Definition at line 75 of file HcalDDDSimConstants.h.

References HcalCellTypes(), hpar, and HcalParameters::rTable.

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

{return hpar->rTable;}

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

private

Definition at line 1070 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 520 of file HcalDDDSimConstants.cc.

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

Referenced by getRTableHF(), 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 560 of file HcalDDDSimConstants.cc.

References cell(), egammaForCoreTracking_cff::depth, dlShort, dzVcal, PVValHelper::eta, HcalParameters::etaMax, HcalParameters::etaMin, muonCSCDigis_cfi::gain, getGain(), getPhis(), getShift(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, idHF2QIE, GetRecoTauVFromDQM_MC_cff::kk, checklumidiff::l, maxDepth, maxLayer_, maxLayerHB_, HcalParameters::noff, nzHB, nzHE, HcalCellType::HcalCell::ok, convertSQLiteXML::ok, phi, HcalParameters::phibin, HcalParameters::phioff, HcalParameters::phitable, HcalCellType::setPhi(), edm::shift, csvLumiCalc::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 = maxLayer_+1; 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 749 of file HcalDDDSimConstants.cc.

References depthEta16, depthEta29, depthMaxDf_, depthMaxSp_, depths, dlShort, runTauDisplay::dtype, HcalParameters::dzVcal, dzVcal, etaHO, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaTable, HcalLayerDepthMap::getDepthMax(), getLayerMax(), HcalParameters::gparHF, HcalForward, hpar, mps_fire::i, idHF2QIE, HcalLayerDepthMap::initialize(), isBH_, gen::k, GetRecoTauVFromDQM_MC_cff::kk, checklumidiff::l, layBHB, layBHE, HcalParameters::Layer0Wt, layerGroup(), layerGroupSize(), layFHB, layFHE, ldmap_, HcalParameters::maxDepth, maxDepth, maxLayerHB_, HcalParameters::modHB, HcalParameters::modHE, nEta, nmodHB, nmodHE, HcalParameters::noff, 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 << " 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;
  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 79 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 1235 of file HcalDDDSimConstants.cc.

References hpar, gen::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 1247 of file HcalDDDSimConstants.cc.

References egammaForCoreTracking_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 1223 of file HcalDDDSimConstants.cc.

References hpar, gen::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 80 of file HcalDDDSimConstants.h.

References ldmap_, maxHFDepth(), numberOfCells(), phiNumber(), printTiles(), csvLumiCalc::unit, and unitPhi().

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

References idHF2QIE, gen::k, maxDepth, and convertSQLiteXML::ok.

Referenced by getMaxDepth(), ldMap(), 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 674 of file HcalDDDSimConstants.cc.

References HcalCellTypes(), createfilelist::int, and pileupDistInMC::num.

Referenced by ldMap(), 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;
}

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

Definition at line 689 of file HcalDDDSimConstants.cc.

References phi.

Referenced by HcalDDDRecConstants::getHCID(), getPhis(), HcalDDDRecConstants::getPhis(), ldMap(), 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 1112 of file HcalDDDSimConstants.cc.

References jets_cff::area, funct::cos(), egammaForCoreTracking_cff::depth, HcalParameters::dxHB, HcalParameters::dyHB, PVValHelper::dz, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, gen::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHB, SiStripPI::mean, 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 " << thetaH/CLHEP::deg 
                   << ":" << thetaL/CLHEP::deg << " Layer " 
                   << layL << ":" << layH-1;
  for (int lay=layL; 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 1153 of file HcalDDDSimConstants.cc.

References jets_cff::area, HcalParameters::dx1HE, HcalParameters::dyHE, PVValHelper::eta, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, gen::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHE, SiStripPI::max, SiStripPI::mean, min(), 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*CLHEP::deg) nphi = 1;
  edm::LogVerbatim("HcalGeom") << "\ntileHE:: Eta/depth " << zside*eta << "|" 
                   << depth << " theta " << thetaH/CLHEP::deg 
                   << ":" << thetaL/CLHEP::deg << " Layer " 
                   << layL << ":" << layH-1 << " phi " << nphi;
  for (int lay=layL; 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.*CLHEP::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 698 of file HcalDDDSimConstants.cc.

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

Referenced by ldMap().

                                           {
 
  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 735 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(), ldMap(), 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 741 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

int HcalDDDSimConstants::depthEta16[2]

private

Definition at line 126 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::depthEta29[2]

private

Definition at line 127 of file HcalDDDSimConstants.h.

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

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

private

Definition at line 134 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

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

private

Definition at line 134 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

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

private

Definition at line 118 of file HcalDDDSimConstants.h.

Referenced by getLayerBack(), and initialize().

double HcalDDDSimConstants::dlShort

private

Definition at line 125 of file HcalDDDSimConstants.h.

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

double HcalDDDSimConstants::dzVcal

private

Definition at line 124 of file HcalDDDSimConstants.h.

Referenced by HcalCellTypes(), and initialize().

double HcalDDDSimConstants::etaHO[4]

private

Definition at line 122 of file HcalDDDSimConstants.h.

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

const HcalParameters* HcalDDDSimConstants::hpar

private

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

std::vector<HcalDetId> HcalDDDSimConstants::idHF2QIE

private

Definition at line 133 of file HcalDDDSimConstants.h.

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

bool HcalDDDSimConstants::isBH_

private

Definition at line 132 of file HcalDDDSimConstants.h.

Referenced by initialize(), and isBH().

int HcalDDDSimConstants::layBHB[3]

private

Definition at line 130 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

int HcalDDDSimConstants::layBHE[4]

private

Definition at line 131 of file HcalDDDSimConstants.h.

Referenced by getLastLayer(), and initialize().

int HcalDDDSimConstants::layFHB[2]

private

Definition at line 128 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

int HcalDDDSimConstants::layFHE[3]

private

Definition at line 129 of file HcalDDDSimConstants.h.

Referenced by getFrontLayer(), and initialize().

HcalLayerDepthMap HcalDDDSimConstants::ldmap_

private

Definition at line 112 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 114 of file HcalDDDSimConstants.h.

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

const int HcalDDDSimConstants::maxLayer_ =18

staticprivate

Definition at line 91 of file HcalDDDSimConstants.h.

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

const int HcalDDDSimConstants::maxLayerHB_ =16

staticprivate

Definition at line 92 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nDepth

private

Definition at line 119 of file HcalDDDSimConstants.h.

const int HcalDDDSimConstants::nDepthMax =9

staticprivate

Definition at line 90 of file HcalDDDSimConstants.h.

int HcalDDDSimConstants::nEta

private

Definition at line 115 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nmodHB

private

Definition at line 120 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

int HcalDDDSimConstants::nmodHE

private

Definition at line 121 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

int HcalDDDSimConstants::nPhiF

private

Definition at line 117 of file HcalDDDSimConstants.h.

Referenced by initialize().

int HcalDDDSimConstants::nR

private

Definition at line 116 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nzHB

private

Definition at line 120 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nzHE

private

Definition at line 121 of file HcalDDDSimConstants.h.

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

double HcalDDDSimConstants::rminHO

private

Definition at line 122 of file HcalDDDSimConstants.h.

Referenced by getEtaHO(), and initialize().

double HcalDDDSimConstants::zVcal

private

Definition at line 123 of file HcalDDDSimConstants.h.

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