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
std::vector< HcalEtaBin >	getEtaBins (const int itype) const
std::pair< double, double >	getEtaLimit (const int i) const
std::pair< double, double >	getEtaPhi (int subdet, int ieta, 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 (int det, int phi, int zside) const
void	getLayerDepth (int ieta, std::map< int, int > &layers) const
int	getLayerFront (int det, int eta, int phi, 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 (int subdet, 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 (int subdet, int ieta, int depth) const
double	getRZ (int subdet, int ieta, int iphi, int depth) const
double	getRZ (int subdet, int layer) 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
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 (const 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 13 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, and initialize().

                                                    : 
  hpar(hp), hcons(hc) {

#ifdef EDM_ML_DEBUG
  std::cout << "HcalDDDRecConstants::HcalDDDRecConstants (const HcalParameters* hp) constructor" << std::endl;
#endif
  initialize();
}

HcalDDDRecConstants::~HcalDDDRecConstants

(

)

Definition at line 23 of file HcalDDDRecConstants.cc.

References gather_cfg::cout.

                                          { 
#ifdef EDM_ML_DEBUG
  std::cout << "HcalDDDRecConstants::destructed!!!" << std::endl;
#endif
}

Member Function Documentation

const HcalDDDSimConstants* HcalDDDRecConstants::dddConstants ( ) const

inline

Definition at line 124 of file HcalDDDRecConstants.h.

References create_public_pileup_plots::bins, getOneEtaBin(), hcons, mps_fire::i, initialize(), layerGroup(), and layerGroupSize().

{return &hcons;}

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

inline

Definition at line 60 of file HcalDDDRecConstants.h.

References HcalDDDRecConstants::HcalID::eta, gconsHB, gconsHE, getDepth(), HcalDDDRecConstants::HcalID::phi, and HcalDigiParam_cfi::zside.

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 50 of file HcalDDDRecConstants.cc.

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

Referenced by getConstHBHE(), and HcalTopology::HcalTopology().

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

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

Definition at line 29 of file HcalDDDRecConstants.cc.

References HcalLayerDepthMap::getLayerDepth(), hcons, hpar, plotBeamSpotDB::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.push_back(layers[lay+1]);
    return depths;
  }
}

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

inline

Definition at line 68 of file HcalDDDRecConstants.h.

References HcalDDDSimConstants::getDepthEta16(), getEtaBins(), getEtaPhi(), hcons, and HcalDDDRecConstants::HcalID::subdet.

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

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

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

(

const int

itype

)

const

Definition at line 64 of file HcalDDDRecConstants.cc.

References create_public_pileup_plots::bins, gather_cfg::cout, spr::find(), plotBeamSpotDB::first, HcalLayerDepthMap::getLayerDepth(), getLayerDepth(), getOneEtaBin(), getPhis(), hcons, mps_fire::i, iEtaMax, iEtaMin, gen::k, LayerTriplets::layers(), HcalDDDSimConstants::ldMap(), edm::second(), HcalLayerDepthMap::validDet(), and HcalDigiParam_cfi::zside.

Referenced by getDepthEta16(), 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 (unsigned k=0; k<phis.size(); ++k) {
      if (std::find(phiSp.begin(),phiSp.end(),(zside*phis[k].first)) ==
          phiSp.end()){
        phiUse.push_back(phis[k]);
      }
    }
      } else {
    phiUse.insert(phiUse.end(),phis.begin(),phis.end());
      }
      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 (unsigned i=0; i<phiSp.size(); ++i) {
      for (unsigned k=0; k<phis.size(); ++k) {
        if (phiSp[i] == zside*phis[k].first) {
          phiUse.push_back(phis[k]);
          break;
        }
      }
    }
    if (phiUse.size() > 0) {
      hcons.ldMap()->getLayerDepth(subdet,ieta,phiUse[0].first,zside,layers);
      getOneEtaBin(subdet,ieta,zside,phiUse,layers,true,bins);
    }
      }
    }
  }
#ifdef EDM_ML_DEBUG
  std::cout << "Prepares " << bins.size() << " eta bins for type " << type 
        << std::endl;
  for (unsigned int i=0; i<bins.size(); ++i) {
    std::cout << "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)
      std::cout << " [" << k << "] " << bins[i].layer[k].first << ":"
        << bins[i].layer[k].second;
    std::cout << std::endl << " and Phi sets";
    for (unsigned int k=0; k<bins[i].phis.size(); ++k)
      std::cout << " " << bins[i].phis[k].first << ":" <<bins[i].phis[k].second;
    std::cout << std::endl;
  }
#endif
  return bins;
}

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

inline

Definition at line 76 of file HcalDDDRecConstants.h.

References HcalDDDRecConstants::HcalID::depth, etaTable, getHCID(), getHFCellParameters(), getLayerDepth(), getLayerFront(), mps_fire::i, and LayerTriplets::layers().

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

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

Definition at line 134 of file HcalDDDRecConstants.cc.

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

Referenced by HGCalSimHitValidation::analyzeHits(), and getDepthEta16().

                                                                   {
  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
  std::cout << "getEtaPhi: subdet|ieta|iphi " << subdet << "|" << ieta << "|"
        << iphi << " eta|phi " << eta << "|" << phi << std::endl;
#endif
  return std::pair<double,double>(eta,phi);
}

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

inline

Definition at line 72 of file HcalDDDRecConstants.h.

References mps_fire::i, iEtaMax, and iEtaMin.

Referenced by 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 74 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 75 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 163 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, stringResolutionProvider_cfi::eta, HcalDDDSimConstants::findDepth(), HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta29(), HcalBarrel, HcalEndcap, HcalOuter, hcons, hpar, ietaMap, iEtaMin, createfilelist::int, layerGroup(), maxLayerHB_, HcalParameters::noff, phi, phibin, HcalParameters::phigroup, HcalDDDSimConstants::phiNumber(), HcalParameters::phioff, phiUnitS, csvLumiCalc::unit, HcalDDDSimConstants::unitPhi(), and HcalDigiParam_cfi::zside.

Referenced by HGCalHitValidation::analyze(), CaloTruthAccumulator::fillSimHits(), getEtaLimit(), 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
  std::cout << "getHCID: input " << subdet << ":" << ieta << ":" << iphi
        << ":" << idepth << ":" << lay << " output " << eta << ":" << phi
        << ":" << depth << std::endl;
#endif
  return HcalDDDRecConstants::HcalID(subdet,eta,phi,depth);
}

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

(

)

const

Definition at line 216 of file HcalDDDRecConstants.cc.

References funct::abs(), gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, HcalDDDSimConstants::getIdHF2QIE(), HcalDDDSimConstants::getPhiTableHF(), HcalDDDSimConstants::getRTableHF(), hcons, iEtaMin, createfilelist::int, gen::k, maxDepth, DetIdAssociatorESProducer_cff::nEta, DetIdAssociatorESProducer_cff::nPhi, conversionPostprocessing_cfi::rMax, and conversionPostprocessing_cfi::rMin.

Referenced by getEtaLimit(), and 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.0*CLHEP::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.push_back(cell1);
      HcalDDDRecConstants::HFCellParameters cell2(-ieta,1,iphi,dphi,nphi,rMin,rMax);
      cells.push_back(cell2);
    }
  }
  if (maxDepth[2] > 2) {
    if (hcons.getIdHF2QIE().size() > 0) {
      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.0*CLHEP::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.push_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.0*CLHEP::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.push_back(cell1);
    HcalDDDRecConstants::HFCellParameters cell2(-ieta,3,iphi,dphi,nphi,rMin,rMax);
    cells.push_back(cell2);
      }
    }
  }
#ifdef EDM_ML_DEBUG
  std::cout << "HcalDDDRecConstants returns " << cells.size() 
        << " HF cell parameters" << std::endl;
  for (unsigned int k=0; k<cells.size(); ++k)
    std::cout << "Cell[" << k <<"] : (" << cells[k].ieta <<", "<< cells[k].depth
          << ", " << cells[k].firstPhi << ", " << cells[k].stepPhi << ", "
          << cells[k].nPhi << ", " << cells[k].rMin << ", "
          << cells[k].rMax << ")" << std::endl;
#endif
  return cells;
}

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

inline

Definition at line 83 of file HcalDDDRecConstants.h.

References HcalDDDSimConstants::getLayer0Wt(), and hcons.

Referenced by HcalHitRelabeller::energyWt().

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

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

Definition at line 273 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, checklumidiff::l, layerGroup(), and layerGroupSize().

Referenced by getEtaBins(), and getEtaLimit().

                                                   {

  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
  std::cout << "getLayerDepth::Input " << ieta << " Output " 
        << layers.size() << " entries" << std::endl;
  for (std::map<int,int>::iterator itr=layers.begin(); itr != layers.end();
       ++itr) std::cout << " [" << itr->first << "] " << itr->second;
  std::cout << std::endl;
#endif
}

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

Definition at line 290 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, stringResolutionProvider_cfi::eta, ALCARECOTkAlBeamHalo_cff::etaMax, HcalDDDSimConstants::getFrontLayer(), HcalLayerDepthMap::getLayerFront(), hcons, gen::k, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), and HcalDigiParam_cfi::zside.

Referenced by getEtaLimit(), and getRZ().

                                                              {
  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);
  if (layFront < 0) {
    int laymin  = hcons.getFrontLayer(subdet, ieta);
    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;
      }
    }
      }
    }
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getLayerFront::Input " << idet << ":" << ieta << ":"
        << iphi << ":" << depth << " Output " << layFront << std::endl;
#endif
  return layFront;
}

int HcalDDDRecConstants::getMaxDepth ( const int  type ) const

inline

Definition at line 84 of file HcalDDDRecConstants.h.

References getMinDepth(), and maxDepth.

Referenced by HGCalDigiClient::beginRun(), HGCalRecHitsClient::beginRun(), HGCalSimHitsClient::beginRun(), HcalSimHitsClient::beginRun(), HcalRecHitsDQMClient::beginRun(), HcalSimHitStudy::bookHistograms(), SimHitsValidationHcal::bookHistograms(), HcalSimHitsValidation::bookHistograms(), HcalTopology::depthBinInformation(), HGCalDigiValidation::dqmBeginRun(), HcalRecHitsAnalyzer::dqmBeginRun(), HGCalRecHitValidation::dqmBeginRun(), HcalDigisValidation::dqmBeginRun(), HGCalSimHitValidation::dqmBeginRun(), HcalTopology::HcalTopology(), 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 316 of file HcalDDDRecConstants.cc.

References HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29M(), HcalDDDSimConstants::getMaxDepth(), hcons, hpar, iEtaMax, createfilelist::int, 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);
    }
  }
  return lmax;
}

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

Definition at line 334 of file HcalDDDRecConstants.cc.

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

Referenced by HcalTopology::depthBinInformation(), HcalGeometry::getClosestCell(), and getMaxDepth().

                                                             {

  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 89 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 90 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 91 of file HcalDDDRecConstants.h.

References nPhiBins.

Referenced by 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 648 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, HcalDDDRecConstants::HcalEtaBin::depthStart, etaTable, plotBeamSpotDB::first, HcalDDDSimConstants::getDepthEta16(), HcalDDDSimConstants::getDepthEta16M(), HcalDDDSimConstants::getDepthEta29(), HcalDDDSimConstants::getDepthEta29M(), HcalBarrel, HcalEndcap, hcons, hpar, iEtaMax, iEtaMin, gen::k, checklumidiff::l, HcalDDDRecConstants::HcalEtaBin::layer, maxLayer_, maxLayerHB_, gen::n, HcalParameters::noff, phibin, HcalDDDRecConstants::HcalEtaBin::phis, and edm::second().

Referenced by dddConstants(), and 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.size() > 0) {
    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
    std::cout << "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)
      std::cout << " [" << l << "] " << itr->first << ":" << itr->second;
    std::cout  << std::endl << " with " << phis.size() << " phis";
    for (unsigned int l=0; l<phis.size(); ++l)
      std::cout << " " << phis[l].first << ":" << phis[l].second;
    std::cout << std::endl;
#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.push_back(std::pair<int,int>(lmin,lmax0));
      } else {
        etabin.layer.push_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.push_back(std::pair<int,int>(lmin,lmax));
    etabin0.phis.insert(etabin0.phis.end(),phis.begin(),phis.end());
    bins.push_back(etabin0);
#ifdef EDM_ML_DEBUG
    std::cout << "etabin0: dStatrt " << etabin0.depthStart << " layers "
          << etabin0.layer.size() << ":" << lmin << ":" << lmax
          << " phis " << phis.size() << std::endl;
    for (unsigned int k=0; k<etabin0.layer.size(); ++k)
      std::cout << " [" << k << "] " << etabin0.layer[k].first << ":"
            << etabin0.layer[k].second;
    std::cout << std::endl;
#endif
      } else if (ieta == hpar->noff[1]) {
      } else {
    etabin.layer.push_back(std::pair<int,int>(lmin,lmax));
    if (dstart < 0) dstart = dep;
      }
    }
  }
  etabin.depthStart = dstart;
  bins.push_back(etabin);
#ifdef EDM_ML_DEBUG
  std::cout << "etabin: dStatrt " << etabin.depthStart << " layers "
        << etabin.layer.size() << ":" << lmin << ":" << lmax
        << " phis " << etabin.phis.size() << std::endl;
  for (unsigned int k=0; k<etabin.layer.size(); ++k)
    std::cout << " [" << k << "] " << etabin.layer[k].first << ":"
          << etabin.layer[k].second;
  std::cout << std::endl;
#endif
}

double HcalDDDRecConstants::getPhiBin ( const int  i ) const

inline

Definition at line 92 of file HcalDDDRecConstants.h.

References mps_fire::i, and phibin.

{return phibin[i];}

double HcalDDDRecConstants::getPhiOff ( const int  i ) const

inline

Definition at line 93 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 94 of file HcalDDDRecConstants.h.

References getPhis(), hpar, and HcalParameters::phioff.

Referenced by HcalTopology::HcalTopology().

{return hpar->phioff;}

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

Definition at line 357 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, etaSimValu, plotBeamSpotDB::first, HcalDDDSimConstants::getPhiCons(), HcalForward, hcons, gen::k, phi, phibin, HcalDDDSimConstants::phiNumber(), edm::second(), HcalDDDSimConstants::unitPhi(), and units().

Referenced by getEtaBins(), and getPhiOffs().

                                                       {

  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((CLHEP::twopi+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.push_back(std::pair<int,double>(iphi,phi));
  }
#ifdef EDM_ML_DEBUG
  std::cout << "getEtaPhi: 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>& HcalDDDRecConstants::getPhiTable ( ) const

inline

Definition at line 96 of file HcalDDDRecConstants.h.

References phibin.

Referenced by HcalTopology::HcalTopology().

{return phibin;}

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

inline

Definition at line 97 of file HcalDDDRecConstants.h.

References getPhiZOne(), getRZ(), getThickActive(), 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 382 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, plotBeamSpotDB::first, HcalLayerDepthMap::getPhis(), HcalLayerDepthMap::getSubdet(), hcons, gen::k, HcalDDDSimConstants::ldMap(), phi, edm::second(), and HcalDigiParam_cfi::zside.

Referenced by getPhiTableHF(), and HcalTopology::getPhiZOne().

                                                                           {

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

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

Definition at line 404 of file HcalDDDRecConstants.cc.

Referenced by HGCalSimHitValidation::analyzeHits(), and getPhiTableHF().

                                                                       {

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

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

Definition at line 409 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, plotBeamSpotDB::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
  std::cout << "getRZ: subdet|ieta|ipho|depth " << subdet << "|" << ieta << "|"
        << iphi << "|" << depth << " lay|rz " << layf << "|" << rz 
        << std::endl;
#endif
  return rz;
}

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

Definition at line 423 of file HcalDDDRecConstants.cc.

References gather_cfg::cout, plotBeamSpotDB::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
  std::cout << "getRZ: subdet|layer " << subdet << "|" << layer << " rz "
        << rz << std::endl;
#endif
  return rz;
}

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

(

const int

type

)

const

Definition at line 438 of file HcalDDDRecConstants.cc.

References create_public_pileup_plots::bins, funct::cos(), gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, HcalDDDRecConstants::HcalActiveLength::depth, stringResolutionProvider_cfi::eta, HcalDDDRecConstants::HcalActiveLength::eta, JetChargeProducer_cfi::exp, plotBeamSpotDB::first, gconsHB, gconsHE, getEtaBins(), mps_fire::i, HcalDDDRecConstants::HcalActiveLength::ieta, gen::k, Scenarios_cff::scale, edm::second(), funct::sin(), HcalDDDRecConstants::HcalActiveLength::stype, lumiQTWidget::t, theta(), HcalDDDRecConstants::HcalActiveLength::thick, HcalDigiParam_cfi::zside, and HcalDDDRecConstants::HcalActiveLength::zside.

Referenced by HcalHBHEMuonSimAnalyzer::beginRun(), HcalRaddamMuon::beginRun(), HcalHBHEMuonAnalyzer::beginRun(), and getPhiTableHF().

                                                        {

  std::vector<HcalDDDRecConstants::HcalActiveLength> actives;
  std::vector<HcalDDDRecConstants::HcalEtaBin> bins = getEtaBins(type);
#ifdef EDM_ML_DEBUG
  unsigned int kount(0);
#endif
  for (unsigned int k=0; k<bins.size(); ++k) {
    int    ieta  = bins[k].ieta;
    int    zside = bins[k].zside;
    int    stype = (bins[k].phis.size() > 4) ? 0 : 1;
    double eta   = 0.5*(bins[k].etaMin+bins[k].etaMax);
    double theta = 2*atan(exp(-eta));
    double scale = 1.0/((type == 0) ? sin(theta) : cos(theta));
    int    depth = bins[k].depthStart;
    for (unsigned int i = 0; i < bins[k].layer.size(); ++i) {
      double thick(0);
      for (int j = bins[k].layer[i].first; j <= bins[k].layer[i].second; ++j) {
    if (type == 0 || j > 1) {
      double t = ((type == 0) ? gconsHB[j-1].second : gconsHE[j-1].second);
      if (t > 0) thick += t;
    }
      }
      thick *= (2.*scale);
      HcalDDDRecConstants::HcalActiveLength active(ieta,depth,zside,stype,zside*eta,thick);
      actives.push_back(active);
      ++depth;
#ifdef EDM_ML_DEBUG
      kount++;
      std::cout << "getThickActive: [" << kount << "] eta:" << active.ieta 
        << ":" << active.eta << " zside " << active.zside << " depth " 
        << active.depth << " type " << active.stype << " thick "
        << active.thick << std::endl;
#endif
    }
  }
  return actives;
}

int HcalDDDRecConstants::getTopoMode ( ) const

inline

Definition at line 103 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 104 of file HcalDDDRecConstants.h.

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

Referenced by HcalTopology::HcalTopology().

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

std::vector< HcalCellType > HcalDDDRecConstants::HcalCellTypes

(

HcalSubdetector

subdet

)

const

Definition at line 478 of file HcalDDDRecConstants.cc.

References KineDebug3::count(), gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, etaSimValu, plotBeamSpotDB::first, getEtaBins(), HcalBarrel, HcalDDDSimConstants::HcalCellTypes(), HcalEndcap, hcons, hpar, iEtaMax, createfilelist::int, HcalParameters::phioff, groupFilesInBlocks::temp, csvLumiCalc::unit, and HcalDDDSimConstants::unitPhi().

Referenced by HcalDDDGeometryLoader::fill(), getTriggerMode(), 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.push_back(cell);
    count.push_back(0);
    dmin.push_back(0);
    dmax.push_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 (unsigned int ic=0; ic<cellsm.size(); ++ic) {
        if (cellsm[ic].depthSegment() >= etabin.layer[il].first &&
        cellsm[ic].depthSegment() <= etabin.layer[il].second &&
        cellsm[ic].etaBin() == temp[il].etaBin() && 
        cellsm[ic].zside()  == temp[il].zside()) {
          if (count[il] == 0) {
        temp[il] = cellsm[ic];
        dmin[il] = cellsm[ic].depthMin();
        dmax[il] = cellsm[ic].depthMax();
          }
          ++count[il];
          if (cellsm[ic].depthMin() < dmin[il]) 
        dmin[il] = cellsm[ic].depthMin();
          if (cellsm[ic].depthMax() > dmax[il]) 
        dmax[il] = cellsm[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.push_back(temp[il]);
      }
    }
#ifdef EDM_ML_DEBUG
    std::cout << "HcalDDDRecConstants: found " << cells.size() 
          << " cells for sub-detector type " << isub << std::endl;
    for (unsigned int ic=0; ic<cells.size(); ++ic)
      std::cout << "Cell[" << ic << "] " << cells[ic] << std::endl;
#endif
    return cells;
  } else {
    return hcons.HcalCellTypes(subdet,-1,-1);
  }
}

HcalDetId HcalDDDRecConstants::idBack

(

const HcalDetId &

id

)

const

Definition at line 601 of file HcalDDDRecConstants.cc.

References detIdSpR_.

Referenced by HcalTopology::idBack(), and maxHFDepth().

                                                               {

  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 592 of file HcalDDDRecConstants.cc.

References particleFlowClusterECALTimeSelected_cfi::depth, and detIdSpR_.

Referenced by HcalTopology::idFront(), and maxHFDepth().

                                                                {

  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 743 of file HcalDDDRecConstants.cc.

References KineDebug3::count(), gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, depthMaxDf_, depthMaxSp_, detIdSp_, detIdSpR_, HcalParameters::drHB, HcalParameters::dzHE, stringResolutionProvider_cfi::eta, HcalParameters::etagroup, ALCARECOTkAlBeamHalo_cff::etaMax, HcalParameters::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, HcalParameters::etaMin, EnergyCorrector::etas, etaSimValu, HcalParameters::etaTable, etaTable, HcalParameters::etaTableHF, Exception, plotBeamSpotDB::first, RemoveAddSevLevel::flag, 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, createfilelist::int, gen::k, GetRecoTauVFromDQM_MC_cff::kk, checklumidiff::l, layerGroup(), layerGroupSize(), HcalDDDSimConstants::ldMap(), genParticles_cff::map, HcalParameters::maxDepth, maxDepth, maxLayer_, maxLayerHB_, HcalParameters::modHB, HcalParameters::modHE, DetIdAssociatorESProducer_cff::nEta, nHalves, nModule, HcalParameters::noff, nPhiBins, phi, HcalParameters::phibin, phibin, HcalParameters::phigroup, HcalParameters::phitable, phiUnitS, HcalParameters::rHB, edm::second(), csvLumiCalc::unit, HcalDDDSimConstants::unitPhi(), HcalLayerDepthMap::validDet(), HcalLayerDepthMap::validEta(), HcalParameters::zHE, and HcalDigiParam_cfi::zside.

Referenced by dddConstants(), and 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.push_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") << "Going beyond the array boundary "
                << hpar->etaTable.size() << " at index " << i 
                << " of etaTable from SimConstant";
      throw cms::Exception("DDException") << "Going beyond the array boundary "
                      << hpar->etaTable.size() 
                      << " at index " << i 
                      << " of etaTable from SimConstant";
    } else {
      etaTable.push_back(hpar->etaTable[ieta]);
      etaSimValu.push_back(std::pair<int,int>(ef,ieta));
    }
    for (int k=0; k<(hpar->etagroup[i]); ++k) ietaMap.push_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.push_back(0);
  ieta = 0;
  phibin.clear(); phiUnitS.clear();
  for (int i=0; i<nEta; ++i) {
    double dphi = (hpar->phigroup[i])*(hpar->phibin[ieta]);
    phibin.push_back(dphi);
    int    nphi = (int)((CLHEP::twopi + 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.push_back(unit);
  }
  for (unsigned int i=0; i<hpar->phitable.size(); ++i)  {
    int  nphi = (int)((CLHEP::twopi + 0.001)/hpar->phitable[i]);
    if (nphi > nPhiBins[2]) nPhiBins[2] = nphi;
  }
#ifdef EDM_ML_DEBUG
  std::cout << "Modified eta/deltaphi table for " << nEta << " bins" << std::endl;
  for (int i=0; i<nEta; ++i)
    std::cout << "Eta[" << i << "] = " << etaTable[i] << ":" << etaTable[i+1]
          << ":" << etaSimValu[i].first << ":" << etaSimValu[i].second
          << " PhiBin[" << i << "] = " << phibin[i]/CLHEP::deg <<std::endl;
  std::cout << "PhiUnitS";
  for (unsigned int i=0; i<phiUnitS.size(); ++i)
    std::cout << " [" << i << "] = " << phiUnitS[i];
  std::cout << std::endl;
  std::cout << "nPhiBins";
  for (unsigned int i=0; i<nPhiBins.size(); ++i)
    std::cout << " [" << i << "] = " << nPhiBins[i];
  std::cout << std::endl;
  std::cout << "EtaTableHF";
  for (unsigned int i=0; i<hpar->etaTableHF.size(); ++i)
    std::cout << " [" << i << "] = " << hpar->etaTableHF[i];
  std::cout << std::endl;
  std::cout << "PhiBinHF";
  for (unsigned int i=0; i<hpar->phitable.size(); ++i)
    std::cout << " [" << i << "] = " << hpar->phitable[i];
  std::cout << std::endl;
#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
      std::cout << "HB " << i << " " << imx << " " << laymax << " " << laymax0 << std::endl;
#endif
      if (maxDepth[0] < laymax0) maxDepth[0] = laymax0;
    }
    if (i >= iEtaMin[1]-1 && i < iEtaMax[1]) {
#ifdef EDM_ML_DEBUG
      std::cout << "HE " << i << " " << imx << " " << laymax << std::endl;
#endif
      if (maxDepth[1] < laymax) maxDepth[1] = laymax;
    }
  }
#ifdef EDM_ML_DEBUG
  for (int i=0; i<4; ++i) 
    std::cout << "Detector Type[" << i << "] iEta " << iEtaMin[i] << ":"
              << iEtaMax[i] << " MaxDepth " << maxDepth[i] << std::endl; 
#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.push_back(std::pair<double,double>(hpar->rHB[i]/CLHEP::cm,
                           hpar->drHB[i]/CLHEP::cm));
  }
#ifdef EDM_ML_DEBUG
  std::cout << "HB with " << nModule[0] << " modules and " << nHalves[0]
        <<" halves and " << gconsHB.size() << " layers" << std::endl;
  for (unsigned int i=0; i<gconsHB.size(); ++i) 
    std::cout << "rHB[" << i << "] = " << gconsHB[i].first << " +- "
          << gconsHB[i].second << std::endl; 
#endif
  nModule[1] = hpar->modHE[0];
  nHalves[1] = hpar->modHE[1];
  for (unsigned int i=0; i<hpar->zHE.size(); ++i) {
    gconsHE.push_back(std::pair<double,double>(hpar->zHE[i]/CLHEP::cm,
                           hpar->dzHE[i]/CLHEP::cm));
  }
#ifdef EDM_ML_DEBUG
  std::cout << "HE with " << nModule[1] << " modules and " << nHalves[1] 
        <<" halves and " << gconsHE.size() << " layers" << std::endl;
  for (unsigned int i=0; i<gconsHE.size(); ++i) 
    std::cout << "zHE[" << i << "] = " << gconsHE[i].first << " +- "
          << gconsHE[i].second << std::endl; 
#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
  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

  //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      
      std::cout << "Eta|Phi|Zside " << eta << ":" << phi << ":" << zside
        << " with " << oldDep.size() << " old Depths" << std::endl;
      unsigned int kk(0);
      for (std::map<int,std::pair<int,int> >::const_iterator itr=oldDep.begin(); itr != oldDep.end(); ++itr,++kk)
    std::cout << "[" << kk << "] " << itr->first << " --> " 
          << itr->second.first << ":" << itr->second.second << "\n";
#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      
      std::cout << "New Depth " << ndepth << " old Depth " << odepth 
            << " max " << maxlay << std::endl;
#endif
      for (unsigned int k=0; k<phis.size(); ++k) {
        zside  = (phis[k] > 0) ? 1 : -1;
        phi    = (phis[k] > 0) ? phis[k] : -phis[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.push_back(newId);
        detIdSpR_[oldId] = ids;
      }
    }
      }
    }
#ifdef EDM_ML_DEBUG
    std::cout << "Map for merging new channels to old channel IDs with "
          << detIdSp_.size() << " entries" << std::endl;
    int l(0);
    for (auto itr : detIdSp_) {
      std::cout << "[" << l << "] Special " << itr.first << " Standard "
        << itr.second << std::endl;
      ++l;
    }
    std::cout << "Reverse Map for mapping old to new IDs with "
          << detIdSpR_.size() << " entries" << std::endl;
    l = 0;
    for (auto itr : detIdSpR_) {
      std::cout << "[" << l << "] Standard " << itr.first << " Special";
      for (auto itr1 : itr.second) 
    std::cout << " " << (itr1);
      std::cout << std::endl;
      ++l;
    }
#endif
  }

}

bool HcalDDDRecConstants::isBH ( ) const

inline

Definition at line 106 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 107 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 123 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 122 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 1021 of file HcalDDDRecConstants.cc.

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

Referenced by dddConstants(), getHCID(), 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 1009 of file HcalDDDRecConstants.cc.

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

Referenced by dddConstants(), 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 108 of file HcalDDDRecConstants.h.

References RemoveAddSevLevel::flag, hcons, idBack(), idFront(), HcalDDDSimConstants::maxHFDepth(), mergedDepthDetId(), nCells(), numberOfCells(), specialRBXHBHE(), and unmergeDepthDetId().

Referenced by HcalTopology::validRaw().

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

HcalDetId HcalDDDRecConstants::mergedDepthDetId

(

const HcalDetId &

id

)

const

Definition at line 585 of file HcalDDDRecConstants.cc.

References detIdSp_.

Referenced by maxHFDepth(), and HcalTopology::mergedDepthDetId().

                                                                         {

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

unsigned int HcalDDDRecConstants::nCells

(

HcalSubdetector

subdet

)

const

Definition at line 562 of file HcalDDDRecConstants.cc.

References getEtaBins(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, mps_fire::i, 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 (unsigned int i=0; i<etabins.size(); ++i) {
      ncell += ((etabins[i].phis.size())*(etabins[i].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 581 of file HcalDDDRecConstants.cc.

References HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

Referenced by maxHFDepth().

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

unsigned int HcalDDDRecConstants::numberOfCells

(

HcalSubdetector

subdet

)

const

Definition at line 543 of file HcalDDDRecConstants.cc.

References HcalBarrel, HcalCellTypes(), HcalEndcap, hcons, mps_fire::i, createfilelist::int, nPhiBins, pileupDistInMC::num, and HcalDDDSimConstants::numberOfCells().

Referenced by HcalDDDGeometryLoader::load(), and maxHFDepth().

                                                                            {

  if (subdet == HcalBarrel || subdet == HcalEndcap) {
    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
    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 625 of file HcalDDDRecConstants.cc.

References detIdSp_, and gen::k.

Referenced by maxHFDepth().

                                                         {
  for (unsigned int k=0; k<idsOld.size(); ++k) {
    std::map<HcalDetId,HcalDetId>::const_iterator itr = detIdSp_.find(idsOld[k]);
    if (itr == detIdSp_.end()) idsNew.push_back(idsOld[k]);
    else                       idsNew.push_back(itr->second);
  }
}

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

Definition at line 634 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.push_back(itr->first);
  } else{
    std::map<HcalDetId,std::vector<HcalDetId>>::const_iterator itr;
    for (itr = detIdSpR_.begin(); itr != detIdSpR_.end(); ++itr) 
      ids.push_back(itr->first);
  }
  return (ids.size() > 0);
}

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

Definition at line 610 of file HcalDDDRecConstants.cc.

References detIdSpR_, and gen::k.

Referenced by maxHFDepth(), and HcalTopology::unmergeDepthDetId().

                                                         {

  ids.clear();
  std::map<HcalDetId,std::vector<HcalDetId>>::const_iterator itr = detIdSpR_.find(id);
  if (itr == detIdSpR_.end()) {
    ids.push_back(id);
  } else {
    for (unsigned k=0; k<itr->second.size(); ++k) {
      HcalDetId hid(id.subdet(),id.ieta(),id.iphi(),(itr->second)[k].depth());
      ids.push_back(hid);
    }
  }
}

bool HcalDDDRecConstants::withSpecialRBXHBHE ( ) const

inline

Definition at line 121 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 151 of file HcalDDDRecConstants.h.

Referenced by initialize().

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

private

Definition at line 151 of file HcalDDDRecConstants.h.

Referenced by initialize().

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

private

Definition at line 152 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 153 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 140 of file HcalDDDRecConstants.h.

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

std::vector<double> HcalDDDRecConstants::etaTable

private

Definition at line 141 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 148 of file HcalDDDRecConstants.h.

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

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

private

Definition at line 149 of file HcalDDDRecConstants.h.

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

const HcalDDDSimConstants& HcalDDDRecConstants::hcons

private

Definition at line 139 of file HcalDDDRecConstants.h.

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

const HcalParameters* HcalDDDRecConstants::hpar

private

Definition at line 138 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 142 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().

std::vector<int> HcalDDDRecConstants::iEtaMax

private

Definition at line 143 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::iEtaMin

private

Definition at line 143 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::maxDepth

private

Definition at line 144 of file HcalDDDRecConstants.h.

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

const int HcalDDDRecConstants::maxLayer_ =18

staticprivate

Definition at line 136 of file HcalDDDRecConstants.h.

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

const int HcalDDDRecConstants::maxLayerHB_ =16

staticprivate

Definition at line 137 of file HcalDDDRecConstants.h.

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

int HcalDDDRecConstants::nHalves[2]

private

Definition at line 150 of file HcalDDDRecConstants.h.

Referenced by initialize().

int HcalDDDRecConstants::nModule[2]

private

Definition at line 150 of file HcalDDDRecConstants.h.

Referenced by initialize().

std::vector<int> HcalDDDRecConstants::nPhiBins

private

Definition at line 145 of file HcalDDDRecConstants.h.

Referenced by getNPhi(), initialize(), and numberOfCells().

std::vector<double> HcalDDDRecConstants::phibin

private

Definition at line 146 of file HcalDDDRecConstants.h.

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

std::vector<int> HcalDDDRecConstants::phiUnitS

private

Definition at line 147 of file HcalDDDRecConstants.h.

Referenced by getHCID(), and initialize().