#include <FastTimeGeometry.h>

Inheritance diagram for FastTimeGeometry:

Public Types
enum	{ k_NumberOfParametersPerShape = 12 }

enum	{ k_NumberOfShapes = 50 }

typedef FastTimeGeometryRecord	AlignedRecord

typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< FlatTrd >	CellVec

typedef std::vector< GlobalPoint >	CornersVec

typedef std::set< DetId >	DetIdSet

typedef PFastTimeRcd	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

	FastTimeGeometry (const FastTimeTopology &topology)

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 override

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

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 FastTimeTopology &	topology () const

virtual	~FastTimeGeometry ()

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

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

ForwardSubdetector	m_subdet

const FastTimeTopology &	m_topology

int	m_Type

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 FastTimeGeometry.h.

Member Typedef Documentation

typedef FastTimeGeometryRecord FastTimeGeometry::AlignedRecord

Definition at line 39 of file FastTimeGeometry.h.

typedef CaloCellGeometry::CCGFloat FastTimeGeometry::CCGFloat

Definition at line 32 of file FastTimeGeometry.h.

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

Definition at line 30 of file FastTimeGeometry.h.

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

Definition at line 37 of file FastTimeGeometry.h.

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

Definition at line 36 of file FastTimeGeometry.h.

typedef PFastTimeRcd FastTimeGeometry::PGeometryRecord

Definition at line 40 of file FastTimeGeometry.h.

typedef CaloCellGeometry::Pt3D FastTimeGeometry::Pt3D

Definition at line 33 of file FastTimeGeometry.h.

typedef CaloCellGeometry::Pt3DVec FastTimeGeometry::Pt3DVec

Definition at line 34 of file FastTimeGeometry.h.

Member Enumeration Documentation

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 42 of file FastTimeGeometry.h.

42 { k_NumberOfParametersPerShape = 12 } ; // FlatTrd

FastTimeGeometry::k_NumberOfParametersPerShape

Definition: FastTimeGeometry.h:42

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 43 of file FastTimeGeometry.h.

43 { k_NumberOfShapes = 50 } ;

FastTimeGeometry::k_NumberOfShapes

Definition: FastTimeGeometry.h:43

Constructor & Destructor Documentation

FastTimeGeometry::FastTimeGeometry ( const FastTimeTopology & topology )

Definition at line 17 of file FastTimeGeometry.cc.

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

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

FastTimeGeometry::~FastTimeGeometry ( )

virtual

Definition at line 31 of file FastTimeGeometry.cc.

31 { }

Member Function Documentation

void FastTimeGeometry::addValidID ( const DetId & id )

protected

Definition at line 168 of file FastTimeGeometry.cc.

                                                  {
   edm::LogError("FastTimeGeom") << "FastTimeGeometry::addValidID is not implemented";
 }

std::string FastTimeGeometry::cellElement ( ) const

Definition at line 134 of file FastTimeGeometry.cc.

References m_Type.

                                               {
   if      (m_Type == 1) return "FastTimeBarrel";
   else if (m_Type == 2) return "FastTimeEndcap";
   else                  return "Unknown";
 }

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

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 157 of file FastTimeGeometry.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 EDM_ML_DEBUG
   //  std::cout << "cellGeomPtr " << m_cellVec[index];
 #endif
   if (0 == cell->param()) return 0;
   return cell;
 }

static std::string FastTimeGeometry::dbString ( )

inlinestatic

Definition at line 45 of file FastTimeGeometry.h.

45 { return "PFastTimeRcd" ; }

void FastTimeGeometry::fillNamedParams ( DDFilteredView fv )

virtual

Definition at line 33 of file FastTimeGeometry.cc.

33 {}

FastTimeGeometry::DetIdSet FastTimeGeometry::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 129 of file FastTimeGeometry.cc.

                                                                                           {
    FastTimeGeometry::DetIdSet dss;
    return dss;
 }

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 109 of file FastTimeGeometry.cc.

References gather_cfg::cout, FastTimeTopology::dddConstants(), FastTimeDDDConstants::getEtaPhi(), FastTimeDDDConstants::getZPhi(), m_Type, PV3DBase< T, PVType, FrameType >::perp(), phi, PV3DBase< T, PVType, FrameType >::phi(), topology(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and ecaldqm::zside().

                                                                  {
   int zside = (r.z() > 0) ? 1 : -1;
   std::pair<int,int> etaZPhi;
   if (m_Type == 1) {
     double zz = (zside > 0) ? r.z() : -r.z();
     etaZPhi = topology().dddConstants().getZPhi(zz,r.phi());
   } else {
     double phi = (zside > 0) ? r.phi() : atan2(r.y(),-r.x());
     etaZPhi = topology().dddConstants().getEtaPhi(r.perp(),phi);
   }
   FastTimeDetId id = FastTimeDetId(m_Type,etaZPhi.first,etaZPhi.second,zside);
 #ifdef EDM_ML_DEBUG
   std::cout << "getClosestCell: for (" << r.x() << ", " << r.y() << ", "
         << r.z() << ")  Id " << id.type() << ":" << id.zside() << ":" 
         << id.ieta() << ":" << id.iphi() << std::endl;
 #endif
 
   return (topology().valid(id) ? DetId(id) : DetId());
 }

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

Returns the corner points of this cell's volume.

Definition at line 103 of file FastTimeGeometry.cc.

References FastTimeTopology::dddConstants(), FastTimeDDDConstants::getCorners(), FastTimeDetId::ieta(), FastTimeDetId::iphi(), m_Type, topology(), and FastTimeDetId::zside().

                                                                              {
 
   FastTimeDetId id_ = FastTimeDetId(id);
   return topology().dddConstants().getCorners(m_Type,id_.ieta(),id_.iphi(),id_.zside());
 }

const CaloCellGeometry * FastTimeGeometry::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 89 of file FastTimeGeometry.cc.

References cellGeomPtr(), FastTimeDetId::geometryCell(), and topology().

                                                                            {
 
   if (id == DetId()) return 0; // nothing to get
   DetId geoId = (DetId)(FastTimeDetId(id).geometryCell());
   const uint32_t cellIndex (topology().detId2denseGeomId(geoId));
   return cellGeomPtr (cellIndex);
 }

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

Definition at line 97 of file FastTimeGeometry.cc.

References FastTimeTopology::dddConstants(), FastTimeDDDConstants::getPosition(), FastTimeDetId::ieta(), FastTimeDetId::iphi(), m_Type, topology(), and FastTimeDetId::zside().

                                                                {
 
   FastTimeDetId id_ = FastTimeDetId(id);
   return topology().dddConstants().getPosition(m_Type,id_.ieta(),id_.iphi(),id_.zside());
 }

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 187 of file FastTimeGeometry.cc.

References cellGeomPtr(), FastTimeTopology::dddConstants(), CaloCellGeometry::getPosition(), FastTimeDDDConstants::getRin(), FastTimeDDDConstants::getRout(), CaloCellGeometry::getTransform(), FastTimeDDDConstants::getZHalf(), runTauDisplay::gp, i, relval_2017::k, k_NumberOfParametersPerShape, k_NumberOfShapes, m_Type, m_validGeomIds, CaloSubdetectorGeometry::numberOfParametersPerShape(), CaloSubdetectorGeometry::numberOfShapes(), CaloSubdetectorGeometry::numberOfTransformParms(), findQualityFiles::rr, topology(), FastTimeTopology::totalGeomModules(), and groupFilesInBlocks::tt.

                                                                   {
 
   unsigned int numberOfCells = topology().totalGeomModules(); // total Geom Modules both sides
   unsigned int numberOfShapes = FastTimeGeometry::k_NumberOfShapes;
   unsigned int numberOfParametersPerShape = FastTimeGeometry::k_NumberOfParametersPerShape;
 
   trVector.reserve( numberOfCells * numberOfTransformParms());
   iVector.reserve( numberOfCells );
   dimVector.reserve( numberOfShapes * numberOfParametersPerShape );
   dinsVector.reserve( numberOfCells );
   
   for (unsigned int k=0; k <topology().totalGeomModules(); ++k) {
     ParmVec params( FastTimeGeometry::k_NumberOfParametersPerShape, 0 );
     params[0] = topology().dddConstants().getZHalf(m_Type);
     params[1] = params[2] = 0;
     params[3] = params[7] = topology().dddConstants().getRin(m_Type);
     params[4] = params[8] = topology().dddConstants().getRout(m_Type);
     params[5] = params[9] = topology().dddConstants().getRout(m_Type);
     params[6] = params[10]= 0;
     params[11]= (k == 0) ? 1.0 : -1.0;
     dimVector.insert( dimVector.end(), params.begin(), params.end());
   }
   
   for (unsigned int i( 0 ); i < numberOfCells; ++i) {
     DetId detId = m_validGeomIds[i];
     dinsVector.push_back( topology().detId2denseGeomId( detId ));
     iVector.push_back(1);
     
     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>& FastTimeGeometry::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 FastTimeGeometry.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>& FastTimeGeometry::getValidGeomDetIds ( void ) const

inline

Definition at line 77 of file FastTimeGeometry.h.

References m_validGeomIds.

77 { return m_validGeomIds; }

FastTimeGeometry::m_validGeomIds

std::vector< DetId > m_validGeomIds

Definition: FastTimeGeometry.h:113

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

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 140 of file FastTimeGeometry.cc.

References gather_cfg::cout, TauDecayModes::dec, FastTimeTopology::detId2denseGeomId(), FastTimeDetId::geometryCell(), m_cellVec, DetId::rawId(), and topology().

                                                              {
   unsigned int cellIndex =  m_cellVec.size();
   if (id != DetId()) {
     DetId geoId = (DetId)(FastTimeDetId(id).geometryCell());
     cellIndex = topology().detId2denseGeomId(geoId);
 #ifdef EDM_ML_DEBUG
     std::cout << "indexFor " << std::hex << id.rawId() << ":" << geoId.rawId()
           << std::dec << " index " << cellIndex << std::endl;
 #endif
   }
   return cellIndex;
 }

void FastTimeGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 35 of file FastTimeGeometry.cc.

35 {

36 }

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

Definition at line 38 of file FastTimeGeometry.cc.

References FlatTrd::localCorners().

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

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

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 45 of file FastTimeGeometry.cc.

References CaloSubdetectorGeometry::cornersMgr(), gather_cfg::cout, FastTimeTopology::dddConstants(), TauDecayModes::dec, FastTimeTopology::decode(), FastTimeTopology::encode(), FastTimeDetId::geometryCell(), m_cellVec, m_Type, m_validGeomIds, CaloSubdetectorGeometry::m_validIds, HLT_25ns10e33_v2_cff::nPhi, FastTimeDDDConstants::numberEtaZ(), FastTimeDDDConstants::numberPhi(), phi, DetId::rawId(), topology(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by FastTimeGeometryLoader::buildGeom().

                                            {
 
   FastTimeTopology::DecodedDetId id = topology().decode(detId);
   DetId geomId = (DetId)(FastTimeDetId(detId).geometryCell());
   int nEtaZ = topology().dddConstants().numberEtaZ(m_Type);
   int nPhi  = topology().dddConstants().numberPhi(m_Type);
 
   const uint32_t cellIndex (topology().detId2denseGeomId(detId));
 
   m_cellVec.at( cellIndex ) = FlatTrd( cornersMgr(), f1, f2, f3, parm ) ;
   m_validGeomIds.at( cellIndex ) = geomId ;
 
 #ifdef EDM_ML_DEBUG
   unsigned int nOld = m_validIds.size();
 #endif
   for (int etaZ = 1; etaZ <= nEtaZ; ++etaZ) {
     id.iEtaZ = etaZ;
     for (int phi = 1; phi <= nPhi; ++phi) {
       id.iPhi  = phi;
       DetId idc = topology().encode(id);
       if (topology().valid(idc)) {
     m_validIds.push_back(idc);
       }
     }
   }
 
 #ifdef EDM_ML_DEBUG
   std::cout << "FastTimeGeometry::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:" << FastTimeDetId(detId)
         << " with valid DetId from " << nOld << " to " << m_validIds.size()
         << std::endl; 
   std::cout << "Cell[" << cellIndex << "] " << std::hex << geomId.rawId() 
         << ":"  << m_validGeomIds[cellIndex].rawId() << std::dec
         << std::endl;
 #endif
 }

static std::string FastTimeGeometry::producerTag ( )

inlinestatic

Definition at line 93 of file FastTimeGeometry.h.

93 { return "FastTime" ; }

unsigned int FastTimeGeometry::sizeForDenseIndex ( ) const

protected

Definition at line 153 of file FastTimeGeometry.cc.

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

                                                        {
   return topology().totalGeomModules();
 }

void FastTimeGeometry::sortDetIds ( void )

Definition at line 182 of file FastTimeGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

Referenced by FastTimeGeometryLoader::build().

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

const FastTimeTopology& FastTimeGeometry::topology ( ) const

inline

Definition at line 96 of file FastTimeGeometry.h.

References m_topology.

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

96 {return m_topology;}

FastTimeGeometry::m_topology

const FastTimeTopology & m_topology

Definition: FastTimeGeometry.h:110

Member Data Documentation

CellVec FastTimeGeometry::m_cellVec

private

Definition at line 112 of file FastTimeGeometry.h.

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

ForwardSubdetector FastTimeGeometry::m_subdet

private

Definition at line 115 of file FastTimeGeometry.h.

const FastTimeTopology& FastTimeGeometry::m_topology

private

Definition at line 110 of file FastTimeGeometry.h.

Referenced by topology().

int FastTimeGeometry::m_Type

private

Definition at line 114 of file FastTimeGeometry.h.

Referenced by cellElement(), getClosestCell(), getCorners(), getPosition(), getSummary(), and newCell().

std::vector<DetId> FastTimeGeometry::m_validGeomIds

private

Definition at line 113 of file FastTimeGeometry.h.

Referenced by cellGeomPtr(), 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