HGCalGeometry Class Reference

#include <HGCalGeometry.h>

Inheritance diagram for HGCalGeometry:

Public Types
enum	{ k_NumberOfParametersPerShape = 12 }
enum	{ k_NumberOfShapes = 50 }
typedef CaloCellGeometry::CCGFloat	CCGFloat
typedef std::vector< FlatTrd >	CellVec
typedef std::vector< GlobalPoint >	CornersVec
typedef std::set< DetId >	DetIdSet
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
	Get a list of all cells within a dR of the given cell. More...
virtual DetId	getClosestCell (const GlobalPoint &r) const
CornersVec	getCorners (const DetId &id) const
	Returns the corner points of this cell's volume. More...
virtual const CaloCellGeometry *	getGeometry (const DetId &id) const
	Get the cell geometry of a given detector id. Should return false if not found. More...
GlobalPoint	getPosition (const DetId &id) const
virtual const std::vector < DetId > &	getValidDetIds (DetId::Detector det=DetId::Detector(0), int subdet=0) const
	Get a list of valid detector ids (for the given subdetector) More...
	HGCalGeometry (const HGCalTopology &topology)
virtual void	initializeParms ()
virtual void	newCell (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId)
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 void	getSummary (TrVec &trVector, IVec &iVector, DimVec &dimVector, IVec &dinsVector) 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	producerTag ()

Protected Member Functions
void	addValidID (const DetId &id)
virtual const CaloCellGeometry *	cellGeomPtr (uint32_t index) const
unsigned int	getClosestCellIndex (const GlobalPoint &r) const
virtual unsigned int	indexFor (const DetId &id) const
virtual 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
std::vector< DetId >	m_validGeomIds
const HGCalTopology &	mTopology

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 25 of file HGCalGeometry.h.

Member Typedef Documentation

typedef CaloCellGeometry::CCGFloat HGCalGeometry::CCGFloat

Definition at line 31 of file HGCalGeometry.h.

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

Definition at line 29 of file HGCalGeometry.h.

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

Definition at line 36 of file HGCalGeometry.h.

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

Definition at line 35 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3D HGCalGeometry::Pt3D

Definition at line 32 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3DVec HGCalGeometry::Pt3DVec

Definition at line 33 of file HGCalGeometry.h.

Member Enumeration Documentation

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 38 of file HGCalGeometry.h.

{ k_NumberOfParametersPerShape = 12 } ; // FlatTrd

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 39 of file HGCalGeometry.h.

{ k_NumberOfShapes = 50 } ; 

Constructor & Destructor Documentation

HGCalGeometry::HGCalGeometry

(

const HGCalTopology &

topology

)

Definition at line 18 of file HGCalGeometry.cc.

References gather_cfg::cout, HGCalTopology::detectorType(), m_halfType, m_subdet, HGCalTopology::subDetector(), topology(), and HGCalTopology::totalGeomModules().

  : mTopology (topology_) {
  m_halfType = topology().detectorType();
  m_subdet   = topology().subDetector();
#ifdef DebugLog
  std::cout << "Expected total # of Geometry Modules " 
        << topology().totalGeomModules() << std::endl;
#endif
}

HGCalGeometry::~HGCalGeometry ( )

virtual

Definition at line 28 of file HGCalGeometry.cc.

{ }

Member Function Documentation

void HGCalGeometry::addValidID ( const DetId &  id )

protected

Definition at line 201 of file HGCalGeometry.cc.

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

std::string HGCalGeometry::cellElement

(

)

const

Definition at line 164 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

protectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 190 of file HGCalGeometry.cc.

References gather_cfg::cout, cmsHarvester::index, m_cellVec, m_validGeomIds, and CaloCellGeometry::param().

Referenced by getGeometry().

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

void HGCalGeometry::fillNamedParams ( DDFilteredView  fv )

virtual

Definition at line 30 of file HGCalGeometry.cc.

{}

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

virtual

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 160 of file HGCalGeometry.cc.

                                                                                     {
   HGCalGeometry::DetIdSet dss;
   return dss;
}

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

virtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 136 of file HGCalGeometry.cc.

References HGCalDDDConstants::assignCell(), gather_cfg::cout, HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::encode(), getClosestCellIndex(), HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, m_cellVec, m_validGeomIds, topology(), and HGCalTopology::DecodedDetId::zside.

                                                                {
  unsigned int cellIndex = getClosestCellIndex(r);
  if (cellIndex < m_cellVec.size()) {
    const HepGeom::Point3D<float> local = m_cellVec[cellIndex].getLocal(r);
    HGCalTopology::DecodedDetId id_ = topology().decode(m_validGeomIds[cellIndex]);
    std::pair<int,int> kxy = 
      topology().dddConstants().assignCell(local.x(),local.y(),id_.iLay,
                       id_.iSubSec,true);
    id_.iCell   = kxy.second;
    id_.iSubSec = kxy.first;
#ifdef DebugLog
    std::cout << "getClosestCell: local " << local << " Id " << id_.zside 
          << ":" << id_.iLay << ":" << id_.iSec << ":" << id_.iSubSec
          << ":" << id_.iCell << " Cell " << m_cellVec[cellIndex];
#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 206 of file HGCalGeometry.cc.

References gather_cfg::cout, getPosition(), relval_steps::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 116 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), i, HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSubSec, HGCalDDDConstants::locateCell(), m_cellVec, 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 = topology().dddConstants().locateCell(id_.iCell,id_.iLay,id_.iSubSec,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

virtual

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

Reimplemented from CaloSubdetectorGeometry.

Definition at line 79 of file HGCalGeometry.cc.

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

                                                                        {
  if (id == DetId()) return 0; // nothing to get
  DetId geoId = (id.subdetId() == HGCEE ? 
         (DetId)(HGCEEDetId(id).geometryCell()) : 
         (DetId)(HGCHEDetId(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 98 of file HGCalGeometry.cc.

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

Referenced by getClosestCellIndex().

                                                              {

  unsigned int cellIndex =  indexFor(id);
  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);
#ifdef DebugLog
    std::cout << "getPosition:: index " << cellIndex << " Local " << xy.first
          << ":" << xy.second << " ID " << id_.iCell << ":" << id_.iLay 
          << " Global " << m_cellVec[cellIndex].getPosition(lcoord)
          << " Cell" << m_cellVec[cellIndex];
#endif
    return m_cellVec[cellIndex].getPosition(lcoord);
  } 
  return GlobalPoint();
}

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

inlinevirtual

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 60 of file HGCalGeometry.h.

References CaloSubdetectorGeometry::m_validIds.

{return m_validIds;}

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

protectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 171 of file HGCalGeometry.cc.

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

Referenced by getCorners(), and getPosition().

                                                          {
  unsigned int cellIndex =  m_cellVec.size();
  if (id != DetId()) {
    DetId geoId = (id.subdetId() == HGCEE ? 
           (DetId)(HGCEEDetId(id).geometryCell()) : 
           (DetId)(HGCHEDetId(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 ( )

virtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 32 of file HGCalGeometry.cc.

                                    {
}

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

virtual

Implements CaloSubdetectorGeometry.

Definition at line 35 of file HGCalGeometry.cc.

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

Referenced by HGCalGeometryLoader::build().

                                            {
  const uint32_t cellIndex (topology().detId2denseGeomId(detId));
  DetId geomId = (detId.subdetId() == HGCEE ? 
          (DetId)(HGCEEDetId(detId).geometryCell()) :
          (DetId)(HGCHEDetId(detId).geometryCell()));
  if (cellIndex >= m_cellVec.size())      m_cellVec.resize (cellIndex+1);
  if (cellIndex >= m_validGeomIds.size()) m_validGeomIds.resize (cellIndex+1);
  m_cellVec     [cellIndex] = FlatTrd( cornersMgr(), f1, f2, f3, parm ) ;
  m_validGeomIds[cellIndex] = geomId ;

  HGCalTopology::DecodedDetId id = topology().decode(detId);
  int cells = topology().dddConstants().maxCells(id.iLay,true);
  unsigned int nOld = m_validIds.size();
  unsigned int nNew = nOld + (unsigned int)(m_halfType ? cells : 2*cells);
  m_validIds.resize(nNew);
  for (int cell = 0; cell < cells; ++cell) {
    id.iCell = cell;
    m_validIds[nOld+cell] = topology().encode(id);
    if (!m_halfType) {
      id.iSubSec = -id.iSubSec;
      m_validIds[nOld+cells+cell] = 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 (m_subdet == HGCEE) std::cout << HGCEEDetId(detId);
  else                   std::cout << HGCHEDetId(detId);
  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 << " "
        << m_cellVec[cellIndex];
#endif
}

static std::string HGCalGeometry::producerTag ( )

inlinestatic

Definition at line 76 of file HGCalGeometry.h.

{ return "HGCal" ; }

unsigned int HGCalGeometry::sizeForDenseIndex ( ) const

protectedvirtual

Definition at line 186 of file HGCalGeometry.cc.

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

                                                    {
  return topology().totalGeomModules();
}

const HGCalTopology& HGCalGeometry::topology ( ) const

inline

Definition at line 79 of file HGCalGeometry.h.

References mTopology.

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

{return mTopology;}

Member Data Documentation

CellVec HGCalGeometry::m_cellVec

private

Definition at line 94 of file HGCalGeometry.h.

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

bool HGCalGeometry::m_halfType

private

Definition at line 96 of file HGCalGeometry.h.

Referenced by HGCalGeometry(), and newCell().

ForwardSubdetector HGCalGeometry::m_subdet

private

Definition at line 97 of file HGCalGeometry.h.

Referenced by cellElement(), HGCalGeometry(), and newCell().

std::vector<DetId> HGCalGeometry::m_validGeomIds

private

Definition at line 95 of file HGCalGeometry.h.

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

const HGCalTopology& HGCalGeometry::mTopology

private

Definition at line 92 of file HGCalGeometry.h.

Referenced by topology().