HcalTopology Class Reference

#include <HcalTopology.h>

Inheritance diagram for HcalTopology:

Public Member Functions
const HcalDDDRecConstants *	dddConstants () const
int	decIEta (const HcalDetId &id, HcalDetId neighbors[2]) const
bool	decIPhi (const HcalDetId &id, HcalDetId &neighbor) const
bool	decrementDepth (HcalDetId &id) const
DetId	denseId2detId (unsigned int) const override
	return a linear packed id More...
HcalCalibDetId	denseId2detIdCALIB (const unsigned int &id) const
	return a Calib DetId from linear packed id More...
void	depthBinInformation (HcalSubdetector subdet, int etaRing, int iphi, int zside, int &nDepthBins, int &startingBin) const
	finds the number of depth bins and which is the number to start with More...
unsigned int	detId2denseId (const DetId &id) const override
	return a linear packed id More...
unsigned int	detId2denseIdCALIB (const DetId &id) const
	return a linear packed id from CALIB More...
unsigned int	detId2denseIdHB (const DetId &id) const
	return a linear packed id from HB More...
unsigned int	detId2denseIdHE (const DetId &id) const
	return a linear packed id from HE More...
unsigned int	detId2denseIdHF (const DetId &id) const
	return a linear packed id from HF More...
unsigned int	detId2denseIdHO (const DetId &id) const
	return a linear packed id from HO More...
unsigned int	detId2denseIdHT (const DetId &id) const
	return a linear packed id from HT More...
int	doublePhiBins () const
std::vector< DetId >	down (const DetId &id) const override
std::vector< DetId >	east (const DetId &id) const override
double	etaMax (HcalSubdetector subdet) const
std::pair< double, double >	etaRange (HcalSubdetector subdet, int ieta) const
int	etaRing (HcalSubdetector subdet, double eta) const
	eta and phi index from eta, phi values More...
void	exclude (const HcalDetId &id)
int	exclude (HcalSubdetector subdet, int ieta1, int ieta2, int iphi1, int iphi2, int depth1=1, int depth2=4)
void	excludeSubdetector (HcalSubdetector subdet)
int	firstHBRing () const
int	firstHEDoublePhiRing () const
int	firstHEQuadPhiRing () const
int	firstHERing () const
int	firstHETripleDepthRing () const
int	firstHFQuadPhiRing () const
int	firstHFRing () const
int	firstHORing () const
unsigned int	getCALIBSize () const
void	getDepthSegmentation (const unsigned ring, std::vector< int > &readoutDepths, const bool flag=false) const
unsigned int	getHBSize () const
unsigned int	getHESize () const
unsigned int	getHFSize () const
unsigned int	getHOSize () const
unsigned int	getHTSize () const
bool	getMergePositionFlag () const
unsigned int	getNumberOfShapes () const
int	getPhiZOne (std::vector< std::pair< int, int > > &phiz) const
	HcalTopology (const HcalDDDRecConstants *hcons, const bool mergePosition=false)
	HcalTopology (HcalTopologyMode::Mode mode, int maxDepthHB, int maxDepthHE, HcalTopologyMode::TriggerMode tmode=HcalTopologyMode::TriggerMode_2009)
HcalDetId	idBack (const HcalDetId &id) const
HcalDetId	idFront (const HcalDetId &id) const
int	incIEta (const HcalDetId &id, HcalDetId neighbors[2]) const
bool	incIPhi (const HcalDetId &id, HcalDetId &neighbor) const
bool	incrementDepth (HcalDetId &id) const
bool	isBH () const
int	lastHBHERing () const
int	lastHBRing () const
int	lastHERing () const
int	lastHFRing () const
int	lastHORing () const
int	maxDepth (void) const
int	maxDepth (HcalSubdetector subdet) const
int	maxDepthHB () const
int	maxDepthHE () const
bool	mergedDepth29 (HcalDetId id) const
HcalDetId	mergedDepthDetId (const HcalDetId &id) const
std::vector< int >	mergedDepthList29 (HcalDetId id) const
HcalTopologyMode::Mode	mode () const
unsigned int	ncells () const override
	return a count of valid cells (for dense indexing use) More...
std::vector< DetId >	north (const DetId &id) const override
int	nPhiBins (int etaRing) const
	how many phi segments in this ring More...
int	nPhiBins (HcalSubdetector subdet, int etaRing) const
int	phiBin (HcalSubdetector subdet, int etaRing, double phi) const
std::pair< int, int >	segmentBoundaries (const unsigned ring, const unsigned depth, const bool flag=false) const
void	setDepthSegmentation (const unsigned ring, const std::vector< int > &readoutDepths, const bool flag)
int	singlePhiBins () const
std::vector< DetId >	south (const DetId &id) const override
int	topoVersion () const override
	return a version which identifies the given topology More...
HcalTopologyMode::TriggerMode	triggerMode () const
void	unmergeDepthDetId (const HcalDetId &id, std::vector< HcalDetId > &ids) const
std::vector< DetId >	up (const DetId &id) const override
bool	valid (const DetId &id) const override
bool	validCalib (const HcalCalibDetId &id) const
bool	validDetId (HcalSubdetector subdet, int ieta, int iphi, int depth) const
bool	validHcal (const HcalDetId &id) const
bool	validHcal (const HcalDetId &id, const unsigned int flag) const
bool	validHT (const HcalTrigTowerDetId &id) const
std::vector< DetId >	west (const DetId &id) const override
bool	withSpecialRBXHBHE () const
Public Member Functions inherited from CaloSubdetectorTopology
	CaloSubdetectorTopology ()
	standard constructor More...
virtual bool	denseIdConsistent (int topoVer) const
	return whether this topology is consistent with the numbering in the given topology More...
virtual std::vector< DetId >	getAllNeighbours (const DetId &id) const
virtual std::vector< DetId >	getNeighbours (const DetId &id, const CaloDirection &dir) const
virtual std::vector< DetId >	getWindow (const DetId &id, const int &northSouthSize, const int &eastWestSize) const
virtual DetId	goDown (const DetId &id) const
virtual DetId	goEast (const DetId &id) const
virtual DetId	goNorth (const DetId &id) const
virtual DetId	goSouth (const DetId &id) const
virtual DetId	goUp (const DetId &id) const
virtual DetId	goWest (const DetId &id) const
virtual	~CaloSubdetectorTopology ()
	virtual destructor More...

Static Public Member Functions
static std::string	producerTag ()

Private Types
enum	{   kHBhalf = 1296, kHEhalf = 1296, kHOhalf = 1080, kHFhalf = 864,   kHThalf = 2088, kZDChalf = 11, kCASTORhalf = 224, kCALIBhalf = 693,   kHThalfPhase1 = 2520, kHcalhalf = kHBhalf + kHEhalf + kHOhalf + kHFhalf }
enum	{ kSizeForDenseIndexingPreLS1 = 2 * kHcalhalf }
enum	{ kHBSizePreLS1 = 2 * kHBhalf }
enum	{ kHESizePreLS1 = 2 * kHEhalf }
enum	{ kHOSizePreLS1 = 2 * kHOhalf }
enum	{ kHFSizePreLS1 = 2 * kHFhalf }
enum	{ kHTSizePreLS1 = 2 * kHThalf }
enum	{ kHTSizePhase1 = 2 * kHThalfPhase1 }
enum	{ kCALIBSizePreLS1 = 2 * kCALIBhalf }
typedef std::map< unsigned, std::vector< int > >	SegmentationMap

Private Member Functions
int	decAIEta (const HcalDetId &id, HcalDetId neighbors[2]) const
unsigned int	detId2denseIdPreLS1 (const DetId &id) const
int	incAIEta (const HcalDetId &id, HcalDetId neighbors[2]) const
bool	isExcluded (const HcalDetId &id) const
bool	validDetIdPreLS1 (const HcalDetId &id) const
bool	validRaw (const HcalDetId &id) const

Private Attributes
unsigned int	CALIBSize_
SegmentationMap	depthSegmentation_
SegmentationMap	depthSegmentationOne_
int	doublePhiBins_
std::vector< double >	dPhiTable
std::vector< double >	dPhiTableHF
std::vector < HcalDDDRecConstants::HcalEtaBin >	etaBinsHB_
std::vector < HcalDDDRecConstants::HcalEtaBin >	etaBinsHE_
int	etaHE2HF_
int	etaHF2HE_
std::vector< double >	etaTable
std::vector< double >	etaTableHF
bool	excludeHB_
bool	excludeHE_
bool	excludeHF_
bool	excludeHO_
std::vector< HcalDetId >	exclusionList_
int	firstHBRing_
int	firstHEDoublePhiRing_
int	firstHEQuadPhiRing_
int	firstHERing_
int	firstHETripleDepthRing_
int	firstHFQuadPhiRing_
int	firstHFRing_
int	firstHORing_
unsigned int	HBSize_
const HcalDDDRecConstants *	hcons_
unsigned int	HESize_
unsigned int	HFSize_
unsigned int	HOSize_
unsigned int	HTSize_
int	lastHBRing_
int	lastHERing_
int	lastHFRing_
int	lastHORing_
int	maxDepthHB_
int	maxDepthHE_
int	maxDepthHF_
int	maxEta_
int	maxPhiHE_
bool	mergePosition_
HcalTopologyMode::Mode	mode_
int	nEtaHB_
int	nEtaHE_
unsigned int	numberOfShapes_
std::vector< double >	phioff
int	singlePhiBins_
int	topoVersion_
HcalTopologyMode::TriggerMode	triggerMode_
std::vector< int >	unitPhi
std::vector< int >	unitPhiHF

Static Private Attributes
static constexpr int	chanCalibHB_ [nchanCalibHB_] = {0, 1, 2}
static constexpr int	chanCalibHE1_ [nchanCalibHE1_] = {0, 1, 3, 4, 5, 6}
static constexpr int	chanCalibHE2_ = 2
static constexpr int	chanCalibHF1_ [nchanCalibHF1_] = {0, 1, 8}
static constexpr int	chanCalibHF2_ = 9
static constexpr int	chanCalibHO_ [nchanCalibHO_] = {0, 1}
static constexpr int	chanCalibHOs_ = 7
static constexpr int	etaCalibHB_ [nEtaCalibHB_] = {-1, 1}
static constexpr int	etaCalibHBX_ = 16
static constexpr int	etaCalibHE_ [nEtaCalibHE_] = {-1, 1}
static constexpr int	etaCalibHEX_ [nEtaCalibHEX_] = {25, 27}
static constexpr int	etaCalibHF_ [nEtaCalibHF_] = {-1, 1}
static constexpr int	etaCalibHO_ [nEtaCalibHO_] = {-2, -1, 0, 1, 2}
static constexpr int	etaCalibHOX_ [nEtaCalibHOX_] = {4, 15}
static constexpr unsigned int	kchanCalibHB_ = nchanCalibHB_ * kPhiCalibHB_
static constexpr unsigned int	kchanCalibHE1_ = nchanCalibHE1_ * kPhiCalibHE_
static constexpr unsigned int	kchanCalibHE2_ = kPhiCalibHE_
static constexpr unsigned int	kchanCalibHF1_ = nchanCalibHF1_ * kPhiCalibHF1_
static constexpr unsigned int	kchanCalibHF2_ = 1
static constexpr unsigned int	kOffCalibHB_ = 0
static constexpr unsigned int	kOffCalibHBX_ = kOffCalibHOX_ + nCalibHOX_
static constexpr unsigned int	kOffCalibHE1_ = kOffCalibHB_ + nCalibHB_
static constexpr unsigned int	kOffCalibHE2_ = kOffCalibHO2_ + nCalibHO2_
static constexpr unsigned int	kOffCalibHEX_ = kOffCalibHBX_ + nCalibHBX_
static constexpr unsigned int	kOffCalibHF1_ = kOffCalibHE1_ + nCalibHE1_
static constexpr unsigned int	kOffCalibHF2_ = kOffCalibHE2_ + nCalibHE2_
static constexpr unsigned int	kOffCalibHFX_ = kOffCalibHEX_ + nCalibHEX_
static constexpr unsigned int	kOffCalibHO1_ = kOffCalibHF1_ + nCalibHF1_
static constexpr unsigned int	kOffCalibHO2_ = kOffCalibHO1_ + nCalibHO1_
static constexpr unsigned int	kOffCalibHOX_ = kOffCalibHF2_ + nCalibHF2_
static constexpr unsigned int	kPhiCalibHB_ = maxPhi_ / mPhiCalibHB_
static constexpr unsigned int	kPhiCalibHBX_ = maxPhi_ / mPhiCalibHBX_
static constexpr unsigned int	kPhiCalibHE_ = maxPhi_ / mPhiCalibHE_
static constexpr unsigned int	kPhiCalibHEX_ = maxPhi_ / mPhiCalibHEX_
static constexpr unsigned int	kPhiCalibHF1_ = maxPhi_ / mPhiCalibHF1_
static constexpr unsigned int	kPhiCalibHO0_ = maxPhi_ / mPhiCalibHO0_
static constexpr unsigned int	kPhiCalibHO1_ = maxPhi_ / mPhiCalibHO1_
static constexpr unsigned int	kPhiCalibHO2_ = 4 * kPhiCalibHO1_ + kPhiCalibHO0_
static constexpr unsigned int	maxPhi_ = 72
static constexpr int	minMaxDepth_ = 4
static constexpr unsigned int	minPhi_ = 1
static constexpr unsigned int	mPhiCalibHB_ = 4
static constexpr unsigned int	mPhiCalibHBX_ = 1
static constexpr unsigned int	mPhiCalibHE_ = 4
static constexpr unsigned int	mPhiCalibHEX_ = 2
static constexpr unsigned int	mPhiCalibHF1_ = 18
static constexpr unsigned int	mPhiCalibHO0_ = 6
static constexpr unsigned int	mPhiCalibHO1_ = 12
static constexpr unsigned int	mPhiCalibHOX_ [nEtaCalibHOX_] = {2, 1}
static constexpr unsigned int	nCalibHB_ = kchanCalibHB_ * nEtaCalibHB_
static constexpr unsigned int	nCalibHBX_ = 2 * kPhiCalibHBX_ * nEtaCalibHBX_
static constexpr unsigned int	nCalibHE1_ = kchanCalibHE1_ * nEtaCalibHE_
static constexpr unsigned int	nCalibHE2_ = kchanCalibHE2_ * nEtaCalibHE_
static constexpr unsigned int	nCalibHEX_ = 2 * kPhiCalibHEX_ * nEtaCalibHEX_
static constexpr unsigned int	nCalibHF1_ = kchanCalibHF1_ * nEtaCalibHF_
static constexpr unsigned int	nCalibHF2_ = nEtaCalibHF_
static constexpr unsigned int	nCalibHO1_ = kPhiCalibHO2_ * nchanCalibHO_
static constexpr unsigned int	nCalibHO2_ = nEtaCalibHO_
static constexpr int	nCalibHOX_ = 2 * (nPhiCalibHOX_[0] + nPhiCalibHOX_[1])
static constexpr unsigned int	nchanCalibHB_ = 3
static constexpr unsigned int	nchanCalibHE1_ = 6
static constexpr unsigned int	nchanCalibHF1_ = 3
static constexpr unsigned int	nchanCalibHO_ = 2
static constexpr unsigned int	nEtaCalibHB_ = 2
static constexpr unsigned int	nEtaCalibHBX_ = 1
static constexpr unsigned int	nEtaCalibHE_ = 2
static constexpr unsigned int	nEtaCalibHEX_ = 2
static constexpr unsigned int	nEtaCalibHF_ = 2
static constexpr unsigned int	nEtaCalibHO_ = 5
static constexpr unsigned int	nEtaCalibHOX_ = 2
static constexpr unsigned int	nPhiCalibHOX_ [nEtaCalibHOX_] = {36, 72}
static constexpr int	phiCalibHF2_ = 1
static constexpr int	phiCalibHO_ [nEtaCalibHO_] = {59, 47, 53, 47, 47}
static constexpr int	phiCalibHOX1_ = 4
static constexpr int	phiCalibHOX2_ = 12
static constexpr int	phiCalibHOX3_ = 6

Additional Inherited Members
Protected Types inherited from CaloSubdetectorTopology
typedef std::pair< int, int >	Coordinate
Protected Member Functions inherited from CaloSubdetectorTopology
Coordinate	getNeighbourIndex (const Coordinate &coord, const CaloDirection &dir) const

Detailed Description

The HcalTopology class contains a set of constants which represent the topology (tower relationship) of the CMS HCAL as built. These constants can be used to determine neighbor relationships and existence of cells.

For use with limited setups (testbeam, cosmic stands, etc), the topology can be limited by creating a rejection list – a list of cells which would normally exist in the full CMS HCAL, but are not present for the specified topology.

Definition at line 26 of file HcalTopology.h.

Member Typedef Documentation

typedef std::map<unsigned, std::vector<int> > HcalTopology::SegmentationMap

private

Definition at line 223 of file HcalTopology.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
kHBhalf
kHEhalf
kHOhalf
kHFhalf
kHThalf
kZDChalf
kCASTORhalf
kCALIBhalf
kHThalfPhase1
kHcalhalf

Definition at line 227 of file HcalTopology.h.

       {
    kHBhalf = 1296,
    kHEhalf = 1296,
    kHOhalf = 1080,
    kHFhalf = 864,
    kHThalf = 2088,
    kZDChalf = 11,
    kCASTORhalf = 224,
    kCALIBhalf = 693,
    kHThalfPhase1 = 2520,
    kHcalhalf = kHBhalf + kHEhalf + kHOhalf + kHFhalf
  };

anonymous enum

private

Enumerator
kSizeForDenseIndexingPreLS1

Definition at line 239 of file HcalTopology.h.

{ kSizeForDenseIndexingPreLS1 = 2 * kHcalhalf };

anonymous enum

private

Enumerator
kHBSizePreLS1

Definition at line 240 of file HcalTopology.h.

{ kHBSizePreLS1 = 2 * kHBhalf };

anonymous enum

private

Enumerator
kHESizePreLS1

Definition at line 241 of file HcalTopology.h.

{ kHESizePreLS1 = 2 * kHEhalf };

anonymous enum

private

Enumerator
kHOSizePreLS1

Definition at line 242 of file HcalTopology.h.

{ kHOSizePreLS1 = 2 * kHOhalf };

anonymous enum

private

Enumerator
kHFSizePreLS1

Definition at line 243 of file HcalTopology.h.

{ kHFSizePreLS1 = 2 * kHFhalf };

anonymous enum

private

Enumerator
kHTSizePreLS1

Definition at line 244 of file HcalTopology.h.

{ kHTSizePreLS1 = 2 * kHThalf };

anonymous enum

private

Enumerator
kHTSizePhase1

Definition at line 245 of file HcalTopology.h.

{ kHTSizePhase1 = 2 * kHThalfPhase1 };

anonymous enum

private

Enumerator
kCALIBSizePreLS1

Definition at line 246 of file HcalTopology.h.

{ kCALIBSizePreLS1 = 2 * kCALIBhalf };

Constructor & Destructor Documentation

HcalTopology::HcalTopology	(	const HcalDDDRecConstants *	hcons,
		const bool	mergePosition = false
	)

Definition at line 17 of file HcalTopology.cc.

References CALIBSize_, dPhiTable, dPhiTableHF, PVValHelper::eta, etaBinsHB_, etaBinsHE_, etaHE2HF_, etaHF2HE_, etaTable, etaTableHF, first, firstHBRing_, firstHEDoublePhiRing_, firstHEQuadPhiRing_, firstHERing_, firstHETripleDepthRing_, firstHFQuadPhiRing_, firstHFRing_, firstHORing_, HcalDDDRecConstants::getDepth(), HcalDDDRecConstants::getEtaBins(), HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getEtaTable(), HcalDDDRecConstants::getEtaTableHF(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNEta(), HcalDDDRecConstants::getNPhi(), HcalDDDRecConstants::getPhiOffs(), HcalDDDRecConstants::getPhiTable(), HcalDDDRecConstants::getPhiTableHF(), HcalDDDRecConstants::getTopoMode(), HcalDDDRecConstants::getTriggerMode(), HBSize_, hcons_, HESize_, HFSize_, HOSize_, HTSize_, mps_fire::i, IPHI_MAX, isotrackApplyRegressor::k, kCALIBSizePreLS1, kHBSizePreLS1, kHESizePreLS1, kHFSizePreLS1, kHOSizePreLS1, kHTSizePhase1, kHTSizePreLS1, kOffCalibHFX_, lastHBRing_, lastHERing_, lastHFRing_, lastHORing_, HcalTopologyMode::LHC, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxEta_, maxPhiHE_, mode_, HLT_FULL_cff::nEta, nEtaHB_, nEtaHE_, numberOfShapes_, phioff, relativeConstraints::ring, edm::second(), setDepthSegmentation(), HcalTopologyMode::SLHC, topoVersion_, triggerMode_, HcalTopologyMode::TriggerMode_2009, unit(), unitPhi, unitPhiHF, and units().

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

    : hcons_(hcons),
      mergePosition_(mergePosition),
      excludeHB_(false),
      excludeHE_(false),
      excludeHO_(false),
      excludeHF_(false),
      firstHBRing_(1),
      firstHERing_(999),
      lastHERing_(0),
      firstHFRing_(29),
      lastHFRing_(41),
      firstHORing_(1),
      lastHORing_(15),
      firstHEDoublePhiRing_(999),
      firstHEQuadPhiRing_(999),
      firstHFQuadPhiRing_(40),
      firstHETripleDepthRing_(999),
      singlePhiBins_(IPHI_MAX),
      doublePhiBins_(36),
      maxPhiHE_(IPHI_MAX) {
  mode_ = (HcalTopologyMode::Mode)(hcons_->getTopoMode());
  triggerMode_ = (HcalTopologyMode::TriggerMode)(hcons_->getTriggerMode());
  maxDepthHB_ = hcons_->getMaxDepth(0);
  maxDepthHE_ = hcons_->getMaxDepth(1);
  maxDepthHF_ = hcons_->getMaxDepth(2);
  etaBinsHB_ = hcons_->getEtaBins(0);
  etaBinsHE_ = hcons_->getEtaBins(1);
  nEtaHB_ = (hcons_->getEtaRange(0)).second - (hcons_->getEtaRange(0)).first + 1;
  lastHBRing_ = firstHBRing_ + nEtaHB_ - 1;
  if (hcons_->getNPhi(1) > maxPhiHE_)
    maxPhiHE_ = hcons_->getNPhi(1);
  for (auto& i : etaBinsHE_) {
    if (firstHERing_ > i.ieta)
      firstHERing_ = i.ieta;
    if (lastHERing_ < i.ieta)
      lastHERing_ = i.ieta;
    int unit = static_cast<int>((i.dphi / 5.0_deg) + 0.01);
    if (unit == 2 && firstHEDoublePhiRing_ > i.ieta)
      firstHEDoublePhiRing_ = i.ieta;
    if (unit == 4 && firstHEQuadPhiRing_ > i.ieta)
      firstHEQuadPhiRing_ = i.ieta;
    if (i.layer.size() > 2 && firstHETripleDepthRing_ > i.ieta)
      firstHETripleDepthRing_ = i.ieta;
  }
  if (firstHERing_ > lastHERing_) {
    firstHERing_ = lastHERing_ = firstHEDoublePhiRing_ = firstHEQuadPhiRing_ = firstHETripleDepthRing_ = nEtaHE_ = 0;
  } else {
    nEtaHE_ = (lastHERing_ - firstHERing_ + 1);
  }
  if (mode_ == HcalTopologyMode::LHC) {
    topoVersion_ = 0;         //DL
    HBSize_ = kHBSizePreLS1;  // qie-per-fiber * fiber/rm * rm/rbx * rbx/barrel * barrel/hcal
    HESize_ = kHESizePreLS1;  // qie-per-fiber * fiber/rm * rm/rbx * rbx/endcap * endcap/hcal
    HOSize_ = kHOSizePreLS1;  // ieta * iphi * 2
    HFSize_ = kHFSizePreLS1;  // ieta * iphi * depth * 2
    CALIBSize_ = kCALIBSizePreLS1;
    numberOfShapes_ = 87;
  } else if (mode_ == HcalTopologyMode::SLHC) {  // need to know more eventually
    topoVersion_ = 10;
    HBSize_ = nEtaHB_ * IPHI_MAX * maxDepthHB_ * 2;
    HESize_ = nEtaHE_ * maxPhiHE_ * maxDepthHE_ * 2;
    HOSize_ = (lastHORing_ - firstHORing_ + 1) * IPHI_MAX * 2;                // ieta * iphi * 2
    HFSize_ = (lastHFRing_ - firstHFRing_ + 1) * IPHI_MAX * maxDepthHF_ * 2;  // ieta * iphi * depth * 2
    CALIBSize_ = kOffCalibHFX_;
    numberOfShapes_ = (maxPhiHE_ > 72) ? 1200 : 500;
  }
  maxEta_ = (lastHERing_ > lastHFRing_) ? lastHERing_ : lastHFRing_;
  if (triggerMode_ == HcalTopologyMode::TriggerMode_2009) {
    HTSize_ = kHTSizePreLS1;
  } else {
    HTSize_ = kHTSizePhase1;
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "Topo sizes " << HBSize_ << ":" << HESize_ << ":" << HOSize_ << ":" << HFSize_ << ":"
                               << HTSize_ << ":" << CALIBSize_ << " for mode " << mode_ << ":" << triggerMode_;
#endif

  //The transition between HE/HF in eta
  etaTableHF = hcons_->getEtaTableHF();
  etaTable = hcons_->getEtaTable();
  dPhiTableHF = hcons_->getPhiTableHF();
  dPhiTable = hcons_->getPhiTable();
  phioff = hcons_->getPhiOffs();
  std::pair<int, int> ietaHF = hcons_->getEtaRange(2);
  etaHE2HF_ = firstHFRing_;
  etaHF2HE_ = lastHERing_;
  if (etaBinsHE_.size() > 1) {
    double eta = etaBinsHE_[etaBinsHE_.size() - 1].etaMax;
    for (unsigned int i = 1; i < etaTableHF.size(); ++i) {
      if (eta < etaTableHF[i]) {
        etaHE2HF_ = ietaHF.first + i - 1;
        break;
      }
    }
    eta = etaTableHF[0];
    for (auto& i : etaBinsHE_) {
      if (eta < i.etaMax) {
        etaHF2HE_ = i.ieta;
        break;
      }
    }
  }
  const double fiveDegInRad = 5.0_deg;
  for (double k : dPhiTable) {
    int units = (int)(k / fiveDegInRad + 0.5);
    unitPhi.emplace_back(units);
  }
  for (double k : dPhiTableHF) {
    int units = (int)(k / fiveDegInRad + 0.5);
    unitPhiHF.emplace_back(units);
  }
  int nEta = hcons_->getNEta();
  for (int ring = 1; ring <= nEta; ++ring) {
    std::vector<int> segmentation = hcons_->getDepth(ring - 1, false);
    setDepthSegmentation(ring, segmentation, false);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "Set segmentation for ring " << ring << " with " << segmentation.size()
                                 << " elements:";
    for (unsigned int k = 0; k < segmentation.size(); ++k)
      edm::LogVerbatim("HCalGeom") << "[" << k << "] " << segmentation[k];
#endif
    segmentation = hcons_->getDepth(ring - 1, true);
    setDepthSegmentation(ring, segmentation, true);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "Set Plan-1 segmentation for ring " << ring << " with " << segmentation.size()
                                 << " elements:";
    for (unsigned int k = 0; k < segmentation.size(); ++k)
      edm::LogVerbatim("HCalGeom") << "[" << k << "] " << segmentation[k];
#endif
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "Constants in HcalTopology " << firstHBRing_ << ":" << lastHBRing_ << " "
                               << firstHERing_ << ":" << lastHERing_ << ":" << firstHEDoublePhiRing_ << ":"
                               << firstHEQuadPhiRing_ << ":" << firstHETripleDepthRing_ << " " << firstHFRing_ << ":"
                               << lastHFRing_ << ":" << firstHFQuadPhiRing_ << " " << firstHORing_ << ":" << lastHORing_
                               << " " << maxDepthHB_ << ":" << maxDepthHE_ << " " << nEtaHB_ << ":" << nEtaHE_ << " "
                               << etaHE2HF_ << ":" << etaHF2HE_ << " " << maxPhiHE_;
#endif
}

HcalTopology::HcalTopology	(	HcalTopologyMode::Mode	mode,
		int	maxDepthHB,
		int	maxDepthHE,
		HcalTopologyMode::TriggerMode	tmode = HcalTopologyMode::TriggerMode_2009
	)

Definition at line 160 of file HcalTopology.cc.

References CALIBSize_, firstHBRing_, firstHERing_, HBSize_, HESize_, HFSize_, HOSize_, HTSize_, IPHI_MAX, kHBSizePreLS1, kHESizePreLS1, kHFSizePreLS1, kHOSizePreLS1, kHTSizePhase1, kHTSizePreLS1, kOffCalibHFX_, lastHBRing_, lastHERing_, HcalTopologyMode::LHC, maxDepthHF_, maxPhiHE_, mode_, nEtaHB_, nEtaHE_, HcalTopologyMode::SLHC, topoVersion_, triggerMode_, and HcalTopologyMode::TriggerMode_2009.

    : hcons_(nullptr),
      mergePosition_(false),
      excludeHB_(false),
      excludeHE_(false),
      excludeHO_(false),
      excludeHF_(false),
      mode_(mode),
      triggerMode_(tmode),
      firstHBRing_(1),
      lastHBRing_(16),
      firstHERing_(16),
      lastHERing_(29),
      firstHFRing_(29),
      lastHFRing_(41),
      firstHORing_(1),
      lastHORing_(15),
      firstHEDoublePhiRing_((mode == HcalTopologyMode::H2 || mode == HcalTopologyMode::H2HE) ? (22) : (21)),
      firstHEQuadPhiRing_(999),
      firstHFQuadPhiRing_(40),
      firstHETripleDepthRing_((mode == HcalTopologyMode::H2 || mode == HcalTopologyMode::H2HE) ? (24) : (27)),
      singlePhiBins_(IPHI_MAX),
      doublePhiBins_(36),
      maxDepthHB_(maxDepthHB),
      maxDepthHE_(maxDepthHE),
      maxDepthHF_(2),
      etaHE2HF_(30),
      etaHF2HE_(29),
      maxPhiHE_(IPHI_MAX),
      HBSize_(kHBSizePreLS1),
      HESize_(kHESizePreLS1),
      HOSize_(kHOSizePreLS1),
      HFSize_(kHFSizePreLS1),
      HTSize_(kHTSizePreLS1),
      CALIBSize_(kCALIBSizePreLS1),
      numberOfShapes_((mode == HcalTopologyMode::SLHC) ? 500 : 87) {
  if (mode_ == HcalTopologyMode::LHC) {
    topoVersion_ = 0;                            //DL
    HBSize_ = kHBSizePreLS1;                     // qie-per-fiber * fiber/rm * rm/rbx * rbx/barrel * barrel/hcal
    HESize_ = kHESizePreLS1;                     // qie-per-fiber * fiber/rm * rm/rbx * rbx/endcap * endcap/hcal
    HOSize_ = kHOSizePreLS1;                     // ieta * iphi * 2
    HFSize_ = kHFSizePreLS1;                     // phi * eta * depth * pm
  } else if (mode_ == HcalTopologyMode::SLHC) {  // need to know more eventually
    HBSize_ = maxDepthHB * 16 * IPHI_MAX * 2;
    HESize_ = maxDepthHE * (29 - 16 + 1) * maxPhiHE_ * 2;
    HOSize_ = 15 * IPHI_MAX * 2;                // ieta * iphi * 2
    HFSize_ = IPHI_MAX * 13 * maxDepthHF_ * 2;  // phi * eta * depth * pm
    CALIBSize_ = kOffCalibHFX_;
    topoVersion_ = 10;
  }
  nEtaHB_ = (lastHBRing_ - firstHBRing_ + 1);
  nEtaHE_ = (lastHERing_ - firstHERing_ + 1);
  if (triggerMode_ == HcalTopologyMode::TriggerMode_2009) {
    HTSize_ = kHTSizePreLS1;
  } else {
    HTSize_ = kHTSizePhase1;
  }

  edm::LogWarning("HCalGeom") << "This is an incomplete constructor of HcalTopology - be warned that many "
                              << "functionalities will not be there - revert from this - get from EventSetup";
}

Member Function Documentation

const HcalDDDRecConstants* HcalTopology::dddConstants ( ) const

inline

Definition at line 164 of file HcalTopology.h.

References hcons_.

Referenced by HcalRaddamMuon::analyze(), HcalHBHEMuonAnalyzer::depth16HE(), HcalHBHEMuonHighEtaAnalyzer::depth16HE(), AlCaHcalHBHEMuonProducer::depth16HE(), CaloParticleDebugger::fillSimHits(), HcalGeometry::getBackPosition(), HcalGeometry::getClosestCell(), HcalGeometry::getPosition(), CaloTruthAccumulator::initializeEvent(), CaloTPGTranscoderULUT::setup(), and HcaluLUTTPGCoder::update().

{ return hcons_; }

int HcalTopology::decAIEta ( const HcalDetId &  id, HcalDetId  neighbors[2]  ) const

private

Get the neighbors of the given cell with lower absolute ieta

Decreasing in |ieta|, there are two neighbors of 40 and 21

Definition at line 770 of file HcalTopology.cc.

References HLT_FULL_cff::depth, etaHF2HE_, firstHEDoublePhiRing(), firstHEQuadPhiRing(), firstHERing(), firstHFQuadPhiRing(), firstHFRing(), HcalBarrel, HcalEndcap, HcalForward, IPHI_MAX, dqmiodumpmetadata::n, valid(), and ecaldqm::zside().

Referenced by decIEta(), and incIEta().

                                                                            {
  int n = 1;
  int aieta = id.ietaAbs();

  if (aieta == firstHEDoublePhiRing()) {
    n = 2;
    neighbors[0] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi(), id.depth());
    neighbors[1] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi() + 1, id.depth());
  } else if (aieta == firstHFQuadPhiRing()) {
    n = 2;
    neighbors[0] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi(), id.depth());
    if (id.iphi() == IPHI_MAX - 1)
      neighbors[1] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), 1, id.depth());
    else
      neighbors[1] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi() + 2, id.depth());
  } else if (aieta == firstHEQuadPhiRing()) {
    n = 2;
    neighbors[0] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi(), id.depth());
    if (id.iphi() == IPHI_MAX - 1)
      neighbors[1] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), 1, id.depth());
    else
      neighbors[1] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi() + 2, id.depth());
  } else if (aieta == 1) {
    neighbors[0] = HcalDetId(id.subdet(), -aieta * id.zside(), id.iphi(), id.depth());
  } else if (aieta == firstHERing() && id.subdet() == HcalEndcap) {
    neighbors[0] = HcalDetId(HcalBarrel, (aieta - 1) * id.zside(), id.iphi(), 1);
  } else if (aieta == firstHFRing() && id.subdet() == HcalForward) {
    neighbors[0] = HcalDetId(HcalEndcap, etaHF2HE_ * id.zside(), id.iphi(), 1);
  } else
    neighbors[0] = HcalDetId(id.subdet(), (aieta - 1) * id.zside(), id.iphi(), id.depth());

  if (!valid(neighbors[0]) && n == 2) {
    if (!valid(neighbors[1]))
      n = 0;
    else {
      n = 1;
      neighbors[0] = neighbors[1];
    }
  }
  if (n == 2 && !valid(neighbors[1]))
    n = 1;
  if (n == 1 && !valid(neighbors[0]))
    n = 0;

  return n;
}

int HcalTopology::decIEta	(	const HcalDetId &	id,
		HcalDetId	neighbors[2]
	)		const

Get the neighbors of the given cell with lower (signed) ieta

Definition at line 737 of file HcalTopology.cc.

References decAIEta(), incAIEta(), and ecaldqm::zside().

Referenced by east().

                                                                           {
  if (id.zside() == 1)
    return decAIEta(id, neighbors);
  else
    return incAIEta(id, neighbors);
}

bool HcalTopology::decIPhi	(	const HcalDetId &	id,
		HcalDetId &	neighbor
	)		const

Get the neighbor (if present) of the given cell with lower iphi

Definition at line 679 of file HcalTopology.cc.

References HLT_FULL_cff::depth, firstHEDoublePhiRing(), firstHEQuadPhiRing(), firstHFQuadPhiRing(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, IPHI_MAX, maxPhiHE_, convertSQLiteXML::ok, valid(), and validRaw().

Referenced by HFStripFilter::runFilter(), and south().

                                                                         {
  bool ok = valid(id);
  if (ok) {
    switch (id.subdet()) {
      case (HcalBarrel):
      case (HcalOuter):
        if (id.iphi() == 1)
          neighbor = HcalDetId(id.subdet(), id.ieta(), IPHI_MAX, id.depth());
        else
          neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() - 1, id.depth());
        break;
      case (HcalEndcap):
        if (id.ietaAbs() >= firstHEQuadPhiRing()) {
          if (id.iphi() == 3)
            neighbor = HcalDetId(id.subdet(), id.ieta(), IPHI_MAX - 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() - 4, id.depth());
        } else if (id.ietaAbs() >= firstHEDoublePhiRing()) {
          if (id.iphi() == 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), IPHI_MAX - 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() - 2, id.depth());
        } else {
          if (id.iphi() == 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), maxPhiHE_, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() - 1, id.depth());
        }
        break;
      case (HcalForward):
        if (id.ietaAbs() >= firstHFQuadPhiRing()) {
          if (id.iphi() == 3)
            neighbor = HcalDetId(id.subdet(), id.ieta(), IPHI_MAX - 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() - 4, id.depth());
        } else {
          if (id.iphi() == 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), IPHI_MAX - 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() - 2, id.depth());
        }
        if (!validRaw(neighbor))
          ok = false;
        break;
      default:
        ok = false;
    }
  }
  return ok;
}

bool HcalTopology::decrementDepth

(

HcalDetId &

id

)

const

Get the detector in front of this one

Definition at line 911 of file HcalTopology.cc.

References HcalDetId::depth(), HLT_FULL_cff::depth, depthBinInformation(), etaBinsHB_, etaBinsHE_, etaHF2HE_, etaRing(), firstHERing(), firstHFRing(), HcalDDDRecConstants::getDepthEta29(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons_, mps_fire::i, HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), lastHERing(), mode_, nEtaHB_, HcalTopologyMode::SLHC, HcalDetId::subdet(), validRaw(), HcalDetId::zside(), and ecaldqm::zside().

Referenced by down().

                                                        {
  HcalSubdetector subdet = detId.subdet();
  int ieta = detId.ieta();
  int etaRing = detId.ietaAbs();
  int depth = detId.depth();
  int iphi = detId.iphi();
  int zside = detId.zside();
  int nDepthBins, startingBin;
  depthBinInformation(subdet, etaRing, iphi, zside, nDepthBins, startingBin);

  // see if the new depth bin exists
  --depth;
  if ((subdet == HcalOuter) || (subdet == HcalEndcap && etaRing == firstHERing())) {
    subdet = HcalBarrel;
    for (int i = 0; i < nEtaHB_; ++i) {
      if (etaRing == etaBinsHB_[i].ieta) {
        depth = etaBinsHB_[i].depthStart + etaBinsHB_[i].layer.size() - 1;
        break;
      }
    }
  } else if (subdet == HcalEndcap && etaRing == lastHERing() && depth == hcons_->getDepthEta29(iphi, zside, 0) &&
             mode_ != HcalTopologyMode::SLHC) {
    (ieta > 0) ? --ieta : ++ieta;
  } else if (depth <= 0) {
    if (subdet == HcalForward && etaRing == firstHFRing()) {
      // overlap
      subdet = HcalEndcap;
      etaRing = etaHF2HE_;
      ieta = (ieta > 0) ? etaRing : -etaRing;
      for (const auto& i : etaBinsHE_) {
        if (etaRing == i.ieta) {
          depth = i.depthStart + i.layer.size() - 1;
          break;
        }
      }
    } else {
      // no more chances
      detId = HcalDetId();
      return false;
    }
  }
  detId = HcalDetId(subdet, ieta, detId.iphi(), depth);
  return validRaw(detId);
}

DetId HcalTopology::denseId2detId ( unsigned int  denseid ) const

overridevirtual

return a linear packed id

Reimplemented from CaloSubdetectorTopology.

Definition at line 1562 of file HcalTopology.cc.

References TauDecayModes::dec, firstHBRing(), firstHERing(), HBSize_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HESize_, HOSize_, recoMuon::in, IPHI_MAX, kHBhalf, kHEhalf, kHFhalf, kHOhalf, kSizeForDenseIndexingPreLS1, lastHBRing(), lastHERing(), maxDepthHB_, maxDepthHE_, maxDepthHF_, maxPhiHE_, ncells(), nEtaHB_, nEtaHE_, DetId::rawId(), sd, and topoVersion_.

Referenced by HcalGeometry::localCorners().

                                                            {
  HcalSubdetector sd(HcalBarrel);
  int ie(0);
  int ip(0);
  int dp(0);
  int in(denseid);
  int iz(1);
  if (topoVersion_ == 0) {  //DL// pre-LS1
    if (denseid < kSizeForDenseIndexingPreLS1) {
      if (in > 2 * (kHBhalf + kHEhalf + kHOhalf) - 1) {  // HF
        sd = HcalForward;
        in -= 2 * (kHBhalf + kHEhalf + kHOhalf);
        iz = (in < kHFhalf ? 1 : -1);
        in %= kHFhalf;
        ip = 4 * (in / 48);
        in %= 48;
        ip += 1 + (in > 21 ? 2 : 0);
        if (3 == ip % 4)
          in -= 22;
        ie = 29 + in / 2;
        dp = 1 + in % 2;
      } else if (in > 2 * (kHBhalf + kHEhalf) - 1) {  // HO
        sd = HcalOuter;
        in -= 2 * (kHBhalf + kHEhalf);
        iz = (in < kHOhalf ? 1 : -1);
        in %= kHOhalf;
        dp = 4;
        ip = 1 + in / 15;
        ie = 1 + (in - 15 * (ip - 1));
      } else if (in > 2 * kHBhalf - 1) {  // Endcap
        sd = HcalEndcap;
        in -= 2 * kHBhalf;
        iz = (in < kHEhalf ? 1 : -1);
        in %= kHEhalf;
        ip = 2 * (in / 36);
        in %= 36;
        ip += 1 + in / 28;
        if (0 == ip % 2)
          in %= 28;
        ie = 15 + (in < 2 ? 1 + in : 2 + (in < 20 ? 1 + (in - 2) / 2 : 9 + (in < 26 ? 1 + (in - 20) / 3 : 3)));
        dp = (in < 1
                  ? 3
                  : (in < 2 ? 1 : (in < 20 ? 1 + (in - 2) % 2 : (in < 26 ? 1 + (in - 20) % 3 : (1 + (in - 26) % 2)))));
      } else {  // barrel
        iz = (in < kHBhalf ? 1 : -1);
        in %= kHBhalf;
        ip = in / 18 + 1;
        in %= 18;
        if (in < 14) {
          dp = 1;
          ie = in + 1;
        } else {
          in %= 14;
          dp = 1 + in % 2;
          ie = 15 + in / 2;
        }
      }
    }
  } else if (topoVersion_ == 10) {
    if (denseid < ncells()) {
      if (denseid >= (HBSize_ + HESize_ + HOSize_)) {
        sd = HcalForward;
        in -= (HBSize_ + HESize_ + HOSize_);
        dp = (in % maxDepthHF_) + 1;
        ip = (in - dp + 1) % (maxDepthHF_ * IPHI_MAX);
        ip = (ip / maxDepthHF_) + 1;
        ie = (in - dp + 1 - maxDepthHF_ * (ip - 1)) / (IPHI_MAX * maxDepthHF_);
        if (ie > 12) {
          ie = 54 - ie;
          iz = -1;
        } else {
          ie += 29;
          iz = 1;
        }
      } else if (denseid >= (HBSize_ + HESize_)) {
        sd = HcalOuter;
        in -= (HBSize_ + HESize_);
        dp = 4;
        ip = (in % IPHI_MAX) + 1;
        ie = (in - ip + 1) / IPHI_MAX;
        if (ie > 14) {
          ie = 30 - ie;
          iz = -1;
        } else {
          ie += 1;
          iz = 1;
        }
      } else if (denseid >= (HBSize_)) {
        sd = HcalEndcap;
        in -= (HBSize_);
        dp = (in % maxDepthHE_) + 1;
        ip = (in - dp + 1) % (maxDepthHE_ * maxPhiHE_);
        ip = (ip / maxDepthHE_) + 1;
        ie = (in - dp + 1 - maxDepthHE_ * (ip - 1)) / (maxPhiHE_ * maxDepthHE_);
        if (ie >= nEtaHE_) {
          ie = lastHERing() + nEtaHE_ - ie;
          iz = -1;
        } else {
          ie = firstHERing() + ie;
          iz = 1;
        }
      } else {
        sd = HcalBarrel;
        dp = (in % maxDepthHB_) + 1;
        ip = (in - dp + 1) % (maxDepthHB_ * IPHI_MAX);
        ip = (ip / maxDepthHB_) + 1;
        ie = (in - dp + 1 - maxDepthHB_ * (ip - 1)) / (IPHI_MAX * maxDepthHB_);
        if (ie >= nEtaHB_) {
          ie = lastHBRing() + nEtaHB_ - ie;
          iz = -1;
        } else {
          ie = firstHBRing() + ie;
          iz = 1;
        }
      }
    }
  }
  HcalDetId hid(sd, iz * int(ie), ip, dp);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalGeom") << "Dens2Det " << topoVersion_ << " i/p " << std::hex << denseid << " : " << hid.rawId()
                               << std::dec << " | " << hid;
#endif
  return hid;
}

HcalCalibDetId HcalTopology::denseId2detIdCALIB

(

const unsigned int &

id

)

const

return a Calib DetId from linear packed id

Definition at line 1386 of file HcalTopology.cc.

References chanCalibHB_, chanCalibHE1_, chanCalibHE2_, chanCalibHF1_, chanCalibHF2_, chanCalibHO_, chanCalibHOs_, etaCalibHB_, etaCalibHBX_, etaCalibHE_, etaCalibHEX_, etaCalibHF_, etaCalibHO_, etaCalibHOX_, HcalCalibDetId::HBX, HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalOuter, HcalCalibDetId::HEX, HcalCalibDetId::HOCrosstalk, gpuClustering::id, kchanCalibHB_, kchanCalibHE1_, kchanCalibHE2_, kchanCalibHF1_, kchanCalibHF2_, kOffCalibHBX_, kOffCalibHE1_, kOffCalibHE2_, kOffCalibHEX_, kOffCalibHF1_, kOffCalibHF2_, kOffCalibHFX_, kOffCalibHO1_, kOffCalibHO2_, kOffCalibHOX_, kPhiCalibHB_, kPhiCalibHBX_, kPhiCalibHE_, kPhiCalibHEX_, kPhiCalibHF1_, kPhiCalibHO0_, kPhiCalibHO1_, kPhiCalibHO2_, mPhiCalibHB_, mPhiCalibHE_, mPhiCalibHF1_, mPhiCalibHO0_, mPhiCalibHO1_, nPhiCalibHOX_, phiCalibHF2_, phiCalibHO_, phiCalibHOX1_, phiCalibHOX2_, and phiCalibHOX3_.

                                                                             {
  HcalCalibDetId id;
  unsigned int hid0(hid);
  if (hid0 < kOffCalibHOX_) {
    HcalSubdetector subdet(HcalEmpty);
    int id0, keta, ieta, iphi, ichan, ctype;
    if (hid0 < kOffCalibHE1_) {
      id0 = static_cast<int>(hid0);
      subdet = HcalBarrel;
      iphi = hid0 % kPhiCalibHB_;
      keta = (hid0 < kchanCalibHB_) ? 0 : 1;
      ieta = etaCalibHB_[keta];
      ichan = (id0 - iphi - keta * kchanCalibHB_) / kPhiCalibHB_;
      iphi = mPhiCalibHB_ * (iphi + 1) - 1;
      ctype = chanCalibHB_[ichan];
    } else if (hid0 < kOffCalibHF1_) {
      hid0 -= kOffCalibHE1_;
      id0 = static_cast<int>(hid0);
      subdet = HcalEndcap;
      iphi = hid0 % kPhiCalibHE_;
      keta = (hid0 < kchanCalibHE1_) ? 0 : 1;
      ieta = etaCalibHE_[keta];
      ichan = (id0 - iphi - keta * kchanCalibHE1_) / kPhiCalibHE_;
      iphi = mPhiCalibHE_ * (iphi + 1) - 1;
      ctype = chanCalibHE1_[ichan];
    } else if (hid0 < kOffCalibHO1_) {
      hid0 -= kOffCalibHF1_;
      id0 = static_cast<int>(hid0);
      subdet = HcalForward;
      iphi = hid0 % kPhiCalibHF1_;
      keta = (hid0 < kchanCalibHF1_) ? 0 : 1;
      ieta = etaCalibHF_[keta];
      ichan = (id0 - iphi - keta * kchanCalibHF1_) / kPhiCalibHF1_;
      iphi = mPhiCalibHF1_ * iphi + 1;
      ctype = chanCalibHF1_[ichan];
    } else if (hid0 < kOffCalibHO2_) {
      hid0 -= kOffCalibHO1_;
      id0 = static_cast<int>(hid0);
      subdet = HcalOuter;
      unsigned int kphi = hid0 % kPhiCalibHO2_;
      if (kphi < 2 * kPhiCalibHO1_) {
        ieta = (kphi >= kPhiCalibHO1_) ? etaCalibHO_[1] : etaCalibHO_[0];
        iphi = kphi % kPhiCalibHO1_;
        ichan = (id0 - iphi - (ieta + 2) * kPhiCalibHO1_) / kPhiCalibHO2_;
        iphi = (iphi + 1) * mPhiCalibHO1_ - 1;
      } else if (kphi < (2 * kPhiCalibHO1_ + kPhiCalibHO0_)) {
        ieta = etaCalibHO_[2];
        iphi = kphi % kPhiCalibHO0_;
        ichan = (id0 - iphi - (ieta + 2) * kPhiCalibHO1_) / kPhiCalibHO2_;
        iphi = (iphi + 1) * mPhiCalibHO0_ - 1;
      } else {
        ieta = (kphi >= 3 * kPhiCalibHO1_ + kPhiCalibHO0_) ? etaCalibHO_[4] : etaCalibHO_[3];
        iphi = kphi % kPhiCalibHO1_;
        ichan = (id0 - iphi - (ieta + 1) * kPhiCalibHO1_ - kPhiCalibHO0_) / kPhiCalibHO2_;
        iphi = (iphi + 1) * mPhiCalibHO1_ - 1;
      }
      ctype = chanCalibHO_[ichan];
    } else if (hid0 < kOffCalibHE2_) {
      hid0 -= kOffCalibHO2_;
      subdet = HcalOuter;
      iphi = phiCalibHO_[hid0];
      ctype = static_cast<int>(chanCalibHOs_);
      ieta = etaCalibHO_[hid0];
    } else if (hid0 < kOffCalibHF2_) {
      hid0 -= kOffCalibHE2_;
      id0 = static_cast<int>(hid0);
      subdet = HcalEndcap;
      iphi = hid0 % kPhiCalibHE_;
      keta = (hid0 < kchanCalibHE2_) ? 0 : 1;
      ieta = etaCalibHE_[keta];
      iphi = mPhiCalibHE_ * (iphi + 1) - 1;
      ctype = chanCalibHE2_;
    } else {
      hid0 -= kOffCalibHF2_;
      id0 = static_cast<int>(hid0);
      subdet = HcalForward;
      keta = (hid0 < kchanCalibHF2_) ? 0 : 1;
      ieta = etaCalibHF_[keta];
      iphi = phiCalibHF2_;
      ctype = chanCalibHF2_;
    }
    id = HcalCalibDetId(subdet, ieta, iphi, ctype);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "CalibrationBox: " << hid << " o/p " << ieta << ":" << iphi << ":" << ichan << ":"
                                 << ctype << " " << id;
#endif
  } else if (hid < kOffCalibHBX_) {
    hid0 -= kOffCalibHOX_;
    int iphi, ieta;
    if (hid0 < nPhiCalibHOX_[1]) {
      iphi = static_cast<int>(hid0) + 1;
      ieta = -etaCalibHOX_[1];
    } else if (hid0 < (nPhiCalibHOX_[1] + nPhiCalibHOX_[0])) {
      hid0 -= nPhiCalibHOX_[1];
      iphi = ((hid0 + phiCalibHOX1_) % phiCalibHOX2_ < phiCalibHOX3_) ? 2 * hid0 + 1 : 2 * hid0 + 2;
      ieta = -etaCalibHOX_[0];
    } else if (hid0 < (nPhiCalibHOX_[1] + 2 * nPhiCalibHOX_[0])) {
      hid0 -= (nPhiCalibHOX_[1] + nPhiCalibHOX_[0]);
      iphi = ((hid0 + phiCalibHOX1_) % phiCalibHOX2_ < phiCalibHOX3_) ? 2 * hid0 + 1 : 2 * hid0 + 2;
      ieta = etaCalibHOX_[0];
    } else {
      hid0 -= (nPhiCalibHOX_[1] + 2 * nPhiCalibHOX_[0]);
      iphi = static_cast<int>(hid0) + 1;
      ieta = etaCalibHOX_[1];
    }
    id = HcalCalibDetId(HcalCalibDetId::HOCrosstalk, ieta, iphi);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "HOCrossTalk: " << hid << ":" << hid0 << " o/p " << ieta << ":" << iphi << " "
                                 << id;
#endif
  } else if (hid < kOffCalibHEX_) {
    hid0 -= kOffCalibHBX_;
    int ieta = (hid0 >= kPhiCalibHBX_) ? etaCalibHBX_ : -etaCalibHBX_;
    int iphi = (hid0 % kPhiCalibHBX_) + 1;
    id = HcalCalibDetId(HcalCalibDetId::HBX, ieta, iphi);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "HBX: " << hid << ":" << hid0 << " o/p " << ieta << ":" << iphi << " " << id;
#endif
  } else if (hid < kOffCalibHFX_) {
    hid0 -= kOffCalibHEX_;
    int iphi = 2 * (hid0 % kPhiCalibHEX_) + 1;
    int ieta = ((hid0 < kPhiCalibHEX_)
                    ? -etaCalibHEX_[0]
                    : ((hid0 < 2 * kPhiCalibHEX_) ? etaCalibHEX_[0]
                                                  : ((hid0 < 3 * kPhiCalibHEX_) ? -etaCalibHEX_[1] : etaCalibHEX_[1])));
    id = HcalCalibDetId(HcalCalibDetId::HEX, ieta, iphi);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "HEX: " << hid << ":" << hid0 << " o/p " << ieta << ":" << iphi << " " << id;
#endif
  }
  return id;
}

void HcalTopology::depthBinInformation	(	HcalSubdetector	subdet,
		int	etaRing,
		int	iphi,
		int	zside,
		int &	nDepthBins,
		int &	startingBin
	)		const

finds the number of depth bins and which is the number to start with

Definition at line 817 of file HcalTopology.cc.

References firstHERing(), firstHETripleDepthRing(), HcalDDDRecConstants::getDepthEta16(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getMinDepth(), HcalTopologyMode::H2HE, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons_, lastHBRing(), lastHERing(), maxDepthHF_, mode_, and HcalTopologyMode::SLHC.

Referenced by CaloTowerConstituentsMap::constituentsOf(), decrementDepth(), HcalTrigTowerGeometry::detIds(), and incrementDepth().

                                                                                                       {
  if (subdet == HcalBarrel) {
    if (mode_ == HcalTopologyMode::SLHC || mode_ == HcalTopologyMode::H2HE) {
      startingBin = hcons_->getMinDepth(0, etaRing, iphi, zside);
      if (etaRing == lastHBRing()) {
        nDepthBins = hcons_->getDepthEta16(1, iphi, zside) - startingBin + 1;
      } else {
        nDepthBins = hcons_->getMaxDepth(0, etaRing, iphi, zside) - startingBin + 1;
      }
    } else {
      if (etaRing <= 14) {
        nDepthBins = 1;
        startingBin = 1;
      } else {
        nDepthBins = 2;
        startingBin = 1;
      }
    }
  } else if (subdet == HcalEndcap) {
    if (mode_ == HcalTopologyMode::SLHC || mode_ == HcalTopologyMode::H2HE) {
      if (etaRing == firstHERing()) {
        startingBin = hcons_->getDepthEta16(2, iphi, zside);
      } else {
        startingBin = hcons_->getMinDepth(1, etaRing, iphi, zside);
      }
      nDepthBins = hcons_->getMaxDepth(1, etaRing, iphi, zside) - startingBin + 1;
    } else {
      if (etaRing == firstHERing()) {
        nDepthBins = 1;
        startingBin = 3;
      } else if (etaRing == 17) {
        nDepthBins = 1;
        startingBin = 1;
      } else if (etaRing == lastHERing()) {
        nDepthBins = 2;
        startingBin = 1;
      } else {
        nDepthBins = (etaRing >= firstHETripleDepthRing()) ? 3 : 2;
        startingBin = 1;
      }
    }
  } else if (subdet == HcalForward) {
    nDepthBins = maxDepthHF_;
    startingBin = 1;
  } else if (subdet == HcalOuter) {
    nDepthBins = 1;
    startingBin = 4;
  } else {
    edm::LogWarning("HCalGeom") << "Bad HCAL subdetector " << subdet;
  }
}

unsigned int HcalTopology::detId2denseId ( const DetId &  id ) const

overridevirtual

return a linear packed id

Reimplemented from CaloSubdetectorTopology.

Definition at line 1519 of file HcalTopology.cc.

References TauDecayModes::dec, HcalDetId::depth(), detId2denseIdPreLS1(), firstHBRing(), firstHERing(), HBSize_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HESize_, HOSize_, HcalDetId::ieta(), HcalDetId::iphi(), IPHI_MAX, lastHBRing(), lastHERing(), maxDepthHB_, maxDepthHE_, maxDepthHF_, maxPhiHE_, nEtaHB_, nEtaHE_, HcalDetId::subdet(), and topoVersion_.

Referenced by HcalGeometry::getGeometry(), HcalGeometry::getGeometryBase(), HcalGeometry::indexFor(), HcalDDDGeometry::newCellImpl(), HcalGeometry::newCellImpl(), and HBHEPhase1Reconstructor::processData().

                                                              {
  unsigned int retval(0);
  if (topoVersion_ == 0) {  // pre-LS1
    retval = detId2denseIdPreLS1(id);
  } else if (topoVersion_ == 10) {
    HcalDetId hid(id);
    if (hid.subdet() == HcalBarrel) {
      retval = (hid.depth() - 1) + maxDepthHB_ * (hid.iphi() - 1);
      if (hid.ieta() > 0)
        retval += maxDepthHB_ * IPHI_MAX * (hid.ieta() - firstHBRing());
      else
        retval += maxDepthHB_ * IPHI_MAX * (hid.ieta() + lastHBRing() + nEtaHB_);
    } else if (hid.subdet() == HcalEndcap) {
      retval = HBSize_;
      retval += (hid.depth() - 1) + maxDepthHE_ * (hid.iphi() - 1);
      if (hid.ieta() > 0)
        retval += maxDepthHE_ * maxPhiHE_ * (hid.ieta() - firstHERing());
      else
        retval += maxDepthHE_ * maxPhiHE_ * (hid.ieta() + lastHERing() + nEtaHE_);
    } else if (hid.subdet() == HcalOuter) {
      retval = HBSize_ + HESize_;
      if (hid.ieta() > 0)
        retval += (hid.iphi() - 1) + IPHI_MAX * (hid.ieta() - 1);
      else
        retval += (hid.iphi() - 1) + IPHI_MAX * (30 + hid.ieta());
    } else if (hid.subdet() == HcalForward) {
      retval = HBSize_ + HESize_ + HOSize_;
      retval += (hid.depth() - 1) + maxDepthHF_ * (hid.iphi() - 1);
      if (hid.ieta() > 0)
        retval += maxDepthHF_ * IPHI_MAX * (hid.ieta() - 29);
      else
        retval += maxDepthHF_ * IPHI_MAX * ((41 + 13) + hid.ieta());
    } else {
      retval = 0xFFFFFFFu;
    }
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalGeom") << "DetId2Dense " << topoVersion_ << " ID " << std::hex << id.rawId() << std::dec
                               << " | " << HcalDetId(id) << " : " << std::hex << retval << std::dec;
#endif
  return retval;
}

unsigned int HcalTopology::detId2denseIdCALIB

(

const DetId &

id

)

const

return a linear packed id from CALIB

Definition at line 1279 of file HcalTopology.cc.

References funct::abs(), HcalCalibDetId::calibFlavor(), HcalCalibDetId::CalibrationBox, HcalCalibDetId::cboxChannel(), chanCalibHB_, chanCalibHE1_, chanCalibHE2_, chanCalibHF1_, chanCalibHF2_, chanCalibHO_, chanCalibHOs_, etaCalibHB_, etaCalibHBX_, etaCalibHE_, etaCalibHEX_, etaCalibHF_, etaCalibHO_, etaCalibHOX_, spr::find(), HcalCalibDetId::HBX, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalCalibDetId::hcalSubdet(), HcalCalibDetId::HEX, HcalCalibDetId::HOCrosstalk, HcalCalibDetId::ieta(), HcalCalibDetId::iphi(), kchanCalibHB_, kchanCalibHE1_, kchanCalibHE2_, kchanCalibHF1_, kchanCalibHF2_, kOffCalibHB_, kOffCalibHBX_, kOffCalibHE1_, kOffCalibHE2_, kOffCalibHEX_, kOffCalibHF1_, kOffCalibHF2_, kOffCalibHO1_, kOffCalibHO2_, kOffCalibHOX_, kPhiCalibHB_, kPhiCalibHBX_, kPhiCalibHE_, kPhiCalibHEX_, kPhiCalibHF1_, kPhiCalibHO0_, kPhiCalibHO1_, kPhiCalibHO2_, mPhiCalibHB_, mPhiCalibHE_, mPhiCalibHEX_, mPhiCalibHF1_, mPhiCalibHO0_, mPhiCalibHO1_, nchanCalibHB_, nchanCalibHE1_, nchanCalibHF1_, nchanCalibHO_, nEtaCalibHB_, nEtaCalibHE_, nEtaCalibHF_, nEtaCalibHO_, nPhiCalibHOX_, HcalCalibDetId::zside(), and ecaldqm::zside().

Referenced by HcalLogicalMapGenerator::buildCALIBMap(), HcalLogicalMapGenerator::buildHOXMap(), HcalLogicalMap::checkHashIds(), HcalLogicalMap::getHcalFrontEndId(), and HcalCondObjectContainerBase::indexFor().

                                                                   {
  HcalCalibDetId tid(id);
  int channel = tid.cboxChannel();
  int ieta = tid.ieta();
  int iphi = tid.iphi();
  int zside = tid.zside();
  unsigned int index = 0xFFFFFFFFu;

  if (tid.calibFlavor() == HcalCalibDetId::CalibrationBox) {
    HcalSubdetector subDet = tid.hcalSubdet();

    if (subDet == HcalBarrel) {
      //std::cout<<"CALIB_HB:  ";
      //dphi = 4 (18 phi values), 3 channel types (0,1,2), eta = -1 or 1
      //total of 18*3*2=108 channels
      auto indx = std::find(etaCalibHB_, etaCalibHB_ + nEtaCalibHB_, ieta);
      auto kndx = std::find(chanCalibHB_, chanCalibHB_ + nchanCalibHB_, channel);
      if (indx != etaCalibHB_ + nEtaCalibHB_ && kndx != chanCalibHB_ + nchanCalibHB_) {
        int keta = static_cast<int>(indx - etaCalibHB_);
        int kchn = static_cast<int>(kndx - chanCalibHB_);
        index = ((iphi + 1) / mPhiCalibHB_ - 1) + kPhiCalibHB_ * kchn + keta * kchanCalibHB_ + kOffCalibHB_;
      }
    } else if (subDet == HcalEndcap) {
      //std::cout<<"CALIB_HE:  ";
      //dphi = 4 (18 phi values), 7 channel types (0,1,2,3,4,5,6), eta = -1/+1
      //total of 18*7*2=252 channels
      //    if (channel > 2) channel -= 1;
      auto indx = std::find(etaCalibHE_, etaCalibHE_ + nEtaCalibHE_, ieta);
      if (indx != etaCalibHE_ + nEtaCalibHE_) {
        int keta = static_cast<int>(indx - etaCalibHE_);
        auto kndx = std::find(chanCalibHE1_, chanCalibHE1_ + nchanCalibHE1_, channel);
        if (kndx != chanCalibHE1_ + nchanCalibHE1_) {
          int kchn = static_cast<int>(kndx - chanCalibHE1_);
          index = ((iphi + 1) / mPhiCalibHE_ - 1) + kPhiCalibHE_ * kchn + keta * kchanCalibHE1_ + kOffCalibHE1_;
        } else if (channel == chanCalibHE2_) {
          index = ((iphi + 1) / mPhiCalibHE_ - 1) + keta * kchanCalibHE2_ + kOffCalibHE2_;
        }
      }
    } else if (subDet == HcalForward) {
      //std::cout<<"CALIB_HF:  ";
      //dphi = 18 (4 phi values), 3 channel types (0,1,8), eta = -1 or +1
      // or iphi = 1; channel = 9, eta = -1 or +1
      //total channels 4*3*2=24 + 2
      auto indx = std::find(etaCalibHF_, etaCalibHF_ + nEtaCalibHF_, ieta);
      if (indx != etaCalibHF_ + nEtaCalibHF_) {
        int keta = static_cast<int>(indx - etaCalibHF_);
        auto kndx = std::find(chanCalibHF1_, chanCalibHF1_ + nchanCalibHF1_, channel);
        if (kndx != chanCalibHF1_ + nchanCalibHF1_) {
          int kchn = static_cast<int>(kndx - chanCalibHF1_);
          index = ((iphi - 1) / mPhiCalibHF1_) + kPhiCalibHF1_ * kchn + keta * kchanCalibHF1_ + kOffCalibHF1_;
        } else if (channel == chanCalibHF2_) {
          index = keta * kchanCalibHF2_ + kOffCalibHF2_;
        }
      }
    } else if (subDet == HcalOuter) {
      //std::cout<<"CALIB_HO:  ";
      auto indx = std::find(etaCalibHO_, etaCalibHO_ + nEtaCalibHO_, ieta);
      if (indx != etaCalibHO_ + nEtaCalibHO_) {
        int keta = static_cast<int>(indx - etaCalibHO_);
        auto kndx = std::find(chanCalibHO_, chanCalibHO_ + nchanCalibHO_, channel);
        //there are 5 special calib crosstalk channels, one in each ring
        if (channel == chanCalibHOs_) {
          index = keta + kOffCalibHO2_;
        } else if (kndx != chanCalibHO_ + nchanCalibHO_) {
          //for HOM/HOP dphi = 6 (12 phi values),  2 channel types (0,1), eta = -2,-1 or 1,2
          //for HO0/YB0 dphi = 12 (6 phi values),  2 channel types (0,1), eta = 0
          int kchn = static_cast<int>(kndx - chanCalibHO_);
          int kphi = (ieta == 0) ? ((iphi + 1) / mPhiCalibHO0_ - 1) : ((iphi + 1) / mPhiCalibHO1_ - 1);
          if (ieta < 0) {
            index = kphi + kPhiCalibHO2_ * kchn + kPhiCalibHO1_ * keta + kOffCalibHO1_;
          } else if (ieta > 0) {
            index = kphi + kPhiCalibHO2_ * kchn + kPhiCalibHO1_ * (keta - 1) + kPhiCalibHO0_ + kOffCalibHO1_;
          } else {
            index = kphi + kPhiCalibHO2_ * kchn + kPhiCalibHO1_ * keta + kOffCalibHO1_;
          }
        }
      }
    } else {
      edm::LogWarning("HGCalGeom") << "HCAL Det Id not valid!";
    }
  } else if (tid.calibFlavor() == HcalCalibDetId::HOCrosstalk) {
    //std::cout<<"HX:  ";
    //for YB0/HO0 phi is grouped in 6 groups of 6 with dphi=2 but the transitions are 1 or 3
    // in such a way that the %36 operation yeilds unique values for every iphi
    //ieta = 0 for HO2M/HO1M ieta=2 for HO1P/HO2P;
    //ieta = 1 YB0/HO0;
    int kphi = ((std::abs(ieta) == etaCalibHOX_[0]) ? ((iphi % 2 == 0) ? (iphi / 2 - 1) : (iphi - 1) / 2) : (iphi - 1));
    if (std::abs(ieta) == etaCalibHOX_[0]) {
      index = kphi + (((zside + 1) * nPhiCalibHOX_[0]) / 2) + nPhiCalibHOX_[1] + kOffCalibHOX_;
    } else if (std::abs(ieta) == etaCalibHOX_[1]) {
      index = kphi + ((zside + 1) * nPhiCalibHOX_[1]) + kOffCalibHOX_;
    }
  } else if (tid.calibFlavor() == HcalCalibDetId::HBX) {
    if (std::abs(ieta) == etaCalibHBX_) {
      index = kOffCalibHBX_ + (iphi - 1) + (zside + 1) * kPhiCalibHBX_ / 2;
    }
  } else if (tid.calibFlavor() == HcalCalibDetId::HEX) {
    if (std::abs(ieta) == etaCalibHEX_[0]) {
      index = kOffCalibHEX_ + (iphi - 1) / mPhiCalibHEX_ + (zside + 1) * kPhiCalibHEX_ / 2;
    } else if (std::abs(ieta) == etaCalibHEX_[1]) {
      index = kOffCalibHEX_ + (iphi - 1) / mPhiCalibHEX_ + (zside + 1) * kPhiCalibHEX_ / 2 + 2 * kPhiCalibHEX_;
    }
  }
  //std::cout << "  " << ieta << "  " << zside << "  " << iphi << "  " << depth << "  " << index << std::endl;
  return index;
}

unsigned int HcalTopology::detId2denseIdHB

(

const DetId &

id

)

const

return a linear packed id from HB

Definition at line 1170 of file HcalTopology.cc.

References HcalDetId::depth(), firstHBRing(), HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), IPHI_MAX, kHBhalf, lastHBRing(), maxDepthHB_, nEtaHB_, topoVersion_, and HcalDetId::zside().

Referenced by HcalLogicalMapGenerator::buildHBEFTMap(), HcalLogicalMap::checkHashIds(), HcalLogicalMap::getHcalFrontEndId(), and HcalCondObjectContainerBase::indexFor().

                                                                {
  HcalDetId hid(id);
  const int ip(hid.iphi());
  const int ie(hid.ietaAbs());
  const int dp(hid.depth());
  const int zn(hid.zside() < 0 ? 1 : 0);
  unsigned int retval = 0xFFFFFFFFu;
  if (topoVersion_ == 0) {
    retval = (ip - 1) * 18 + dp - 1 + ie - (ie < 16 ? 1 : 0) + zn * kHBhalf;
  } else if (topoVersion_ == 10) {
    retval = (dp - 1) + maxDepthHB_ * (ip - 1);
    if (hid.ieta() > 0)
      retval += maxDepthHB_ * IPHI_MAX * (hid.ieta() - firstHBRing());
    else
      retval += maxDepthHB_ * IPHI_MAX * (hid.ieta() + lastHBRing() + nEtaHB_);
  }
  return retval;
}

unsigned int HcalTopology::detId2denseIdHE

(

const DetId &

id

)

const

return a linear packed id from HE

Definition at line 1189 of file HcalTopology.cc.

References HcalDetId::depth(), firstHERing(), HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), kHEhalf, lastHERing(), maxDepthHE_, maxPhiHE_, nEtaHE_, topoVersion_, and HcalDetId::zside().

Referenced by HcalLogicalMapGenerator::buildHBEFTMap(), HcalLogicalMap::checkHashIds(), HcalLogicalMap::getHcalFrontEndId(), and HcalCondObjectContainerBase::indexFor().

                                                                {
  HcalDetId hid(id);
  const int ip(hid.iphi());
  const int ie(hid.ietaAbs());
  const int dp(hid.depth());
  const int zn(hid.zside() < 0 ? 1 : 0);
  unsigned int retval = 0xFFFFFFFFu;
  if (topoVersion_ == 0) {
    retval = (ip - 1) * 8 + (ip / 2) * 20 +
             ((ie == 16 || ie == 17)
                  ? ie - 16
                  : ((ie >= 18 && ie <= 20)
                         ? 2 + 2 * (ie - 18) + dp - 1
                         : ((ie >= 21 && ie <= 26) ? 8 + 2 * (ie - 21) + dp - 1
                                                   : ((ie >= 27 && ie <= 28) ? 20 + 3 * (ie - 27) + dp - 1
                                                                             : 26 + 2 * (ie - 29) + dp - 1)))) +
             zn * kHEhalf;
  } else if (topoVersion_ == 10) {
    retval = (dp - 1) + maxDepthHE_ * (ip - 1);
    if (hid.ieta() > 0)
      retval += maxDepthHE_ * maxPhiHE_ * (hid.ieta() - firstHERing());
    else
      retval += maxDepthHE_ * maxPhiHE_ * (hid.ieta() + lastHERing() + nEtaHE_);
  }
  return retval;
}

unsigned int HcalTopology::detId2denseIdHF

(

const DetId &

id

)

const

return a linear packed id from HF

Definition at line 1234 of file HcalTopology.cc.

References HcalDetId::depth(), HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), IPHI_MAX, kHFhalf, maxDepthHF_, topoVersion_, and HcalDetId::zside().

Referenced by HcalLogicalMapGenerator::buildHBEFTMap(), HcalLogicalMap::checkHashIds(), HcalLogicalMap::getHcalFrontEndId(), and HcalCondObjectContainerBase::indexFor().

                                                                {
  HcalDetId hid(id);
  const int ip(hid.iphi());
  const int ie(hid.ietaAbs());
  const int dp(hid.depth());
  const int zn(hid.zside() < 0 ? 1 : 0);

  unsigned int retval = 0xFFFFFFFFu;
  if (topoVersion_ == 0) {
    retval = ((ip - 1) / 4) * 4 + ((ip - 1) / 2) * 22 + 2 * (ie - 29) + (dp - 1) + zn * kHFhalf;
  } else if (topoVersion_ == 10) {
    retval = dp - 1 + 2 * (ip - 1);
    if (hid.ieta() > 0)
      retval += maxDepthHF_ * IPHI_MAX * (hid.ieta() - 29);
    else
      retval += maxDepthHF_ * IPHI_MAX * ((41 + 13) + hid.ieta());
  }
  return retval;
}

unsigned int HcalTopology::detId2denseIdHO

(

const DetId &

id

)

const

return a linear packed id from HO

Definition at line 1216 of file HcalTopology.cc.

References HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), IPHI_MAX, kHOhalf, topoVersion_, and HcalDetId::zside().

Referenced by HcalLogicalMapGenerator::buildHOXMap(), HcalLogicalMap::checkHashIds(), HcalLogicalMap::getHcalFrontEndId(), and HcalCondObjectContainerBase::indexFor().

                                                                {
  HcalDetId hid(id);
  const int ip(hid.iphi());
  const int ie(hid.ietaAbs());
  const int zn(hid.zside() < 0 ? 1 : 0);

  unsigned int retval = 0xFFFFFFFFu;
  if (topoVersion_ == 0) {
    retval = (ip - 1) * 15 + (ie - 1) + zn * kHOhalf;
  } else if (topoVersion_ == 10) {
    if (hid.ieta() > 0)
      retval = (ip - 1) + IPHI_MAX * (hid.ieta() - 1);
    else
      retval = (ip - 1) + IPHI_MAX * (30 + hid.ieta());
  }
  return retval;
}

unsigned int HcalTopology::detId2denseIdHT

(

const DetId &

id

)

const

return a linear packed id from HT

Definition at line 1254 of file HcalTopology.cc.

References HcalTrigTowerDetId::ietaAbs(), HcalTrigTowerDetId::iphi(), kHThalf, kHTSizePhase1, kHTSizePreLS1, HcalTrigTowerDetId::version(), HcalTrigTowerDetId::zside(), and ecaldqm::zside().

Referenced by HcalLogicalMap::checkHashIds(), CaloTPGTranscoderULUT::getOutputLUTId(), and HcalCondObjectContainerBase::indexFor().

                                                                {
  HcalTrigTowerDetId tid(id);
  int zside = tid.zside();
  unsigned int ietaAbs = tid.ietaAbs();
  unsigned int iphi = tid.iphi();
  unsigned int ivers = tid.version();

  unsigned int index;
  if (ivers == 0) {
    if ((iphi - 1) % 4 == 0)
      index = (iphi - 1) * 32 + (ietaAbs - 1) - (12 * ((iphi - 1) / 4));
    else
      index = (iphi - 1) * 28 + (ietaAbs - 1) + (4 * (((iphi - 1) / 4) + 1));
    if (zside == -1)
      index += kHThalf;
  } else {
    index = kHTSizePreLS1;
    if (zside == -1)
      index += ((kHTSizePhase1 - kHTSizePreLS1) / 2);
    index += (36 * (ietaAbs - 30) + ((iphi - 1) / 2));
  }

  return index;
}

unsigned int HcalTopology::detId2denseIdPreLS1 ( const DetId &  id ) const

private

Definition at line 1141 of file HcalTopology.cc.

References HcalDetId::depth(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalDetId::ietaAbs(), HcalDetId::iphi(), kHBhalf, kHEhalf, kHFhalf, kHOhalf, sd, HcalDetId::subdet(), and HcalDetId::zside().

Referenced by detId2denseId().

                                                                    {
  HcalDetId hid(id);
  const HcalSubdetector sd(hid.subdet());
  const int ip(hid.iphi());
  const int ie(hid.ietaAbs());
  const int dp(hid.depth());
  const int zn(hid.zside() < 0 ? 1 : 0);
  unsigned int retval =
      ((sd == HcalBarrel)
           ? (ip - 1) * 18 + dp - 1 + ie - (ie < 16 ? 1 : 0) + zn * kHBhalf
           : ((sd == HcalEndcap)
                  ? 2 * kHBhalf + (ip - 1) * 8 + (ip / 2) * 20 +
                        ((ie == 16 || ie == 17)
                             ? ie - 16
                             : ((ie >= 18 && ie <= 20)
                                    ? 2 + 2 * (ie - 18) + dp - 1
                                    : ((ie >= 21 && ie <= 26)
                                           ? 8 + 2 * (ie - 21) + dp - 1
                                           : ((ie >= 27 && ie <= 28) ? 20 + 3 * (ie - 27) + dp - 1
                                                                     : 26 + 2 * (ie - 29) + dp - 1)))) +
                        zn * kHEhalf
                  : ((sd == HcalOuter)
                         ? 2 * kHBhalf + 2 * kHEhalf + (ip - 1) * 15 + (ie - 1) + zn * kHOhalf
                         : ((sd == HcalForward) ? 2 * kHBhalf + 2 * kHEhalf + 2 * kHOhalf + ((ip - 1) / 4) * 4 +
                                                      ((ip - 1) / 2) * 22 + 2 * (ie - 29) + (dp - 1) + zn * kHFhalf
                                                : 0xFFFFFFFFu))));
  return retval;
}

int HcalTopology::doublePhiBins ( ) const

inline

Definition at line 106 of file HcalTopology.h.

References doublePhiBins_.

{ return doublePhiBins_; }

std::vector< DetId > HcalTopology::down ( const DetId &  id ) const

overridevirtual

Get the neighbors of the given cell in down direction (inward)

Implements CaloSubdetectorTopology.

Definition at line 458 of file HcalTopology.cc.

References HcalDetId::changeForm(), decrementDepth(), gpuClustering::id, and HcalDetId::oldFormat().

                                                         {
  HcalDetId neighbor = id;
  std::vector<DetId> vNeighborsDetId;
  if (decrementDepth(neighbor)) {
    if (neighbor.oldFormat())
      neighbor.changeForm();
    vNeighborsDetId.emplace_back(neighbor);
  }
  return vNeighborsDetId;
}

std::vector< DetId > HcalTopology::east ( const DetId &  id ) const

overridevirtual

Flag=0 for unmerged Id's; =1 for for merged Id's; =2 for either Get the neighbors of the given cell in east direction

Implements CaloSubdetectorTopology.

Definition at line 403 of file HcalTopology.cc.

References decIEta(), and mps_fire::i.

Referenced by spr::newHCALIdEW().

                                                         {
  std::vector<DetId> vNeighborsDetId;
  HcalDetId neighbors[2];
  for (int i = 0; i < decIEta(HcalDetId(id), neighbors); i++) {
    if (neighbors[i].oldFormat())
      neighbors[i].changeForm();
    vNeighborsDetId.emplace_back(DetId(neighbors[i].rawId()));
  }
  return vNeighborsDetId;
}

double HcalTopology::etaMax

(

HcalSubdetector

subdet

)

const

Definition at line 1101 of file HcalTopology.cc.

References PVValHelper::eta, etaTable, etaTableHF, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, lastHBRing_, lastHERing_, lastHORing_, and nEtaHE_.

Referenced by HcalGeometry::getCells(), and HcalGeometry::getClosestCell().

                                                        {
  double eta(0);
  switch (subdet) {
    case (HcalBarrel):
      if (lastHBRing_ < (int)(etaTable.size()))
        eta = etaTable[lastHBRing_];
      break;
    case (HcalEndcap):
      if (lastHERing_ < (int)(etaTable.size()) && nEtaHE_ > 0)
        eta = etaTable[lastHERing_];
      break;
    case (HcalOuter):
      if (lastHORing_ < (int)(etaTable.size()))
        eta = etaTable[lastHORing_];
      break;
    case (HcalForward):
      if (!etaTableHF.empty())
        eta = etaTableHF[etaTableHF.size() - 1];
      break;
    default:
      eta = 0;
  }
  return eta;
}

std::pair< double, double > HcalTopology::etaRange	(	HcalSubdetector	subdet,
		int	ieta
	)		const

Definition at line 1125 of file HcalTopology.cc.

References etaTable, etaTableHF, firstHFRing_, HcalForward, cuy::ii, lastHERing_, HcalTopologyMode::LHC, and mode_.

Referenced by HcalHF_PETalgorithm::HFSetFlagFromPET(), HcalHFStatusBitFromRecHits::hfSetFlagFromRecHits(), HcalHF_S9S1algorithm::HFSetFlagFromS9S1(), HcaluLUTTPGCoder::make_cosh_ieta_map(), and HcalTrigTowerGeometry::towerEtaBounds().

                                                                                     {
  int ieta = (keta > 0) ? keta : -keta;
  if (subdet == HcalForward) {
    if (ieta >= firstHFRing_) {
      unsigned int ii = (unsigned int)(ieta - firstHFRing_);
      if (ii + 1 < etaTableHF.size())
        return std::pair<double, double>(etaTableHF[ii], etaTableHF[ii + 1]);
    }
  } else {
    int ietal = (mode_ == HcalTopologyMode::LHC && ieta == lastHERing_ - 1) ? (ieta + 1) : ieta;
    if ((ietal < (int)(etaTable.size())) && (ieta > 0))
      return std::pair<double, double>(etaTable[ieta - 1], etaTable[ietal]);
  }
  return std::pair<double, double>(0, 0);
}

int HcalTopology::etaRing	(	HcalSubdetector	subdet,
		double	eta
	)		const

eta and phi index from eta, phi values

Definition at line 998 of file HcalTopology.cc.

References etaTable, etaTableHF, firstHBRing_, firstHFRing_, HcalForward, isotrackApplyRegressor::k, and lastHERing_.

Referenced by decrementDepth(), HcalGeometry::etaRing(), and incrementDepth().

                                                                 {
  int etaring = firstHBRing_;
  if (bc == HcalForward) {
    etaring = firstHFRing_;
    for (unsigned int k = 0; k < etaTableHF.size() - 1; ++k) {
      if (abseta < etaTableHF[k + 1]) {
        etaring += k;
        break;
      }
    }
  } else {
    for (unsigned int k = 0; k < etaTable.size() - 1; ++k) {
      if (abseta < etaTable[k + 1]) {
        etaring += k;
        break;
      }
    }
    if (abseta >= etaTable[etaTable.size() - 1])
      etaring = lastHERing_;
  }
  return etaring;
}

void HcalTopology::exclude

(

const HcalDetId &

id

)

Add a cell to exclusion list

Definition at line 377 of file HcalTopology.cc.

References exclusionList_, mps_fire::i, gpuClustering::id, and cuda_std::lower_bound().

Referenced by exclude(), and HcalTopologyRestrictionParser::parse().

                                              {
  std::vector<HcalDetId>::iterator i = std::lower_bound(exclusionList_.begin(), exclusionList_.end(), id);
  if (i == exclusionList_.end() || *i != id) {
    exclusionList_.insert(i, id);
  }
}

int HcalTopology::exclude	(	HcalSubdetector	subdet,
		int	ieta1,
		int	ieta2,
		int	iphi1,
		int	iphi2,
		int	depth1 = 1,
		int	depth2 = 4
	)

Exclude an eta/phi/depth range for a given subdetector

Definition at line 469 of file HcalTopology.cc.

References HLT_FULL_cff::depth, exclude(), excludeHB_, excludeHE_, excludeHF_, excludeHO_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, gpuClustering::id, SiStripPI::max, SiStripPI::min, dqmiodumpmetadata::n, and validRaw().

                                                                                                                    {
  bool exed = false;
  // first, check the full detector exclusions...  (fast)
  switch (subdet) {
    case (HcalBarrel):
      exed = excludeHB_;
      break;
    case (HcalEndcap):
      exed = excludeHE_;
      break;
    case (HcalOuter):
      exed = excludeHO_;
      break;
    case (HcalForward):
      exed = excludeHF_;
      break;
    default:
      exed = false;
  }
  if (exed)
    return 0;  // if the whole detector is excluded...

  int ieta_l = std::min(ieta1, ieta2);
  int ieta_h = std::max(ieta1, ieta2);
  int iphi_l = std::min(iphi1, iphi2);
  int iphi_h = std::max(iphi1, iphi2);
  int depth_l = std::min(depth1, depth2);
  int depth_h = std::max(depth1, depth2);

  int n = 0;
  for (int ieta = ieta_l; ieta <= ieta_h; ieta++)
    for (int iphi = iphi_l; iphi <= iphi_h; iphi++)
      for (int depth = depth_l; depth <= depth_h; depth++) {
        HcalDetId id(subdet, ieta, iphi, depth);
        if (validRaw(id)) {  // use 'validRaw' to include check validity in "uncut" detector
          exclude(id);
          n++;
        }
      }
  return n;
}

void HcalTopology::excludeSubdetector

(

HcalSubdetector

subdet

)

Exclude an entire subdetector

Definition at line 384 of file HcalTopology.cc.

References excludeHB_, excludeHE_, excludeHF_, excludeHO_, HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

Referenced by HcalTopologyRestrictionParser::parse().

                                                            {
  switch (subdet) {
    case (HcalBarrel):
      excludeHB_ = true;
      break;
    case (HcalEndcap):
      excludeHE_ = true;
      break;
    case (HcalOuter):
      excludeHO_ = true;
      break;
    case (HcalForward):
      excludeHF_ = true;
      break;
    default:
      break;
  }
}

int HcalTopology::firstHBRing ( ) const

inline

Definition at line 91 of file HcalTopology.h.

References firstHBRing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertCTtoHcal(), CaloTowerTopology::convertHcaltoCT(), denseId2detId(), detId2denseId(), detId2denseIdHB(), HcaluLUTTPGCoder::init(), and CaloTowerHardcodeGeometryLoader::load().

{ return firstHBRing_; }

int HcalTopology::firstHEDoublePhiRing ( ) const

inline

Definition at line 101 of file HcalTopology.h.

References firstHEDoublePhiRing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), decAIEta(), decIPhi(), HcalTrigTowerGeometry::detIds(), HcalLutManager::getHEFineGrainLUTs(), CaloTowersCreationAlgo::getThresholdAndWeight(), incAIEta(), incIPhi(), HcalHardcodeGeometryLoader::makeHECells(), nPhiBins(), HcalTrigTowerGeometry::towerIds(), and validRaw().

{ return firstHEDoublePhiRing_; }

int HcalTopology::firstHEQuadPhiRing ( ) const

inline

Definition at line 102 of file HcalTopology.h.

References firstHEQuadPhiRing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), decAIEta(), decIPhi(), incAIEta(), incIPhi(), nPhiBins(), and validRaw().

{ return firstHEQuadPhiRing_; }

int HcalTopology::firstHERing ( ) const

inline

Definition at line 93 of file HcalTopology.h.

References firstHERing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertCTtoHcal(), CaloTowerTopology::convertHcaltoCT(), decAIEta(), decrementDepth(), denseId2detId(), depthBinInformation(), detId2denseId(), detId2denseIdHE(), HcalTrigTowerGeometry::detIds(), HcaluLUTTPGCoder::init(), CaloTowerHardcodeGeometryLoader::load(), validDetIdPreLS1(), and validRaw().

{ return firstHERing_; }

int HcalTopology::firstHETripleDepthRing ( ) const

inline

Definition at line 104 of file HcalTopology.h.

References firstHETripleDepthRing_.

Referenced by depthBinInformation(), and HcalHardcodeGeometryLoader::makeHECells().

{ return firstHETripleDepthRing_; }

int HcalTopology::firstHFQuadPhiRing ( ) const

inline

Definition at line 103 of file HcalTopology.h.

References firstHFQuadPhiRing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), decAIEta(), decIPhi(), incAIEta(), incIPhi(), nPhiBins(), validDetIdPreLS1(), and validRaw().

{ return firstHFQuadPhiRing_; }

int HcalTopology::firstHFRing ( ) const

inline

Definition at line 96 of file HcalTopology.h.

References firstHFRing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertCTtoHcal(), CaloTowerTopology::convertHcaltoCT(), decAIEta(), decrementDepth(), HcalTrigTowerGeometry::detIds(), HcalTrigTowerGeometry::firstHFRingInTower(), HcaluLUTTPGCoder::init(), CaloTowerHardcodeGeometryLoader::load(), CaloTowerHardcodeGeometryLoader::makeCell(), validDetIdPreLS1(), and validRaw().

{ return firstHFRing_; }

int HcalTopology::firstHORing ( ) const

inline

Definition at line 98 of file HcalTopology.h.

References firstHORing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), and CaloTowerTopology::convertHcaltoCT().

{ return firstHORing_; }

unsigned int HcalTopology::getCALIBSize ( ) const

inline

Definition at line 137 of file HcalTopology.h.

References CALIBSize_.

Referenced by HcalCondObjectContainerBase::sizeFor().

{ return CALIBSize_; }

void HcalTopology::getDepthSegmentation	(	const unsigned	ring,
		std::vector< int > &	readoutDepths,
		const bool	flag = false
	)		const

for each of the ~17 depth segments, specify which readout bin they belong to if the ring is not found, the first one with a lower ring will be returned.

Definition at line 1065 of file HcalTopology.cc.

References depthSegmentation_, depthSegmentationOne_, Exception, and relativeConstraints::ring.

Referenced by HcalHardcodeGeometryLoader::load(), and segmentBoundaries().

                                                                                                                {
  // if it doesn't exist, return the first entry with a lower index.  So if we only
  // have entries for 1 and 17, any input from 1-16 should return the entry for ring 1
  SegmentationMap::const_iterator pos;
  if (!one) {
    pos = depthSegmentation_.upper_bound(ring);
    if (pos == depthSegmentation_.begin()) {
      throw cms::Exception("HcalTopology") << "No depth segmentation found for ring" << ring;
    }
  } else {
    pos = depthSegmentationOne_.upper_bound(ring);
    if (pos == depthSegmentationOne_.begin()) {
      throw cms::Exception("HcalTopology") << "No depth segmentation found for ring" << ring;
    }
  }
  --pos;
  // pos now refers to the last element with key <= ring.
  readoutDepths = pos->second;
}

unsigned int HcalTopology::getHBSize ( ) const

inline

Definition at line 132 of file HcalTopology.h.

References HBSize_.

Referenced by HcalGeometry::init(), and HcalCondObjectContainerBase::sizeFor().

{ return HBSize_; }

unsigned int HcalTopology::getHESize ( ) const

inline

Definition at line 133 of file HcalTopology.h.

References HESize_.

Referenced by HcalGeometry::init(), and HcalCondObjectContainerBase::sizeFor().

{ return HESize_; }

unsigned int HcalTopology::getHFSize ( ) const

inline

Definition at line 135 of file HcalTopology.h.

References HFSize_.

Referenced by HcalGeometry::init(), HcalHardcodeGeometryLoader::load(), HcalFlexiHardcodeGeometryLoader::load(), and HcalCondObjectContainerBase::sizeFor().

{ return HFSize_; }

unsigned int HcalTopology::getHOSize ( ) const

inline

Definition at line 134 of file HcalTopology.h.

References HOSize_.

Referenced by HcalGeometry::init(), and HcalCondObjectContainerBase::sizeFor().

{ return HOSize_; }

unsigned int HcalTopology::getHTSize ( ) const

inline

Definition at line 136 of file HcalTopology.h.

References HTSize_.

Referenced by CaloTPGTranscoderULUT::setup(), and HcalCondObjectContainerBase::sizeFor().

{ return HTSize_; }

bool HcalTopology::getMergePositionFlag ( ) const

inline

Definition at line 167 of file HcalTopology.h.

References mergePosition_.

Referenced by CaloTowersCreationAlgo::assignHitHcal(), CaloTowerConstituentsMap::constituentsOf(), CaloTowersCreationAlgo::hadShwPosFromCells(), HcalRecHitsAnalyzer::hcalSevLvl(), PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::importRecHits(), HcalGeometry::init(), HcalHitSelection::skim(), and PFRecHitQTestHCALChannel::test().

{ return mergePosition_; }

unsigned int HcalTopology::getNumberOfShapes ( ) const

inline

Definition at line 161 of file HcalTopology.h.

References numberOfShapes_.

Referenced by HcalGeometry::numberOfShapes().

{ return numberOfShapes_; }

int HcalTopology::getPhiZOne ( std::vector< std::pair< int, int > > &  phiz ) const

inline

Definition at line 130 of file HcalTopology.h.

References HcalDDDRecConstants::getPhiZOne(), and hcons_.

{ return hcons_->getPhiZOne(phiz); }

HcalDetId HcalTopology::idBack ( const HcalDetId &  id ) const

inline

Definition at line 171 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::idBack().

Referenced by CaloTowersCreationAlgo::hadShwPosFromCells(), and CaloTowersCreationAlgo::hadShwrPos().

{ return hcons_->idBack(id); }

HcalDetId HcalTopology::idFront ( const HcalDetId &  id ) const

inline

Definition at line 170 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::idFront().

Referenced by CaloTowersCreationAlgo::assignHitHcal(), HcalGeometry::getGeometry(), HcalGeometry::getPosition(), EgammaHcalIsolation::goodHitEnergy(), CaloTowersCreationAlgo::hadShwPosFromCells(), CaloTowersCreationAlgo::hadShwrPos(), egamma::hasActiveHcal(), CaloTowersCreationAlgo::hcalChanStatusForCaloTower(), HcalRecHitsAnalyzer::hcalSevLvl(), PFHcalRecHitCreator< Digi, Geometry, Layer, Detector >::importRecHits(), HcalHitSelection::skim(), and PFRecHitQTestHCALChannel::test().

{ return hcons_->idFront(id); }

int HcalTopology::incAIEta ( const HcalDetId &  id, HcalDetId  neighbors[2]  ) const

private

Get the neighbors of the given cell with higher absolute ieta

Increasing in |ieta|, there is always at most one neighbor

Definition at line 745 of file HcalTopology.cc.

References HLT_FULL_cff::depth, etaHE2HF_, firstHEDoublePhiRing(), firstHEQuadPhiRing(), firstHFQuadPhiRing(), HcalBarrel, HcalEndcap, HcalForward, lastHBRing(), lastHERing(), dqmiodumpmetadata::n, valid(), and ecaldqm::zside().

Referenced by decIEta(), and incIEta().

                                                                            {
  int n = 1;
  int aieta = id.ietaAbs();

  if (aieta == firstHEDoublePhiRing() - 1 && (id.iphi() % 2) == 0)
    neighbors[0] = HcalDetId(id.subdet(), (aieta + 1) * id.zside(), id.iphi() - 1, id.depth());
  else if (aieta == firstHFQuadPhiRing() - 1 && ((id.iphi() + 1) % 4) != 0)
    neighbors[0] =
        HcalDetId(id.subdet(), (aieta + 1) * id.zside(), ((id.iphi() == 1) ? (71) : (id.iphi() - 2)), id.depth());
  else if (aieta == firstHEQuadPhiRing() - 1 && ((id.iphi() + 1) % 4) != 0)
    neighbors[0] =
        HcalDetId(id.subdet(), (aieta + 1) * id.zside(), ((id.iphi() == 1) ? (71) : (id.iphi() - 2)), id.depth());
  else if (aieta == lastHBRing() && id.subdet() == HcalBarrel)
    neighbors[0] = HcalDetId(HcalEndcap, (aieta + 1) * id.zside(), id.iphi(), 1);
  else if (aieta == lastHERing() && id.subdet() == HcalEndcap)
    neighbors[0] = HcalDetId(HcalForward, etaHE2HF_ * id.zside(), id.iphi(), 1);
  else
    neighbors[0] = HcalDetId(id.subdet(), (aieta + 1) * id.zside(), id.iphi(), id.depth());

  if (!valid(neighbors[0]))
    n = 0;
  return n;
}

int HcalTopology::incIEta	(	const HcalDetId &	id,
		HcalDetId	neighbors[2]
	)		const

Get the neighbors of the given cell with higher (signed) ieta

Definition at line 730 of file HcalTopology.cc.

References decAIEta(), incAIEta(), and ecaldqm::zside().

Referenced by west().

                                                                           {
  if (id.zside() == 1)
    return incAIEta(id, neighbors);
  else
    return decAIEta(id, neighbors);
}

bool HcalTopology::incIPhi	(	const HcalDetId &	id,
		HcalDetId &	neighbor
	)		const

Get the neighbor (if present) of the given cell with higher iphi

Definition at line 627 of file HcalTopology.cc.

References HLT_FULL_cff::depth, firstHEDoublePhiRing(), firstHEQuadPhiRing(), firstHFQuadPhiRing(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, IPHI_MAX, maxPhiHE_, convertSQLiteXML::ok, valid(), and validRaw().

Referenced by north(), and HFStripFilter::runFilter().

                                                                         {
  bool ok = valid(id);
  if (ok) {
    switch (id.subdet()) {
      case (HcalBarrel):
      case (HcalOuter):
        if (id.iphi() == IPHI_MAX)
          neighbor = HcalDetId(id.subdet(), id.ieta(), 1, id.depth());
        else
          neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() + 1, id.depth());
        break;
      case (HcalEndcap):
        if (id.ietaAbs() >= firstHEQuadPhiRing()) {
          if (id.iphi() == IPHI_MAX - 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), 3, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() + 4, id.depth());
        } else if (id.ietaAbs() >= firstHEDoublePhiRing()) {
          if (id.iphi() == IPHI_MAX - 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() + 2, id.depth());
        } else {
          if (id.iphi() == maxPhiHE_)
            neighbor = HcalDetId(id.subdet(), id.ieta(), 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() + 1, id.depth());
        }
        break;
      case (HcalForward):
        if (id.ietaAbs() >= firstHFQuadPhiRing()) {
          if (id.iphi() == IPHI_MAX - 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), 3, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() + 4, id.depth());
        } else {
          if (id.iphi() == IPHI_MAX - 1)
            neighbor = HcalDetId(id.subdet(), id.ieta(), 1, id.depth());
          else
            neighbor = HcalDetId(id.subdet(), id.ieta(), id.iphi() + 2, id.depth());
        }
        if (!validRaw(neighbor))
          ok = false;
        break;
      default:
        ok = false;
    }
  }
  return ok;
}

bool HcalTopology::incrementDepth

(

HcalDetId &

id

)

const

Get the detector behind this one

Definition at line 870 of file HcalTopology.cc.

References HcalDetId::depth(), HLT_FULL_cff::depth, depthBinInformation(), etaRing(), HcalDDDRecConstants::getDepthEta16(), HcalTopologyMode::H2HE, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons_, HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), lastHBRing(), lastHERing(), lastHORing(), mode_, HcalTopologyMode::SLHC, HcalDetId::subdet(), validRaw(), HcalDetId::zside(), and ecaldqm::zside().

Referenced by HcalGeometry::getClosestCell(), and up().

                                                        {
  HcalSubdetector subdet = detId.subdet();
  int ieta = detId.ieta();
  int etaRing = detId.ietaAbs();
  int depth = detId.depth();
  int iphi = detId.iphi();
  int zside = detId.zside();
  int nDepthBins(0), startingBin(0);
  depthBinInformation(subdet, etaRing, iphi, zside, nDepthBins, startingBin);

  // see if the new depth bin exists
  ++depth;
  if (depth >= (startingBin + nDepthBins)) {
    // handle on a case-by-case basis
    if (subdet == HcalBarrel && etaRing < lastHORing()) {
      // HO
      subdet = HcalOuter;
      depth = 4;
    } else if (subdet == HcalBarrel && etaRing == lastHBRing()) {
      // overlap
      subdet = HcalEndcap;
      if (mode_ == HcalTopologyMode::SLHC || mode_ == HcalTopologyMode::H2HE)
        depth = hcons_->getDepthEta16(2, iphi, zside);
    } else if (subdet == HcalEndcap && etaRing == lastHERing() - 1 && mode_ != HcalTopologyMode::SLHC) {
      // guard ring HF29 is behind HE 28
      subdet = HcalForward;
      (ieta > 0) ? ++ieta : --ieta;
      depth = 1;
    } else if (subdet == HcalEndcap && etaRing == lastHERing() && mode_ != HcalTopologyMode::SLHC) {
      // split cells go to bigger granularity.  Ring 29 -> 28
      (ieta > 0) ? --ieta : ++ieta;
    } else {
      // no more chances
      detId = HcalDetId();
      return false;
    }
  }
  detId = HcalDetId(subdet, ieta, iphi, depth);
  return validRaw(detId);
}

bool HcalTopology::isBH ( ) const

inline

Definition at line 162 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::isBH().

Referenced by CaloTowerTopology::CaloTowerTopology().

{ return ((hcons_ == nullptr) ? false : hcons_->isBH()); }

bool HcalTopology::isExcluded ( const HcalDetId &  id ) const

private

Definition at line 349 of file HcalTopology.cc.

References excludeHB_, excludeHE_, excludeHF_, excludeHO_, exclusionList_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, gpuClustering::id, and cuda_std::lower_bound().

Referenced by validHcal().

                                                       {
  bool exed = false;
  // first, check the full detector exclusions...  (fast)
  switch (id.subdet()) {
    case (HcalBarrel):
      exed = excludeHB_;
      break;
    case (HcalEndcap):
      exed = excludeHE_;
      break;
    case (HcalOuter):
      exed = excludeHO_;
      break;
    case (HcalForward):
      exed = excludeHF_;
      break;
    default:
      exed = false;
  }
  // next, check the list (slower)
  if (!exed && !exclusionList_.empty()) {
    std::vector<HcalDetId>::const_iterator i = std::lower_bound(exclusionList_.begin(), exclusionList_.end(), id);
    if (i != exclusionList_.end() && *i == id)
      exed = true;
  }
  return exed;
}

int HcalTopology::lastHBHERing ( ) const

inline

Definition at line 95 of file HcalTopology.h.

References lastHBRing_, lastHERing_, and SiStripPI::max.

Referenced by CaloTPGTranscoderULUT::getOutputLUTSize().

{ return std::max(lastHBRing_, lastHERing_); }

int HcalTopology::lastHBRing ( ) const

inline

Definition at line 92 of file HcalTopology.h.

References lastHBRing_.

Referenced by HcalTriggerPrimitiveAlgo::analyzeQIE11(), HcalTrigPrimDigiProducer::beginRun(), CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertHcaltoCT(), denseId2detId(), depthBinInformation(), detId2denseId(), detId2denseIdHB(), HcalTrigTowerGeometry::detIds(), CaloTPGTranscoderULUT::getOutputLUTSize(), reco::hcalSubdetector(), incAIEta(), incrementDepth(), HcaluLUTTPGCoder::init(), CaloTPGTranscoderULUT::isOnlyQIE11(), CaloTowerHardcodeGeometryLoader::load(), HcalTrigTowerGeometry::towerIds(), HcaluLUTTPGCoder::update(), and validRaw().

{ return lastHBRing_; }

int HcalTopology::lastHERing ( ) const

inline

Definition at line 94 of file HcalTopology.h.

References lastHERing_.

Referenced by CaloTowersCreationAlgo::assignHitHcal(), HcalTrigPrimDigiProducer::beginRun(), CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertHcaltoCT(), decrementDepth(), denseId2detId(), depthBinInformation(), detId2denseId(), detId2denseIdHE(), HcalTrigTowerGeometry::detIds(), HcalLutManager::getHEFineGrainLUTs(), CaloTPGTranscoderULUT::getOutputLUTSize(), reco::hcalSubdetector(), incAIEta(), incrementDepth(), HcaluLUTTPGCoder::init(), CaloTowerHardcodeGeometryLoader::load(), HcalHardcodeGeometryLoader::load(), CaloTowersCreationAlgo::makeHcalDropChMap(), HcalHardcodeGeometryLoader::makeHECells(), HcalTrigTowerGeometry::towerIds(), HcaluLUTTPGCoder::update(), and validRaw().

{ return lastHERing_; }

int HcalTopology::lastHFRing ( ) const

inline

Definition at line 97 of file HcalTopology.h.

References lastHFRing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertHcaltoCT(), reco::hcalSubdetector(), HcaluLUTTPGCoder::init(), CaloTowerHardcodeGeometryLoader::load(), validDetIdPreLS1(), and validRaw().

{ return lastHFRing_; }

int HcalTopology::lastHORing ( ) const

inline

Definition at line 99 of file HcalTopology.h.

References lastHORing_.

Referenced by CaloTowerTopology::CaloTowerTopology(), CaloTowerTopology::convertHcaltoCT(), HcalTrigTowerGeometry::detIds(), incrementDepth(), and validRaw().

{ return lastHORing_; }

int HcalTopology::maxDepth

(

void

)

const

Definition at line 981 of file HcalTopology.cc.

References SiStripPI::max, maxDepthHB_, maxDepthHE_, maxDepthHF_, and minMaxDepth_.

Referenced by HcalGeometry::getCells(), and HcaluLUTTPGCoder::init().

                                 {
  int maxd1 = std::max(maxDepthHB_, maxDepthHE_);
  int maxd2 = std::max(maxDepthHF_, minMaxDepth_);
  return std::max(maxd1, maxd2);
}

int HcalTopology::maxDepth

(

HcalSubdetector

subdet

)

const

Definition at line 987 of file HcalTopology.cc.

References HcalBarrel, HcalEndcap, HcalForward, maxDepthHB_, maxDepthHE_, and maxDepthHF_.

                                                   {
  if (bc == HcalBarrel)
    return maxDepthHB_;
  else if (bc == HcalEndcap)
    return maxDepthHE_;
  else if (bc == HcalForward)
    return maxDepthHF_;
  else
    return 4;
}

int HcalTopology::maxDepthHB ( ) const

inline

Definition at line 139 of file HcalTopology.h.

References maxDepthHB_.

Referenced by RecAnalyzerMinbias::beginRun().

{ return maxDepthHB_; }

int HcalTopology::maxDepthHE ( ) const

inline

Definition at line 140 of file HcalTopology.h.

References maxDepthHE_.

Referenced by RecAnalyzerMinbias::beginRun(), and HcalTrigTowerGeometry::towerIds().

{ return maxDepthHE_; }

bool HcalTopology::mergedDepth29 ( HcalDetId  id ) const

inline

Definition at line 111 of file HcalTopology.h.

References HLT_FULL_cff::depth, hcons_, and HcalDDDRecConstants::mergedDepthList29().

Referenced by CaloTowersCreationAlgo::assignHitHcal(), and CaloTowersCreationAlgo::makeHcalDropChMap().

{ return hcons_->mergedDepthList29(id.ieta(), id.iphi(), id.depth()); }

HcalDetId HcalTopology::mergedDepthDetId ( const HcalDetId &  id ) const

inline

Definition at line 166 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::mergedDepthDetId().

Referenced by CaloTowerConstituentsMap::constituentsOf(), HcalGeometry::correctId(), CaloTowersCreationAlgo::makeHcalDropChMap(), and SimplePlan1RechitCombiner::mapRechit().

{ return hcons_->mergedDepthDetId(id); }

std::vector<int> HcalTopology::mergedDepthList29 ( HcalDetId  id ) const

inline

Definition at line 112 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::mergedDepthList29().

{ return hcons_->mergedDepthList29(id.ieta(), id.iphi()); }

HcalTopologyMode::Mode HcalTopology::mode ( void  ) const

inline

Definition at line 34 of file HcalTopology.h.

References mode_.

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

{ return mode_; }

unsigned int HcalTopology::ncells ( ) const

overridevirtual

return a count of valid cells (for dense indexing use)

Reimplemented from CaloSubdetectorTopology.

Definition at line 1687 of file HcalTopology.cc.

References HBSize_, HESize_, HFSize_, and HOSize_.

Referenced by denseId2detId(), HcalGeometry::getSummary(), HcalHardcodeGeometryLoader::load(), HcalFlexiHardcodeGeometryLoader::load(), and HcalGeometry::sizeForDenseIndex().

{ return HBSize_ + HESize_ + HOSize_ + HFSize_; }

std::vector< DetId > HcalTopology::north ( const DetId &  id ) const

overridevirtual

Get the neighbors of the given cell in north direction

Implements CaloSubdetectorTopology.

Definition at line 425 of file HcalTopology.cc.

References HcalDetId::changeForm(), incIPhi(), HcalDetId::oldFormat(), and DetId::rawId().

Referenced by spr::newHCALIdNS().

                                                          {
  std::vector<DetId> vNeighborsDetId;
  HcalDetId neighbor;
  if (incIPhi(HcalDetId(id), neighbor)) {
    if (neighbor.oldFormat())
      neighbor.changeForm();
    vNeighborsDetId.emplace_back(DetId(neighbor.rawId()));
  }
  return vNeighborsDetId;
}

int HcalTopology::nPhiBins

(

int

etaRing

)

const

how many phi segments in this ring

Definition at line 956 of file HcalTopology.cc.

References doublePhiBins_, firstHEDoublePhiRing(), firstHEQuadPhiRing(), firstHFQuadPhiRing(), HcalDDDRecConstants::getEtaRange(), HcalBarrel, hcons_, and singlePhiBins_.

Referenced by CaloTowerHardcodeGeometryLoader::makeCell().

                                            {
  int lastPhiBin = singlePhiBins_;
  if (etaRing >= firstHFQuadPhiRing() || etaRing >= firstHEQuadPhiRing())
    lastPhiBin = doublePhiBins_ / 2;
  else if (etaRing >= firstHEDoublePhiRing())
    lastPhiBin = doublePhiBins_;
  if (hcons_ && etaRing >= hcons_->getEtaRange(1).first && etaRing <= hcons_->getEtaRange(1).second) {
    return nPhiBins(HcalBarrel, etaRing);
  }
  return lastPhiBin;
}

int HcalTopology::nPhiBins	(	HcalSubdetector	subdet,
		int	etaRing
	)		const

Definition at line 968 of file HcalTopology.cc.

References dPhiTable, dPhiTableHF, firstHBRing_, firstHFRing_, and HcalForward.

                                                                {
  double phiTableVal;
  if (bc == HcalForward) {
    phiTableVal = dPhiTableHF[etaRing - firstHFRing_];
  } else {
    phiTableVal = dPhiTable[etaRing - firstHBRing_];
  }
  int lastPhiBin = 0;
  if (phiTableVal != 0.0)
    lastPhiBin = static_cast<int>((2._pi / phiTableVal) + 0.001);
  return lastPhiBin;
}

int HcalTopology::phiBin	(	HcalSubdetector	subdet,
		int	etaRing,
		double	phi
	)		const

Definition at line 1021 of file HcalTopology.cc.

References dPhiTable, dPhiTableHF, firstHBRing_, firstHFRing_, HcalBarrel, HcalEndcap, HcalForward, phioff, unit(), unitPhi, and unitPhiHF.

Referenced by HcalGeometry::phiBin().

                                                                          {
  //put phi in correct range (0->2pi)
  int index(0);
  if (bc == HcalBarrel) {
    index = (etaring - firstHBRing_);
    phi -= phioff[0];
  } else if (bc == HcalEndcap) {
    index = (etaring - firstHBRing_);
    phi -= phioff[1];
  } else if (bc == HcalForward) {
    index = (etaring - firstHFRing_);
    if (index < static_cast<int>(dPhiTableHF.size())) {
      if (index >= 0 && unitPhiHF[index] > 2)
        phi -= phioff[4];
      else
        phi -= phioff[2];
    }
  }
  if (index < 0)
    index = 0;
  if (phi < 0.0)
    phi += 2._pi;
  else if (phi > 2._pi)
    phi -= 2._pi;
  int phibin(1), unit(1);
  if (bc == HcalForward) {
    if (index < (int)(dPhiTableHF.size())) {
      unit = unitPhiHF[index];
      phibin = static_cast<int>(phi / dPhiTableHF[index]) + 1;
    }
  } else {
    if (index < (int)(dPhiTable.size())) {
      phibin = static_cast<int>(phi / dPhiTable[index]) + 1;
      unit = unitPhi[index];
    }
  }
  int iphi(phibin);
  if (unit == 2)
    iphi = (phibin - 1) * 2 + 1;
  else if (unit == 4)
    iphi = (phibin - 1) * 4 + 3;
  return iphi;
}

static std::string HcalTopology::producerTag ( )

inlinestatic

Definition at line 43 of file HcalTopology.h.

{ return "HCAL"; }

std::pair< int, int > HcalTopology::segmentBoundaries	(	const unsigned	ring,
		const unsigned	depth,
		const bool	flag = false
	)		const

returns the boundaries of the depth segmentation, so that the first result is the first segment, and the second result is the first one of the next segment. Used for calculating physical bounds.

Definition at line 1093 of file HcalTopology.cc.

References d1, HLT_FULL_cff::depth, getDepthSegmentation(), cuda_std::lower_bound(), and cuda_std::upper_bound().

                                                                                                                 {
  std::vector<int> readoutDepths;
  getDepthSegmentation(ring, readoutDepths, one);
  int d1 = std::lower_bound(readoutDepths.begin(), readoutDepths.end(), depth) - readoutDepths.begin();
  int d2 = std::upper_bound(readoutDepths.begin(), readoutDepths.end(), depth) - readoutDepths.begin();
  return std::pair<int, int>(d1, d2);
}

void HcalTopology::setDepthSegmentation	(	const unsigned	ring,
		const std::vector< int > &	readoutDepths,
		const bool	flag
	)

Definition at line 1085 of file HcalTopology.cc.

References depthSegmentation_, depthSegmentationOne_, and relativeConstraints::ring.

Referenced by HcalTopology().

                                                                                                                {
  if (one) {
    depthSegmentationOne_[ring] = readoutDepths;
  } else {
    depthSegmentation_[ring] = readoutDepths;
  }
}

int HcalTopology::singlePhiBins ( ) const

inline

Definition at line 105 of file HcalTopology.h.

References singlePhiBins_.

{ return singlePhiBins_; }

std::vector< DetId > HcalTopology::south ( const DetId &  id ) const

overridevirtual

Get the neighbors of the given cell in south direction

Implements CaloSubdetectorTopology.

Definition at line 436 of file HcalTopology.cc.

References HcalDetId::changeForm(), decIPhi(), HcalDetId::oldFormat(), and DetId::rawId().

Referenced by spr::newHCALIdNS().

                                                          {
  std::vector<DetId> vNeighborsDetId;
  HcalDetId neighbor;
  if (decIPhi(HcalDetId(id), neighbor)) {
    if (neighbor.oldFormat())
      neighbor.changeForm();
    vNeighborsDetId.emplace_back(DetId(neighbor.rawId()));
  }
  return vNeighborsDetId;
}

int HcalTopology::topoVersion ( ) const

overridevirtual

return a version which identifies the given topology

Reimplemented from CaloSubdetectorTopology.

Definition at line 1689 of file HcalTopology.cc.

References topoVersion_.

Referenced by HcalCondObjectContainerBase::HcalCondObjectContainerBase(), and HcalCondObjectContainerBase::setTopo().

{ return topoVersion_; }

HcalTopologyMode::TriggerMode HcalTopology::triggerMode ( ) const

inline

Definition at line 35 of file HcalTopology.h.

References triggerMode_.

Referenced by CaloTPGTranscoderULUT::getOutputLUTSize(), HcalTrigTowerGeometry::HcalTrigTowerGeometry(), CaloTPGTranscoderULUT::isOnlyQIE11(), and HcaluLUTTPGCoder::update().

{ return triggerMode_; }

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

inline

Definition at line 168 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::unmergeDepthDetId().

Referenced by HcalCollapseAnalyzer::analyze(), HcalGeometry::getBackPosition(), and HcalGeometry::getCorners().

{ hcons_->unmergeDepthDetId(id, ids); }

std::vector< DetId > HcalTopology::up ( const DetId &  id ) const

overridevirtual

Get the neighbors of the given cell in up direction (outward)

Implements CaloSubdetectorTopology.

Definition at line 447 of file HcalTopology.cc.

References HcalDetId::changeForm(), gpuClustering::id, incrementDepth(), and HcalDetId::oldFormat().

Referenced by spr::matrixHCALIdsDepth().

                                                       {
  HcalDetId neighbor = id;
  std::vector<DetId> vNeighborsDetId;
  if (incrementDepth(neighbor)) {
    if (neighbor.oldFormat())
      neighbor.changeForm();
    vNeighborsDetId.emplace_back(neighbor);
  }
  return vNeighborsDetId;
}

bool HcalTopology::valid ( const DetId &  id ) const

overridevirtual

Is this a valid cell id?

Reimplemented from CaloSubdetectorTopology.

Definition at line 225 of file HcalTopology.cc.

References cms::cuda::assert(), DetId::Hcal, and validHcal().

Referenced by RecAnalyzerMinbias::beginRun(), decAIEta(), decIPhi(), HcalGeometry::getCells(), incAIEta(), incIPhi(), hcalCalib::makeTextFile(), CaloMiscalibMapHcal::prefillMap(), hcalCalib::ReadPhiSymCor(), CaloTPGTranscoderULUT::setup(), HcaluLUTTPGCoder::update(), and HcaluLUTTPGCoder::updateXML().

                                              {
  assert(id.det() == DetId::Hcal);
  return validHcal(id);
}

bool HcalTopology::validCalib

(

const HcalCalibDetId &

id

)

const

Definition at line 277 of file HcalTopology.cc.

References funct::abs(), HcalCalibDetId::calibFlavor(), HcalCalibDetId::CalibrationBox, HcalCalibDetId::CastorRadFacility, HcalCalibDetId::cboxChannel(), officialStyle::chan, chanCalibHB_, chanCalibHE1_, chanCalibHE2_, chanCalibHF1_, chanCalibHF2_, chanCalibHO_, chanCalibHOs_, etaCalibHB_, etaCalibHBX_, etaCalibHE_, etaCalibHEX_, etaCalibHF_, etaCalibHO_, etaCalibHOX_, spr::find(), HcalCalibDetId::HBX, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalCalibDetId::hcalSubdet(), HcalCalibDetId::HEX, HcalCalibDetId::HOCrosstalk, HcalCalibDetId::ieta(), HcalCalibDetId::iphi(), HcalCalibDetId::LASERMON, maxPhi_, minPhi_, nchanCalibHB_, nchanCalibHE1_, nchanCalibHF1_, nchanCalibHO_, nEtaCalibHB_, nEtaCalibHE_, nEtaCalibHEX_, nEtaCalibHF_, nEtaCalibHO_, nEtaCalibHOX_, convertSQLiteXML::ok, and HcalCalibDetId::uMNqie.

Referenced by validRaw().

                                                             {
  bool ok(false);
  if (tid.calibFlavor() == HcalCalibDetId::CalibrationBox) {
    HcalSubdetector subdet = tid.hcalSubdet();
    int ieta = tid.ieta();
    int chan = tid.cboxChannel();
    unsigned int iphi = static_cast<unsigned int>(tid.iphi());
    if (subdet == HcalBarrel) {
      if ((std::find(etaCalibHB_, etaCalibHB_ + nEtaCalibHB_, ieta) != (etaCalibHB_ + nEtaCalibHB_)) &&
          (std::find(chanCalibHB_, chanCalibHB_ + nchanCalibHB_, chan) != (chanCalibHB_ + nchanCalibHB_)) &&
          (iphi >= minPhi_) && (iphi <= maxPhi_))
        ok = true;
    } else if (subdet == HcalEndcap) {
      if ((std::find(etaCalibHE_, etaCalibHE_ + nEtaCalibHE_, ieta) != (etaCalibHE_ + nEtaCalibHE_)) &&
          (std::find(chanCalibHE1_, chanCalibHE1_ + nchanCalibHE1_, chan) != (chanCalibHE1_ + nchanCalibHE1_) ||
           (chan == chanCalibHE2_)) &&
          (iphi >= minPhi_) && (iphi <= maxPhi_))
        ok = true;
    } else if (subdet == HcalForward) {
      if ((std::find(etaCalibHF_, etaCalibHF_ + nEtaCalibHF_, ieta) != (etaCalibHF_ + nEtaCalibHF_)) &&
          (std::find(chanCalibHF1_, chanCalibHF1_ + nchanCalibHF1_, chan) != (chanCalibHF1_ + nchanCalibHF1_) ||
           (chan == chanCalibHF2_)) &&
          (iphi >= minPhi_) && (iphi <= maxPhi_))
        ok = true;
    } else if (subdet == HcalOuter) {
      if ((std::find(etaCalibHO_, etaCalibHO_ + nEtaCalibHO_, ieta) != (etaCalibHO_ + nEtaCalibHO_)) &&
          (std::find(chanCalibHO_, chanCalibHO_ + nchanCalibHO_, chan) != (chanCalibHO_ + nchanCalibHO_) ||
           (chan == chanCalibHOs_)) &&
          (iphi >= minPhi_) && (iphi <= maxPhi_))
        ok = true;
    }
  } else if (tid.calibFlavor() == HcalCalibDetId::HOCrosstalk) {
    int ieta = std::abs(tid.ieta());
    unsigned int iphi = static_cast<unsigned int>(tid.iphi());
    if ((std::find(etaCalibHOX_, etaCalibHOX_ + nEtaCalibHOX_, ieta) != (etaCalibHOX_ + nEtaCalibHOX_)) &&
        (iphi >= minPhi_) && (iphi <= maxPhi_))
      ok = true;
  } else if (tid.calibFlavor() == HcalCalibDetId::HBX) {
    int ieta = std::abs(tid.ieta());
    unsigned int iphi = static_cast<unsigned int>(tid.iphi());
    if ((ieta == etaCalibHBX_) && (iphi >= minPhi_) && (iphi <= maxPhi_))
      ok = true;
  } else if (tid.calibFlavor() == HcalCalibDetId::HEX) {
    int ieta = std::abs(tid.ieta());
    unsigned int iphi = static_cast<unsigned int>(tid.iphi());
    if ((std::find(etaCalibHEX_, etaCalibHEX_ + nEtaCalibHEX_, ieta) != (etaCalibHEX_ + nEtaCalibHEX_)) &&
        (iphi >= minPhi_) && (iphi <= maxPhi_))
      ok = true;
  } else if ((tid.calibFlavor() == HcalCalibDetId::uMNqie) || (tid.calibFlavor() == HcalCalibDetId::LASERMON) ||
             (tid.calibFlavor() == HcalCalibDetId::CastorRadFacility)) {
    ok = true;
  }
  return ok;
}

bool HcalTopology::validDetId	(	HcalSubdetector	subdet,
		int	ieta,
		int	iphi,
		int	depth
	)		const

Definition at line 235 of file HcalTopology.cc.

References validHcal().

                                                                                         {
  return validHcal(HcalDetId(subdet, ieta, iphi, depth));
}

bool HcalTopology::validDetIdPreLS1 ( const HcalDetId &  id ) const

private

Is this a valid cell id, ignoring the exclusion list

Basic rules used to derive this code:

HB has 72 towers in iphi. Ieta 1-14 have depth=1, Ieta 15-16 have depth=1 or 2.

HE ieta=16-20 have 72 towers in iphi ieta=21-29 have 36 towers in iphi ieta=16 is depth 3 only ieta=17 is depth 1 only ieta=18-26 & 29 have depth 1 and 2 ieta=27-28 has depth 1-3

HF ieta=29-39 have 36 in iphi ieta=40-41 have 18 in iphi (71,3,7,11...) all have two depths

HO has 15 towers in ieta and 72 in iphi and depth = 4 (one value)

At H2:

HE ieta 17 is two depths HE ieta 22- have 36 towers in iphi (starts one higher) HE ieta 24- has three depths

Definition at line 537 of file HcalTopology.cc.

References HLT_FULL_cff::depth, firstHERing(), firstHFQuadPhiRing(), firstHFRing(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, IPHI_MAX, lastHFRing(), and sd.

                                                             {
  const HcalSubdetector sd(id.subdet());
  const int ie(id.ietaAbs());
  const int ip(id.iphi());
  const int dp(id.depth());

  return ((ip >= 1) && (ip <= IPHI_MAX) && (dp >= 1) && (ie >= 1) &&
          (((sd == HcalBarrel) && (((ie <= 14) && (dp == 1)) || (((ie == 15) || (ie == 16)) && (dp <= 2)))) ||
           ((sd == HcalEndcap) &&
            (((ie == firstHERing()) && (dp == 3)) || ((ie == 17) && (dp == 1)) ||
             ((ie >= 18) && (ie <= 20) && (dp <= 2)) || ((ie >= 21) && (ie <= 26) && (dp <= 2) && (ip % 2 == 1)) ||
             ((ie >= 27) && (ie <= 28) && (dp <= 3) && (ip % 2 == 1)) || ((ie == 29) && (dp <= 2) && (ip % 2 == 1)))) ||
           ((sd == HcalOuter) && (ie <= 15) && (dp == 4)) ||
           ((sd == HcalForward) && (dp <= 2) &&
            (((ie >= firstHFRing()) && (ie < firstHFQuadPhiRing()) && (ip % 2 == 1)) ||
             ((ie >= firstHFQuadPhiRing()) && (ie <= lastHFRing()) && (ip % 4 == 3))))));
}

bool HcalTopology::validHcal

(

const HcalDetId &

id

)

const

Is this a valid cell id?

Definition at line 230 of file HcalTopology.cc.

References isExcluded(), and validRaw().

Referenced by HcalDigitizer::accumulateCaloHits(), HcalTestSimHitID::analyze(), HcalHBHEMuonAnalyzer::analyze(), HcalHBHEMuonHighEtaAnalyzer::analyzeTracks(), CaloTowersCreationAlgo::hadShwrPos(), AlCaHcalHBHEMuonProducer::produce(), valid(), validDetId(), and validHcal().

                                                      {
  // check the raw rules
  return validRaw(id) && !isExcluded(id);
}

bool HcalTopology::validHcal	(	const HcalDetId &	id,
		const unsigned int	flag
	)		const

Is this a valid cell in context of Plan1

Definition at line 332 of file HcalTopology.cc.

References hcons_, HcalDDDRecConstants::isPlan1MergedId(), HcalDDDRecConstants::isPlan1ToBeMergedId(), or, and validHcal().

                                                                               {
  /* original logic show here because condensed form below is rather terse
  // check the raw rules
  bool ok = validHcal(id);
  if (flag == 0) { // This is all what is needed
  } else if (flag == 1) { // See if it is in the to be merged list and merged list
    if (hcons_->isPlan1MergedId(id))          ok = true;
    else if (hcons_->isPlan1ToBeMergedId(id)) ok = false;
  } else if (!ok) {
    ok = hcons_->isPlan1MergedId(id);
  }
  return ok;
  */
  return (flag > 0 and hcons_->isPlan1MergedId(id)) or
         ((flag != 1 or !hcons_->isPlan1ToBeMergedId(id)) and validHcal(id));
}

bool HcalTopology::validHT

(

const HcalTrigTowerDetId &

id

)

const

Definition at line 239 of file HcalTopology.cc.

References HLT_FULL_cff::depth, firstHFRing_, IPHI_MAX, lastHBRing_, maxDepthHE_, or, triggerMode_, HcalTopologyMode::TriggerMode_2009, HcalTopologyMode::TriggerMode_2017, HcalTopologyMode::TriggerMode_2018legacy, and relval_steps::version.

Referenced by CaloTPGTranscoderULUT::HTvalid(), CaloTPGTranscoderULUT::loadHCALCompress(), L1EGCrystalClusterEmulatorProducer::produce(), and validRaw().

                                                             {
  if (id.iphi() < 1 || id.iphi() > IPHI_MAX || id.ieta() == 0)
    return false;
  if (id.depth() != 0)
    return false;
  if (maxDepthHE_ == 0) {
    if (id.ietaAbs() > lastHBRing_ && id.ietaAbs() < firstHFRing_)
      return false;
  }
  // Version 2 TPs should be for HBHE when using 1TS filter scheme
  if (id.version() == 0 or id.version() == 2) {
    if (id.ietaAbs() > 28) {
      if (triggerMode_ >= HcalTopologyMode::TriggerMode_2017)
        return false;
      if (triggerMode_ == HcalTopologyMode::TriggerMode_2018legacy)
        return false;
      if ((id.iphi() % 4) != 1)
        return false;
      if (id.ietaAbs() > 32)
        return false;
    }
  } else if (id.version() == 1) {
    if (triggerMode_ == HcalTopologyMode::TriggerMode_2009)
      return false;
    if (id.ietaAbs() < 30 || id.ietaAbs() > 41)
      return false;
    if (id.ietaAbs() > 29 && ((id.iphi() % 2) == 0))
      return false;
    if (id.ietaAbs() > 39 && ((id.iphi() % 4) != 3))
      return false;
  } else if (id.version() > 2) {
    // only versions 0, 1, and 2 are supported
    return false;
  }

  return true;
}

bool HcalTopology::validRaw ( const HcalDetId &  id ) const

private

Is this a valid cell id?

Definition at line 556 of file HcalTopology.cc.

References HLT_FULL_cff::depth, etaBinsHE_, firstHEDoublePhiRing(), firstHEQuadPhiRing(), firstHERing(), firstHFQuadPhiRing(), firstHFRing(), HcalDDDRecConstants::getDepthEta16(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNoff(), HcalTopologyMode::H2, HcalTopologyMode::H2HE, HcalBarrel, HcalEndcap, HcalForward, HcalOther, HcalOuter, HcalTriggerTower, hcons_, mps_fire::i, IPHI_MAX, lastHBRing(), lastHERing(), lastHFRing(), lastHORing(), maxEta_, HcalDDDRecConstants::maxHFDepth(), HLT_FULL_cff::maxPhi, maxPhiHE_, mode_, convertSQLiteXML::ok, HcalTopologyMode::SLHC, validCalib(), validHT(), and ecaldqm::zside().

Referenced by decIPhi(), decrementDepth(), exclude(), incIPhi(), incrementDepth(), and validHcal().

                                                     {
  bool ok = true;
  int ieta = id.ieta();
  int aieta = id.ietaAbs();
  int depth = id.depth();
  int iphi = id.iphi();
  int zside = id.zside();
  HcalSubdetector subdet = id.subdet();
  int maxPhi = (subdet == HcalEndcap) ? maxPhiHE_ : IPHI_MAX;
  if ((ieta == 0 || iphi <= 0 || iphi > maxPhi) || aieta > maxEta_)
    ok = false;  // outer limits

  if (ok) {
    if (subdet == HcalBarrel) {
      if (mode_ == HcalTopologyMode::SLHC || mode_ == HcalTopologyMode::H2HE) {
        if ((aieta > lastHBRing()) || (depth > hcons_->getMaxDepth(0, aieta, iphi, zside)) ||
            (depth < hcons_->getMinDepth(0, aieta, iphi, zside)))
          ok = false;
      } else {
        if (aieta > lastHBRing() || depth > 2 || (aieta <= 14 && depth > 1))
          ok = false;
      }
    } else if (subdet == HcalEndcap) {
      if (mode_ == HcalTopologyMode::SLHC || mode_ == HcalTopologyMode::H2HE) {
        if ((depth > hcons_->getMaxDepth(1, aieta, iphi, zside)) ||
            (depth < hcons_->getMinDepth(1, aieta, iphi, zside)) || (aieta < firstHERing()) || (aieta > lastHERing())) {
          ok = false;
        } else {
          for (const auto& i : etaBinsHE_) {
            if (aieta == i.ieta) {
              if (aieta >= firstHEDoublePhiRing() && (iphi % 2) == 0)
                ok = false;
              if (aieta >= firstHEQuadPhiRing() && (iphi % 4) != 3)
                ok = false;
              if (aieta + 1 == hcons_->getNoff(1)) {
                if (depth < 1)
                  ok = false;
              } else {
                if (depth < i.depthStart)
                  ok = false;
              }
              break;
            }
          }
        }
      } else {
        if (depth > hcons_->getMaxDepth(1, aieta, iphi, zside) || aieta < firstHERing() || aieta > lastHERing() ||
            (aieta == firstHERing() && depth != hcons_->getDepthEta16(2, iphi, zside)) ||
            (aieta == 17 && depth != 1 && mode_ != HcalTopologyMode::H2) ||  // special case at H2
            (((aieta >= 17 && aieta < firstHETripleDepthRing()) || aieta == lastHERing()) && depth > 2) ||
            (aieta >= firstHEDoublePhiRing() && (iphi % 2) == 0))
          ok = false;
      }
    } else if (subdet == HcalOuter) {
      if (aieta > lastHORing() || iphi > IPHI_MAX || depth != 4)
        ok = false;
    } else if (subdet == HcalForward) {
      if (aieta < firstHFRing() || aieta > lastHFRing() || ((iphi % 2) == 0) ||
          (depth > hcons_->maxHFDepth(ieta, iphi)) || (aieta >= firstHFQuadPhiRing() && ((iphi + 1) % 4) != 0))
        ok = false;
    } else if (subdet == HcalTriggerTower) {
      ok = validHT(HcalTrigTowerDetId(id.rawId()));
    } else if (subdet == HcalOther) {
      ok = validCalib(HcalCalibDetId(id.rawId()));
    } else {
      ok = false;
    }
  }
  return ok;
}

std::vector< DetId > HcalTopology::west ( const DetId &  id ) const

overridevirtual

Get the neighbors of the given cell in west direction

Implements CaloSubdetectorTopology.

Definition at line 414 of file HcalTopology.cc.

References mps_fire::i, and incIEta().

Referenced by spr::newHCALIdEW().

                                                         {
  std::vector<DetId> vNeighborsDetId;
  HcalDetId neighbors[2];
  for (int i = 0; i < incIEta(HcalDetId(id), neighbors); i++) {
    if (neighbors[i].oldFormat())
      neighbors[i].changeForm();
    vNeighborsDetId.emplace_back(DetId(neighbors[i].rawId()));
  }
  return vNeighborsDetId;
}

bool HcalTopology::withSpecialRBXHBHE ( ) const

inline

Definition at line 165 of file HcalTopology.h.

References hcons_, and HcalDDDRecConstants::withSpecialRBXHBHE().

{ return hcons_->withSpecialRBXHBHE(); }

Member Data Documentation

unsigned int HcalTopology::CALIBSize_

private

Definition at line 213 of file HcalTopology.h.

Referenced by getCALIBSize(), and HcalTopology().

constexpr int HcalTopology::chanCalibHB_[nchanCalibHB_] = {0, 1, 2}

staticprivate

Definition at line 252 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::chanCalibHE1_[nchanCalibHE1_] = {0, 1, 3, 4, 5, 6}

staticprivate

Definition at line 259 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::chanCalibHE2_ = 2

staticprivate

Definition at line 284 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::chanCalibHF1_[nchanCalibHF1_] = {0, 1, 8}

staticprivate

Definition at line 266 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::chanCalibHF2_ = 9

staticprivate

Definition at line 288 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::chanCalibHO_[nchanCalibHO_] = {0, 1}

staticprivate

Definition at line 273 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::chanCalibHOs_ = 7

staticprivate

Definition at line 272 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

SegmentationMap HcalTopology::depthSegmentation_

private

Definition at line 224 of file HcalTopology.h.

Referenced by getDepthSegmentation(), and setDepthSegmentation().

SegmentationMap HcalTopology::depthSegmentationOne_

private

Definition at line 225 of file HcalTopology.h.

Referenced by getDepthSegmentation(), and setDepthSegmentation().

int HcalTopology::doublePhiBins_

private

Definition at line 203 of file HcalTopology.h.

Referenced by doublePhiBins(), and nPhiBins().

std::vector<double> HcalTopology::dPhiTable

private

Definition at line 216 of file HcalTopology.h.

Referenced by HcalTopology(), nPhiBins(), and phiBin().

std::vector<double> HcalTopology::dPhiTableHF

private

Definition at line 216 of file HcalTopology.h.

Referenced by HcalTopology(), nPhiBins(), and phiBin().

std::vector<HcalDDDRecConstants::HcalEtaBin> HcalTopology::etaBinsHB_

private

Definition at line 198 of file HcalTopology.h.

Referenced by decrementDepth(), and HcalTopology().

std::vector<HcalDDDRecConstants::HcalEtaBin> HcalTopology::etaBinsHE_

private

Definition at line 198 of file HcalTopology.h.

Referenced by decrementDepth(), HcalTopology(), and validRaw().

constexpr int HcalTopology::etaCalibHB_[nEtaCalibHB_] = {-1, 1}

staticprivate

Definition at line 253 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::etaCalibHBX_ = 16

staticprivate

Definition at line 300 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::etaCalibHE_[nEtaCalibHE_] = {-1, 1}

staticprivate

Definition at line 260 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::etaCalibHEX_[nEtaCalibHEX_] = {25, 27}

staticprivate

Definition at line 305 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::etaCalibHF_[nEtaCalibHF_] = {-1, 1}

staticprivate

Definition at line 267 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::etaCalibHO_[nEtaCalibHO_] = {-2, -1, 0, 1, 2}

staticprivate

Definition at line 274 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

constexpr int HcalTopology::etaCalibHOX_[nEtaCalibHOX_] = {4, 15}

staticprivate

Definition at line 293 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), detId2denseIdCALIB(), and validCalib().

int HcalTopology::etaHE2HF_

private

Definition at line 205 of file HcalTopology.h.

Referenced by HcalTopology(), and incAIEta().

int HcalTopology::etaHF2HE_

private

Definition at line 205 of file HcalTopology.h.

Referenced by decAIEta(), decrementDepth(), and HcalTopology().

std::vector<double> HcalTopology::etaTable

private

Definition at line 216 of file HcalTopology.h.

Referenced by etaMax(), etaRange(), etaRing(), and HcalTopology().

std::vector<double> HcalTopology::etaTableHF

private

Definition at line 216 of file HcalTopology.h.

Referenced by etaMax(), etaRange(), etaRing(), and HcalTopology().

bool HcalTopology::excludeHB_

private

Definition at line 188 of file HcalTopology.h.

Referenced by exclude(), excludeSubdetector(), and isExcluded().

bool HcalTopology::excludeHE_

private

Definition at line 188 of file HcalTopology.h.

Referenced by exclude(), excludeSubdetector(), and isExcluded().

bool HcalTopology::excludeHF_

private

Definition at line 188 of file HcalTopology.h.

Referenced by exclude(), excludeSubdetector(), and isExcluded().

bool HcalTopology::excludeHO_

private

Definition at line 188 of file HcalTopology.h.

Referenced by exclude(), excludeSubdetector(), and isExcluded().

std::vector<HcalDetId> HcalTopology::exclusionList_

private

Definition at line 187 of file HcalTopology.h.

Referenced by exclude(), and isExcluded().

int HcalTopology::firstHBRing_

private

Definition at line 193 of file HcalTopology.h.

Referenced by etaRing(), firstHBRing(), HcalTopology(), nPhiBins(), and phiBin().

int HcalTopology::firstHEDoublePhiRing_

private

Definition at line 201 of file HcalTopology.h.

Referenced by firstHEDoublePhiRing(), and HcalTopology().

int HcalTopology::firstHEQuadPhiRing_

private

Definition at line 201 of file HcalTopology.h.

Referenced by firstHEQuadPhiRing(), and HcalTopology().

int HcalTopology::firstHERing_

private

Definition at line 194 of file HcalTopology.h.

Referenced by firstHERing(), and HcalTopology().

int HcalTopology::firstHETripleDepthRing_

private

Definition at line 202 of file HcalTopology.h.

Referenced by firstHETripleDepthRing(), and HcalTopology().

int HcalTopology::firstHFQuadPhiRing_

private

Definition at line 201 of file HcalTopology.h.

Referenced by firstHFQuadPhiRing(), and HcalTopology().

int HcalTopology::firstHFRing_

private

Definition at line 195 of file HcalTopology.h.

Referenced by etaRange(), etaRing(), firstHFRing(), HcalTopology(), nPhiBins(), phiBin(), and validHT().

int HcalTopology::firstHORing_

private

Definition at line 196 of file HcalTopology.h.

Referenced by firstHORing(), and HcalTopology().

unsigned int HcalTopology::HBSize_

private

Definition at line 208 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), getHBSize(), HcalTopology(), and ncells().

const HcalDDDRecConstants* HcalTopology::hcons_

private

Definition at line 185 of file HcalTopology.h.

Referenced by dddConstants(), decrementDepth(), depthBinInformation(), getPhiZOne(), HcalTopology(), idBack(), idFront(), incrementDepth(), isBH(), mergedDepth29(), mergedDepthDetId(), mergedDepthList29(), nPhiBins(), unmergeDepthDetId(), validHcal(), validRaw(), and withSpecialRBXHBHE().

unsigned int HcalTopology::HESize_

private

Definition at line 209 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), getHESize(), HcalTopology(), and ncells().

unsigned int HcalTopology::HFSize_

private

Definition at line 211 of file HcalTopology.h.

Referenced by getHFSize(), HcalTopology(), and ncells().

unsigned int HcalTopology::HOSize_

private

Definition at line 210 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), getHOSize(), HcalTopology(), and ncells().

unsigned int HcalTopology::HTSize_

private

Definition at line 212 of file HcalTopology.h.

Referenced by getHTSize(), and HcalTopology().

constexpr unsigned int HcalTopology::kchanCalibHB_ = nchanCalibHB_ * kPhiCalibHB_

staticprivate

Definition at line 254 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kchanCalibHE1_ = nchanCalibHE1_ * kPhiCalibHE_

staticprivate

Definition at line 261 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kchanCalibHE2_ = kPhiCalibHE_

staticprivate

Definition at line 285 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kchanCalibHF1_ = nchanCalibHF1_ * kPhiCalibHF1_

staticprivate

Definition at line 268 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kchanCalibHF2_ = 1

staticprivate

Definition at line 289 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHB_ = 0

staticprivate

Definition at line 249 of file HcalTopology.h.

Referenced by detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHBX_ = kOffCalibHOX_ + nCalibHOX_

staticprivate

Definition at line 298 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHE1_ = kOffCalibHB_ + nCalibHB_

staticprivate

Definition at line 256 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHE2_ = kOffCalibHO2_ + nCalibHO2_

staticprivate

Definition at line 283 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHEX_ = kOffCalibHBX_ + nCalibHBX_

staticprivate

Definition at line 303 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHF1_ = kOffCalibHE1_ + nCalibHE1_

staticprivate

Definition at line 263 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHF2_ = kOffCalibHE2_ + nCalibHE2_

staticprivate

Definition at line 287 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHFX_ = kOffCalibHEX_ + nCalibHEX_

staticprivate

Definition at line 308 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and HcalTopology().

constexpr unsigned int HcalTopology::kOffCalibHO1_ = kOffCalibHF1_ + nCalibHF1_

staticprivate

Definition at line 270 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHO2_ = kOffCalibHO1_ + nCalibHO1_

staticprivate

Definition at line 281 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kOffCalibHOX_ = kOffCalibHF2_ + nCalibHF2_

staticprivate

Definition at line 291 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHB_ = maxPhi_ / mPhiCalibHB_

staticprivate

Definition at line 251 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHBX_ = maxPhi_ / mPhiCalibHBX_

staticprivate

Definition at line 301 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHE_ = maxPhi_ / mPhiCalibHE_

staticprivate

Definition at line 258 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHEX_ = maxPhi_ / mPhiCalibHEX_

staticprivate

Definition at line 306 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHF1_ = maxPhi_ / mPhiCalibHF1_

staticprivate

Definition at line 265 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHO0_ = maxPhi_ / mPhiCalibHO0_

staticprivate

Definition at line 277 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHO1_ = maxPhi_ / mPhiCalibHO1_

staticprivate

Definition at line 278 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::kPhiCalibHO2_ = 4 * kPhiCalibHO1_ + kPhiCalibHO0_

staticprivate

Definition at line 279 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

int HcalTopology::lastHBRing_

private

Definition at line 193 of file HcalTopology.h.

Referenced by etaMax(), HcalTopology(), lastHBHERing(), lastHBRing(), and validHT().

int HcalTopology::lastHERing_

private

Definition at line 194 of file HcalTopology.h.

Referenced by etaMax(), etaRange(), etaRing(), HcalTopology(), lastHBHERing(), and lastHERing().

int HcalTopology::lastHFRing_

private

Definition at line 195 of file HcalTopology.h.

Referenced by HcalTopology(), and lastHFRing().

int HcalTopology::lastHORing_

private

Definition at line 196 of file HcalTopology.h.

Referenced by etaMax(), HcalTopology(), and lastHORing().

int HcalTopology::maxDepthHB_

private

Definition at line 204 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), detId2denseIdHB(), HcalTopology(), maxDepth(), and maxDepthHB().

int HcalTopology::maxDepthHE_

private

Definition at line 204 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), detId2denseIdHE(), HcalTopology(), maxDepth(), maxDepthHE(), and validHT().

int HcalTopology::maxDepthHF_

private

Definition at line 204 of file HcalTopology.h.

Referenced by denseId2detId(), depthBinInformation(), detId2denseId(), detId2denseIdHF(), HcalTopology(), and maxDepth().

int HcalTopology::maxEta_

private

Definition at line 206 of file HcalTopology.h.

Referenced by HcalTopology(), and validRaw().

constexpr unsigned int HcalTopology::maxPhi_ = 72

staticprivate

Definition at line 248 of file HcalTopology.h.

Referenced by validCalib().

int HcalTopology::maxPhiHE_

private

Definition at line 206 of file HcalTopology.h.

Referenced by decIPhi(), denseId2detId(), detId2denseId(), detId2denseIdHE(), HcalTopology(), incIPhi(), and validRaw().

bool HcalTopology::mergePosition_

private

Definition at line 186 of file HcalTopology.h.

Referenced by getMergePositionFlag().

constexpr int HcalTopology::minMaxDepth_ = 4

staticprivate

Definition at line 247 of file HcalTopology.h.

Referenced by maxDepth().

constexpr unsigned int HcalTopology::minPhi_ = 1

staticprivate

Definition at line 248 of file HcalTopology.h.

Referenced by validCalib().

HcalTopologyMode::Mode HcalTopology::mode_

private

Definition at line 190 of file HcalTopology.h.

Referenced by decrementDepth(), depthBinInformation(), etaRange(), HcalTopology(), incrementDepth(), mode(), and validRaw().

constexpr unsigned int HcalTopology::mPhiCalibHB_ = 4

staticprivate

Definition at line 251 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::mPhiCalibHBX_ = 1

staticprivate

Definition at line 299 of file HcalTopology.h.

constexpr unsigned int HcalTopology::mPhiCalibHE_ = 4

staticprivate

Definition at line 258 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::mPhiCalibHEX_ = 2

staticprivate

Definition at line 304 of file HcalTopology.h.

Referenced by detId2denseIdCALIB().

constexpr unsigned int HcalTopology::mPhiCalibHF1_ = 18

staticprivate

Definition at line 265 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::mPhiCalibHO0_ = 6

staticprivate

Definition at line 276 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::mPhiCalibHO1_ = 12

staticprivate

Definition at line 276 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

constexpr unsigned int HcalTopology::mPhiCalibHOX_[nEtaCalibHOX_] = {2, 1}

staticprivate

Definition at line 294 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHB_ = kchanCalibHB_ * nEtaCalibHB_

staticprivate

Definition at line 255 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHBX_ = 2 * kPhiCalibHBX_ * nEtaCalibHBX_

staticprivate

Definition at line 302 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHE1_ = kchanCalibHE1_ * nEtaCalibHE_

staticprivate

Definition at line 262 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHE2_ = kchanCalibHE2_ * nEtaCalibHE_

staticprivate

Definition at line 286 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHEX_ = 2 * kPhiCalibHEX_ * nEtaCalibHEX_

staticprivate

Definition at line 307 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHF1_ = kchanCalibHF1_ * nEtaCalibHF_

staticprivate

Definition at line 269 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHF2_ = nEtaCalibHF_

staticprivate

Definition at line 290 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHO1_ = kPhiCalibHO2_ * nchanCalibHO_

staticprivate

Definition at line 280 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nCalibHO2_ = nEtaCalibHO_

staticprivate

Definition at line 282 of file HcalTopology.h.

constexpr int HcalTopology::nCalibHOX_ = 2 * (nPhiCalibHOX_[0] + nPhiCalibHOX_[1])

staticprivate

Definition at line 296 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nchanCalibHB_ = 3

staticprivate

Definition at line 250 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nchanCalibHE1_ = 6

staticprivate

Definition at line 257 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nchanCalibHF1_ = 3

staticprivate

Definition at line 264 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nchanCalibHO_ = 2

staticprivate

Definition at line 271 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nEtaCalibHB_ = 2

staticprivate

Definition at line 250 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nEtaCalibHBX_ = 1

staticprivate

Definition at line 299 of file HcalTopology.h.

constexpr unsigned int HcalTopology::nEtaCalibHE_ = 2

staticprivate

Definition at line 257 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nEtaCalibHEX_ = 2

staticprivate

Definition at line 304 of file HcalTopology.h.

Referenced by validCalib().

constexpr unsigned int HcalTopology::nEtaCalibHF_ = 2

staticprivate

Definition at line 264 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nEtaCalibHO_ = 5

staticprivate

Definition at line 271 of file HcalTopology.h.

Referenced by detId2denseIdCALIB(), and validCalib().

constexpr unsigned int HcalTopology::nEtaCalibHOX_ = 2

staticprivate

Definition at line 292 of file HcalTopology.h.

Referenced by validCalib().

int HcalTopology::nEtaHB_

private

Definition at line 199 of file HcalTopology.h.

Referenced by decrementDepth(), denseId2detId(), detId2denseId(), detId2denseIdHB(), and HcalTopology().

int HcalTopology::nEtaHE_

private

Definition at line 199 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), detId2denseIdHE(), etaMax(), and HcalTopology().

constexpr unsigned int HcalTopology::nPhiCalibHOX_[nEtaCalibHOX_] = {36, 72}

staticprivate

Definition at line 295 of file HcalTopology.h.

Referenced by denseId2detIdCALIB(), and detId2denseIdCALIB().

unsigned int HcalTopology::numberOfShapes_

private

Definition at line 214 of file HcalTopology.h.

Referenced by getNumberOfShapes(), and HcalTopology().

constexpr int HcalTopology::phiCalibHF2_ = 1

staticprivate

Definition at line 288 of file HcalTopology.h.

Referenced by denseId2detIdCALIB().

constexpr int HcalTopology::phiCalibHO_[nEtaCalibHO_] = {59, 47, 53, 47, 47}

staticprivate

Definition at line 275 of file HcalTopology.h.

Referenced by denseId2detIdCALIB().

constexpr int HcalTopology::phiCalibHOX1_ = 4

staticprivate

Definition at line 297 of file HcalTopology.h.

Referenced by denseId2detIdCALIB().

constexpr int HcalTopology::phiCalibHOX2_ = 12

staticprivate

Definition at line 297 of file HcalTopology.h.

Referenced by denseId2detIdCALIB().

constexpr int HcalTopology::phiCalibHOX3_ = 6

staticprivate

Definition at line 297 of file HcalTopology.h.

Referenced by denseId2detIdCALIB().

std::vector<double> HcalTopology::phioff

private

Definition at line 217 of file HcalTopology.h.

Referenced by HcalTopology(), and phiBin().

int HcalTopology::singlePhiBins_

private

Definition at line 203 of file HcalTopology.h.

Referenced by nPhiBins(), and singlePhiBins().

int HcalTopology::topoVersion_

private

Definition at line 220 of file HcalTopology.h.

Referenced by denseId2detId(), detId2denseId(), detId2denseIdHB(), detId2denseIdHE(), detId2denseIdHF(), detId2denseIdHO(), HcalTopology(), and topoVersion().

HcalTopologyMode::TriggerMode HcalTopology::triggerMode_

private

Definition at line 191 of file HcalTopology.h.

Referenced by HcalTopology(), triggerMode(), and validHT().

std::vector<int> HcalTopology::unitPhi

private

Definition at line 218 of file HcalTopology.h.

Referenced by HcalTopology(), and phiBin().

std::vector<int> HcalTopology::unitPhiHF

private

Definition at line 218 of file HcalTopology.h.

Referenced by HcalTopology(), and phiBin().