#include <HGCalGeometry.h>

Inheritance diagram for HGCalGeometry:

Public Types
enum	{ k_NumberOfParametersPerShape = 12 }

enum	{ k_NumberOfShapes = 50 }

typedef HGCalGeometryRecord	AlignedRecord

typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< FlatTrd >	CellVec

typedef std::vector< GlobalPoint >	CornersVec

typedef std::set< DetId >	DetIdSet

typedef PHGCalRcd	PGeometryRecord

typedef CaloCellGeometry::Pt3D	Pt3D

typedef CaloCellGeometry::Pt3DVec	Pt3DVec

Public Types inherited from CaloSubdetectorGeometry
typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector < CaloCellGeometry const * >	CellSet

typedef std::set< DetId >	DetIdSet

typedef std::vector< CCGFloat >	DimVec

typedef std::vector< unsigned int >	IVec

typedef CaloCellGeometry::ParMgr	ParMgr

typedef CaloCellGeometry::ParVec	ParVec

typedef CaloCellGeometry::ParVecVec	ParVecVec

typedef std::vector< CCGFloat >	TrVec

Public Member Functions
std::string	cellElement () const

virtual void	fillNamedParams (DDFilteredView fv)

virtual DetIdSet	getCells (const GlobalPoint &r, double dR) const override
	Get a list of all cells within a dR of the given cell. More...

virtual DetId	getClosestCell (const GlobalPoint &r) const override

CornersVec	getCorners (const DetId &id) const
	Returns the corner points of this cell's volume. More...

virtual const CaloCellGeometry *	getGeometry (const DetId &id) const override
	Get the cell geometry of a given detector id. Should return false if not found. More...

GlobalPoint	getPosition (const DetId &id) const

virtual void	getSummary (CaloSubdetectorGeometry::TrVec &trVector, CaloSubdetectorGeometry::IVec &iVector, CaloSubdetectorGeometry::DimVec &dimVector, CaloSubdetectorGeometry::IVec &dinsVector) const

virtual const std::vector < DetId > &	getValidDetIds (DetId::Detector det=DetId::Detector(0), int subdet=0) const override
	Get a list of valid detector ids (for the given subdetector) More...

const std::vector< DetId > &	getValidGeomDetIds (void) const

	HGCalGeometry (const HGCalTopology &topology)

virtual void	initializeParms () override

void	localCorners (Pt3DVec &lc, const CCGFloat *pv, unsigned int i, Pt3D &ref)

virtual void	newCell (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId) override

void	sortDetIds ()

const HGCalTopology &	topology () const

virtual	~HGCalGeometry ()

Public Member Functions inherited from CaloSubdetectorGeometry
void	allocateCorners (CaloCellGeometry::CornersVec::size_type n)

void	allocatePar (ParVec::size_type n, unsigned int m)

	CaloSubdetectorGeometry ()

CaloCellGeometry::CornersMgr *	cornersMgr ()

CCGFloat	deltaEta (const DetId &detId) const

CCGFloat	deltaPhi (const DetId &detId) const

virtual void	fillDefaultNamedParameters () const

virtual CellSet	getCellSet (const GlobalPoint &r, double dR) const

virtual unsigned int	numberOfParametersPerShape () const

virtual unsigned int	numberOfShapes () const

virtual unsigned int	numberOfTransformParms () const

ParMgr *	parMgr ()

const ParMgr *	parMgrConst () const

ParVecVec &	parVecVec ()

const ParVecVec &	parVecVec () const

virtual bool	present (const DetId &id) const
	is this detid present in the geometry? More...

virtual	~CaloSubdetectorGeometry ()
	The base class DOES assume that it owns the CaloCellGeometry objects. More...

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

static std::string	producerTag ()

Protected Member Functions
void	addValidID (const DetId &id)

virtual const CaloCellGeometry *	cellGeomPtr (uint32_t index) const override

unsigned int	getClosestCellIndex (const GlobalPoint &r) const

virtual unsigned int	indexFor (const DetId &id) const override

unsigned int	sizeForDenseIndex () const

Protected Member Functions inherited from CaloSubdetectorGeometry
void	addValidID (const DetId &id)

virtual unsigned int	sizeForDenseIndex (const DetId &id) const

Private Attributes
CellVec	m_cellVec

bool	m_halfType

ForwardSubdetector	m_subdet

const HGCalTopology &	m_topology

std::vector< DetId >	m_validGeomIds

Additional Inherited Members
Static Protected Member Functions inherited from CaloSubdetectorGeometry
static CCGFloat	deltaR (const GlobalPoint &p1, const GlobalPoint &p2)

Protected Attributes inherited from CaloSubdetectorGeometry
ParVecVec	m_parVecVec

std::vector< DetId >	m_validIds

Detailed Description

Definition at line 26 of file HGCalGeometry.h.

Member Typedef Documentation

typedef HGCalGeometryRecord HGCalGeometry::AlignedRecord

Definition at line 39 of file HGCalGeometry.h.

typedef CaloCellGeometry::CCGFloat HGCalGeometry::CCGFloat

Definition at line 32 of file HGCalGeometry.h.

typedef std::vector<FlatTrd> HGCalGeometry::CellVec

Definition at line 30 of file HGCalGeometry.h.

typedef std::vector<GlobalPoint> HGCalGeometry::CornersVec

Definition at line 37 of file HGCalGeometry.h.

typedef std::set<DetId> HGCalGeometry::DetIdSet

Definition at line 36 of file HGCalGeometry.h.

typedef PHGCalRcd HGCalGeometry::PGeometryRecord

Definition at line 40 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3D HGCalGeometry::Pt3D

Definition at line 33 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3DVec HGCalGeometry::Pt3DVec

Definition at line 34 of file HGCalGeometry.h.

Member Enumeration Documentation

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 42 of file HGCalGeometry.h.

42 { k_NumberOfParametersPerShape = 12 } ; // FlatTrd

HGCalGeometry::k_NumberOfParametersPerShape

Definition: HGCalGeometry.h:42

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 43 of file HGCalGeometry.h.

43 { k_NumberOfShapes = 50 } ;

HGCalGeometry::k_NumberOfShapes

Definition: HGCalGeometry.h:43

Constructor & Destructor Documentation

HGCalGeometry::HGCalGeometry ( const HGCalTopology & topology )

Definition at line 24 of file HGCalGeometry.cc.

References gather_cfg::cout, CaloSubdetectorGeometry::m_validIds, topology(), and HGCalTopology::totalGeomModules().

   : m_topology( topology_ ),
     m_cellVec( topology_.totalGeomModules()),
     m_validGeomIds( topology_.totalGeomModules()),
     m_halfType( topology_.detectorType()),
     m_subdet( topology_.subDetector()) {
   
   m_validIds.reserve( topology().totalModules());
 #ifdef DebugLog
   std::cout << "Expected total # of Geometry Modules " 
         << topology().totalGeomModules() << std::endl;
 #endif
 }

HGCalGeometry::~HGCalGeometry ( )

virtual

Definition at line 38 of file HGCalGeometry.cc.

38 { }

Member Function Documentation

void HGCalGeometry::addValidID ( const DetId & id )

protected

Definition at line 265 of file HGCalGeometry.cc.

                                               {
   edm::LogError("HGCalGeom") << "HGCalGeometry::addValidID is not implemented";
 }

std::string HGCalGeometry::cellElement ( ) const

Definition at line 223 of file HGCalGeometry.cc.

References HGCEE, HGCHEB, HGCHEF, and m_subdet.

                                            {
   if      (m_subdet == HGCEE)  return "HGCalEE";
   else if (m_subdet == HGCHEF) return "HGCalHEFront";
   else if (m_subdet == HGCHEB) return "HGCalHEBack";
   else                         return "Unknown";
 }

const CaloCellGeometry * HGCalGeometry::cellGeomPtr ( uint32_t index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 254 of file HGCalGeometry.cc.

References m_cellVec, m_validGeomIds, and CaloCellGeometry::param().

Referenced by getGeometry(), and getSummary().

                                                                        {
   if ((index >= m_cellVec.size()) || (m_validGeomIds[index].rawId() == 0)) 
     return 0;
   const CaloCellGeometry* cell ( &m_cellVec[ index ] ) ;
 #ifdef DebugLog
   //  std::cout << "cellGeomPtr " << m_cellVec[index];
 #endif
   if (0 == cell->param()) return 0;
   return cell;
 }

static std::string HGCalGeometry::dbString ( )

inlinestatic

Definition at line 45 of file HGCalGeometry.h.

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

45 { return "PHGCalRcd" ; }

void HGCalGeometry::fillNamedParams ( DDFilteredView fv )

virtual

Definition at line 40 of file HGCalGeometry.cc.

40 {}

HGCalGeometry::DetIdSet HGCalGeometry::getCells	(	const GlobalPoint &	r,
		double	dR
	)		const

overridevirtual

Get a list of all cells within a dR of the given cell.

The default implementation makes a loop over all cell geometries. Cleverer implementations are suggested to use rough conversions between eta/phi and ieta/iphi and test on the boundaries.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 219 of file HGCalGeometry.cc.

                                                                                      {
    HGCalGeometry::DetIdSet dss;
    return dss;
 }

DetId HGCalGeometry::getClosestCell ( const GlobalPoint & r ) const

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 183 of file HGCalGeometry.cc.

References HGCalDDDConstants::assignCell(), gather_cfg::cout, HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::encode(), HGCalDDDConstants::geomMode(), getClosestCellIndex(), HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, m_cellVec, m_validGeomIds, Square, topology(), HGCalDDDConstants::waferTypeT(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and HGCalTopology::DecodedDetId::zside.

                                                               {
   unsigned int cellIndex = getClosestCellIndex(r);
   if (cellIndex < m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(m_validGeomIds[cellIndex]);
     HepGeom::Point3D<float> local;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       local = m_cellVec[cellIndex].getLocal(r);
     } else {
       if (r.z() > 0) local = HepGeom::Point3D<float>(r.x(),r.y(),0);
       else           local = HepGeom::Point3D<float>(-r.x(),r.y(),0);
     }
     std::pair<int,int> kxy = 
       topology().dddConstants().assignCell(local.x(),local.y(),id_.iLay,
                        id_.iSubSec,true);
     id_.iCell   = kxy.second;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       id_.iSubSec = kxy.first;
     } else {
       id_.iSec    = kxy.first;
       id_.iSubSec = topology().dddConstants().waferTypeT(kxy.first);
       if (id_.iSubSec != 1) id_.iSubSec = -1;
     }
 #ifdef DebugLog
     std::cout << "getClosestCell: local " << local << " Id " << id_.zside 
           << ":" << id_.iLay << ":" << id_.iSec << ":" << id_.iSubSec
           << ":" << id_.iCell << std::endl;
 #endif
 
     //check if returned cell is valid
     if (id_.iCell>=0) return topology().encode(id_);
   }
 
   //if not valid or out of bounds return a null DetId
   return DetId();
 }

unsigned int HGCalGeometry::getClosestCellIndex ( const GlobalPoint & r ) const

protected

Definition at line 270 of file HGCalGeometry.cc.

References gather_cfg::cout, getPosition(), relval_2017::k, m_cellVec, M_PI, PV3DBase< T, PVType, FrameType >::phi(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by getClosestCell().

                                                                            {
 
   float phip = r.phi();
   float zp   = r.z();
   unsigned int cellIndex =  m_cellVec.size();
   float dzmin(9999), dphimin(9999), dphi10(0.175);
   for (unsigned int k=0; k<m_cellVec.size(); ++k) {
     float dphi = phip-m_cellVec[k].phiPos();
     while (dphi >   M_PI) dphi -= 2*M_PI;
     while (dphi <= -M_PI) dphi += 2*M_PI;
     if (fabs(dphi) < dphi10) {
       float dz = fabs(zp - m_cellVec[k].getPosition().z());
       if (dz < (dzmin+0.001)) {
     dzmin     = dz;
     if (fabs(dphi) < (dphimin+0.01)) {
       cellIndex = k;
       dphimin   = fabs(dphi);
     } else {
       if (cellIndex >= m_cellVec.size()) cellIndex = k;
     }
       }
     }
   }
 #ifdef DebugLog
   std::cout << "getClosestCellIndex::Input " << zp << ":" << phip << " Index "
         << cellIndex;
   if (cellIndex < m_cellVec.size()) 
     std::cout << " Cell z " << m_cellVec[cellIndex].getPosition().z() << ":"
           << dzmin << " phi " << m_cellVec[cellIndex].phiPos() << ":"
           << dphimin;
   std::cout << std::endl;
 
 #endif
   return cellIndex;
 }

HGCalGeometry::CornersVec HGCalGeometry::getCorners ( const DetId & id ) const

Returns the corner points of this cell's volume.

Definition at line 158 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), i, HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, Square, topology(), and create_public_lumi_plots::xy.

                                                                        {
 
   HGCalGeometry::CornersVec co (8, GlobalPoint(0,0,0));
   unsigned int cellIndex =  indexFor(id);
   if (cellIndex <  m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(id);
     std::pair<float,float> xy;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       xy = topology().dddConstants().locateCell(id_.iCell,id_.iLay,id_.iSubSec,true);
     } else {
       xy = topology().dddConstants().locateCellHex(id_.iCell,id_.iSec,true);
     }
     float dz = m_cellVec[cellIndex].param()[0];
     float dx = 0.5*m_cellVec[cellIndex].param()[11];
     static const int signx[] = {-1,-1,1,1,-1,-1,1,1};
     static const int signy[] = {-1,1,1,-1,-1,1,1,-1};
     static const int signz[] = {-1,-1,-1,-1,1,1,1,1};
     for (unsigned int i = 0; i != 8; ++i) {
       const HepGeom::Point3D<float> lcoord(xy.first+signx[i]*dx,xy.second+signy[i]*dx,signz[i]*dz);
       co[i] = m_cellVec[cellIndex].getPosition(lcoord);
     }
   }
   return co;
 }

const CaloCellGeometry * HGCalGeometry::getGeometry ( const DetId & id ) const

overridevirtual

Get the cell geometry of a given detector id. Should return false if not found.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 111 of file HGCalGeometry.cc.

References cellGeomPtr(), HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), HGCEE, Square, and topology().

Referenced by PFHGCalRecHitCreator< DET, Layer, subdet >::importRecHits().

                                                                         {
   if (id == DetId()) return 0; // nothing to get
   DetId geoId;
   if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
     geoId = (id.subdetId() == HGCEE ? 
          (DetId)(HGCEEDetId(id).geometryCell()) : 
          (DetId)(HGCHEDetId(id).geometryCell()));
   } else {
     geoId = (DetId)(HGCalDetId(id).geometryCell());
   }
   const uint32_t cellIndex (topology().detId2denseGeomId(geoId));
   /*
   if (cellIndex <  m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(id);
     std::pair<float,float> xy = topology().dddConstants().locateCell(id_iCell,id_iLay,id_.iSubSec,true);
     const HepGeom::Point3D<float> lcoord(xy.first,xy.second,0);
     std::auto_ptr<FlatTrd> cellGeom(new FlatTrd(m_cellVec[cellIndex],lcoord));
     return cellGeom.release();
   }
   */
   return cellGeomPtr (cellIndex);
 
 }

GlobalPoint HGCalGeometry::getPosition ( const DetId & id ) const

Definition at line 135 of file HGCalGeometry.cc.

References gather_cfg::cout, HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, Square, topology(), and create_public_lumi_plots::xy.

Referenced by getClosestCellIndex().

                                                             {
 
   unsigned int cellIndex =  indexFor(id);
   GlobalPoint glob;
   if (cellIndex <  m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(id);
     std::pair<float,float> xy;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       xy = topology().dddConstants().locateCell(id_.iCell,id_.iLay,id_.iSubSec,true);
     } else {
       xy = topology().dddConstants().locateCellHex(id_.iCell,id_.iSec,true);
     }
     const HepGeom::Point3D<float> lcoord(xy.first,xy.second,0);
     glob = m_cellVec[cellIndex].getPosition(lcoord);
 #ifdef DebugLog
     std::cout << "getPosition:: index " << cellIndex << " Local " << lcoord.x()
           << ":" << lcoord.y() << " ID " << id_.iCell << ":" << id_.iLay 
           << " Global " << glob << std::endl;
 #endif
   } 
   return glob;
 }

void HGCalGeometry::getSummary	(	CaloSubdetectorGeometry::TrVec &	trVector,
		CaloSubdetectorGeometry::IVec &	iVector,
		CaloSubdetectorGeometry::DimVec &	dimVector,
		CaloSubdetectorGeometry::IVec &	dinsVector
	)		const

virtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 324 of file HGCalGeometry.cc.

References HGCalParameters::hgtrap::alpha, HGCalParameters::hgtrap::bl, cellGeomPtr(), HGCalParameters::hgtrap::cellSize, HGCalTopology::dddConstants(), HGCalParameters::hgtrap::dz, HGCalDDDConstants::getModule(), CaloCellGeometry::getPosition(), CaloCellGeometry::getTransform(), HGCalDDDConstants::getTrForm(), HGCalDDDConstants::getTrFormN(), HGCalParameters::hgtrap::h, HGCEE, i, relval_2017::k, k_NumberOfParametersPerShape, k_NumberOfShapes, HGCalParameters::hgtrform::lay, HGCEEDetId::layer(), HGCHEDetId::layer(), m_validGeomIds, CaloSubdetectorGeometry::numberOfParametersPerShape(), CaloSubdetectorGeometry::numberOfShapes(), CaloSubdetectorGeometry::numberOfTransformParms(), findQualityFiles::rr, HGCalDDDConstants::sectors(), Square, DetId::subdetId(), HGCalParameters::hgtrap::tl, topology(), HGCalTopology::totalGeomModules(), groupFilesInBlocks::tt, HGCalDDDConstants::volumes(), and HGCalDDDConstants::waferInLayer().

                                                                {
 
   unsigned int numberOfCells = topology().totalGeomModules(); // total Geom Modules both sides
   unsigned int numberOfShapes = HGCalGeometry::k_NumberOfShapes;
   unsigned int numberOfParametersPerShape = HGCalGeometry::k_NumberOfParametersPerShape;
 
   trVector.reserve( numberOfCells * numberOfTransformParms());
   iVector.reserve( numberOfCells );
   dimVector.reserve( numberOfShapes * numberOfParametersPerShape );
   dinsVector.reserve( numberOfCells );
   
   if (topology().geomMode() == HGCalGeometryMode::Square) {
     for (unsigned int k=0; k <topology().dddConstants().volumes(); ++k) {
       HGCalParameters::hgtrap vol = topology().dddConstants().getModule(k,false,true);
       ParmVec params( HGCalGeometry::k_NumberOfParametersPerShape, 0 );
       params[0] = vol.dz;
       params[1] = params[2] = 0;
       params[3] = params[7] = vol.h;
       params[4] = params[8] = vol.bl;
       params[5] = params[9] = vol.tl;
       params[6] = params[10]= vol.alpha;
       params[11]= vol.cellSize;
       dimVector.insert( dimVector.end(), params.begin(), params.end());
     }
   } else {
     for (unsigned itr=0; itr<topology().dddConstants().getTrFormN(); ++itr) {
       HGCalParameters::hgtrform mytr = topology().dddConstants().getTrForm(itr);
       int layer  = mytr.lay;
       for (int wafer=0; wafer<topology().dddConstants().sectors(); ++wafer) {
         if (topology().dddConstants().waferInLayer(wafer,layer,true)) {
       HGCalParameters::hgtrap vol = topology().dddConstants().getModule(wafer, true, true);
       ParmVec params( HGCalGeometry::k_NumberOfParametersPerShape, 0 );
       params[0] = vol.dz;
       params[1] = params[2] = 0;
       params[3] = params[7] = vol.h;
       params[4] = params[8] = vol.bl;
       params[5] = params[9] = vol.tl;
       params[6] = params[10]= vol.alpha;
       params[11]= vol.cellSize;
       dimVector.insert( dimVector.end(), params.begin(), params.end());
     }
       }
     }
   }
   
   for( unsigned int i( 0 ); i < numberOfCells; ++i )
   {
     DetId detId = m_validGeomIds[i];
     int layer = ((detId.subdetId() ==  ForwardSubdetector::HGCEE) ?
          (HGCEEDetId(detId).layer()) :
          (HGCHEDetId(detId).layer()));
     dinsVector.push_back( topology().detId2denseGeomId( detId ));
     iVector.push_back( layer );
     
     Tr3D tr;
     const CaloCellGeometry* ptr( cellGeomPtr( i ));
     if( 0 != ptr )
     {
       ptr->getTransform( tr, ( Pt3DVec* ) 0 );
 
       if( Tr3D() == tr ) // there is no rotation
       {
     const GlobalPoint& gp( ptr->getPosition()); 
     tr = HepGeom::Translate3D( gp.x(), gp.y(), gp.z());
       }
 
       const CLHEP::Hep3Vector tt( tr.getTranslation());
       trVector.push_back( tt.x());
       trVector.push_back( tt.y());
       trVector.push_back( tt.z());
       if( 6 == numberOfTransformParms())
       {
     const CLHEP::HepRotation rr( tr.getRotation());
     const ROOT::Math::Transform3D rtr( rr.xx(), rr.xy(), rr.xz(), tt.x(),
                        rr.yx(), rr.yy(), rr.yz(), tt.y(),
                        rr.zx(), rr.zy(), rr.zz(), tt.z());
     ROOT::Math::EulerAngles ea;
     rtr.GetRotation( ea );
     trVector.push_back( ea.Phi());
     trVector.push_back( ea.Theta());
     trVector.push_back( ea.Psi());
       }
     }
   }
 }

virtual const std::vector<DetId>& HGCalGeometry::getValidDetIds	(	DetId::Detector	det = `DetId::Detector(0)`,
		int	subdet = `0`
	)		const

inlineoverridevirtual

Get a list of valid detector ids (for the given subdetector)

Note: The implementation in this class is relevant for SubdetectorGeometries which handle only a single subdetector at a time. It does not look at the det and subdet arguments.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 76 of file HGCalGeometry.h.

References CaloSubdetectorGeometry::m_validIds.

76 { return m_validIds; }

CaloSubdetectorGeometry::m_validIds

std::vector< DetId > m_validIds

Definition: CaloSubdetectorGeometry.h:120

const std::vector<DetId>& HGCalGeometry::getValidGeomDetIds ( void ) const

inline

Definition at line 77 of file HGCalGeometry.h.

References m_validGeomIds.

77 { return m_validGeomIds; }

HGCalGeometry::m_validGeomIds

std::vector< DetId > m_validGeomIds

Definition: HGCalGeometry.h:113

unsigned int HGCalGeometry::indexFor ( const DetId & id ) const

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 230 of file HGCalGeometry.cc.

References gather_cfg::cout, TauDecayModes::dec, HGCalTopology::detId2denseGeomId(), HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), HGCEE, m_cellVec, DetId::rawId(), Square, and topology().

Referenced by getCorners(), and getPosition().

                                                           {
   unsigned int cellIndex =  m_cellVec.size();
   if (id != DetId()) {
     DetId geoId;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       geoId = (id.subdetId() == HGCEE ? 
            (DetId)(HGCEEDetId(id).geometryCell()) : 
            (DetId)(HGCHEDetId(id).geometryCell()));
     } else {
       geoId = (DetId)(HGCalDetId(id).geometryCell());
     }
     cellIndex = topology().detId2denseGeomId(geoId);
 #ifdef DebugLog
     std::cout << "indexFor " << std::hex << id.rawId() << ":" << geoId.rawId()
           << std::dec << " index " << cellIndex << std::endl;
 #endif
   }
   return cellIndex;
 }

void HGCalGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 42 of file HGCalGeometry.cc.

42 {

43 }

void HGCalGeometry::localCorners	(	Pt3DVec &	lc,
		const CCGFloat *	pv,
		unsigned int	i,
		Pt3D &	ref
	)

Definition at line 45 of file HGCalGeometry.cc.

References FlatTrd::localCorners().

                                       {
   FlatTrd::localCorners( lc, pv, ref ) ;
 }

void HGCalGeometry::newCell	(	const GlobalPoint &	f1,
		const GlobalPoint &	f2,
		const GlobalPoint &	f3,
		const CCGFloat *	parm,
		const DetId &	detId
	)

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 52 of file HGCalGeometry.cc.

References CaloSubdetectorGeometry::cornersMgr(), gather_cfg::cout, HGCalTopology::dddConstants(), TauDecayModes::dec, HGCalTopology::decode(), alcaDQMUpload::encode(), HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), HGCEE, m_cellVec, m_halfType, m_subdet, m_validGeomIds, CaloSubdetectorGeometry::m_validIds, HGCalDDDConstants::maxCells(), HGCalDDDConstants::numberCellsHexagon(), DetId::rawId(), Square, DetId::subdetId(), topology(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by HGCalGeometryLoader::buildGeom().

                                             {
 
   DetId geomId;
   int   cells;
   HGCalTopology::DecodedDetId id = topology().decode(detId);
   if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
     geomId = (detId.subdetId() == HGCEE ? 
           (DetId)(HGCEEDetId(detId).geometryCell()) :
           (DetId)(HGCHEDetId(detId).geometryCell()));
     cells = topology().dddConstants().maxCells(id.iLay,true);
   } else {
     geomId = (DetId)(HGCalDetId(detId).geometryCell());
     cells = topology().dddConstants().numberCellsHexagon(id.iSec);
   }
   const uint32_t cellIndex (topology().detId2denseGeomId(detId));
 
   m_cellVec.at( cellIndex ) = FlatTrd( cornersMgr(), f1, f2, f3, parm ) ;
   m_validGeomIds.at( cellIndex ) = geomId ;
 
 #ifdef DebugLog
   unsigned int nOld = m_validIds.size();
 #endif
   for (int cell = 0; cell < cells; ++cell) {
     id.iCell = cell;
     m_validIds.push_back( topology().encode(id));
     if ((topology().dddConstants().geomMode() == HGCalGeometryMode::Square) &&
     (!m_halfType)) {
       id.iSubSec = -id.iSubSec;
       m_validIds.push_back( topology().encode(id));
       id.iSubSec = -id.iSubSec;
     }
   }
 
 #ifdef DebugLog
   std::cout << "HGCalGeometry::newCell-> [" << cellIndex << "]"
         << " front:" << f1.x() << '/' << f1.y() << '/' << f1.z() 
             << " back:" <<  f2.x() << '/' << f2.y() << '/' << f2.z()
         << " eta|phi " << m_cellVec[cellIndex].etaPos() << ":"
         << m_cellVec[cellIndex].phiPos() << " id:";
   if (topology().dddConstants().geomMode() != HGCalGeometryMode::Square) {
     std::cout << HGCalDetId(detId);
   } else if (m_subdet == HGCEE) {
     std::cout << HGCEEDetId(detId);
   } else {
     std::cout << HGCHEDetId(detId);
   }
   unsigned int nNew = m_validIds.size();
   std::cout << " with valid DetId from " << nOld << " to " << nNew
         << std::endl; 
   std::cout << "Cell[" << cellIndex << "] " << std::hex << geomId.rawId() 
         << ":"  << m_validGeomIds[cellIndex].rawId() << std::dec
         << std::endl;
 #endif
 }

static std::string HGCalGeometry::producerTag ( )

inlinestatic

Definition at line 93 of file HGCalGeometry.h.

Referenced by PCaloGeometryBuilder::beginRun().

93 { return "HGCal" ; }

unsigned int HGCalGeometry::sizeForDenseIndex ( ) const

protected

Definition at line 250 of file HGCalGeometry.cc.

References topology(), and HGCalTopology::totalGeomModules().

                                                     {
   return topology().totalGeomModules();
 }

void HGCalGeometry::sortDetIds ( void )

Definition at line 318 of file HGCalGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

Referenced by HGCalGeometryLoader::build().

                                      {
   m_validIds.shrink_to_fit();
   std::sort( m_validIds.begin(), m_validIds.end(), rawIdSort());
 }

const HGCalTopology& HGCalGeometry::topology ( ) const

inline

Definition at line 96 of file HGCalGeometry.h.

References m_topology.

Referenced by getClosestCell(), getCorners(), getGeometry(), getPosition(), getSummary(), HGCalGeometry(), indexFor(), newCell(), and sizeForDenseIndex().

96 {return m_topology;}

HGCalGeometry::m_topology

const HGCalTopology & m_topology

Definition: HGCalGeometry.h:110

Member Data Documentation

CellVec HGCalGeometry::m_cellVec

private

Definition at line 112 of file HGCalGeometry.h.

Referenced by cellGeomPtr(), getClosestCell(), getClosestCellIndex(), getCorners(), getPosition(), indexFor(), and newCell().

bool HGCalGeometry::m_halfType

private

Definition at line 114 of file HGCalGeometry.h.

Referenced by newCell().

ForwardSubdetector HGCalGeometry::m_subdet

private

Definition at line 115 of file HGCalGeometry.h.

Referenced by cellElement(), and newCell().

const HGCalTopology& HGCalGeometry::m_topology

private

Definition at line 110 of file HGCalGeometry.h.

Referenced by topology().

std::vector<DetId> HGCalGeometry::m_validGeomIds

private

Definition at line 113 of file HGCalGeometry.h.

Referenced by cellGeomPtr(), getClosestCell(), getSummary(), getValidGeomDetIds(), and newCell().

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation