HcalDDDRecConstants Class Reference

#include <HcalDDDRecConstants.h>

Classes
struct	HcalActiveLength
struct	HcalEtaBin
struct	HcalID
struct	HFCellParameters

Public Member Functions
const HcalDDDSimConstants *	dddConstants () const
std::vector< std::pair< double, double > >	getConstHBHE (const int &type) const
std::vector< int >	getDepth (const int &det, const int &phi, const int &zside, const unsigned int &eta) const
std::vector< int >	getDepth (const unsigned int &eta, const bool &extra) const
int	getDepthEta16 (const int &det, const int &iphi, const int &zside) const
int	getDepthEta29 (const int &iphi, const int &zside, const int &type) const
std::vector< HcalEtaBin >	getEtaBins (const int &itype) const
std::pair< double, double >	getEtaLimit (const int &i) const
std::pair< double, double >	getEtaPhi (const int &subdet, const int &ieta, const int &iphi) const
std::pair< int, int >	getEtaRange (const int &i) const
const std::vector< double > &	getEtaTable () const
const std::vector< double > &	getEtaTableHF () const
HcalID	getHCID (int subdet, int ieta, int iphi, int lay, int idepth) const
std::vector< HFCellParameters >	getHFCellParameters () const
double	getLayer0Wt (const int &det, const int &phi, const int &zside) const
int	getLayerBack (const int &det, const int &eta, const int &phi, const int &depth) const
void	getLayerDepth (const int &ieta, std::map< int, int > &layers) const
int	getLayerFront (const int &det, const int &eta, const int &phi, const int &depth) const
int	getMaxDepth (const int &type) const
int	getMaxDepth (const int &itype, const int &ieta, const int &iphi, const int &zside) const
int	getMinDepth (const int &itype, const int &ieta, const int &iphi, const int &zside) const
int	getNEta () const
int	getNoff (const int &i) const
int	getNPhi (const int &type) const
double	getPhiBin (const int &i) const
double	getPhiOff (const int &i) const
const std::vector< double > &	getPhiOffs () const
std::vector< std::pair< int, double > >	getPhis (const int &subdet, const int &ieta) const
const std::vector< double > &	getPhiTable () const
const std::vector< double > &	getPhiTableHF () const
int	getPhiZOne (std::vector< std::pair< int, int >> &phiz) const
double	getRZ (const int &subdet, const int &ieta, const int &depth) const
double	getRZ (const int &subdet, const int &ieta, const int &iphi, const int &depth) const
double	getRZ (const int &subdet, const int &layer) const
std::pair< double, double >	getRZ (const HcalDetId &id) const
std::vector< HcalActiveLength >	getThickActive (const int &type) const
int	getTopoMode () const
int	getTriggerMode () const
std::vector< HcalCellType >	HcalCellTypes (HcalSubdetector) const
	HcalDDDRecConstants (const HcalParameters *hp, const HcalDDDSimConstants &hc)
HcalDetId	idBack (const HcalDetId &id) const
HcalDetId	idFront (const HcalDetId &id) const
bool	isBH () const
bool	isPlan1 (const HcalDetId &id) const
bool	isPlan1MergedId (const HcalDetId &id) const
bool	isPlan1ToBeMergedId (const HcalDetId &id) const
int	maxHFDepth (int ieta, int iphi) const
HcalDetId	mergedDepthDetId (const HcalDetId &id) const
bool	mergedDepthList29 (int ieta, int iphi, int depth) const
std::vector< int >	mergedDepthList29 (int ieta, int iphi) const
unsigned int	nCells (HcalSubdetector) const
unsigned int	nCells () const
unsigned int	numberOfCells (HcalSubdetector) const
void	specialRBXHBHE (const std::vector< HcalDetId > &, std::vector< HcalDetId > &) const
bool	specialRBXHBHE (bool flag, std::vector< HcalDetId > &) const
void	unmergeDepthDetId (const HcalDetId &id, std::vector< HcalDetId > &ids) const
bool	withSpecialRBXHBHE () const
	~HcalDDDRecConstants ()

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

Private Attributes
std::pair< int, int >	depthMaxDf_
std::pair< int, int >	depthMaxSp_
std::map< HcalDetId, HcalDetId >	detIdSp_
std::map< HcalDetId, std::vector< HcalDetId > >	detIdSpR_
std::vector< std::pair< int, int > >	etaSimValu
std::vector< double >	etaTable
std::vector< std::pair< double, double > >	gconsHB
std::vector< std::pair< double, double > >	gconsHE
const HcalDDDSimConstants &	hcons
const HcalParameters *	hpar
std::vector< int >	ietaMap
std::vector< int >	iEtaMax
std::vector< int >	iEtaMin
std::vector< int >	maxDepth
int	nHalves [2]
int	nModule [2]
std::vector< int >	nPhiBins
std::vector< double >	phibin
std::vector< int >	phiUnitS

Static Private Attributes
static const int	maxLayer_ = 18
static const int	maxLayerHB_ = 16

Detailed Description

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

Author

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

Definition at line 23 of file HcalDDDRecConstants.h.

Constructor & Destructor Documentation

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

Definition at line 16 of file HcalDDDRecConstants.cc.

References initialize().

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

HcalDDDRecConstants::~HcalDDDRecConstants

(

)

Definition at line 24 of file HcalDDDRecConstants.cc.

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

Member Function Documentation

const HcalDDDSimConstants* HcalDDDRecConstants::dddConstants ( ) const

inline

Definition at line 126 of file HcalDDDRecConstants.h.

References hcons.

{ return &hcons; }

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

inline

Definition at line 54 of file HcalDDDRecConstants.h.

References gconsHB, and gconsHE.

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

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

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

Definition at line 53 of file HcalDDDRecConstants.cc.

References HcalLayerDepthMap::getLayerDepth(), hcons, LayerTriplets::layers(), and HcalDDDSimConstants::ldMap().

Referenced by HcalTopology::HcalTopology().

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

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

Definition at line 30 of file HcalDDDRecConstants.cc.

References HcalLayerDepthMap::getLayerDepth(), hcons, hpar, dqmdumpme::last, HcalParameters::layerGroupEtaRec, LayerTriplets::layers(), and HcalDDDSimConstants::ldMap().

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

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

inline

Definition at line 66 of file HcalDDDRecConstants.h.

References HcalDDDSimConstants::getDepthEta16(), and hcons.

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

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

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

inline

Definition at line 69 of file HcalDDDRecConstants.h.

References HcalDDDSimConstants::getDepthEta29(), and hcons.

Referenced by HcalTopology::decrementDepth().

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

std::vector< HcalDDDRecConstants::HcalEtaBin > HcalDDDRecConstants::getEtaBins

(

const int &

itype

)

const

Definition at line 69 of file HcalDDDRecConstants.cc.

References spr::find(), first, HcalLayerDepthMap::getLayerDepth(), getLayerDepth(), getOneEtaBin(), getPhis(), hcons, mps_fire::i, iEtaMax, iEtaMin, isotrackApplyRegressor::k, LayerTriplets::layers(), HcalDDDSimConstants::ldMap(), phi, HcalLayerDepthMap::validDet(), and ecaldqm::zside().

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

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

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

inline

Definition at line 77 of file HcalDDDRecConstants.h.

References etaTable, and mps_fire::i.

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

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

Definition at line 133 of file HcalDDDRecConstants.cc.

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

Referenced by HcalGeomCheck::analyze().

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

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

inline

Definition at line 74 of file HcalDDDRecConstants.h.

References mps_fire::i, iEtaMax, and iEtaMin.

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

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

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

inline

Definition at line 75 of file HcalDDDRecConstants.h.

References etaTable.

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

{ return etaTable; }

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

inline

Definition at line 76 of file HcalDDDRecConstants.h.

References HcalParameters::etaTableHF, and hpar.

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

{ return hpar->etaTableHF; }

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

Definition at line 162 of file HcalDDDRecConstants.cc.

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

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

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

std::vector< HcalDDDRecConstants::HFCellParameters > HcalDDDRecConstants::getHFCellParameters

(

)

const

Definition at line 211 of file HcalDDDRecConstants.cc.

References funct::abs(), cells, HcalDDDSimConstants::getIdHF2QIE(), HcalDDDSimConstants::getPhiTableHF(), HcalDDDSimConstants::getRTableHF(), hcons, iEtaMin, isotrackApplyRegressor::k, maxDepth, HLT_FULL_cff::nEta, and nphi.

Referenced by HcalFlexiHardcodeGeometryLoader::makeHFCells().

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

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

inline

Definition at line 85 of file HcalDDDRecConstants.h.

References HcalDDDSimConstants::getLayer0Wt(), and hcons.

Referenced by HcalHitRelabeller::energyWt().

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

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

Definition at line 278 of file HcalDDDRecConstants.cc.

References PVValHelper::eta, Puppi_cff::etaMax, HcalDDDSimConstants::getLastLayer(), HcalLayerDepthMap::getLayerBack(), hcons, isotrackApplyRegressor::k, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), and ecaldqm::zside().

Referenced by getRZ().

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

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

Definition at line 265 of file HcalDDDRecConstants.cc.

References cmsLHEtoEOSManager::l, layerGroup(), and layerGroupSize().

Referenced by getEtaBins().

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

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

Definition at line 301 of file HcalDDDRecConstants.cc.

References PVValHelper::eta, Puppi_cff::etaMax, HcalDDDSimConstants::getFrontLayer(), HcalLayerDepthMap::getLayerFront(), HcalEndcap, hcons, isotrackApplyRegressor::k, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), and ecaldqm::zside().

Referenced by getRZ(), and getThickActive().

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

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

inline

Definition at line 88 of file HcalDDDRecConstants.h.

References maxDepth.

Referenced by HcalSimHitCheck::beginRun(), HcalSimHitsClient::beginRun(), HcalRecHitsDQMClient::beginRun(), HcalSimHitStudy::bookHistograms(), SimHitsValidationHcal::bookHistograms(), HcalSimHitsValidation::bookHistograms(), HcalTopology::depthBinInformation(), HcalDigisValidation::dqmBeginRun(), HcalRecHitsAnalyzer::dqmBeginRun(), HcalTopology::HcalTopology(), mergedDepthList29(), HcalTopologyIdealEP::produce(), and HcalTopology::validRaw().

{ return maxDepth[type]; }

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

Definition at line 331 of file HcalDDDRecConstants.cc.

References HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29M(), HcalDDDSimConstants::getMaxDepth(), hcons, hpar, iEtaMax, layerGroup(), layerGroupSize(), maxLayer_, maxLayerHB_, and HcalParameters::noff.

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

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

Definition at line 354 of file HcalDDDRecConstants.cc.

References HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getMinDepth(), hcons, iEtaMin, layerGroup(), layerGroupSize(), and maxDepth.

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

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

int HcalDDDRecConstants::getNEta ( ) const

inline

Definition at line 91 of file HcalDDDRecConstants.h.

References HcalParameters::etagroup, and hpar.

Referenced by HcalTopology::HcalTopology().

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

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

inline

Definition at line 92 of file HcalDDDRecConstants.h.

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

Referenced by HcalTopology::validRaw().

{ return hpar->noff[i]; }

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

inline

Definition at line 93 of file HcalDDDRecConstants.h.

References nPhiBins.

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

{ return nPhiBins[type]; }

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

private

Definition at line 708 of file HcalDDDRecConstants.cc.

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

Referenced by getEtaBins().

                                                                                             {
  unsigned int lymax = (subdet == HcalBarrel) ? maxLayerHB_ + 1 : maxLayer_ + 1;
  int type = (subdet == HcalBarrel) ? 0 : 1;
  double dphi = phibin[ieta - 1];
  HcalDDDRecConstants::HcalEtaBin etabin =
      HcalDDDRecConstants::HcalEtaBin(ieta, zside, dphi, etaTable[ieta - 1], etaTable[ieta]);
  etabin.phis.insert(etabin.phis.end(), phis.begin(), phis.end());
  int n = (ieta == iEtaMax[type]) ? 0 : 1;
  HcalDDDRecConstants::HcalEtaBin etabin0 =
      HcalDDDRecConstants::HcalEtaBin(ieta, zside, dphi, etaTable[ieta - 1], etaTable[ieta + n]);
  etabin0.depthStart = hcons.getDepthEta29(phis[0].first, zside, 0) + 1;
  int dstart = -1;
  int lmin(0), lmax(0);

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

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

inline

Definition at line 94 of file HcalDDDRecConstants.h.

References mps_fire::i, and phibin.

{ return phibin[i]; }

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

inline

Definition at line 95 of file HcalDDDRecConstants.h.

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

{ return hpar->phioff[i]; }

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

inline

Definition at line 96 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::phioff.

Referenced by HcalTopology::HcalTopology().

{ return hpar->phioff; }

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

Definition at line 374 of file HcalDDDRecConstants.cc.

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

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

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

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

inline

Definition at line 98 of file HcalDDDRecConstants.h.

References phibin.

Referenced by HcalTopology::HcalTopology().

{ return phibin; }

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

inline

Definition at line 99 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::phitable.

Referenced by HcalTopology::HcalTopology().

{ return hpar->phitable; }

int HcalDDDRecConstants::getPhiZOne

(

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

phiz

)

const

Definition at line 398 of file HcalDDDRecConstants.cc.

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

Referenced by HcalTopology::getPhiZOne().

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

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

Definition at line 417 of file HcalDDDRecConstants.cc.

Referenced by HcalGeomCheck::analyze(), HcalHBHEMuonAnalyzer::goodCell(), HcalHBHEMuonHighEtaAnalyzer::goodCell(), and AlCaHcalHBHEMuonProducer::goodCell().

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

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

Definition at line 421 of file HcalDDDRecConstants.cc.

References first, gconsHB, gconsHE, getLayerFront(), and HcalBarrel.

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

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

Definition at line 432 of file HcalDDDRecConstants.cc.

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

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

std::pair< double, double > HcalDDDRecConstants::getRZ

(

const HcalDetId &

id

)

const

Definition at line 442 of file HcalDDDRecConstants.cc.

References HLT_FULL_cff::depth, first, gconsHB, gconsHE, getLayerBack(), getLayerFront(), HcalBarrel, edm::second(), and ecaldqm::zside().

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

std::vector< HcalDDDRecConstants::HcalActiveLength > HcalDDDRecConstants::getThickActive

(

const int &

type

)

const

Definition at line 466 of file HcalDDDRecConstants.cc.

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

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

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

int HcalDDDRecConstants::getTopoMode ( ) const

inline

Definition at line 106 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::topologyMode.

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

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

int HcalDDDRecConstants::getTriggerMode ( ) const

inline

Definition at line 107 of file HcalDDDRecConstants.h.

References hpar, and HcalParameters::topologyMode.

Referenced by HcalTopology::HcalTopology().

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

std::vector< HcalCellType > HcalDDDRecConstants::HcalCellTypes

(

HcalSubdetector

subdet

)

const

Definition at line 521 of file HcalDDDRecConstants.cc.

References cells, submitPVResolutionJobs::count, HLT_FULL_cff::depth, defaults_cfi::etabins, etaSimValu, first, getEtaBins(), HcalBarrel, HcalDDDSimConstants::HcalCellTypes(), HcalEndcap, hcons, hpar, iEtaMax, HcalParameters::phioff, groupFilesInBlocks::temp, unit(), and HcalDDDSimConstants::unitPhi().

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

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

HcalDetId HcalDDDRecConstants::idBack

(

const HcalDetId &

id

)

const

Definition at line 664 of file HcalDDDRecConstants.cc.

References detIdSpR_.

Referenced by HcalTopology::idBack().

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

HcalDetId HcalDDDRecConstants::idFront

(

const HcalDetId &

id

)

const

Definition at line 656 of file HcalDDDRecConstants.cc.

References HLT_FULL_cff::depth, and detIdSpR_.

Referenced by HcalTopology::idFront().

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

void HcalDDDRecConstants::initialize ( void  )

private

Definition at line 802 of file HcalDDDRecConstants.cc.

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

Referenced by HcalDDDRecConstants().

                                         {
  //Eta grouping
  int nEta = (int)(hpar->etagroup.size());
  if (nEta != (int)(hpar->phigroup.size())) {
    edm::LogError("HCalGeom") << "HcalDDDRecConstants: sizes of the vectors "
                              << " etaGroup (" << nEta << ") and phiGroup (" << hpar->phigroup.size()
                              << ") do not match";
    throw cms::Exception("DDException") << "HcalDDDRecConstants: inconsistent array sizes" << nEta << ":"
                                        << hpar->phigroup.size();
  }

  // First eta table
  iEtaMin = hpar->etaMin;
  iEtaMax = hpar->etaMax;
  etaTable.clear();
  ietaMap.clear();
  etaSimValu.clear();
  int ieta(0), ietaHB(0), ietaHE(0), ietaHEM(0);
  etaTable.emplace_back(hpar->etaTable[ieta]);
  for (int i = 0; i < nEta; ++i) {
    int ef = ieta + 1;
    ieta += (hpar->etagroup[i]);
    if (ieta >= (int)(hpar->etaTable.size())) {
      edm::LogError("HCalGeom") << "HcalDDDRecConstants: Going beyond the array boundary " << hpar->etaTable.size()
                                << " at index " << i << " of etaTable from SimConstant";
      throw cms::Exception("DDException")
          << "HcalDDDRecConstants: Going beyond the array boundary " << hpar->etaTable.size() << " at index " << i
          << " of etaTable from SimConstant";
    } else {
      etaTable.emplace_back(hpar->etaTable[ieta]);
      etaSimValu.emplace_back(std::pair<int, int>(ef, ieta));
    }
    for (int k = 0; k < (hpar->etagroup[i]); ++k)
      ietaMap.emplace_back(i + 1);
    if (ieta <= hpar->etaMax[0])
      ietaHB = i + 1;
    if (ieta <= hpar->etaMin[1])
      ietaHE = i + 1;
    if (ieta <= hpar->etaMax[1])
      ietaHEM = i + 1;
  }
  iEtaMin[1] = ietaHE;
  iEtaMax[0] = ietaHB;
  iEtaMax[1] = ietaHEM;

  // Then Phi bins
  nPhiBins.clear();
  for (unsigned int k = 0; k < 4; ++k)
    nPhiBins.emplace_back(0);
  ieta = 0;
  phibin.clear();
  phiUnitS.clear();
  for (int i = 0; i < nEta; ++i) {
    double dphi = (hpar->phigroup[i]) * (hpar->phibin[ieta]);
    phibin.emplace_back(dphi);
    int nphi = (int)((2._pi + 0.001) / dphi);
    if (ieta <= iEtaMax[0]) {
      if (nphi > nPhiBins[0])
        nPhiBins[3] = nPhiBins[0] = nphi;
    }
    if (ieta >= iEtaMin[1]) {
      if (nphi > nPhiBins[1])
        nPhiBins[1] = nphi;
    }
    ieta += (hpar->etagroup[i]);
  }
  for (unsigned int i = 1; i < hpar->etaTable.size(); ++i) {
    int unit = hcons.unitPhi(hpar->phibin[i - 1]);
    phiUnitS.emplace_back(unit);
  }
  for (double i : hpar->phitable) {
    int nphi = (int)((2._pi + 0.001) / i);
    if (nphi > nPhiBins[2])
      nPhiBins[2] = nphi;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants: Modified eta/deltaphi table for " << nEta << " bins";
  for (int i = 0; i < nEta; ++i)
    edm::LogVerbatim("HCalGeom") << "Eta[" << i << "] = " << etaTable[i] << ":" << etaTable[i + 1] << ":"
                                 << etaSimValu[i].first << ":" << etaSimValu[i].second << " PhiBin[" << i
                                 << "] = " << convertRadToDeg(phibin[i]);
  for (unsigned int i = 0; i < phiUnitS.size(); ++i)
    edm::LogVerbatim("HCalGeom") << " PhiUnitS[" << i << "] = " << phiUnitS[i];
  for (unsigned int i = 0; i < nPhiBins.size(); ++i)
    edm::LogVerbatim("HCalGeom") << " nPhiBins[" << i << "] = " << nPhiBins[i];
  for (unsigned int i = 0; i < hpar->etaTableHF.size(); ++i)
    edm::LogVerbatim("HCalGeom") << " EtaTableHF[" << i << "] = " << hpar->etaTableHF[i];
  for (unsigned int i = 0; i < hpar->phitable.size(); ++i)
    edm::LogVerbatim("HCalGeom") << " PhiBinHF[" << i << "] = " << hpar->phitable[i];
#endif

  //Now the depths
  maxDepth = hpar->maxDepth;
  maxDepth[0] = maxDepth[1] = 0;
  for (int i = 0; i < nEta; ++i) {
    unsigned int imx = layerGroupSize(i);
    int laymax = (imx > 0) ? layerGroup(i, imx - 1) : 0;
    if (i < iEtaMax[0]) {
      int laymax0 = (imx > 16) ? layerGroup(i, 16) : laymax;
      if (i + 1 == iEtaMax[0])
        laymax0 = hcons.getDepthEta16M(1);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:HB " << i << " " << imx << " " << laymax << " " << laymax0;
#endif
      if (maxDepth[0] < laymax0)
        maxDepth[0] = laymax0;
    }
    if (i >= iEtaMin[1] - 1 && i < iEtaMax[1]) {
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:HE " << i << " " << imx << " " << laymax;
#endif
      if (maxDepth[1] < laymax)
        maxDepth[1] = laymax;
    }
  }
#ifdef EDM_ML_DEBUG
  for (int i = 0; i < 4; ++i)
    edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:Detector Type[" << i << "] iEta " << iEtaMin[i] << ":"
                                 << iEtaMax[i] << " MaxDepth " << maxDepth[i];
#endif

  //Now the geometry constants
  nModule[0] = hpar->modHB[0];
  nHalves[0] = hpar->modHB[1];
  for (unsigned int i = 0; i < hpar->rHB.size(); ++i) {
    gconsHB.emplace_back(std::pair<double, double>(hpar->rHB[i] / CLHEP::cm, hpar->drHB[i] / CLHEP::cm));
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:HB with " << nModule[0] << " modules and " << nHalves[0]
                               << " halves and " << gconsHB.size() << " layers";
  for (unsigned int i = 0; i < gconsHB.size(); ++i)
    edm::LogVerbatim("HCalGeom") << "rHB[" << i << "] = " << gconsHB[i].first << " +- " << gconsHB[i].second;
#endif
  nModule[1] = hpar->modHE[0];
  nHalves[1] = hpar->modHE[1];
  for (unsigned int i = 0; i < hpar->zHE.size(); ++i) {
    gconsHE.emplace_back(std::pair<double, double>(hpar->zHE[i] / CLHEP::cm, hpar->dzHE[i] / CLHEP::cm));
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:HE with " << nModule[1] << " modules and " << nHalves[1]
                               << " halves and " << gconsHE.size() << " layers";
  for (unsigned int i = 0; i < gconsHE.size(); ++i)
    edm::LogVerbatim("HCalGeom") << "zHE[" << i << "] = " << gconsHE[i].first << " +- " << gconsHE[i].second;
#endif

  //Special RBX
  depthMaxSp_ = hcons.getMaxDepthDet(0);
  if (depthMaxSp_.first == 0) {
    depthMaxSp_ = depthMaxDf_ = std::pair<int, int>(2, maxDepth[1]);
  } else if (depthMaxSp_.first == 1) {
    depthMaxDf_ = std::pair<int, int>(1, maxDepth[0]);
    if (depthMaxSp_.second > maxDepth[0])
      maxDepth[0] = depthMaxSp_.second;
  } else {
    depthMaxDf_ = std::pair<int, int>(2, maxDepth[1]);
    if (depthMaxSp_.second > maxDepth[1])
      maxDepth[1] = depthMaxSp_.second;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:Detector type and maximum depth for all RBX "
                               << depthMaxDf_.first << ":" << depthMaxDf_.second << " and for special RBX "
                               << depthMaxSp_.first << ":" << depthMaxSp_.second;
#endif

  //Map of special DetId's
  std::vector<int> phis;
  HcalSubdetector subdet = HcalSubdetector(hcons.ldMap()->validDet(phis));
  detIdSp_.clear();
  detIdSpR_.clear();
  if ((subdet == HcalBarrel) || (subdet == HcalEndcap)) {
    int phi = (phis[0] > 0) ? phis[0] : -phis[0];
    int zside = (phis[0] > 0) ? 1 : -1;
    int lymax = (subdet == HcalBarrel) ? maxLayerHB_ + 1 : maxLayer_ + 1;
    std::pair<int, int> etas = hcons.ldMap()->validEta();
    for (int eta = etas.first; eta <= etas.second; ++eta) {
      std::map<int, std::pair<int, int>> oldDep;
      int depth(0);
      int lmin = layerGroup(eta - 1, 0);
      for (int lay = 0; lay < lymax; ++lay) {
        int depc = layerGroup(eta - 1, lay);
        if (depth != depc) {
          if (depth != 0)
            oldDep[depth] = std::pair<int, int>(lmin, lay - 1);
          depth = depc;
          lmin = lay;
        }
      }
      if (depth != 0)
        oldDep[depth] = std::pair<int, int>(lmin, lymax - 1);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:Eta|Phi|Zside " << eta << ":" << phi << ":" << zside
                                   << " with " << oldDep.size() << " old Depths";
      unsigned int kk(0);
      for (std::map<int, std::pair<int, int>>::const_iterator itr = oldDep.begin(); itr != oldDep.end(); ++itr, ++kk)
        edm::LogVerbatim("HCalGeom") << "[" << kk << "] " << itr->first << " --> " << itr->second.first << ":"
                                     << itr->second.second;
#endif
      std::pair<int, int> depths = hcons.ldMap()->getDepths(eta);
      for (int ndepth = depths.first; ndepth <= depths.second; ++ndepth) {
        bool flag = ((subdet == HcalBarrel && eta == iEtaMax[0] && ndepth > hcons.getDepthEta16(subdet, phi, zside)) ||
                     (subdet == HcalEndcap && eta == iEtaMin[1] && ndepth < hcons.getDepthEta16(subdet, phi, zside)));
        if (!flag) {
          std::vector<int> count(oldDep.size(), 0);
          int layFront = hcons.ldMap()->getLayerFront(subdet, eta, phi, zside, ndepth);
          int layBack = hcons.ldMap()->getLayerBack(subdet, eta, phi, zside, ndepth);
          for (int lay = layFront; lay <= layBack; ++lay) {
            unsigned int l(0);
            for (std::map<int, std::pair<int, int>>::iterator itr = oldDep.begin(); itr != oldDep.end(); ++itr, ++l) {
              if (lay >= (itr->second).first && lay <= (itr->second).second) {
                ++count[l];
                break;
              }
            }
          }
          int odepth(0), maxlay(0);
          unsigned int l(0);
          for (std::map<int, std::pair<int, int>>::iterator itr = oldDep.begin(); itr != oldDep.end(); ++itr, ++l) {
            if (count[l] > maxlay) {
              odepth = itr->first;
              maxlay = count[l];
            }
          }
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:New Depth " << ndepth << " old Depth " << odepth
                                       << " max " << maxlay;
#endif
          for (int k : phis) {
            zside = (k > 0) ? 1 : -1;
            phi = (k > 0) ? k : -k;
            if (subdet == HcalEndcap && eta == hpar->noff[1] && ndepth > hcons.getDepthEta29M(0, true))
              break;
            HcalDetId newId(subdet, zside * eta, phi, ndepth);
            HcalDetId oldId(subdet, zside * eta, phi, odepth);
            detIdSp_[newId] = oldId;
            std::vector<HcalDetId> ids;
            std::map<HcalDetId, std::vector<HcalDetId>>::iterator itr = detIdSpR_.find(oldId);
            if (itr != detIdSpR_.end())
              ids = itr->second;
            ids.emplace_back(newId);
            detIdSpR_[oldId] = ids;
          }
        }
      }
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:Map for merging new channels to old channel"
                                 << " IDs with " << detIdSp_.size() << " entries";
    int l(0);
    for (auto itr : detIdSp_) {
      edm::LogVerbatim("HCalGeom") << "[" << l << "] Special " << itr.first << " Standard " << itr.second;
      ++l;
    }
    edm::LogVerbatim("HCalGeom") << "HcalDDDRecConstants:Reverse Map for mapping old to new IDs with "
                                 << detIdSpR_.size() << " entries";
    l = 0;
    for (auto itr : detIdSpR_) {
      edm::LogVerbatim("HCalGeom") << "[" << l << "] Standard " << itr.first << " Special";
      for (auto itr1 : itr.second)
        edm::LogVerbatim("HCalGeom") << "ID " << (itr1);
      ++l;
    }
#endif
  }
}

bool HcalDDDRecConstants::isBH ( ) const

inline

Definition at line 109 of file HcalDDDRecConstants.h.

References hcons, and HcalDDDSimConstants::isBH().

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

{ return hcons.isBH(); }

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

inline

Definition at line 110 of file HcalDDDRecConstants.h.

References detIdSp_.

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

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

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

inline

Definition at line 125 of file HcalDDDRecConstants.h.

References detIdSpR_.

Referenced by HcalTopology::validHcal().

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

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

inline

Definition at line 124 of file HcalDDDRecConstants.h.

References detIdSp_.

Referenced by HcalTopology::validHcal().

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

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

private

Definition at line 1080 of file HcalDDDRecConstants.cc.

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

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

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

unsigned int HcalDDDRecConstants::layerGroupSize ( int  eta ) const

private

Definition at line 1067 of file HcalDDDRecConstants.cc.

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

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

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

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

inline

Definition at line 111 of file HcalDDDRecConstants.h.

References hcons, and HcalDDDSimConstants::maxHFDepth().

Referenced by HcalTopology::validRaw().

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

HcalDetId HcalDDDRecConstants::mergedDepthDetId

(

const HcalDetId &

id

)

const

Definition at line 648 of file HcalDDDRecConstants.cc.

References detIdSp_.

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

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

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

Definition at line 581 of file HcalDDDRecConstants.cc.

References funct::abs(), PVValHelper::eta, HcalDDDSimConstants::getDepthEta29(), getMaxDepth(), hcons, iEtaMax, and ecaldqm::zside().

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

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

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

Definition at line 594 of file HcalDDDRecConstants.cc.

References funct::abs(), HLT_FULL_cff::depth, PVValHelper::eta, HcalDDDSimConstants::getDepthEta29(), getMaxDepth(), hcons, iEtaMax, and ecaldqm::zside().

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

unsigned int HcalDDDRecConstants::nCells

(

HcalSubdetector

subdet

)

const

Definition at line 626 of file HcalDDDRecConstants.cc.

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

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

unsigned int HcalDDDRecConstants::nCells

(

)

const

Definition at line 644 of file HcalDDDRecConstants.cc.

References HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

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

unsigned int HcalDDDRecConstants::numberOfCells

(

HcalSubdetector

subdet

)

const

Definition at line 609 of file HcalDDDRecConstants.cc.

References HcalBarrel, HcalCellTypes(), HcalEndcap, hcons, pileupDistInMC::num, and HcalDDDSimConstants::numberOfCells().

Referenced by HcalDDDGeometryLoader::load().

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

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

Definition at line 685 of file HcalDDDRecConstants.cc.

References detIdSp_, and isotrackApplyRegressor::k.

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

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

Definition at line 695 of file HcalDDDRecConstants.cc.

References detIdSp_, and detIdSpR_.

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

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

Definition at line 672 of file HcalDDDRecConstants.cc.

References detIdSpR_, and isotrackApplyRegressor::k.

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

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

bool HcalDDDRecConstants::withSpecialRBXHBHE ( ) const

inline

Definition at line 123 of file HcalDDDRecConstants.h.

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

Referenced by HcalTopology::withSpecialRBXHBHE().

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

Member Data Documentation

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

private

Definition at line 155 of file HcalDDDRecConstants.h.

Referenced by initialize().

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

private

Definition at line 155 of file HcalDDDRecConstants.h.

Referenced by initialize().

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

private

Definition at line 156 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 157 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 144 of file HcalDDDRecConstants.h.

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

std::vector<double> HcalDDDRecConstants::etaTable

private

Definition at line 145 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 152 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 153 of file HcalDDDRecConstants.h.

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

const HcalDDDSimConstants& HcalDDDRecConstants::hcons

private

Definition at line 143 of file HcalDDDRecConstants.h.

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

const HcalParameters* HcalDDDRecConstants::hpar

private

Definition at line 142 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::ietaMap

private

Definition at line 146 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().

std::vector<int> HcalDDDRecConstants::iEtaMax

private

Definition at line 147 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::iEtaMin

private

Definition at line 147 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::maxDepth

private

Definition at line 148 of file HcalDDDRecConstants.h.

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

const int HcalDDDRecConstants::maxLayer_ = 18

staticprivate

Definition at line 140 of file HcalDDDRecConstants.h.

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

const int HcalDDDRecConstants::maxLayerHB_ = 16

staticprivate

Definition at line 141 of file HcalDDDRecConstants.h.

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

int HcalDDDRecConstants::nHalves[2]

private

Definition at line 154 of file HcalDDDRecConstants.h.

Referenced by initialize().

int HcalDDDRecConstants::nModule[2]

private

Definition at line 154 of file HcalDDDRecConstants.h.

Referenced by initialize().

std::vector<int> HcalDDDRecConstants::nPhiBins

private

Definition at line 149 of file HcalDDDRecConstants.h.

Referenced by getNPhi(), and initialize().

std::vector<double> HcalDDDRecConstants::phibin

private

Definition at line 150 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::phiUnitS

private

Definition at line 151 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().