HGCalTopology Class Reference

#include <HGCalTopology.h>

Inheritance diagram for HGCalTopology:

Classes
struct	DecodedDetId

Public Member Functions
unsigned int	allGeomModules () const
const HGCalDDDConstants &	dddConstants () const
DecodedDetId	decode (const DetId &id) const
DetId	denseId2detId (uint32_t denseId) const override
DetId::Detector	detector () const
bool	detectorType () const
virtual uint32_t	detId2denseGeomId (const DetId &id) const
uint32_t	detId2denseId (const DetId &id) const override
	Dense indexing. More...
std::vector< DetId >	down (const DetId &id) const override
std::vector< DetId >	east (const DetId &id) const override
DetId	encode (const DecodedDetId &id_) const
DecodedDetId	geomDenseId2decId (const uint32_t &hi) const
HGCalGeometryMode::GeometryMode	geomMode () const
	Geometry mode. More...
DetId	goEast (const DetId &id) const override
	move the Topology east (positive ix) More...
DetId	goNorth (const DetId &id) const override
	move the Topology north (increment iy) More...
DetId	goSouth (const DetId &id) const override
	move the Topology south (decrement iy) More...
DetId	goWest (const DetId &id) const override
	move the Topology west (negative ix) More...
	HGCalTopology (const HGCalDDDConstants &hdcons, int subdet)
	create a new Topology More...
bool	isHFNose () const
bool	maskCell (const DetId &id, int corners=3) const
std::vector< DetId >	neighbors (const DetId &id) const
std::vector< DetId >	north (const DetId &id) const override
DetId	offsetBy (const DetId startId, int nrStepsX, int nrStepsY) const
std::vector< DetId >	south (const DetId &id) const override
ForwardSubdetector	subDetector () const
DetId	switchZSide (const DetId startId) const
unsigned int	totalGeomModules () const
unsigned int	totalModules () const
std::vector< DetId >	up (const DetId &id) const override
bool	valid (const DetId &id) const override
	Is this a valid cell id. More...
bool	validHashIndex (uint32_t ix) const
std::vector< DetId >	west (const DetId &id) const override
	~HGCalTopology () override
	virtual destructor More...
Public Member Functions inherited from CaloSubdetectorTopology
	CaloSubdetectorTopology ()
	standard constructor More...
virtual DetId	denseId2detId (unsigned int) const
	return a linear packed id 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	goUp (const DetId &id) const
virtual unsigned int	ncells () const
	return a count of valid cells (for dense indexing use) More...
virtual int	topoVersion () const
	return a version which identifies the given topology More...
virtual	~CaloSubdetectorTopology ()
	virtual destructor More...

Static Public Attributes
static const int	subSectors_ = 2
	Use subSector in square mode as wafer type in hexagon mode. More...

Private Member Functions
void	addHGCSCintillatorId (std::vector< DetId > &ids, int zside, int type, int lay, int iradius, int iphi) const
	add DetId of Scintillator and Silicon type if valid More...
void	addHGCSiliconId (std::vector< DetId > &ids, int det, int zside, int type, int lay, int waferU, int waferV, int cellU, int cellV) const
DetId	changeXY (const DetId &id, int nrStepsX, int nrStepsY) const
	move the nagivator along x, y More...
DetId	changeZ (const DetId &id, int nrStepsZ) const
	move the nagivator along z More...

Private Attributes
int	cellMax_
int	cells_
DetId::Detector	det_
int	firstLay_
const HGCalDDDConstants &	hdcons_
int	kHGeomHalf_
int	kHGhalf_
unsigned int	kSizeForDenseIndexing
int	layers_
HGCalGeometryMode::GeometryMode	mode_
int	sectors_
ForwardSubdetector	subdet_
int	types_
int	waferMax_
int	waferOff_

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

Definition at line 12 of file HGCalTopology.h.

Constructor & Destructor Documentation

HGCalTopology::HGCalTopology	(	const HGCalDDDConstants &	hdcons,
		int	subdet
	)

create a new Topology

Definition at line 11 of file HGCalTopology.cc.

References cellMax_, cells_, det_, firstLay_, HGCalDDDConstants::firstLayer(), DetId::Forward, ForwardEmpty, HGCalDDDConstants::geomMode(), hdcons_, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, HFNose, createfilelist::int, kHGeomHalf_, kHGhalf_, kSizeForDenseIndexing, HGCalDDDConstants::layers(), layers_, HGCalDDDConstants::maxCells(), HGCalDDDConstants::maxCellUV(), mode_, HGCalDDDConstants::sectors(), sectors_, subdet_, HGCalGeometryMode::Trapezoid, types_, waferMax_, waferOff_, and HGCalDDDConstants::waferUVMax().

                                                                     : hdcons_(hdcons) {
  sectors_ = hdcons_.sectors();
  layers_ = hdcons_.layers(true);
  cells_ = hdcons_.maxCells(true);
  mode_ = hdcons_.geomMode();
  cellMax_ = hdcons_.maxCellUV();
  waferOff_ = hdcons_.waferUVMax();
  waferMax_ = 2 * waferOff_ + 1;
  kHGhalf_ = sectors_ * layers_ * cells_;
  firstLay_ = hdcons_.firstLayer();
  if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
    det_ = DetId::Forward;
    subdet_ = (ForwardSubdetector)(det);
    kHGeomHalf_ = sectors_ * layers_;
    types_ = 2;
  } else if (det == (int)(DetId::Forward)) {
    det_ = DetId::Forward;
    subdet_ = HFNose;
    kHGeomHalf_ = sectors_ * layers_;
    types_ = 3;
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    det_ = (DetId::Detector)(det);
    subdet_ = ForwardEmpty;
    kHGeomHalf_ = sectors_ * layers_ * cellMax_;
    types_ = 2;
  } else {
    det_ = (DetId::Detector)(det);
    subdet_ = ForwardEmpty;
    kHGeomHalf_ = sectors_ * layers_;
    types_ = 3;
  }
  kSizeForDenseIndexing = (unsigned int)(2 * kHGhalf_);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "HGCalTopology initialized for detector " << det_ << ":" << subdet_ << " having "
                                << sectors_ << " Sectors, " << layers_ << " Layers from " << firstLay_ << ", " << cells_
                                << " cells and total channels " << kSizeForDenseIndexing << ":" << (2 * kHGeomHalf_)
                                << std::endl;
#endif
}

HGCalTopology::~HGCalTopology ( )

inlineoverride

virtual destructor

Definition at line 18 of file HGCalTopology.h.

{}

Member Function Documentation

void HGCalTopology::addHGCSCintillatorId ( std::vector< DetId > &  ids, int  zside, int  type, int  lay, int  iradius, int  iphi  ) const

private

add DetId of Scintillator and Silicon type if valid

Definition at line 539 of file HGCalTopology.cc.

References hdcons_, triggerObjects_cff::id, and HGCalDDDConstants::isValidTrap().

Referenced by isHFNose(), and neighbors().

                                                                                      {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "addHGCSCintillatorId " << zside << ":" << type << ":" << lay << ":" << iradius
                                << ":" << iphi << " ==> Validity " << hdcons_.isValidTrap(lay, iradius, iphi);
#endif
  if (hdcons_.isValidTrap(lay, iradius, iphi)) {
    HGCScintillatorDetId id(type, lay, zside * iradius, iphi);
    ids.emplace_back(DetId(id));
  }
}

void HGCalTopology::addHGCSiliconId ( std::vector< DetId > &  ids, int  det, int  zside, int  type, int  lay, int  waferU, int  waferV, int  cellU, int  cellV  ) const

private

Definition at line 551 of file HGCalTopology.cc.

References hdcons_, triggerObjects_cff::id, and HGCalDDDConstants::isValidHex8().

Referenced by isHFNose(), and neighbors().

                                                                                                                      {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "addHGCSiliconId " << det << ":" << zside << ":" << type << ":" << lay << ":"
                                << waferU << ":" << waferV << ":" << cellU << ":" << cellV << " ==> Validity "
                                << hdcons_.isValidHex8(lay, waferU, waferV, cellU, cellV);
#endif
  if (hdcons_.isValidHex8(lay, waferU, waferV, cellU, cellV)) {
    HGCSiliconDetId id((DetId::Detector)(det), zside, type, lay, waferU, waferV, cellU, cellV);
    ids.emplace_back(DetId(id));
  }
}

unsigned int HGCalTopology::allGeomModules

(

)

const

Definition at line 51 of file HGCalTopology.cc.

References hdcons_, mode_, HGCalDDDConstants::numberCells(), HGCalGeometryMode::Trapezoid, and HGCalDDDConstants::wafers().

Referenced by HGCalGeometryLoader::build(), and totalGeomModules().

                                                 {
  return ((mode_ == HGCalGeometryMode::Trapezoid) ? (unsigned int)(2 * hdcons_.numberCells(true))
                                                  : (unsigned int)(2 * hdcons_.wafers()));
}

DetId HGCalTopology::changeXY ( const DetId &  id, int  nrStepsX, int  nrStepsY  ) const

private

move the nagivator along x, y

Definition at line 628 of file HGCalTopology.cc.

Referenced by goEast(), goNorth(), goSouth(), goWest(), isHFNose(), and offsetBy().

{ return DetId(); }

DetId HGCalTopology::changeZ ( const DetId &  id, int  nrStepsZ  ) const

private

move the nagivator along z

Definition at line 630 of file HGCalTopology.cc.

Referenced by down(), isHFNose(), and up().

{ return DetId(); }

const HGCalDDDConstants& HGCalTopology::dddConstants ( ) const

inline

Definition at line 96 of file HGCalTopology.h.

References hdcons_, offsetBy(), and switchZSide().

Referenced by hgc_digi_utils::addCellMetadata(), PFRecHitQTestThresholdInThicknessNormalizedMIPs::beginEvent(), CaloTruthAccumulator::beginLuminosityBlock(), HGCalGeometryLoader::build(), HGCDigitizer::checkPosition(), HGCalTriggerGeometryV9Imp2::detIdWaferType(), HGCalTriggerGeometryHexImp2::detIdWaferType(), HGCalTriggerGeometryHexLayerBasedImp1::detIdWaferType(), HGCalTriggerGeometryV9Imp1::detIdWaferType(), HGCalTriggerTools::eventSetup(), HGCalTriggerGeometryHexImp2::fillInvalidTriggerCells(), HGCalTriggerGeometryHexLayerBasedImp1::fillInvalidTriggerCells(), HGCalTriggerGeometryV9Imp1::fillInvalidTriggerCells(), HGCalTriggerGeometryHexImp1::fillMaps(), CaloParticleDebugger::fillSimHits(), HGCalGeometry::get8Corners(), HGCalTriggerGeometryV9Imp1::getCellsFromTriggerCell(), HGCalGeometry::getClosestCell(), HGCalGeometry::getCorners(), HGCalTriggerTools::getLayerZ(), HGCalTriggerGeometryV9Imp2::getModuleFromTriggerCell(), HGCalTriggerGeometryV9Imp1::getNeighborsFromTriggerCell(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), HGCalGeometry::getSummary(), HGCalTriggerGeometryV9Imp2::getTriggerCellsFromModule(), HGCalTriggerGeometryHexLayerBasedImp1::initialize(), HGCalTriggerGeometryV9Imp1::initialize(), HGCalTriggerGeometryV9Imp2::initialize(), HGCalTriggerGeometryHexImp2::lastTriggerLayer(), HGCalTriggerGeometryGenericMapping::lastTriggerLayer(), hgcal::RecHitTools::maskCell(), maskCell(), HGCalGeometry::neighborZ(), HGCalGeometry::newCell(), HGCalTriggerTools::sensorCellThicknessV8(), HGCalRecHitWorkerSimple::set(), HGCalRadiationMap::setGeometry(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry(), and HGCalTriggerTools::simToReco().

{ return hdcons_; }

HGCalTopology::DecodedDetId HGCalTopology::decode

(

const DetId &

id

)

const

Definition at line 564 of file HGCalTopology.cc.

References HGCalTopology::DecodedDetId::det, det_, DetId::Forward, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, HFNose, HGCalTopology::DecodedDetId::iCell1, HGCalTopology::DecodedDetId::iCell2, triggerObjects_cff::id, training_settings::idx, HGCalTopology::DecodedDetId::iLay, createfilelist::int, HGCalTopology::DecodedDetId::iSec1, HGCalTopology::DecodedDetId::iSec2, HGCalTopology::DecodedDetId::iType, mode_, subdet_, HGCalGeometryMode::Trapezoid, and HGCalTopology::DecodedDetId::zSide.

Referenced by detId2denseGeomId(), detId2denseId(), HGCalGeometry::get8Corners(), HGCalGeometry::getClosestCell(), HGCalGeometry::getCorners(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), neighbors(), HGCalGeometry::neighborZ(), HGCalGeometry::newCell(), switchZSide(), and valid().

                                                                          {
  HGCalTopology::DecodedDetId idx;
  if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
    HGCalDetId id(startId);
    idx.iCell1 = id.cell();
    idx.iCell2 = 0;
    idx.iLay = id.layer();
    idx.iSec1 = id.wafer();
    idx.iSec2 = 0;
    idx.iType = id.waferType();
    idx.zSide = id.zside();
    idx.det = id.subdetId();
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    HGCScintillatorDetId id(startId);
    idx.iCell1 = id.iphi();
    idx.iCell2 = 0;
    idx.iLay = id.layer();
    idx.iSec1 = id.ietaAbs();
    idx.iSec2 = 0;
    idx.iType = id.type();
    idx.zSide = id.zside();
    idx.det = (int)(id.subdet());
  } else if (det_ == DetId::Forward && subdet_ == ForwardSubdetector::HFNose) {
    HFNoseDetId id(startId);
    idx.iCell1 = id.cellU();
    idx.iCell2 = id.cellV();
    idx.iLay = id.layer();
    idx.iSec1 = id.waferU();
    idx.iSec2 = id.waferV();
    idx.iType = id.type();
    idx.zSide = id.zside();
    idx.det = (int)(id.subdet());
  } else {
    HGCSiliconDetId id(startId);
    idx.iCell1 = id.cellU();
    idx.iCell2 = id.cellV();
    idx.iLay = id.layer();
    idx.iSec1 = id.waferU();
    idx.iSec2 = id.waferV();
    idx.iType = id.type();
    idx.zSide = id.zside();
    idx.det = (int)(id.subdet());
  }
  return idx;
}

DetId HGCalTopology::denseId2detId ( uint32_t  denseId ) const

override

Definition at line 387 of file HGCalTopology.cc.

References cellMax_, encode(), firstLay_, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, triggerObjects_cff::id, createfilelist::int, kHGhalf_, layers_, mode_, sectors_, HGCalGeometryMode::Trapezoid, types_, validHashIndex(), waferMax_, and waferOff_.

Referenced by geomMode().

                                                    {
  HGCalTopology::DecodedDetId id;
  if (validHashIndex(hi)) {
    id.zSide = ((int)(hi) < kHGhalf_ ? -1 : 1);
    int di = ((int)(hi) % kHGhalf_);
    if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
      int type = (di % types_);
      id.iType = (type == 0 ? -1 : 1);
      id.iSec1 = (((di - type) / types_) % sectors_);
      id.iLay = (((((di - type) / types_) - id.iSec1 + 1) / sectors_) % layers_ + 1);
      id.iCell1 = (((((di - type) / types_) - id.iSec1 + 1) / sectors_ - id.iLay + 1) / layers_ + 1);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "Input Hex " << hi << " o/p " << id.zSide << ":" << id.iLay << ":" << id.iType
                                    << ":" << id.iSec1 << ":" << id.iCell1;
#endif
    } else if (mode_ == HGCalGeometryMode::Trapezoid) {
      int type = (di % types_);
      id.iType = type;
      id.iSec1 = (((di - type) / types_) % sectors_) + 1;
      id.iLay = (((((di - type) / types_) - id.iSec1 + 1) / sectors_) % layers_ + firstLay_);
      id.iCell1 = (((((di - type) / types_) - id.iSec1 + 1) / sectors_ - id.iLay + firstLay_) / layers_ + 1);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "Input Trap " << hi << " o/p " << id.zSide << ":" << id.iLay << ":" << id.iType
                                    << ":" << id.iSec1 << ":" << id.iCell1;
#endif
    } else {
      int type = (di % types_);
      id.iType = type;
      di = (di - type) / types_;
      id.iSec2 = (di % waferMax_) - waferOff_;
      di = (di - id.iSec2 - waferOff_) / waferMax_;
      id.iSec1 = (di % waferMax_) - waferOff_;
      di = (di - id.iSec1 - waferOff_) / waferMax_;
      id.iLay = (di % layers_) + 1;
      di = (di - id.iLay + 1) / layers_;
      id.iCell2 = (di % cellMax_);
      id.iCell1 = (di - id.iCell2) / cellMax_;
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "Input Hex8 " << hi << " o/p " << id.zSide << ":" << id.iLay << ":" << id.iType
                                    << ":" << id.iSec1 << ":" << id.iSec2 << ":" << id.iCell1 << ":" << id.iCell2;
#endif
    }
  }
  return encode(id);
}

DetId::Detector HGCalTopology::detector ( ) const

inline

Definition at line 115 of file HGCalTopology.h.

References det_.

Referenced by HGCalGeometryLoader::build().

{ return det_; }

bool HGCalTopology::detectorType ( ) const

inline

Definition at line 117 of file HGCalTopology.h.

Referenced by HGCalHitValidation::analyzeHGCalSimHit(), and HGCalTriggerTools::simToReco().

{ return false; }

uint32_t HGCalTopology::detId2denseGeomId ( const DetId &  id ) const

virtual

Definition at line 433 of file HGCalTopology.cc.

References cellMax_, decode(), firstLay_, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, training_settings::idx, kHGeomHalf_, layers_, mode_, sectors_, HGCalGeometryMode::Trapezoid, waferMax_, and waferOff_.

Referenced by geomMode(), HGCalGeometry::getGeometry(), HGCalGeometry::getSummary(), HGCalGeometry::indexFor(), HGCalGeometry::newCell(), and HGCalGeometry::present().

                                                                 {
  HGCalTopology::DecodedDetId id = decode(idin);
  uint32_t idx;
  if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
    idx = (uint32_t)(((id.zSide > 0) ? kHGeomHalf_ : 0) + (id.iLay - 1) * sectors_ + id.iSec1);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Geom Hex I/P " << id.zSide << ":" << id.iLay << ":" << id.iSec1 << ":" << id.iType
                                  << " Constants " << kHGeomHalf_ << ":" << layers_ << ":" << sectors_ << " o/p "
                                  << idx;
#endif
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    idx = (uint32_t)(((id.zSide > 0) ? kHGeomHalf_ : 0) +
                     (((id.iLay - firstLay_) * sectors_ + id.iSec1 - 1) * cellMax_ + id.iCell1 - 1));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Geom Trap I/P " << id.zSide << ":" << id.iLay << ":" << id.iSec1 << ":"
                                  << id.iCell1 << ":" << id.iType << " Constants " << kHGeomHalf_ << ":" << layers_
                                  << ":" << firstLay_ << ":" << sectors_ << ":" << cellMax_ << " o/p " << idx;
#endif
  } else {
    idx = (uint32_t)(((id.zSide > 0) ? kHGeomHalf_ : 0) +
                     (((id.iLay - 1) * waferMax_ + id.iSec1 + waferOff_) * waferMax_ + id.iSec2 + waferOff_));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Geom Hex8 I/P " << id.zSide << ":" << id.iLay << ":" << id.iSec1 << ":"
                                  << id.iSec2 << ":" << id.iType << " Constants " << kHGeomHalf_ << ":" << layers_
                                  << ":" << waferMax_ << ":" << waferOff_ << " o/p " << idx;
#endif
  }
  return idx;
}

uint32_t HGCalTopology::detId2denseId ( const DetId &  id ) const

overridevirtual

Dense indexing.

Reimplemented from CaloSubdetectorTopology.

Definition at line 348 of file HGCalTopology.cc.

References cellMax_, decode(), firstLay_, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, training_settings::idx, kHGeomHalf_, kHGhalf_, layers_, mode_, sectors_, HGCalGeometryMode::Trapezoid, types_, waferMax_, and waferOff_.

Referenced by geomMode().

                                                             {
  HGCalTopology::DecodedDetId id = decode(idin);
  uint32_t idx;
  if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
    int type = (id.iType > 0) ? 1 : 0;
    idx = (uint32_t)(((id.zSide > 0) ? kHGhalf_ : 0) +
                     ((((id.iCell1 - 1) * layers_ + id.iLay - 1) * sectors_ + id.iSec1) * types_ + type));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Input Hex " << id.zSide << ":" << id.iLay << ":" << id.iSec1 << ":" << id.iCell1
                                  << ":" << id.iType << " Constants " << kHGeomHalf_ << ":" << layers_ << ":"
                                  << sectors_ << ":" << types_ << " o/p " << idx;
#endif
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    idx =
        (uint32_t)(((id.zSide > 0) ? kHGhalf_ : 0) +
                   ((((id.iCell1 - 1) * layers_ + id.iLay - firstLay_) * sectors_ + id.iSec1 - 1) * types_ + id.iType));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Input Trap " << id.zSide << ":" << id.iLay << ":" << id.iSec1 << ":" << id.iCell1
                                  << ":" << id.iType << " Constants " << kHGeomHalf_ << ":" << layers_ << ":"
                                  << sectors_ << ":" << types_ << " o/p " << idx;
#endif
  } else {
    idx =
        (uint32_t)(((id.zSide > 0) ? kHGhalf_ : 0) +
                   (((((id.iCell1 * cellMax_ + id.iCell2) * layers_ + id.iLay - 1) * waferMax_ + id.iSec1 + waferOff_) *
                         waferMax_ +
                     id.iSec2 + waferOff_) *
                        types_ +
                    id.iType));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Input Hex8 " << id.zSide << ":" << id.iLay << ":" << id.iSec1 << ":" << id.iSec2
                                  << ":" << id.iCell1 << ":" << id.iCell2 << ":" << id.iType << " Constants "
                                  << kHGeomHalf_ << ":" << cellMax_ << ":" << layers_ << ":" << waferMax_ << ":"
                                  << waferOff_ << ":" << types_ << " o/p " << idx;
#endif
  }
  return idx;
}

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

inlineoverridevirtual

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

Implements CaloSubdetectorTopology.

Definition at line 68 of file HGCalTopology.h.

References changeZ(), and neighbors().

                                                        {
    DetId nextId = changeZ(id, -1);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }

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

inlineoverridevirtual

Get the neighbors of the given cell in east direction

Implements CaloSubdetectorTopology.

Definition at line 42 of file HGCalTopology.h.

References goEast().

                                                        {
    DetId nextId = goEast(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }

DetId HGCalTopology::encode

(

const DecodedDetId &

id_

)

const

Definition at line 610 of file HGCalTopology.cc.

References HGCalTopology::DecodedDetId::det, det_, DetId::Forward, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, HFNose, HGCalTopology::DecodedDetId::iCell1, HGCalTopology::DecodedDetId::iCell2, triggerObjects_cff::id, HGCalTopology::DecodedDetId::iLay, HGCalTopology::DecodedDetId::iSec1, HGCalTopology::DecodedDetId::iSec2, HGCalTopology::DecodedDetId::iType, mode_, DetId::rawId(), subdet_, HGCalGeometryMode::Trapezoid, and HGCalTopology::DecodedDetId::zSide.

Referenced by denseId2detId(), HGCalGeometry::getClosestCell(), HGCalGeometry::newCell(), and switchZSide().

                                                                      {
  DetId id;
  if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
    id =
        HGCalDetId((ForwardSubdetector)(idx.det), idx.zSide, idx.iLay, ((idx.iType > 0) ? 1 : 0), idx.iSec1, idx.iCell1)
            .rawId();
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    id = HGCScintillatorDetId(idx.iType, idx.iLay, idx.zSide * idx.iSec1, idx.iCell1).rawId();
  } else if (det_ == DetId::Forward && subdet_ == ForwardSubdetector::HFNose) {
    id = HFNoseDetId(idx.zSide, idx.iType, idx.iLay, idx.iSec1, idx.iSec2, idx.iCell1, idx.iCell2).rawId();
  } else {
    id = HGCSiliconDetId(
             (DetId::Detector)(idx.det), idx.zSide, idx.iType, idx.iLay, idx.iSec1, idx.iSec2, idx.iCell1, idx.iCell2)
             .rawId();
  }
  return id;
}

HGCalTopology::DecodedDetId HGCalTopology::geomDenseId2decId

(

const uint32_t &

hi

)

const

Definition at line 498 of file HGCalTopology.cc.

References cellMax_, firstLay_, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, triggerObjects_cff::id, createfilelist::int, kHGeomHalf_, layers_, mode_, sectors_, totalGeomModules(), HGCalGeometryMode::Trapezoid, waferMax_, and waferOff_.

                                                                                   {
  HGCalTopology::DecodedDetId id;
  if (hi < totalGeomModules()) {
    id.zSide = ((int)(hi) < kHGeomHalf_ ? -1 : 1);
    int di = ((int)(hi) % kHGeomHalf_);
    if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
      id.iSec1 = (di % sectors_);
      di = (di - id.iSec1) / sectors_;
      id.iLay = (di % layers_) + 1;
      id.iType = ((di - id.iLay + 1) / layers_ == 0) ? -1 : 1;
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "Geom Hex I/P " << hi << " O/P " << id.zSide << ":" << id.iType << ":" << id.iLay
                                    << ":" << id.iSec1;
#endif
    } else if (mode_ == HGCalGeometryMode::Trapezoid) {
      id.iCell1 = (di % cellMax_) + 1;
      di = (di - id.iCell1 + 1) / cellMax_;
      id.iSec1 = (di % sectors_) + 1;
      di = (di - id.iSec1 + 1) / sectors_;
      id.iLay = (di % layers_) + firstLay_;
      id.iType = (di - id.iLay + firstLay_) / layers_;
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "Geom Trap I/P " << hi << " O/P " << id.zSide << ":" << id.iType << ":"
                                    << id.iLay << ":" << id.iSec1 << ":" << id.iCell1;
#endif
    } else {
      id.iSec2 = (di % waferMax_) - waferOff_;
      di = (di - id.iSec2 - waferOff_) / waferMax_;
      id.iSec1 = (di % waferMax_) - waferOff_;
      di = (di - id.iSec1 - waferOff_) / waferMax_;
      id.iLay = (di % layers_) + 1;
      id.iType = (di - id.iLay + 1) / layers_;
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "Geom Hex8 I/P " << hi << " O/P " << id.zSide << ":" << id.iType << ":"
                                    << id.iLay << ":" << id.iSec1 << ":" << id.iSec2;
#endif
    }
  }
  return id;
}

HGCalGeometryMode::GeometryMode HGCalTopology::geomMode ( ) const

inline

Geometry mode.

Definition at line 79 of file HGCalTopology.h.

References denseId2detId(), detId2denseGeomId(), detId2denseId(), mode_, and valid().

Referenced by HGCalRecHitValidation::analyze(), HGCalDigiValidation::analyze(), and HGCalGeometryLoader::build().

{ return mode_; }

DetId HGCalTopology::goEast ( const DetId &  id ) const

inlineoverridevirtual

move the Topology east (positive ix)

Reimplemented from CaloSubdetectorTopology.

Definition at line 41 of file HGCalTopology.h.

References changeXY().

Referenced by east().

{ return changeXY(id, +1, 0); }

DetId HGCalTopology::goNorth ( const DetId &  id ) const

inlineoverridevirtual

move the Topology north (increment iy)

Reimplemented from CaloSubdetectorTopology.

Definition at line 21 of file HGCalTopology.h.

References changeXY().

Referenced by north().

{ return changeXY(id, 0, +1); }

DetId HGCalTopology::goSouth ( const DetId &  id ) const

inlineoverridevirtual

move the Topology south (decrement iy)

Reimplemented from CaloSubdetectorTopology.

Definition at line 31 of file HGCalTopology.h.

References changeXY().

Referenced by south().

{ return changeXY(id, 0, -1); }

DetId HGCalTopology::goWest ( const DetId &  id ) const

inlineoverridevirtual

move the Topology west (negative ix)

Reimplemented from CaloSubdetectorTopology.

Definition at line 51 of file HGCalTopology.h.

References changeXY().

Referenced by west().

{ return changeXY(id, -1, 0); }

bool HGCalTopology::isHFNose ( ) const

inline

Definition at line 118 of file HGCalTopology.h.

References addHGCSCintillatorId(), addHGCSiliconId(), changeXY(), changeZ(), HGCalTopology::DecodedDetId::det, det_, DetId::Forward, HFNose, photons_cff::ids, LEDCalibrationChannels::iphi, subdet_, and ecaldqm::zside().

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getGeometryDetId(), HGCalGeometry::getSummary(), and HGCalGeometry::newCell().

                        {
    return (((det_ == DetId::Forward) && (subdet_ == ForwardSubdetector::HFNose)) ? true : false);
  }

bool HGCalTopology::maskCell ( const DetId &  id, int  corners = 3  ) const

inline

Definition at line 94 of file HGCalTopology.h.

References dddConstants(), and HGCalDDDConstants::maskCell().

{ return dddConstants().maskCell(id, corners); }

std::vector< DetId > HGCalTopology::neighbors

(

const DetId &

id

)

const

Definition at line 56 of file HGCalTopology.cc.

References addHGCSCintillatorId(), addHGCSiliconId(), HGCalDDDConstants::BottomCorner, HGCalDDDConstants::BottomLeftCorner, HGCalDDDConstants::BottomLeftEdge, HGCalDDDConstants::BottomRightCorner, HGCalDDDConstants::BottomRightEdge, HGCalDDDConstants::cellType(), HGCalDDDConstants::CentralType, decode(), HGCalDDDConstants::getTypeHex(), HGCalDDDConstants::getUVMax(), hdcons_, HGCalGeometryMode::Hexagon8, HGCalGeometryMode::Hexagon8Full, photons_cff::ids, HGCalDDDConstants::LeftEdge, mode_, HGCalDDDConstants::modifyUV(), N, HGCalDDDConstants::RightEdge, RandomServiceHelper::t1, RandomServiceHelper::t2, HGCalDDDConstants::TopCorner, HGCalDDDConstants::TopLeftCorner, HGCalDDDConstants::TopLeftEdge, HGCalDDDConstants::TopRightCorner, HGCalDDDConstants::TopRightEdge, HGCalGeometryMode::Trapezoid, testProducerWithPsetDescEmpty_cfi::u1, and MetAnalyzer::u2.

Referenced by down().

                                                                 {
  std::vector<DetId> ids;
  HGCalTopology::DecodedDetId id = decode(idin);
  if ((mode_ == HGCalGeometryMode::Hexagon8) || (mode_ == HGCalGeometryMode::Hexagon8Full)) {
    HGCalDDDConstants::CellType celltype = hdcons_.cellType(id.iType, id.iCell1, id.iCell2);
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Type:WaferU:WaferV " << id.iType << ":" << id.iCell1 << ":" << id.iCell2
                                  << " CellType "
                                  << static_cast<std::underlying_type<HGCalDDDConstants::CellType>::type>(celltype);
#endif
    switch (celltype) {
      case (HGCalDDDConstants::CellType::CentralType): {
        // cell within the wafer
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 0";
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::BottomLeftEdge): {
        // bottom left edge
        int wu1(id.iSec1), wv1(id.iSec2 - 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 1 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 2 * N1 - 1, v1 + N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 2 * N1 - 1, v1 + N1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::LeftEdge): {
        // left edege
        int wu1(id.iSec1 + 1), wv1(id.iSec2);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int u1 = hdcons_.modifyUV(id.iCell1, id.iType, t1);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 2 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << u1;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 + N1, 2 * N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 + N1 - 1, 2 * N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::TopLeftEdge): {
        // top left edge
        int wu1(id.iSec1 + 1), wv1(id.iSec2 + 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 3 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, v1 + 1, v1 + N1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, v1, v1 + N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::TopRightEdge): {
        // top right edge
        int wu1(id.iSec1), wv1(id.iSec2 + 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 4 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 0, v1 - N1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 0, v1 - N1 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::RightEdge): {
        // right edge
        int wu1(id.iSec1 - 1), wv1(id.iSec2);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int u1 = hdcons_.modifyUV(id.iCell1, id.iType, t1);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 5 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << u1;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 - N1, 0);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 - N1 + 1, 0);
        break;
      }
      case (HGCalDDDConstants::CellType::BottomRightEdge): {
        // bottom right edge
        int wu1(id.iSec1 - 1), wv1(id.iSec2 - 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int u1 = hdcons_.modifyUV(id.iCell1, id.iType, t1);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 6 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << u1;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 + N1 - 1, u1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 + N1, u1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::BottomCorner): {
        // bottom corner
        int wu1(id.iSec1), wv1(id.iSec2 - 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
        int wu2(id.iSec1 - 1), wv2(id.iSec2 - 1);
        int t2 = hdcons_.getTypeHex(id.iLay, wu2, wv2);
        int N2 = hdcons_.getUVMax(t2);
        int u2 = hdcons_.modifyUV(id.iCell1, id.iType, t2);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 11 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1 << ":" << t2
                                      << ":" << N2 << ":" << u2;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 2 * N1 - 1, v1 + N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 2 * N1 - 1, v1 + N1);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, u2 + N2, u2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::BottomLeftCorner): {
        // bottom left corner
        int wu1(id.iSec1 + 1), wv1(id.iSec2);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int u1 = hdcons_.modifyUV(id.iCell1, id.iType, t1);
        int wu2(id.iSec1), wv2(id.iSec2 - 1);
        int t2 = hdcons_.getTypeHex(id.iLay, wu2, wv2);
        int N2 = hdcons_.getUVMax(t2);
        int v2 = hdcons_.modifyUV(id.iCell2, id.iType, t2);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 12 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << u1 << ":" << t2
                                      << ":" << N2 << ":" << v2;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 + N1, 2 * N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, 2 * N2 - 1, v2 + N2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, 2 * N2 - 1, v2 + N2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::TopLeftCorner): {
        // top left corner
        int wu1(id.iSec1 + 1), wv1(id.iSec2 + 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
        int wu2(id.iSec1 + 1), wv2(id.iSec2);
        int t2 = hdcons_.getTypeHex(id.iLay, wu2, wv2);
        int N2 = hdcons_.getUVMax(t2);
        int u2 = hdcons_.modifyUV(id.iCell1, id.iType, t2);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 13 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1 << ":" << t2
                                      << ":" << N2 << ":" << u2;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, v1 + 1, N1 + v1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, v1, N1 + v1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, u2 + N2 - 1, 2 * N2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::TopCorner): {
        // top corner
        int wu1(id.iSec1 + 1), wv1(id.iSec2 + 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
        int wu2(id.iSec1), wv2(id.iSec2 + 1);
        int t2 = hdcons_.getTypeHex(id.iLay, wu2, wv2);
        int N2 = hdcons_.getUVMax(t2);
        int v2 = hdcons_.modifyUV(id.iCell2, id.iType, t2);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 14 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1 << ":" << t2
                                      << ":" << N2 << ":" << v2;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, v1 + 1, v1 + N1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, v1, v1 + N1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 + 1);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, 0, v2 - N2 + 1);
        break;
      }
      case (HGCalDDDConstants::CellType::TopRightCorner): {
        // top right corner
        int wu1(id.iSec1), wv1(id.iSec2 + 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int v1 = hdcons_.modifyUV(id.iCell2, id.iType, t1);
        int wu2(id.iSec1 - 1), wv2(id.iSec2);
        int t2 = hdcons_.getTypeHex(id.iLay, wu2, wv2);
        int N2 = hdcons_.getUVMax(t2);
        int u2 = hdcons_.modifyUV(id.iCell1, id.iType, t2);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 15 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << v1 << ":" << t2
                                      << ":" << N2 << ":" << u2;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, 0, v1 - N1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, u2 - N2, 0);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, u2 - N2 + 1, 0);
        break;
      }
      case (HGCalDDDConstants::CellType::BottomRightCorner): {
        // bottom right corner
        int wu1(id.iSec1 - 1), wv1(id.iSec2 - 1);
        int t1 = hdcons_.getTypeHex(id.iLay, wu1, wv1);
        int N1 = hdcons_.getUVMax(t1);
        int u1 = hdcons_.modifyUV(id.iCell1, id.iType, t1);
        int wu2(id.iSec1 - 1), wv2(id.iSec2);
        int t2 = hdcons_.getTypeHex(id.iLay, wu2, wv2);
        int N2 = hdcons_.getUVMax(t2);
        int u2 = hdcons_.modifyUV(id.iCell1, id.iType, t2);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Cell Type 16 "
                                      << ":" << wu1 << ":" << wv1 << ":" << t1 << ":" << N1 << ":" << u1 << ":" << t2
                                      << ":" << N2 << ":" << u2;
#endif
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 + 1, id.iCell2);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, id.iType, id.iLay, id.iSec1, id.iSec2, id.iCell1 - 1, id.iCell2 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t1, id.iLay, wu1, wv1, u1 + N1 - 1, u1 - 1);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, u2 - N2, 0);
        addHGCSiliconId(ids, id.det, id.zSide, t2, id.iLay, wu2, wv2, u2 - N2 + 1, 0);
        break;
      }
      default:
        // Not valid u, v
        int N = hdcons_.getUVMax(id.iType);
        edm::LogWarning("HGCalGeom") << "u:v " << id.iCell1 << ":" << id.iCell2 << " Tests " << (id.iCell1 > 2 * N - 1)
                                     << ":" << (id.iCell2 > 2 * N - 1) << ":" << (id.iCell2 >= (id.iCell1 + N)) << ":"
                                     << (id.iCell1 > (id.iCell2 + N)) << " ERROR";
    }
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    int iphi1 = (id.iCell1 > 1) ? id.iCell1 - 1 : hdcons_.getUVMax(id.iType);
    int iphi2 = (id.iCell1 < hdcons_.getUVMax(id.iType)) ? id.iCell1 + 1 : 1;
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1 - 1, id.iCell1);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1 - 1, iphi1);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1, iphi1);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1 + 1, iphi1);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1 + 1, id.iCell1);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1 + 1, iphi2);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1, iphi2);
    addHGCSCintillatorId(ids, id.zSide, id.iType, id.iLay, id.iSec1 - 1, iphi2);
  }
  return ids;
}

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

inlineoverridevirtual

Get the neighbors of the given cell in north direction

Implements CaloSubdetectorTopology.

Definition at line 22 of file HGCalTopology.h.

References goNorth().

                                                         {
    DetId nextId = goNorth(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }

DetId HGCalTopology::offsetBy	(	const DetId	startId,
		int	nrStepsX,
		int	nrStepsY
	)		const

returns a new DetId offset by nrStepsX and nrStepsY (can be negative), returns DetId(0) if invalid

Definition at line 479 of file HGCalTopology.cc.

References changeXY(), DetId::det(), DetId::Forward, triggerObjects_cff::id, createfilelist::int, subdet_, DetId::subdetId(), and valid().

Referenced by dddConstants().

                                                                                   {
  if (startId.det() == DetId::Forward && startId.subdetId() == (int)(subdet_)) {
    DetId id = changeXY(startId, nrStepsX, nrStepsY);
    if (valid(id))
      return id;
  }
  return DetId(0);
}

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

inlineoverridevirtual

Get the neighbors of the given cell in south direction

Implements CaloSubdetectorTopology.

Definition at line 32 of file HGCalTopology.h.

References goSouth().

                                                         {
    DetId nextId = goSouth(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }

ForwardSubdetector HGCalTopology::subDetector ( ) const

inline

Definition at line 116 of file HGCalTopology.h.

References subdet_.

Referenced by HGCalGeometryLoader::build().

{ return subdet_; }

DetId HGCalTopology::switchZSide

(

const DetId

startId

)

const

Definition at line 488 of file HGCalTopology.cc.

References decode(), encode(), triggerObjects_cff::id, valid(), and HGCalTopology::DecodedDetId::zSide.

Referenced by dddConstants().

                                                          {
  HGCalTopology::DecodedDetId id_ = decode(startId);
  id_.zSide = -id_.zSide;
  DetId id = encode(id_);
  if (valid(id))
    return id;
  else
    return DetId(0);
}

unsigned int HGCalTopology::totalGeomModules ( ) const

inline

Definition at line 91 of file HGCalTopology.h.

References allGeomModules(), and kHGeomHalf_.

Referenced by HGCalGeometryLoader::build(), geomDenseId2decId(), HGCalGeometry::getSummary(), HGCalGeometry::HGCalGeometry(), and HGCalGeometry::sizeForDenseIndex().

{ return (unsigned int)(2 * kHGeomHalf_); }

unsigned int HGCalTopology::totalModules ( ) const

inline

Definition at line 90 of file HGCalTopology.h.

References kSizeForDenseIndexing.

Referenced by HGCalGeometry::HGCalGeometry().

{ return kSizeForDenseIndexing; }

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

inlineoverridevirtual

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

Implements CaloSubdetectorTopology.

Definition at line 60 of file HGCalTopology.h.

References changeZ().

                                                      {
    DetId nextId = changeZ(id, +1);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }

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

overridevirtual

Is this a valid cell id.

Reimplemented from CaloSubdetectorTopology.

Definition at line 463 of file HGCalTopology.cc.

References cells_, decode(), DetId::det(), det_, RemoveAddSevLevel::flag, hdcons_, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, createfilelist::int, HGCalDDDConstants::isValidHex(), HGCalDDDConstants::isValidHex8(), HGCalDDDConstants::isValidTrap(), layers_, mode_, sectors_, subdet_, DetId::subdetId(), and HGCalGeometryMode::Trapezoid.

Referenced by HGCDigitizer::checkPosition(), geomMode(), HGCalGeometry::newCell(), offsetBy(), switchZSide(), HGCalTriggerGeometryV9Imp2::validCellId(), HGCalTriggerGeometryHexImp2::validCellId(), HGCalTriggerGeometryHexLayerBasedImp1::validCellId(), and HGCalTriggerGeometryV9Imp1::validCellId().

                                                 {
  HGCalTopology::DecodedDetId id = decode(idin);
  bool flag;
  if ((mode_ == HGCalGeometryMode::Hexagon) || (mode_ == HGCalGeometryMode::HexagonFull)) {
    flag = (idin.det() == det_ && idin.subdetId() == (int)(subdet_) && id.iCell1 >= 0 && id.iCell1 < cells_ &&
            id.iLay > 0 && id.iLay <= layers_ && id.iSec1 >= 0 && id.iSec1 <= sectors_);
    if (flag)
      flag = hdcons_.isValidHex(id.iLay, id.iSec1, id.iCell1, true);
  } else if (mode_ == HGCalGeometryMode::Trapezoid) {
    flag = ((idin.det() == det_) && hdcons_.isValidTrap(id.iLay, id.iSec1, id.iCell1));
  } else {
    flag = ((idin.det() == det_) && hdcons_.isValidHex8(id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2));
  }
  return flag;
}

bool HGCalTopology::validHashIndex ( uint32_t  ix ) const

inline

Definition at line 88 of file HGCalTopology.h.

References kSizeForDenseIndexing.

Referenced by denseId2detId().

{ return (ix < kSizeForDenseIndexing); }

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

inlineoverridevirtual

Get the neighbors of the given cell in west direction

Implements CaloSubdetectorTopology.

Definition at line 52 of file HGCalTopology.h.

References goWest().

                                                        {
    DetId nextId = goWest(id);
    std::vector<DetId> vNeighborsDetId;
    if (!(nextId == DetId(0)))
      vNeighborsDetId.emplace_back(nextId);
    return vNeighborsDetId;
  }

Member Data Documentation

int HGCalTopology::cellMax_

private

Definition at line 141 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), and HGCalTopology().

int HGCalTopology::cells_

private

Definition at line 140 of file HGCalTopology.h.

Referenced by HGCalTopology(), and valid().

DetId::Detector HGCalTopology::det_

private

Definition at line 138 of file HGCalTopology.h.

Referenced by decode(), detector(), encode(), HGCalTopology(), isHFNose(), and valid().

int HGCalTopology::firstLay_

private

Definition at line 141 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), and HGCalTopology().

const HGCalDDDConstants& HGCalTopology::hdcons_

private

Definition at line 135 of file HGCalTopology.h.

Referenced by addHGCSCintillatorId(), addHGCSiliconId(), allGeomModules(), dddConstants(), HGCalTopology(), neighbors(), and valid().

int HGCalTopology::kHGeomHalf_

private

Definition at line 142 of file HGCalTopology.h.

Referenced by detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), HGCalTopology(), and totalGeomModules().

int HGCalTopology::kHGhalf_

private

Definition at line 142 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseId(), and HGCalTopology().

unsigned int HGCalTopology::kSizeForDenseIndexing

private

Definition at line 143 of file HGCalTopology.h.

Referenced by HGCalTopology(), totalModules(), and validHashIndex().

int HGCalTopology::layers_

private

Definition at line 140 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), HGCalTopology(), and valid().

HGCalGeometryMode::GeometryMode HGCalTopology::mode_

private

Definition at line 136 of file HGCalTopology.h.

Referenced by allGeomModules(), decode(), denseId2detId(), detId2denseGeomId(), detId2denseId(), encode(), geomDenseId2decId(), geomMode(), HGCalTopology(), neighbors(), and valid().

int HGCalTopology::sectors_

private

Definition at line 140 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), HGCalTopology(), and valid().

ForwardSubdetector HGCalTopology::subdet_

private

Definition at line 139 of file HGCalTopology.h.

Referenced by decode(), encode(), HGCalTopology(), isHFNose(), offsetBy(), subDetector(), and valid().

const int HGCalTopology::subSectors_ = 2

static

Use subSector in square mode as wafer type in hexagon mode.

Definition at line 104 of file HGCalTopology.h.

int HGCalTopology::types_

private

Definition at line 140 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseId(), and HGCalTopology().

int HGCalTopology::waferMax_

private

Definition at line 141 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), and HGCalTopology().

int HGCalTopology::waferOff_

private

Definition at line 141 of file HGCalTopology.h.

Referenced by denseId2detId(), detId2denseGeomId(), detId2denseId(), geomDenseId2decId(), and HGCalTopology().