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
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
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_*cells_;
    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 19 of file HGCalTopology.h.

{ }  

Member Function Documentation

unsigned int HGCalTopology::allGeomModules

(

)

const

Definition at line 57 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 358 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 365 of file HGCalTopology.cc.

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

                                                                {

  return DetId();
}

const HGCalDDDConstants& HGCalTopology::dddConstants ( ) const

inline

Definition at line 101 of file HGCalTopology.h.

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

Referenced by hgc_digi_utils::addCellMetadata(), PFRecHitQTestThresholdInThicknessNormalizedMIPs::beginEvent(), CaloTruthAccumulator::beginLuminosityBlock(), HGCalGeometryLoader::build(), HGCDigitizer::checkPosition(), HGCalTriggerGeometryHexImp2::detIdWaferType(), HGCalTriggerGeometryHexLayerBasedImp1::detIdWaferType(), HGCalTriggerGeometryV9Imp1::detIdWaferType(), HGCalTriggerTools::eventSetup(), HGCalTriggerNtupleHGCTriggerCells::fill(), HGCalTriggerGeometryHexImp2::fillInvalidTriggerCells(), HGCalTriggerGeometryHexLayerBasedImp1::fillInvalidTriggerCells(), HGCalTriggerGeometryV9Imp1::fillInvalidTriggerCells(), HGCalTriggerGeometryHexImp1::fillMaps(), CaloParticleDebugger::fillSimHits(), HGCalGeometry::get8Corners(), HGCalTriggerGeometryV9Imp1::getCellsFromTriggerCell(), HGCalGeometry::getClosestCell(), HGCalGeometry::getCorners(), HGCalTriggerTools::getLayerZ(), HGCalTriggerGeometryV9Imp1::getNeighborsFromTriggerCell(), HGCalGeometry::getPosition(), HGCalGeometry::getSummary(), HGCalTriggerGeometryHexLayerBasedImp1::initialize(), HGCalTriggerGeometryV9Imp1::initialize(), HGCalGeometry::newCell(), HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run(), HGCalRecHitWorkerSimple::set(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry(), HGCalTriggerNtupleHGCTriggerCells::simhits(), HGCalTriggerNtupleHGCDigis::simhits(), HGCalTriggerCellBestChoiceCodecImpl::triggerCellSums(), and HGCalTriggerCellThresholdCodecImpl::triggerCellSums().

{return hdcons_;}

HGCalTopology::DecodedDetId HGCalTopology::decode

(

const DetId &

id

)

const

Definition at line 294 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::getPosition(), 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 111 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 121 of file HGCalTopology.h.

References det_.

Referenced by HGCalGeometryLoader::build().

{ return det_;}

bool HGCalTopology::detectorType ( ) const

inline

Definition at line 123 of file HGCalTopology.h.

Referenced by HGCalHitValidation::analyzeHGCalSimHit(), HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run(), HGCalTriggerNtupleHGCTriggerCells::simhits(), and HGCalTriggerNtupleHGCDigis::simhits().

{ return false;}

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

virtual

Definition at line 163 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 63 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 77 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;
  }

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 49 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 342 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(), CaloGeometryDBEP< T, U >::produceAligned(), 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 246 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_.

Referenced by CaloGeometryDBEP< T, U >::produceAligned().

                                                                                   {

  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 86 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 46 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 22 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 34 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 58 of file HGCalTopology.h.

References changeXY().

Referenced by west().

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

bool HGCalTopology::isHFNose ( ) const

inline

Definition at line 124 of file HGCalTopology.h.

References changeXY(), changeZ(), 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);
  }

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 25 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 227 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 37 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 122 of file HGCalTopology.h.

References subdet_.

Referenced by HGCalGeometryLoader::build().

{return subdet_;}

DetId HGCalTopology::switchZSide

(

const DetId

startId

)

const

Definition at line 237 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 98 of file HGCalTopology.h.

References allGeomModules(), and kHGeomHalf_.

Referenced by HGCalGeometryLoader::build(), geomDenseId2decId(), HGCalGeometry::getSummary(), HGCalGeometry::HGCalGeometry(), CaloGeometryDBEP< T, U >::produceAligned(), and HGCalGeometry::sizeForDenseIndex().

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

unsigned int HGCalTopology::totalModules ( ) const

inline

Definition at line 97 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 69 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 207 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 geomMode(), HGCalGeometry::newCell(), offsetBy(), switchZSide(), 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 95 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 61 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 142 of file HGCalTopology.h.

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

int HGCalTopology::cells_

private

Definition at line 141 of file HGCalTopology.h.

Referenced by HGCalTopology(), and valid().

DetId::Detector HGCalTopology::det_

private

Definition at line 139 of file HGCalTopology.h.

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

int HGCalTopology::firstLay_

private

Definition at line 142 of file HGCalTopology.h.

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

const HGCalDDDConstants& HGCalTopology::hdcons_

private

Definition at line 136 of file HGCalTopology.h.

Referenced by allGeomModules(), dddConstants(), HGCalTopology(), and valid().

int HGCalTopology::kHGeomHalf_

private

Definition at line 143 of file HGCalTopology.h.

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

int HGCalTopology::kHGhalf_

private

Definition at line 143 of file HGCalTopology.h.

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

unsigned int HGCalTopology::kSizeForDenseIndexing

private

Definition at line 144 of file HGCalTopology.h.

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

int HGCalTopology::layers_

private

Definition at line 141 of file HGCalTopology.h.

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

HGCalGeometryMode::GeometryMode HGCalTopology::mode_

private

Definition at line 137 of file HGCalTopology.h.

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

int HGCalTopology::sectors_

private

Definition at line 141 of file HGCalTopology.h.

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

ForwardSubdetector HGCalTopology::subdet_

private

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

int HGCalTopology::types_

private

Definition at line 141 of file HGCalTopology.h.

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

int HGCalTopology::waferMax_

private

Definition at line 142 of file HGCalTopology.h.

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

int HGCalTopology::waferOff_

private

Definition at line 142 of file HGCalTopology.h.

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