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
bool	tileTrapezoid () 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	valid (const DetId &id, int cornerMin) const
bool	validHashIndex (uint32_t ix) const
bool	validModule (const DetId &id, int cornerMin) const
bool	waferHexagon6 () const
bool	waferHexagon8 () const
std::vector< DetId >	west (const DetId &id) const override
	~HGCalTopology () override=default
	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_
int	kHGhalfType_
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::HGCalTopology	(	const HGCalDDDConstants &	hdcons,
		int	subdet
	)

create a new Topology

Definition at line 13 of file HGCalTopology.cc.

References cellMax_, cells_, det_, firstLay_, HGCalDDDConstants::firstLayer(), DetId::Forward, ForwardEmpty, HGCalDDDConstants::geomMode(), hdcons_, HFNose, kHGeomHalf_, kHGhalf_, kHGhalfType_, kSizeForDenseIndexing, HGCalDDDConstants::layers(), layers_, HGCalDDDConstants::maxCells(), HGCalDDDConstants::maxCellUV(), mode_, HGCalDDDConstants::sectors(), sectors_, subdet_, tileTrapezoid(), types_, waferHexagon6(), 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 (waferHexagon6()) {
    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 (tileTrapezoid()) {
    det_ = (DetId::Detector)(det);
    subdet_ = ForwardEmpty;
    kHGeomHalf_ = sectors_ * layers_ * cellMax_;
    types_ = 3;
  } else {
    det_ = (DetId::Detector)(det);
    subdet_ = ForwardEmpty;
    kHGeomHalf_ = sectors_ * layers_;
    types_ = 3;
  }
  kHGhalfType_ = sectors_ * layers_ * cells_ * types_;
  kSizeForDenseIndexing = static_cast<unsigned int>(2 * kHGhalf_);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "HGCalTopology initialized for detector " << det << ":" << det_ << ":" << subdet_
                                << " having " << sectors_ << " Sectors, " << layers_ << " Layers from " << firstLay_
                                << ", " << cells_ << " cells and total channels " << kSizeForDenseIndexing << ":"
                                << (2 * kHGeomHalf_);
#endif
}

~HGCalTopology()

HGCalTopology::~HGCalTopology ( )

overridedefault

virtual destructor

Member Function Documentation

addHGCSCintillatorId()

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 565 of file HGCalTopology.cc.

References hdcons_, l1ctLayer2EG_cff::id, LEDCalibrationChannels::iphi, HGCalDDDConstants::isValidTrap(), and ecaldqm::zside().

Referenced by neighbors().

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

addHGCSiliconId()

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 577 of file HGCalTopology.cc.

References hdcons_, l1ctLayer2EG_cff::id, HGCalDDDConstants::isValidHex8(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and ecaldqm::zside().

Referenced by 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, false);
#endif
  if (hdcons_.isValidHex8(lay, waferU, waferV, cellU, cellV, false)) {
    HGCSiliconDetId id((DetId::Detector)(det), zside, type, lay, waferU, waferV, cellU, cellV);
    ids.emplace_back(DetId(id));
  }
}

allGeomModules()

unsigned int HGCalTopology::allGeomModules

(

)

const

Definition at line 54 of file HGCalTopology.cc.

References hdcons_, HGCalDDDConstants::numberCells(), tileTrapezoid(), and HGCalDDDConstants::wafers().

Referenced by HGCalGeometryLoader::build().

                                                 {
  return (tileTrapezoid() ? (unsigned int)(2 * hdcons_.numberCells(true)) : (unsigned int)(2 * hdcons_.wafers()));
}

changeXY()

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

private

move the nagivator along x, y

Definition at line 664 of file HGCalTopology.cc.

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

{ return DetId(); }

changeZ()

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

private

move the nagivator along z

Definition at line 666 of file HGCalTopology.cc.

Referenced by down(), and up().

{ return DetId(); }

dddConstants()

const HGCalDDDConstants& HGCalTopology::dddConstants ( ) const

inline

Definition at line 98 of file HGCalTopology.h.

References hdcons_.

Referenced by PFRecHitQTestThresholdInThicknessNormalizedMIPs::beginEvent(), HGCalGeometryLoader::build(), HGCDigitizer::checkPosition(), HGCalTriggerGeometryV9Imp2::detIdWaferType(), HGCalTriggerGeometryV9Imp3::detIdWaferType(), CaloParticleDebugger::fillSimHits(), HGCalGeometry::get8Corners(), HGCalGeometry::getClosestCell(), HGCalGeometry::getClosestCellHex(), HGCalGeometry::getCorners(), HGCalTriggerTools::getLayerZ(), HGCalTriggerGeometryV9Imp2::getModuleFromTriggerCell(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), HGCalGeometry::getSummary(), HGCalTriggerGeometryV9Imp2::getTriggerCellsFromModule(), HGCalTriggerGeometryV9Imp2::initialize(), HGCalTriggerGeometryV9Imp3::initialize(), CaloTruthAccumulator::initializeEvent(), hgcal::RecHitTools::maskCell(), maskCell(), HGCalGeometry::neighborZ(), HGCalGeometry::newCell(), HGCalRecHitWorkerSimple::set(), HGCalRadiationMap::setGeometry(), HGCalTriggerTools::setGeometry(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry(), and HGCalTriggerTools::simToReco().

{ return hdcons_; }

decode()

HGCalTopology::DecodedDetId HGCalTopology::decode

(

const DetId &

id

)

const

Definition at line 590 of file HGCalTopology.cc.

References det_, DetId::Forward, HFNose, l1ctLayer2EG_cff::id, heavyIonCSV_trainingSettings::idx, createfilelist::int, subdet_, tileTrapezoid(), and waferHexagon6().

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

                                                                          {
  HGCalTopology::DecodedDetId idx;
  if (waferHexagon6()) {
    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 (tileTrapezoid()) {
    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;
}

denseId2detId()

DetId HGCalTopology::denseId2detId ( uint32_t  denseId ) const

override

Definition at line 389 of file HGCalTopology.cc.

References cellMax_, encode(), firstLay_, l1ctLayer2EG_cff::id, createfilelist::int, kHGhalfType_, layers_, sectors_, tileTrapezoid(), types_, validHashIndex(), waferHexagon6(), waferMax_, and waferOff_.

                                                    {
  HGCalTopology::DecodedDetId id;
  if (validHashIndex(hi)) {
    id.zSide = ((int)(hi) < kHGhalfType_ ? -1 : 1);
    int di = ((int)(hi) % kHGhalfType_);
    if (waferHexagon6()) {
      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 (tileTrapezoid()) {
      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);
}

detector()

DetId::Detector HGCalTopology::detector ( ) const

inline

Definition at line 117 of file HGCalTopology.h.

References det_.

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getClosestCell(), and HGCalGeometry::getClosestCellHex().

{ return det_; }

detectorType()

bool HGCalTopology::detectorType ( ) const

inline

Definition at line 119 of file HGCalTopology.h.

{ return false; }

detId2denseGeomId()

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

virtual

Definition at line 435 of file HGCalTopology.cc.

References cellMax_, decode(), firstLay_, heavyIonCSV_trainingSettings::idx, kHGeomHalf_, layers_, sectors_, tileTrapezoid(), waferHexagon6(), waferMax_, and waferOff_.

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

                                                                 {
  HGCalTopology::DecodedDetId id = decode(idin);
  uint32_t idx;
  if (waferHexagon6()) {
    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 (tileTrapezoid()) {
    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;
}

detId2denseId()

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

overridevirtual

Dense indexing.

Reimplemented from CaloSubdetectorTopology.

Definition at line 350 of file HGCalTopology.cc.

References cellMax_, decode(), firstLay_, heavyIonCSV_trainingSettings::idx, kHGeomHalf_, kHGhalfType_, layers_, sectors_, tileTrapezoid(), types_, waferHexagon6(), waferMax_, and waferOff_.

                                                             {
  HGCalTopology::DecodedDetId id = decode(idin);
  uint32_t idx;
  if (waferHexagon6()) {
    int type = (id.iType > 0) ? 1 : 0;
    idx = (uint32_t)(((id.zSide > 0) ? kHGhalfType_ : 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 (tileTrapezoid()) {
    idx =
        (uint32_t)(((id.zSide > 0) ? kHGhalfType_ : 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) ? kHGhalfType_ : 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;
}

down()

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().

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

east()

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;
  }

encode()

DetId HGCalTopology::encode

(

const DecodedDetId &

id_

)

const

Definition at line 636 of file HGCalTopology.cc.

References det_, DetId::Forward, hdcons_, HFNose, l1ctLayer2EG_cff::id, heavyIonCSV_trainingSettings::idx, HGCScintillatorDetId::layer(), DetId::rawId(), HGCScintillatorDetId::ring(), HGCScintillatorDetId::setSiPM(), HGCScintillatorDetId::setType(), subdet_, tileTrapezoid(), HGCalDDDConstants::tileType(), and waferHexagon6().

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

                                                                      {
  DetId id;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeomX") << "Encode " << idx.det << ":" << idx.zSide << ":" << idx.iType << ":" << idx.iLay
                                 << ":" << idx.iSec1 << ":" << idx.iSec2 << ":" << idx.iCell1 << ":" << idx.iCell2;
#endif
  if (waferHexagon6()) {
    id =
        HGCalDetId((ForwardSubdetector)(idx.det), idx.zSide, idx.iLay, ((idx.iType > 0) ? 1 : 0), idx.iSec1, idx.iCell1)
            .rawId();
  } else if (tileTrapezoid()) {
    HGCScintillatorDetId hid(idx.iType, idx.iLay, idx.zSide * idx.iSec1, idx.iCell1);
    std::pair<int, int> typm = hdcons_.tileType(hid.layer(), hid.ring(), 0);
    if (typm.first >= 0) {
      hid.setType(typm.first);
      hid.setSiPM(typm.second);
    }
    id = hid.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;
}

geomDenseId2decId()

HGCalTopology::DecodedDetId HGCalTopology::geomDenseId2decId

(

const uint32_t &

hi

)

const

Definition at line 524 of file HGCalTopology.cc.

References cellMax_, firstLay_, l1ctLayer2EG_cff::id, createfilelist::int, kHGeomHalf_, layers_, sectors_, tileTrapezoid(), totalGeomModules(), waferHexagon6(), waferMax_, and waferOff_.

                                                                                   {
  HGCalTopology::DecodedDetId id;
  if (hi < totalGeomModules()) {
    id.zSide = ((int)(hi) < kHGeomHalf_ ? -1 : 1);
    int di = ((int)(hi) % kHGeomHalf_);
    if (waferHexagon6()) {
      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 (tileTrapezoid()) {
      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;
}

geomMode()

HGCalGeometryMode::GeometryMode HGCalTopology::geomMode ( ) const

inline

Geometry mode.

Definition at line 79 of file HGCalTopology.h.

References mode_.

Referenced by HGCalGeometryLoader::build().

{ return mode_; }

goEast()

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); }

goNorth()

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); }

goSouth()

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); }

goWest()

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); }

isHFNose()

bool HGCalTopology::isHFNose ( ) const

inline

Definition at line 120 of file HGCalTopology.h.

References det_, DetId::Forward, HFNose, and subdet_.

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

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

maskCell()

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

inline

Definition at line 96 of file HGCalTopology.h.

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

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

neighbors()

std::vector< DetId > HGCalTopology::neighbors

(

const DetId &

id

)

const

Definition at line 58 of file HGCalTopology.cc.

References addHGCSCintillatorId(), addHGCSiliconId(), HGCalCell::bottomCorner, HGCalCell::bottomLeftCorner, HGCalCell::bottomLeftEdge, HGCalCell::bottomRightCorner, HGCalCell::bottomRightEdge, HGCalDDDConstants::cellType(), HGCalCell::centralCell, decode(), HGCalDDDConstants::getTypeHex(), HGCalDDDConstants::getUVMax(), hdcons_, HGCalCell::leftEdge, HGCalDDDConstants::modifyUV(), N, HGCalCell::rightEdge, RandomServiceHelper::t1, RandomServiceHelper::t2, tileTrapezoid(), HGCalCell::topCorner, HGCalCell::topLeftCorner, HGCalCell::topLeftEdge, HGCalCell::topRightCorner, HGCalCell::topRightEdge, testProducerWithPsetDescEmpty_cfi::u1, MetAnalyzer::u2, and waferHexagon8().

                                                                 {
  std::vector<DetId> ids;
  HGCalTopology::DecodedDetId id = decode(idin);
  if (waferHexagon8()) {
    int celltype = hdcons_.cellType(
        id.iType, id.iCell1, id.iCell2, id.zSide, 1, -1);  // Temporary fix - later for v17 define fwd back and orient
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HGCalGeom") << "Type:WaferU:WaferV " << id.iType << ":" << id.iCell1 << ":" << id.iCell2
                                  << " CellType " << celltype;
#endif
    switch (celltype) {
      case (HGCalCell::centralCell): {
        // 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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (HGCalCell::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 (tileTrapezoid()) {
    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;
}

north()

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;
  }

offsetBy()

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 505 of file HGCalTopology.cc.

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

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

south()

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;
  }

subDetector()

ForwardSubdetector HGCalTopology::subDetector ( ) const

inline

Definition at line 118 of file HGCalTopology.h.

References subdet_.

Referenced by HGCalGeometryLoader::build().

{ return subdet_; }

switchZSide()

DetId HGCalTopology::switchZSide

(

const DetId

startId

)

const

Definition at line 514 of file HGCalTopology.cc.

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

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

tileTrapezoid()

bool HGCalTopology::tileTrapezoid ( ) const

inline

Definition at line 124 of file HGCalTopology.h.

References hdcons_, and HGCalDDDConstants::tileTrapezoid().

Referenced by allGeomModules(), HGCalGeometryLoader::build(), decode(), denseId2detId(), detId2denseGeomId(), detId2denseId(), encode(), geomDenseId2decId(), HGCalGeometry::getClosestCell(), HGCalGeometry::getGeometryDetId(), HGCalGeometry::getPosition(), HGCalGeometry::getSummary(), HGCalGeometry::getWaferPosition(), HGCalTopology(), neighbors(), HGCalGeometry::newCell(), and valid().

{ return hdcons_.tileTrapezoid(); }

totalGeomModules()

unsigned int HGCalTopology::totalGeomModules ( ) const

inline

Definition at line 93 of file HGCalTopology.h.

References kHGeomHalf_.

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

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

totalModules()

unsigned int HGCalTopology::totalModules ( ) const

inline

Definition at line 92 of file HGCalTopology.h.

References kSizeForDenseIndexing.

Referenced by HGCalGeometry::HGCalGeometry().

{ return kSizeForDenseIndexing; }

up()

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;
  }

valid() [1/2]

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

overridevirtual

Is this a valid cell id.

Reimplemented from CaloSubdetectorTopology.

Definition at line 465 of file HGCalTopology.cc.

References cells_, decode(), DetId::det(), det_, RemoveAddSevLevel::flag, hdcons_, createfilelist::int, HGCalDDDConstants::isValidHex(), HGCalDDDConstants::isValidHex8(), HGCalDDDConstants::isValidTrap(), layers_, sectors_, subdet_, DetId::subdetId(), tileTrapezoid(), and waferHexagon6().

Referenced by HGCDigitizer::checkPosition(), HGCalGeometry::newCell(), offsetBy(), switchZSide(), valid(), HGCalTriggerGeometryV9Imp2::validCell(), HGCalTriggerGeometryV9Imp3::validCell(), HGCalTriggerGeometryV9Imp2::validCellId(), HGCalTriggerGeometryV9Imp3::validCellId(), and validModule().

                                                 {
  HGCalTopology::DecodedDetId id = decode(idin);
  bool flag;
  if (waferHexagon6()) {
    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 (tileTrapezoid()) {
    flag = ((idin.det() == det_) && hdcons_.isValidTrap(id.zSide, id.iLay, id.iSec1, id.iCell1));
  } else {
    flag = ((idin.det() == det_) && hdcons_.isValidHex8(id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2, false));
  }
  return flag;
}

valid() [2/2]

bool HGCalTopology::valid	(	const DetId &	id,
		int	cornerMin
	)		const

Definition at line 481 of file HGCalTopology.cc.

References decode(), DetId::det(), det_, RemoveAddSevLevel::flag, hdcons_, HGCalDDDConstants::isValidHex8(), gpuClustering::pixelStatus::mask, HGCalDDDConstants::maskCell(), valid(), HGCalTypes::WaferCornerMin, and waferHexagon8().

                                                                {
  if (waferHexagon8()) {
    HGCalTopology::DecodedDetId id = decode(idin);
    bool mask = (cornerMin < HGCalTypes::WaferCornerMin) ? false : hdcons_.maskCell(idin, cornerMin);
    bool flag = ((idin.det() == det_) &&
                 hdcons_.isValidHex8(
                     id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2, (cornerMin >= HGCalTypes::WaferCornerMin)));
    return (flag && (!mask));
  } else {
    return valid(idin);
  }
}

validHashIndex()

bool HGCalTopology::validHashIndex ( uint32_t  ix ) const

inline

Definition at line 89 of file HGCalTopology.h.

References kSizeForDenseIndexing.

Referenced by denseId2detId().

{ return (ix < kSizeForDenseIndexing); }

validModule()

bool HGCalTopology::validModule	(	const DetId &	id,
		int	cornerMin
	)		const

Definition at line 494 of file HGCalTopology.cc.

References decode(), DetId::det(), det_, hdcons_, DetId::HGCalEE, DetId::HGCalHSi, HGCalDDDConstants::isValidHex8(), valid(), and HGCalTypes::WaferCornerMin.

                                                                      {
  if (idin.det() != det_) {
    return false;
  } else if ((idin.det() == DetId::HGCalEE) || (idin.det() == DetId::HGCalHSi)) {
    HGCalTopology::DecodedDetId id = decode(idin);
    return hdcons_.isValidHex8(id.iLay, id.iSec1, id.iSec2, (cornerMin >= HGCalTypes::WaferCornerMin));
  } else {
    return valid(idin);
  }
}

waferHexagon6()

bool HGCalTopology::waferHexagon6 ( ) const

inline

Definition at line 125 of file HGCalTopology.h.

References hdcons_, and HGCalDDDConstants::waferHexagon6().

Referenced by HGCalGeometryLoader::build(), decode(), denseId2detId(), detId2denseGeomId(), detId2denseId(), encode(), geomDenseId2decId(), HGCalGeometry::get8Corners(), HGCalGeometry::getClosestCell(), HGCalGeometry::getCorners(), HGCalGeometry::getGeometryDetId(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), HGCalGeometry::getSummary(), HGCalTopology(), HGCalGeometry::newCell(), and valid().

{ return hdcons_.waferHexagon6(); }

waferHexagon8()

bool HGCalTopology::waferHexagon8 ( ) const

inline

Definition at line 126 of file HGCalTopology.h.

References hdcons_, and HGCalDDDConstants::waferHexagon8().

Referenced by HGCalGeometry::getClosestCellHex(), neighbors(), and valid().

{ return hdcons_.waferHexagon8(); }

west()

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

cellMax_

int HGCalTopology::cellMax_

private

Definition at line 147 of file HGCalTopology.h.

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

cells_

int HGCalTopology::cells_

private

Definition at line 146 of file HGCalTopology.h.

Referenced by HGCalTopology(), and valid().

det_

DetId::Detector HGCalTopology::det_

private

Definition at line 144 of file HGCalTopology.h.

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

firstLay_

int HGCalTopology::firstLay_

private

Definition at line 147 of file HGCalTopology.h.

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

hdcons_

const HGCalDDDConstants& HGCalTopology::hdcons_

private

Definition at line 141 of file HGCalTopology.h.

Referenced by addHGCSCintillatorId(), addHGCSiliconId(), allGeomModules(), dddConstants(), encode(), HGCalTopology(), neighbors(), tileTrapezoid(), valid(), validModule(), waferHexagon6(), and waferHexagon8().

kHGeomHalf_

int HGCalTopology::kHGeomHalf_

private

Definition at line 148 of file HGCalTopology.h.

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

kHGhalf_

int HGCalTopology::kHGhalf_

private

Definition at line 148 of file HGCalTopology.h.

Referenced by HGCalTopology().

kHGhalfType_

int HGCalTopology::kHGhalfType_

private

Definition at line 148 of file HGCalTopology.h.

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

kSizeForDenseIndexing

unsigned int HGCalTopology::kSizeForDenseIndexing

private

Definition at line 149 of file HGCalTopology.h.

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

layers_

int HGCalTopology::layers_

private

Definition at line 146 of file HGCalTopology.h.

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

mode_

HGCalGeometryMode::GeometryMode HGCalTopology::mode_

private

Definition at line 142 of file HGCalTopology.h.

Referenced by geomMode(), and HGCalTopology().

sectors_

int HGCalTopology::sectors_

private

Definition at line 146 of file HGCalTopology.h.

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

subdet_

ForwardSubdetector HGCalTopology::subdet_

private

Definition at line 145 of file HGCalTopology.h.

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

subSectors_

const int HGCalTopology::subSectors_ = 2

static

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

Definition at line 106 of file HGCalTopology.h.

types_

int HGCalTopology::types_

private

Definition at line 146 of file HGCalTopology.h.

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

waferMax_

int HGCalTopology::waferMax_

private

Definition at line 147 of file HGCalTopology.h.

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

waferOff_

int HGCalTopology::waferOff_

private

Definition at line 147 of file HGCalTopology.h.

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