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 (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, int zside, 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, int phi, int zside, int depth, int lay)
double	getEtaHO (double &etaR, double &x, double &y, double &z) const
const std::vector< double > &	getEtaTableHF () const
const std::vector< double > &	getGparHF () const
const std::vector< HcalDetId > &	getIdHF2QIE () const
std::pair< int, int >	getiEtaRange (const int i) 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, 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, 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 (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 (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, const 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 (const int eta, const int i) const
unsigned int	layerGroup (const int det, const int eta, const int phi, const int zside, const int i) const
unsigned int	layerGroupSize (const 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_
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 11 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, 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();
  std::cout << "HcalDDDSimConstants: " << cellTypes.size()
        << " cells of type HCal (All)\n";
#endif
}

HcalDDDSimConstants::~HcalDDDSimConstants

(

)

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

References gather_cfg::cout, deltaEta(), HcalCellType::HcalCell::deta, HcalCellType::HcalCell::dphi, HcalParameters::drHB, HcalCellType::HcalCell::drz, HcalParameters::dzHE, stringResolutionProvider_cfi::eta, HcalCellType::HcalCell::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
  std::cout << "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 << std::endl;
#endif
  return tmp;
}

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

private

Definition at line 905 of file HcalDDDSimConstants.cc.

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

Referenced by cell(), and ldMap().

                                            {

  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
  std::cout << "HcalDDDSimConstants::deltaEta " << etaR << " " 
        << depth << " ==> " << tmp << std::endl;
#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 123 of file HcalDDDSimConstants.cc.

References particleFlowClusterECALTimeSelected_cfi::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	(	int	layer,
		const std::vector< HcalParameters::LayerItem > &	layerGroup
	)		const

Definition at line 131 of file HcalDDDSimConstants.cc.

References mps_fire::i, and hcalTTPDigis_cfi::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 143 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.push_back(std::pair<double,double>(hpar->rHB[i],hpar->drHB[i]));
    }
  } else {
    for (unsigned int i=0; i<hpar->zHE.size(); ++i) {
      gcons.push_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 158 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::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
  std::cout << "getDepthEta16: " << det << ":" << depth << std::endl;
#endif
  return depth;
}

int HcalDDDSimConstants::getDepthEta16M

(

const int

det

)

const

Definition at line 169 of file HcalDDDSimConstants.cc.

References particleFlowClusterECALTimeSelected_cfi::depth, depthEta16, HcalLayerDepthMap::getDepth16(), ldmap_, HcalLayerDepthMap::validDet(), and HcalDigiParam_cfi::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,
		int	zside,
		int	i
	)		const

Definition at line 183 of file HcalDDDSimConstants.cc.

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

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

                                          {

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

References particleFlowClusterECALTimeSelected_cfi::depth, depthEta29, HcalLayerDepthMap::getMaxDepthLastHE(), ldmap_, HcalLayerDepthMap::validDet(), and HcalDigiParam_cfi::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 207 of file HcalDDDSimConstants.cc.

References depthEta16, stringResolutionProvider_cfi::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 232 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(), ldMap(), 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, const int  depth = 1  ) const

private

Definition at line 944 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, 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
  std::cout << "HcalDDDSimConstants::getEta " << etaR << " " 
        << zside << " " << depth << " ==> " << tmp << "\n";
#endif
  return tmp;
}

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

private

Definition at line 984 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, cmsBatch::log, funct::tan(), and tmp.

                                                                       {

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

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

Definition at line 256 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::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
  std::cout << "HcalDDDEsimConstants:getEtaDepth: I/P " << det << ":" << etaR
        << ":" << phi << ":" << zside << ":" << depth << ":" << lay << "\n";
#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
  std::cout << "HcalDDDEsimConstants:getEtaDepth: O/P " << etaR << ":" << depth
        << std::endl;
#endif
  return std::pair<int,int>(etaR,depth);
}

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

Definition at line 302 of file HcalDDDSimConstants.cc.

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

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

inline

Definition at line 51 of file HcalDDDSimConstants.h.

References HcalParameters::etaTableHF, and hpar.

{return hpar->etaTableHF;}

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

private

Definition at line 1016 of file HcalDDDSimConstants.cc.

References muonCSCDigis_cfi::gain, HcalParameters::HBGains, HcalEndcap, HcalForward, HcalOuter, HcalParameters::HEGains, HcalParameters::HFGains, and hpar.

Referenced by HcalCellTypes(), and ldMap().

                                           {

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

References HcalParameters::gparHF, and hpar.

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

{return hpar->gparHF;}

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

inline

Definition at line 53 of file HcalDDDSimConstants.h.

References 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 50 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]);}

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

Definition at line 328 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 356 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 336 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 371 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,
		bool	partialOnly
	)		const

Definition at line 377 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,
		bool	partialOnly
	)		const

Definition at line 405 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 433 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 442 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 464 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, plotBeamSpotDB::first, getPhiCons(), createfilelist::int, gen::k, phi, phiNumber(), edm::second(), 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.push_back(std::pair<int,double>(iphi,phi));
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getPhis: subdet|ieta|iphi " << subdet << "|" << ieta 
        << " with " << phis.size() << " phi bins" << std::endl;
  for (unsigned int k=0; k<phis.size(); ++k)
    std::cout << "[" << k << "] iphi " << phis[k].first << " phi "
          << phis[k].second/CLHEP::deg << std::endl;
#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 995 of file HcalDDDSimConstants.cc.

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

Referenced by HcalCellTypes(), and ldMap().

                                         {

  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 486 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	(	HcalSubdetector	subdet,
		int	ieta = -1,
		int	depth = -1
	)		const

Definition at line 526 of file HcalDDDSimConstants.cc.

References cell(), particleFlowClusterECALTimeSelected_cfi::depth, dlShort, dzVcal, stringResolutionProvider_cfi::eta, HcalParameters::etaMax, HcalParameters::etaMin, plotBeamSpotDB::first, muonCSCDigis_cfi::gain, getGain(), getPhis(), getShift(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hpar, idHF2QIE, gen::k, 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 HcalDigiParam_cfi::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.size() > 0) ? 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.push_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.push_back(temp2);
      // For HF look at extra cells
      if ((subdet == HcalForward) && (idHF2QIE.size() > 0)) {
        HcalCellType temp3(subdet, eta, zside+2, depth, temp1,
                   shift, gain, hsize);
        std::vector<int> phiMiss3;
        for (unsigned int k=0; k<phis.size(); ++k) {
          bool ok(false);
          for (unsigned int l=0; l<idHF2QIE.size(); ++l) {
        if ((eta*zside     == idHF2QIE[l].ieta()) && 
            (phis[k].first == idHF2QIE[l].iphi())) {
          ok = true;
          break;
        }
          }
          if (!ok) phiMiss3.push_back(phis[k].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.push_back(temp3);
      }
    }
      }
    }
  }
  return cellTypes;
}

void HcalDDDSimConstants::initialize ( void  )

private

Definition at line 710 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, 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, HcalParameters::Layer0Wt, layerGroup(), layerGroupSize(), 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, HcalDigiParam_cfi::zside, and zVcal.

Referenced by HcalDDDSimConstants(), and ldMap().

                                           {

  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) {
    std::cout << "HcalDDDSimConstants:Read LayerGroup" << i << ":";
    for (unsigned int k=0; k<layerGroupSize( i ); k++) 
      std::cout << " [" << k << "] = " << layerGroup( i, k );
    std::cout << std::endl;
  }
#endif

  // Geometry parameters for HF
  dlShort   = hpar->gparHF[0];
  zVcal     = hpar->gparHF[4];
  dzVcal    = hpar->dzVcal;
#ifdef EDM_ML_DEBUG
  std::cout << "HcalDDDSimConstants: dlShort " << dlShort << " zVcal " << zVcal
        << " and dzVcal " << dzVcal << std::endl;
#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)
    std::cout << "Detector Type [" << i << "] iEta " << hpar->etaMin[i] << ":" 
          << hpar->etaMax[i] << " MaxDepth " << maxDepth[i] << std::endl;
#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].push_back(ll);
    }

#ifdef EDM_ML_DEBUG
    std::cout << "Depth " << i << " with " << depths[i].size() << " etas:";
    for (int k=0; k<nEta-1; ++k) std::cout << " " << depths[i][k];
    std::cout << std::endl;
#endif
  }

  nzHB   = hpar->modHB[1];
  nmodHB = hpar->modHB[0];
  nzHE   = hpar->modHE[1];
  nmodHE = hpar->modHE[0];
#ifdef EDM_ML_DEBUG
  std::cout << "HcalDDDSimConstants:: " << nzHB << ":" << nmodHB
        << " barrel and " << nzHE << ":" << nmodHE
        << " endcap half-sectors" << std::endl;
#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
  std::cout << "HO Eta boundaries " << etaHO[0] << " " << etaHO[1]
        << " " << etaHO[2] << " " << etaHO[3] << std::endl;
  std::cout << "HO Parameters " << rminHO << " " << hpar->zHO.size();
  for (int i=0; i<4; ++i) std::cout << " eta[" << i << "] = " << etaHO[i];
  for (unsigned int i=0; i<hpar->zHO.size(); ++i) 
    std::cout << " zHO[" << i << "] = " << hpar->zHO[i];
  std::cout << std::endl;
#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
  std::cout << "isBH_ " << hpar->noff.size() << ":" << noffsize << ":" 
        << noffl << ":" << isBH_ << std::endl;
  std::cout << "Depth index at ieta = 16 for HB (max) " << depthEta16[0] 
        << " HE (min) " << depthEta16[1] << "; max depth for itea = 29 : ("
        << depthEta29[0] << ":" << depthEta29[1] << ")" << std::endl;
#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.push_back(HcalDetId(HcalForward,hpar->noff[kk+1],hpar->noff[kk+2],1));
      kk += 2;
    }
  }
#ifdef EDM_ML_DEBUG
  std::cout << idHF2QIE.size() << " detector channels having 2 QIE cards:" 
        << std::endl;
  for (unsigned int k=0; k<idHF2QIE.size(); ++k) 
    std::cout << " [" << k << "] " << idHF2QIE[k] << std::endl;
#endif

  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+5) {
      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 (dtype == 1 || dtype == 2) {
    if ((int)(hpar->noff.size()) >= (noffk+7+nphi+3*ndeps)) {
      std::vector<int> ifi, iet, ily, idp;
      for (int i=0; i<nphi; ++i) ifi.push_back(hpar->noff[noffk+7+i]);
      for (int i=0; i<ndeps;++i) {
        iet.push_back(hpar->noff[noffk+7+nphi+3*i]);
        ily.push_back(hpar->noff[noffk+7+nphi+3*i+1]);
        idp.push_back(hpar->noff[noffk+7+nphi+3*i+2]);
      }
#ifdef EDM_ML_DEBUG
      std::cout << "Initialize HcalLayerDepthMap for Detector " << dtype
            << " etaMax " << hpar->etaMax[dtype] << " with " << nphi 
            << " sectors";
      for (int i=0; i<nphi; ++i) std::cout << ":" << ifi[i];
      std::cout << " and " << ndeps << " depth sections" << std::endl;
      for (int i=0; i<ndeps;++i)
        std::cout << " [" << i << "] " << iet[i] << "  " << ily[i] << "  "
              << idp[i] << std::endl;
      std::cout << "Maximum depth for last HE Eta tower " << depthEta29[0]
            << ":" << ndp16 << ":" << ndp29 << " L0 Wt " 
            << hpar->Layer0Wt[dtype-1] << ":" << wt << std::endl;
#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));
    }
      }
    }
  }
  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
  std::cout << "Detector type and maximum depth for all RBX " 
        << depthMaxDf_.first << ":" << depthMaxDf_.second
        << " and for special RBX " << depthMaxSp_.first << ":" 
        << depthMaxSp_.second << std::endl;
#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 ( const int  eta, const int  i  ) const

private

Definition at line 1132 of file HcalDDDSimConstants.cc.

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

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

                                            {
  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 ( const int  det, const int  eta, const int  phi, const int  zside, const int  i  ) const

private

Definition at line 1145 of file HcalDDDSimConstants.cc.

References particleFlowClusterECALTimeSelected_cfi::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 ( const int  eta ) const

private

Definition at line 1120 of file HcalDDDSimConstants.cc.

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

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

                                                                    {
  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 deltaEta(), getEta(), getGain(), getShift(), initialize(), layerGroup(), layerGroupSize(), ldmap_, maxHFDepth(), numberOfCells(), phiNumber(), printTileHB(), printTileHE(), printTiles(), alignCSCRings::r, csvLumiCalc::unit, and unitPhi().

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

{return &ldmap_;}

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

Definition at line 625 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.size() > 0) {
    bool ok(false);
    for (unsigned int k=0; k<idHF2QIE.size(); ++k) {
      if ((eta == idHF2QIE[k].ieta()) && (iphi == idHF2QIE[k].iphi())) {
    ok = true; break;
      }
    }
    if (!ok) mxdepth = 2;
  }
  return mxdepth;
}

unsigned int HcalDDDSimConstants::numberOfCells

(

HcalSubdetector

subdet

)

const

Definition at line 639 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, HcalCellTypes(), mps_fire::i, createfilelist::int, ecaldqm::binning::nPhiBins, and pileupDistInMC::num.

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

                                                                                 {

  unsigned int num = 0;
  std::vector<HcalCellType> cellTypes = HcalCellTypes(subdet);
  for (unsigned int i=0; i<cellTypes.size(); i++) {
    num += (unsigned int)(cellTypes[i].nPhiBins());
  }
#ifdef EDM_ML_DEBUG
  std::cout << "HcalDDDSimConstants:numberOfCells " 
        << cellTypes.size()  << " " << num 
        << " for subdetector " << subdet << std::endl;
#endif
  return num;
}

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

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

References funct::cos(), gather_cfg::cout, HcalParameters::dxHB, HcalParameters::dyHB, allConversions_cfi::dz, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, gen::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHB, min(), HcalParameters::rhoxHB, and funct::sin().

Referenced by ldMap(), and printTiles().

                                                              {
  std::cout << "HcalDDDSimConstants::printTileHB for eta " << eta << " and depth " << depth << "\n";
  
  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);
  std::cout << "\ntileHB:: eta|depth " << zside*eta << "|" << depth << " theta " << thetaH/CLHEP::deg << ":" << thetaL/CLHEP::deg << " Layer " << layL << ":" << layH-1 << "\n";
  for (int lay=layL; lay<layH; ++lay) {
    std::vector<double> area(2,0);
    int kk=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];
      kk++;
    }
      }
    }
    if (area[0] > 0) std::cout << std::setw(2) << lay << " Area " << std::setw(8) << area[0] << " " << std::setw(8) << area[1] << "\n";
  }
}

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

private

Definition at line 1066 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, HcalParameters::dx1HE, HcalParameters::dyHE, stringResolutionProvider_cfi::eta, HcalParameters::etaTable, JetChargeProducer_cfi::exp, getLayerBack(), getLayerFront(), hpar, gen::k, GetRecoTauVFromDQM_MC_cff::kk, HcalParameters::layHE, hpstanc_transforms::max, min(), HcalParameters::phibin, HcalParameters::rhoxHE, funct::tan(), and HcalParameters::zxHE.

Referenced by ldMap(), and 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;
  std::cout << "\ntileHE:: Eta/depth " << zside*eta << "|" << depth << " theta " << thetaH/CLHEP::deg << ":" << thetaL/CLHEP::deg << " Layer " << layL << ":" << layH-1 << " phi " << nphi << "\n";
  for (int lay=layL; lay<layH; ++lay) {
    std::vector<double> area(4,0);
    int kk=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]);
      } 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;
      }
      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) {
    std::cout << std::setw(2) << lay0 << " Area " << std::setw(8) << area[0] << " " << std::setw(8) << area[1] << "\n";
      } else {
    std::cout << 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] << "\n";
      }
    }
  }
}

void HcalDDDSimConstants::printTiles

(

)

const

Definition at line 663 of file HcalDDDSimConstants.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, stringResolutionProvider_cfi::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, HcalParameters::etaMin, getDepthEta16(), getMaxDepth(), hpar, ldmap_, printTileHB(), printTileHE(), HcalLayerDepthMap::validDet(), and HcalDigiParam_cfi::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;
  std::cout << "Tile Information for HB from " << hpar->etaMin[0] << " to " << hpar->etaMax[0] << "\n\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);
    }
  }

  std::cout << "\nTile Information for HE from " << hpar->etaMin[1] << " to " << hpar->etaMax[1] << "\n\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 696 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 702 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 128 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::depthEta29[2]

private

Definition at line 129 of file HcalDDDSimConstants.h.

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

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

private

Definition at line 132 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

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

private

Definition at line 132 of file HcalDDDSimConstants.h.

Referenced by getMaxDepthDet(), and initialize().

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

private

Definition at line 120 of file HcalDDDSimConstants.h.

Referenced by getLayerBack(), and initialize().

double HcalDDDSimConstants::dlShort

private

Definition at line 127 of file HcalDDDSimConstants.h.

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

double HcalDDDSimConstants::dzVcal

private

Definition at line 126 of file HcalDDDSimConstants.h.

Referenced by HcalCellTypes(), and initialize().

double HcalDDDSimConstants::etaHO[4]

private

Definition at line 124 of file HcalDDDSimConstants.h.

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

const HcalParameters* HcalDDDSimConstants::hpar

private

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

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

bool HcalDDDSimConstants::isBH_

private

Definition at line 130 of file HcalDDDSimConstants.h.

Referenced by initialize(), and isBH().

HcalLayerDepthMap HcalDDDSimConstants::ldmap_

private

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

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

const int HcalDDDSimConstants::maxLayer_ =18

staticprivate

Definition at line 113 of file HcalDDDSimConstants.h.

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

const int HcalDDDSimConstants::maxLayerHB_ =16

staticprivate

Definition at line 114 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nDepth

private

Definition at line 121 of file HcalDDDSimConstants.h.

const int HcalDDDSimConstants::nDepthMax =9

staticprivate

Definition at line 112 of file HcalDDDSimConstants.h.

int HcalDDDSimConstants::nEta

private

Definition at line 117 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nmodHB

private

Definition at line 122 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

int HcalDDDSimConstants::nmodHE

private

Definition at line 123 of file HcalDDDSimConstants.h.

Referenced by getModHalfHBHE(), and initialize().

int HcalDDDSimConstants::nPhiF

private

Definition at line 119 of file HcalDDDSimConstants.h.

Referenced by initialize().

int HcalDDDSimConstants::nR

private

Definition at line 118 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nzHB

private

Definition at line 122 of file HcalDDDSimConstants.h.

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

int HcalDDDSimConstants::nzHE

private

Definition at line 123 of file HcalDDDSimConstants.h.

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

double HcalDDDSimConstants::rminHO

private

Definition at line 124 of file HcalDDDSimConstants.h.

Referenced by getEtaHO(), and initialize().

double HcalDDDSimConstants::zVcal

private

Definition at line 125 of file HcalDDDSimConstants.h.

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