#include <EcalEndcapGeometry.h>

Inheritance diagram for EcalEndcapGeometry:

Public Types
typedef EcalEndcapGeometryRecord	AlignedRecord

typedef EEAlignmentRcd	AlignmentRecord

typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< TruncatedPyramid >	CellVec

typedef EEDetId	DetIdType

typedef IdealGeometryRecord	IdealRecord

typedef EcalEndcapNumberingScheme	NumberingScheme

typedef EZArrayFL< EBDetId >	OrderedListOfEBDetId

typedef PEcalEndcapRcd	PGeometryRecord

typedef CaloCellGeometry::Pt3D	Pt3D

typedef CaloCellGeometry::Pt3DVec	Pt3DVec

typedef std::vector< OrderedListOfEBDetId * >	VecOrdListEBDetIdPtr

Public Types inherited from CaloSubdetectorGeometry
typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< std::shared_ptr< const CaloCellGeometry > >	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
CCGFloat	avgAbsZFrontFaceCenter () const

	EcalEndcapGeometry ()

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

const OrderedListOfEBDetId *	getClosestBarrelCells (EEDetId id) const

DetId	getClosestCell (const GlobalPoint &r) const override

int	getNumberOfCrystalPerModule () const

int	getNumberOfModules () const

void	initializeParms () override

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

unsigned int	numberOfParametersPerShape () const override

unsigned int	numberOfShapes () const override

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

void	setNumberOfCrystalPerModule (const int nncrys)

void	setNumberOfModules (const int nnmods)

	~EcalEndcapGeometry () override

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

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

	CaloSubdetectorGeometry ()

	CaloSubdetectorGeometry (const CaloSubdetectorGeometry &)=delete
	avoid copies More...

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 std::shared_ptr< const CaloCellGeometry >	getGeometry (const DetId &id) const
	Get the cell geometry of a given detector id. Should return false if not found. More...

virtual void	getSummary (TrVec &trVector, IVec &iVector, DimVec &dimVector, IVec &dinsVector) 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...

virtual unsigned int	numberOfTransformParms () const

CaloSubdetectorGeometry &	operator= (const CaloSubdetectorGeometry &)=delete

ParMgr *	parMgr ()

const ParMgr *	parMgrConst () const

ParVecVec &	parVecVec ()

const ParVecVec &	parVecVec () const

virtual bool	valid (const DetId &id) const

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

Static Public Member Functions
static unsigned int	alignmentTransformIndexGlobal (const DetId &id)

static unsigned int	alignmentTransformIndexLocal (const DetId &id)

static std::string	dbString ()

static DetId	detIdFromLocalAlignmentIndex (unsigned int iLoc)

static std::string	hitString ()

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

static unsigned int	numberOfAlignments ()

static std::string	producerTag ()

Static Public Attributes
static constexpr int	k_NumberOfCellsForCorners = EEDetId::kSizeForDenseIndexing

static constexpr int	k_NumberOfParametersPerShape = 11

static constexpr int	k_NumberOfShapes = 1

Protected Member Functions
const CaloCellGeometry *	getGeometryRawPtr (uint32_t index) const override

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

virtual std::shared_ptr< const CaloCellGeometry >	cellGeomPtr (uint32_t index) const

virtual unsigned int	indexFor (const DetId &id) const

virtual unsigned int	sizeForDenseIndex (const DetId &id) const

Private Member Functions
EEDetId	gId (float x, float y, float z) const

unsigned int	xindex (CCGFloat x, CCGFloat z) const

unsigned int	yindex (CCGFloat y, CCGFloat z) const

Static Private Member Functions
static int	myPhi (int i)

Private Attributes
int	_nncrys
	number of crystals per module More...

int	_nnmods
	number of modules More...

CCGFloat	m_avgZ

std::atomic< EZMgrFL< EBDetId > * >	m_borderMgr

std::atomic< VecOrdListEBDetIdPtr * >	m_borderPtrVec

CellVec	m_cellVec

std::atomic< bool >	m_check

CCGFloat	m_del

unsigned int	m_nref

CCGFloat	m_wref

CCGFloat	m_xhi [2]

CCGFloat	m_xlo [2]

CCGFloat	m_xoff [2]

CCGFloat	m_yhi [2]

CCGFloat	m_ylo [2]

CCGFloat	m_yoff [2]

CCGFloat	zeN

CCGFloat	zeP

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 20 of file EcalEndcapGeometry.h.

Member Typedef Documentation

◆ AlignedRecord

typedef EcalEndcapGeometryRecord EcalEndcapGeometry::AlignedRecord

Definition at line 29 of file EcalEndcapGeometry.h.

◆ AlignmentRecord

typedef EEAlignmentRcd EcalEndcapGeometry::AlignmentRecord

Definition at line 30 of file EcalEndcapGeometry.h.

◆ CCGFloat

typedef CaloCellGeometry::CCGFloat EcalEndcapGeometry::CCGFloat

Definition at line 24 of file EcalEndcapGeometry.h.

◆ CellVec

typedef std::vector<TruncatedPyramid> EcalEndcapGeometry::CellVec

Definition at line 22 of file EcalEndcapGeometry.h.

◆ DetIdType

typedef EEDetId EcalEndcapGeometry::DetIdType

Definition at line 39 of file EcalEndcapGeometry.h.

◆ IdealRecord

typedef IdealGeometryRecord EcalEndcapGeometry::IdealRecord

Definition at line 28 of file EcalEndcapGeometry.h.

◆ NumberingScheme

typedef EcalEndcapNumberingScheme EcalEndcapGeometry::NumberingScheme

Definition at line 37 of file EcalEndcapGeometry.h.

◆ OrderedListOfEBDetId

typedef EZArrayFL<EBDetId> EcalEndcapGeometry::OrderedListOfEBDetId

Definition at line 33 of file EcalEndcapGeometry.h.

◆ PGeometryRecord

typedef PEcalEndcapRcd EcalEndcapGeometry::PGeometryRecord

Definition at line 31 of file EcalEndcapGeometry.h.

◆ Pt3D

typedef CaloCellGeometry::Pt3D EcalEndcapGeometry::Pt3D

Definition at line 25 of file EcalEndcapGeometry.h.

◆ Pt3DVec

typedef CaloCellGeometry::Pt3DVec EcalEndcapGeometry::Pt3DVec

Definition at line 26 of file EcalEndcapGeometry.h.

◆ VecOrdListEBDetIdPtr

typedef std::vector<OrderedListOfEBDetId*> EcalEndcapGeometry::VecOrdListEBDetIdPtr

Definition at line 35 of file EcalEndcapGeometry.h.

Constructor & Destructor Documentation

◆ EcalEndcapGeometry()

EcalEndcapGeometry::EcalEndcapGeometry ( void )

Definition at line 17 of file EcalEndcapGeometry.cc.

References m_xhi, m_xlo, m_xoff, m_yhi, m_ylo, and m_yoff.

     : _nnmods(316),
       _nncrys(25),
       zeP(0.),
       zeN(0.),
       m_wref(0.),
       m_del(0.),
       m_nref(0),
       m_borderMgr(nullptr),
       m_borderPtrVec(nullptr),
       m_avgZ(-1),
       m_check(false),
       m_cellVec(k_NumberOfCellsForCorners) {
   m_xlo[0] = 999.;
   m_xlo[1] = 999.;
   m_xhi[0] = -999.;
   m_xhi[1] = -999.;
   m_ylo[0] = 999.;
   m_ylo[1] = 999.;
   m_yhi[0] = -999.;
   m_yhi[1] = -999.;
   m_xoff[0] = 0.;
   m_xoff[1] = 0.;
   m_yoff[0] = 0.;
   m_yoff[1] = 0.;
 }

◆ ~EcalEndcapGeometry()

EcalEndcapGeometry::~EcalEndcapGeometry ( )

override

Definition at line 44 of file EcalEndcapGeometry.cc.

References m_borderMgr, m_borderPtrVec, and findQualityFiles::v.

                                         {
   if (m_borderPtrVec) {
     auto ptr = m_borderPtrVec.load(std::memory_order_acquire);
     for (auto& v : (*ptr)) {
       delete v;
       v = nullptr;
     }
     delete m_borderPtrVec.load();
   }
   delete m_borderMgr.load();
 }

Member Function Documentation

◆ alignmentTransformIndexGlobal()

unsigned int EcalEndcapGeometry::alignmentTransformIndexGlobal ( const DetId & id )

static

Definition at line 69 of file EcalEndcapGeometry.cc.

References DetId::Ecal.

                                                                              {
   return (unsigned int)DetId::Ecal - 1;
 }

◆ alignmentTransformIndexLocal()

unsigned int EcalEndcapGeometry::alignmentTransformIndexLocal ( const DetId & id )

static

Definition at line 56 of file EcalEndcapGeometry.cc.

References cms::cuda::assert(), CaloGenericDetId::isEE(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), and ecaldqm::zside().

                                                                              {
   const CaloGenericDetId gid(id);
 
   assert(gid.isEE());
   unsigned int index(EEDetId(id).ix() / 51 + (EEDetId(id).zside() < 0 ? 0 : 2));
 
   return index;
 }

◆ avgAbsZFrontFaceCenter()

CCGFloat EcalEndcapGeometry::avgAbsZFrontFaceCenter ( ) const

Definition at line 435 of file EcalEndcapGeometry.cc.

References CaloSubdetectorGeometry::cellGeomPtr(), mps_fire::i, m_avgZ, m_cellVec, and m_check.

                                                           {
   if (!m_check.load(std::memory_order_acquire)) {
     CCGFloat sum(0);
     for (unsigned int i(0); i != m_cellVec.size(); ++i) {
       auto cell = cellGeomPtr(i);
       if (nullptr != cell) {
         sum += fabs(cell->getPosition().z());
       }
     }
     m_avgZ = sum / m_cellVec.size();
     m_check.store(true, std::memory_order_release);
   }
   return m_avgZ;
 }

◆ dbString()

static std::string EcalEndcapGeometry::dbString ( )

inlinestatic

Definition at line 47 of file EcalEndcapGeometry.h.

47 { return "PEcalEndcapRcd"; }

◆ detIdFromLocalAlignmentIndex()

DetId EcalEndcapGeometry::detIdFromLocalAlignmentIndex ( unsigned int iLoc )

static

Definition at line 65 of file EcalEndcapGeometry.cc.

Referenced by FakeCaloAlignmentEP::produceEEAli(), and TestCaloAlignmentEP::produceEEAli().

                                                                         {
   return EEDetId(20 + 50 * (iLoc % 2), 50, 2 * (iLoc / 2) - 1);
 }

◆ getCells()

CaloSubdetectorGeometry::DetIdSet EcalEndcapGeometry::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 313 of file EcalEndcapGeometry.cc.

References funct::abs(), reco::deltaR2(), HGC3DClusterGenMatchSelector_cfi::dR, ALPAKA_ACCELERATOR_NAMESPACE::dR2(), PVValHelper::eta, CaloCellGeometry::etaPos(), CaloSubdetectorGeometry::getCells(), l1ctLayer2EG_cff::id, m_cellVec, m_nref, M_PI, m_wref, m_xhi, m_xlo, m_yhi, m_ylo, phi, CaloCellGeometry::phiPos(), alignCSCRings::r, corrVsCorr::rx, corrVsCorr::ry, funct::tan(), EEDetId::validDetId(), xindex(), geometryCSVtoXML::xx, yindex(), geometryCSVtoXML::yy, zeP, and geometryCSVtoXML::zz.

                                                                                                   {
   CaloSubdetectorGeometry::DetIdSet dis;  // return object
   if (0.000001 < dR) {
     if (dR > M_PI / 2.)  // this code assumes small dR
     {
       dis = CaloSubdetectorGeometry::getCells(r, dR);  // base class version
     } else {
       const float dR2(dR * dR);
       const float reta(r.eta());
       const float rphi(r.phi());
       const float rx(r.x());
       const float ry(r.y());
       const float rz(r.z());
       const float fac(std::abs(zeP / rz));
       const float xx(rx * fac);  // xyz at endcap z
       const float yy(ry * fac);
       const float zz(rz * fac);
 
       const float xang(std::atan(xx / zz));
       const float lowX(zz > 0 ? zz * std::tan(xang - dR) : zz * std::tan(xang + dR));
       const float highX(zz > 0 ? zz * std::tan(xang + dR) : zz * std::tan(xang - dR));
       const float yang(std::atan(yy / zz));
       const float lowY(zz > 0 ? zz * std::tan(yang - dR) : zz * std::tan(yang + dR));
       const float highY(zz > 0 ? zz * std::tan(yang + dR) : zz * std::tan(yang - dR));
 
       const float refxlo(0 > rz ? m_xlo[0] : m_xlo[1]);
       const float refxhi(0 > rz ? m_xhi[0] : m_xhi[1]);
       const float refylo(0 > rz ? m_ylo[0] : m_ylo[1]);
       const float refyhi(0 > rz ? m_yhi[0] : m_yhi[1]);
 
       if (lowX < refxhi &&  // proceed if any possible overlap with the endcap
           lowY < refyhi && highX > refxlo && highY > refylo) {
         const int ix_ctr(xindex(xx, rz));
         const int iy_ctr(yindex(yy, rz));
         const int iz(rz > 0 ? 1 : -1);
 
         const int ix_hi(ix_ctr + int((highX - xx) / m_wref) + 2);
         const int ix_lo(ix_ctr - int((xx - lowX) / m_wref) - 2);
 
         const int iy_hi(iy_ctr + int((highY - yy) / m_wref) + 2);
         const int iy_lo(iy_ctr - int((yy - lowY) / m_wref) - 2);
 
         for (int kx(ix_lo); kx <= ix_hi; ++kx) {
           if (kx > 0 && kx <= (int)m_nref) {
             for (int ky(iy_lo); ky <= iy_hi; ++ky) {
               if (ky > 0 && ky <= (int)m_nref) {
                 if (EEDetId::validDetId(kx, ky, iz))  // reject invalid ids
                 {
                   const EEDetId id(kx, ky, iz);
                   const CaloCellGeometry* cell(&m_cellVec[id.denseIndex()]);
                   const float eta(cell->etaPos());
                   const float phi(cell->phiPos());
                   if (reco::deltaR2(eta, phi, reta, rphi) < dR2)
                     dis.insert(id);
                 }
               }
             }
           }
         }
       }
     }
   }
   return dis;
 }

◆ getClosestBarrelCells()

const EcalEndcapGeometry::OrderedListOfEBDetId * EcalEndcapGeometry::getClosestBarrelCells ( EEDetId id ) const

Definition at line 378 of file EcalEndcapGeometry.cc.

References degrees(), CaloSubdetectorGeometry::getGeometry(), mps_fire::i, l1tPhase2CaloJetEmulator_cfi::iEta, m_borderMgr, m_borderPtrVec, myPhi(), phi, nano_mu_digi_cff::rawId, and ecaldqm::zside().

Referenced by spr::simpleMove(), and PFECALHashNavigator::stdsimplemove().

                                                                                                         {
   OrderedListOfEBDetId* ptr(nullptr);
   auto ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
   if (!ptrVec) {
     if (0 != id.rawId() && nullptr != getGeometry(id)) {
       const float phi(370. + getGeometry(id)->getPosition().phi().degrees());
       const int iPhi(1 + int(phi) % 360);
       const int iz(id.zside());
       if (!m_borderMgr.load(std::memory_order_acquire)) {
         EZMgrFL<EBDetId>* expect = nullptr;
         auto ptrMgr = new EZMgrFL<EBDetId>(720 * 9, 9);
         bool exchanged = m_borderMgr.compare_exchange_strong(expect, ptrMgr, std::memory_order_acq_rel);
         if (!exchanged)
           delete ptrMgr;
       }
       VecOrdListEBDetIdPtr* expect = nullptr;
       auto ptrVec = new VecOrdListEBDetIdPtr();
       ptrVec->reserve(720);
       for (unsigned int i(0); i != 720; ++i) {
         const int kz(360 > i ? -1 : 1);
         const int iEta(kz * 85);
         const int iEtam1(kz * 84);
         const int iEtam2(kz * 83);
         const int jPhi(i % 360 + 1);
         OrderedListOfEBDetId& olist(*new OrderedListOfEBDetId(m_borderMgr.load(std::memory_order_acquire)));
         olist[0] = EBDetId(iEta, jPhi);
         olist[1] = EBDetId(iEta, myPhi(jPhi + 1));
         olist[2] = EBDetId(iEta, myPhi(jPhi - 1));
         olist[3] = EBDetId(iEtam1, jPhi);
         olist[4] = EBDetId(iEtam1, myPhi(jPhi + 1));
         olist[5] = EBDetId(iEtam1, myPhi(jPhi - 1));
         olist[6] = EBDetId(iEta, myPhi(jPhi + 2));
         olist[7] = EBDetId(iEta, myPhi(jPhi - 2));
         olist[8] = EBDetId(iEtam2, jPhi);
         ptrVec->emplace_back(&olist);
       }
       bool exchanged = m_borderPtrVec.compare_exchange_strong(expect, ptrVec, std::memory_order_acq_rel);
       if (!exchanged)
         delete ptrVec;
       ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
       ptr = (*ptrVec)[(iPhi - 1) + (0 > iz ? 0 : 360)];
     }
   }
   return ptr;
 }

◆ getClosestCell()

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 216 of file EcalEndcapGeometry.cc.

References A, B, correctionTermsCaloMet_cff::C, HLT_2024v11_cff::distance, MillePedeFileConverter_cfg::e, CaloSubdetectorGeometry::getGeometry(), gId(), mps_fire::i, EEDetId::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), EEDetId::iy(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), hltrates_dqm_sourceclient-live_cfg::offset, point, alignCSCRings::r, EEDetId::validDetId(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, PV3DBase< T, PVType, FrameType >::z(), and EEDetId::zside().

Referenced by CaloGeometryHelper::getClosestCell().

                                                                    {
   try {
     EEDetId mycellID(gId(r.x(), r.y(), r.z()));  // educated guess
 
     if (EEDetId::validDetId(mycellID.ix(), mycellID.iy(), mycellID.zside())) {
       // now get points in convenient ordering
 
       Pt3D A;
       Pt3D B;
       Pt3D C;
       Pt3D point(r.x(), r.y(), r.z());
       // D.K. : equation of plane : AA*x+BB*y+CC*z+DD=0;
       // finding equation for each edge
 
       // ================================================================
       CCGFloat x, y, z;
       unsigned offset = 0;
       int zsign = 1;
       //================================================================
 
       // compute the distance of the point with respect of the 4 crystal lateral planes
 
       if (nullptr != getGeometry(mycellID)) {
         const GlobalPoint& myPosition = getGeometry(mycellID)->getPosition();
 
         x = myPosition.x();
         y = myPosition.y();
         z = myPosition.z();
 
         offset = 0;
         // This will disappear when Andre has applied his fix
 
         if (z > 0) {
           if (x > 0 && y > 0)
             offset = 1;
           else if (x < 0 && y > 0)
             offset = 2;
           else if (x > 0 && y < 0)
             offset = 0;
           else if (x < 0 && y < 0)
             offset = 3;
           zsign = 1;
         } else {
           if (x > 0 && y > 0)
             offset = 3;
           else if (x < 0 && y > 0)
             offset = 2;
           else if (x > 0 && y < 0)
             offset = 0;
           else if (x < 0 && y < 0)
             offset = 1;
           zsign = -1;
         }
         GlobalPoint corners[8];
         for (unsigned ic = 0; ic < 4; ++ic) {
           corners[ic] = getGeometry(mycellID)->getCorners()[(unsigned)((zsign * ic + offset) % 4)];
           corners[4 + ic] = getGeometry(mycellID)->getCorners()[(unsigned)(4 + (zsign * ic + offset) % 4)];
         }
 
         CCGFloat SS[4];
         for (short i = 0; i < 4; ++i) {
           A = Pt3D(corners[i % 4].x(), corners[i % 4].y(), corners[i % 4].z());
           B = Pt3D(corners[(i + 1) % 4].x(), corners[(i + 1) % 4].y(), corners[(i + 1) % 4].z());
           C = Pt3D(corners[4 + (i + 1) % 4].x(), corners[4 + (i + 1) % 4].y(), corners[4 + (i + 1) % 4].z());
           Pl3D plane(A, B, C);
           plane.normalize();
           CCGFloat distance = plane.distance(point);
           if (corners[0].z() < 0.)
             distance = -distance;
           SS[i] = distance;
         }
 
         // Only one move in necessary direction
 
         const bool yout(0 > SS[0] * SS[2]);
         const bool xout(0 > SS[1] * SS[3]);
 
         if (yout || xout) {
           const int ydel(!yout ? 0 : (0 < SS[0] ? -1 : 1));
           const int xdel(!xout ? 0 : (0 < SS[1] ? -1 : 1));
           const unsigned int ix(mycellID.ix() + xdel);
           const unsigned int iy(mycellID.iy() + ydel);
           const unsigned int iz(mycellID.zside());
           if (EEDetId::validDetId(ix, iy, iz))
             mycellID = EEDetId(ix, iy, iz);
         }
 
         return mycellID;
       }
       return DetId(0);
     }
   } catch (cms::Exception& e) {
     return DetId(0);
   }
   return DetId(0);
 }

◆ getGeometryRawPtr()

const CaloCellGeometry * EcalEndcapGeometry::getGeometryRawPtr ( uint32_t index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 450 of file EcalEndcapGeometry.cc.

References m_cellVec, and CaloCellGeometry::param().

                                                                                   {
   // Modify the RawPtr class
   const CaloCellGeometry* cell(&m_cellVec[index]);
   return (m_cellVec.size() < index || nullptr == cell->param() ? nullptr : cell);
 }

◆ getNumberOfCrystalPerModule()

int EcalEndcapGeometry::getNumberOfCrystalPerModule ( ) const

inline

Definition at line 58 of file EcalEndcapGeometry.h.

References _nncrys.

58 { return _nncrys; }

EcalEndcapGeometry::_nncrys

int _nncrys

number of crystals per module

Definition: EcalEndcapGeometry.h:111

◆ getNumberOfModules()

int EcalEndcapGeometry::getNumberOfModules ( ) const

inline

Definition at line 56 of file EcalEndcapGeometry.h.

References _nnmods.

56 { return _nnmods; }

EcalEndcapGeometry::_nnmods

int _nnmods

number of modules

Definition: EcalEndcapGeometry.h:108

◆ gId()

EEDetId EcalEndcapGeometry::gId	(	float	x,
		float	y,
		float	z
	)		const

private

Definition at line 192 of file EcalEndcapGeometry.cc.

References mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), dqmdumpme::k, EEDetId::validDetId(), x, xindex(), y, yindex(), z, zeN, and zeP.

Referenced by getClosestCell().

                                                                {
   const CCGFloat fac(fabs((0 > z ? zeN : zeP) / z));
   const unsigned int ix(xindex(x * fac, z));
   const unsigned int iy(yindex(y * fac, z));
   const unsigned int iz(z > 0 ? 1 : -1);
 
   if (EEDetId::validDetId(ix, iy, iz)) {
     return EEDetId(ix, iy, iz);  // first try is on target
   } else                         // try nearby coordinates, spiraling out from center
   {
     for (unsigned int i(1); i != 6; ++i) {
       for (unsigned int k(0); k != 8; ++k) {
         const int jx(0 == k || 4 == k || 5 == k ? +i : (1 == k || 5 < k ? -i : 0));
         const int jy(2 == k || 4 == k || 6 == k ? +i : (3 == k || 5 == k || 7 == k ? -i : 0));
         if (EEDetId::validDetId(ix + jx, iy + jy, iz)) {
           return EEDetId(ix + jx, iy + jy, iz);
         }
       }
     }
   }
   return EEDetId();  // nowhere near any crystal
 }

◆ hitString()

static std::string EcalEndcapGeometry::hitString ( )

inlinestatic

Definition at line 75 of file EcalEndcapGeometry.h.

75 { return "EcalHitsEE"; }

◆ initializeParms()

void EcalEndcapGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 73 of file EcalEndcapGeometry.cc.

References CaloSubdetectorGeometry::cellGeomPtr(), EEDetId::detIdFromDenseIndex(), mps_fire::i, EEDetId::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), EEDetId::iy(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), m_cellVec, m_del, m_nref, m_wref, m_xhi, m_xlo, m_xoff, m_yhi, m_ylo, m_yoff, AlCaHLTBitMon_ParallelJobs::p, x, y, z, zeN, zeP, and geometryCSVtoXML::zz.

                                          {
   zeP = 0.;
   zeN = 0.;
   unsigned nP = 0;
   unsigned nN = 0;
   m_nref = 0;
 
   for (uint32_t i(0); i != m_cellVec.size(); ++i) {
     auto cell = cellGeomPtr(i);
     if (nullptr != cell) {
       const CCGFloat z(cell->getPosition().z());
       if (z > 0.) {
         zeP += z;
         ++nP;
       } else {
         zeN += z;
         ++nN;
       }
       const EEDetId myId(EEDetId::detIdFromDenseIndex(i));
       const unsigned int ix(myId.ix());
       const unsigned int iy(myId.iy());
       if (ix > m_nref)
         m_nref = ix;
       if (iy > m_nref)
         m_nref = iy;
     }
   }
   if (0 < nP)
     zeP /= (CCGFloat)nP;
   if (0 < nN)
     zeN /= (CCGFloat)nN;
 
   m_xlo[0] = 999;
   m_xhi[0] = -999;
   m_ylo[0] = 999;
   m_yhi[0] = -999;
   m_xlo[1] = 999;
   m_xhi[1] = -999;
   m_ylo[1] = 999;
   m_yhi[1] = -999;
   for (uint32_t i(0); i != m_cellVec.size(); ++i) {
     auto cell = cellGeomPtr(i);
     if (nullptr != cell) {
       const GlobalPoint& p(cell->getPosition());
       const CCGFloat z(p.z());
       const CCGFloat zz(0 > z ? zeN : zeP);
       const CCGFloat x(p.x() * zz / z);
       const CCGFloat y(p.y() * zz / z);
 
       if (0 > z && x < m_xlo[0])
         m_xlo[0] = x;
       if (0 < z && x < m_xlo[1])
         m_xlo[1] = x;
       if (0 > z && y < m_ylo[0])
         m_ylo[0] = y;
       if (0 < z && y < m_ylo[1])
         m_ylo[1] = y;
 
       if (0 > z && x > m_xhi[0])
         m_xhi[0] = x;
       if (0 < z && x > m_xhi[1])
         m_xhi[1] = x;
       if (0 > z && y > m_yhi[0])
         m_yhi[0] = y;
       if (0 < z && y > m_yhi[1])
         m_yhi[1] = y;
     }
   }
 
   m_xoff[0] = (m_xhi[0] + m_xlo[0]) / 2.;
   m_xoff[1] = (m_xhi[1] + m_xlo[1]) / 2.;
   m_yoff[0] = (m_yhi[0] + m_ylo[0]) / 2.;
   m_yoff[1] = (m_yhi[1] + m_ylo[1]) / 2.;
 
   m_del = (m_xhi[0] - m_xlo[0] + m_xhi[1] - m_xlo[1] + m_yhi[0] - m_ylo[0] + m_yhi[1] - m_ylo[1]);
 
   if (1 != m_nref)
     m_wref = m_del / (4. * (m_nref - 1));
 
   m_xlo[0] -= m_wref / 2;
   m_xlo[1] -= m_wref / 2;
   m_xhi[0] += m_wref / 2;
   m_xhi[1] += m_wref / 2;
 
   m_ylo[0] -= m_wref / 2;
   m_ylo[1] -= m_wref / 2;
   m_yhi[0] += m_wref / 2;
   m_yhi[1] += m_wref / 2;
 
   m_del += m_wref;
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("EcalGeom") << "zeP=" << zeP << ", zeN=" << zeN << ", nP=" << nP << ", nN=" << nN;
 
   edm::LogVerbatim("EcalGeom") << "xlo[0]=" << m_xlo[0] << ", xlo[1]=" << m_xlo[1] << ", xhi[0]=" << m_xhi[0]
                                << ", xhi[1]=" << m_xhi[1] << "\nylo[0]=" << m_ylo[0] << ", ylo[1]=" << m_ylo[1]
                                << ", yhi[0]=" << m_yhi[0] << ", yhi[1]=" << m_yhi[1];
 
   edm::LogVerbatim("EcalGeom") << "xoff[0]=" << m_xoff[0] << ", xoff[1]" << m_xoff[1] << ", yoff[0]=" << m_yoff[0]
                                << ", yoff[1]" << m_yoff[1];
 
   edm::LogVerbatim("EcalGeom") << "nref=" << m_nref << ", m_wref=" << m_wref;
 #endif
 }

◆ localCorners()

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

static

Definition at line 424 of file EcalEndcapGeometry.cc.

References TruncatedPyramid::localCorners().

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

◆ myPhi()

static int EcalEndcapGeometry::myPhi ( int i )

inlinestaticprivate

Definition at line 102 of file EcalEndcapGeometry.h.

References mps_fire::i.

Referenced by getClosestBarrelCells().

                           {
     i += 720;
     return (1 + (i - 1) % 360);
   }

◆ newCell()

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

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 428 of file EcalEndcapGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), CaloSubdetectorGeometry::cornersMgr(), hcalRecHitTable_cff::detId, DeadROC_duringRun::f1, DeadROC_duringRun::f2, and m_cellVec.

                                                                                                                    {
   const unsigned int cellIndex(EEDetId(detId).denseIndex());
   m_cellVec[cellIndex] = TruncatedPyramid(cornersMgr(), f1, f2, f3, parm);
   addValidID(detId);
 }

◆ numberOfAlignments()

static unsigned int EcalEndcapGeometry::numberOfAlignments ( )

inlinestatic

Definition at line 79 of file EcalEndcapGeometry.h.

Referenced by FakeCaloAlignmentEP::produceEEAli(), and TestCaloAlignmentEP::produceEEAli().

79 { return 4; }

◆ numberOfParametersPerShape()

unsigned int EcalEndcapGeometry::numberOfParametersPerShape ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 50 of file EcalEndcapGeometry.h.

References k_NumberOfParametersPerShape.

50 { return k_NumberOfParametersPerShape; }

EcalEndcapGeometry::k_NumberOfParametersPerShape

static constexpr int k_NumberOfParametersPerShape

Definition: EcalEndcapGeometry.h:45

◆ numberOfShapes()

unsigned int EcalEndcapGeometry::numberOfShapes ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 49 of file EcalEndcapGeometry.h.

References k_NumberOfShapes.

49 { return k_NumberOfShapes; }

EcalEndcapGeometry::k_NumberOfShapes

static constexpr int k_NumberOfShapes

Definition: EcalEndcapGeometry.h:43

◆ present()

bool EcalEndcapGeometry::present ( const DetId & id ) const

overridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 456 of file EcalEndcapGeometry.cc.

References DetId::Ecal, EcalEndcap, EEDetId::ix(), EEDetId::iy(), EEDetId::validDetId(), and EEDetId::zside().

Referenced by spr::extraIds().

                                                       {
   if (id.det() == DetId::Ecal && id.subdetId() == EcalEndcap) {
     EEDetId eeId(id);
     if (EEDetId::validDetId(eeId.ix(), eeId.iy(), eeId.zside()))
       return true;
   }
   return false;
 }

◆ producerTag()

static std::string EcalEndcapGeometry::producerTag ( )

inlinestatic

Definition at line 77 of file EcalEndcapGeometry.h.

Referenced by PCaloGeometryBuilder::beginRun(), CaloGeometryBuilder::CaloGeometryBuilder(), and PCaloGeometryBuilder::PCaloGeometryBuilder().

77 { return "EcalEndcap"; }

◆ setNumberOfCrystalPerModule()

void EcalEndcapGeometry::setNumberOfCrystalPerModule ( const int nncrys )

inline

Definition at line 62 of file EcalEndcapGeometry.h.

References _nncrys.

62 { _nncrys = nncrys; }

EcalEndcapGeometry::_nncrys

int _nncrys

number of crystals per module

Definition: EcalEndcapGeometry.h:111

◆ setNumberOfModules()

void EcalEndcapGeometry::setNumberOfModules ( const int nnmods )

inline

Definition at line 60 of file EcalEndcapGeometry.h.

References _nnmods.

60 { _nnmods = nnmods; }

EcalEndcapGeometry::_nnmods

int _nnmods

number of modules

Definition: EcalEndcapGeometry.h:108

◆ xindex()

unsigned int EcalEndcapGeometry::xindex	(	CCGFloat	x,
		CCGFloat	z
	)		const

private

Definition at line 178 of file EcalEndcapGeometry.cc.

References mps_fire::i, m_nref, m_wref, m_xlo, x, and z.

Referenced by getCells(), and gId().

                                                                     {
   const CCGFloat xlo(0 > z ? m_xlo[0] : m_xlo[1]);
   const int i(1 + int((x - xlo) / m_wref));
 
   return (1 > i ? 1 : (m_nref < (unsigned int)i ? m_nref : (unsigned int)i));
 }

◆ yindex()

unsigned int EcalEndcapGeometry::yindex	(	CCGFloat	y,
		CCGFloat	z
	)		const

private

Definition at line 185 of file EcalEndcapGeometry.cc.

References mps_fire::i, m_nref, m_wref, m_ylo, y, and z.

Referenced by getCells(), and gId().

                                                                     {
   const CCGFloat ylo(0 > z ? m_ylo[0] : m_ylo[1]);
   const int i(1 + int((y - ylo) / m_wref));
 
   return (1 > i ? 1 : (m_nref < (unsigned int)i ? m_nref : (unsigned int)i));
 }

Member Data Documentation

◆ _nncrys

int EcalEndcapGeometry::_nncrys

private

number of crystals per module

Definition at line 111 of file EcalEndcapGeometry.h.

Referenced by getNumberOfCrystalPerModule(), and setNumberOfCrystalPerModule().

◆ _nnmods

int EcalEndcapGeometry::_nnmods

private

number of modules

Definition at line 108 of file EcalEndcapGeometry.h.

Referenced by getNumberOfModules(), and setNumberOfModules().

◆ k_NumberOfCellsForCorners

constexpr int EcalEndcapGeometry::k_NumberOfCellsForCorners = EEDetId::kSizeForDenseIndexing

static

Definition at line 41 of file EcalEndcapGeometry.h.

◆ k_NumberOfParametersPerShape

constexpr int EcalEndcapGeometry::k_NumberOfParametersPerShape = 11

static

Definition at line 45 of file EcalEndcapGeometry.h.

Referenced by numberOfParametersPerShape().

◆ k_NumberOfShapes

constexpr int EcalEndcapGeometry::k_NumberOfShapes = 1

static

Definition at line 43 of file EcalEndcapGeometry.h.

Referenced by numberOfShapes().

◆ m_avgZ

CCGFloat EcalEndcapGeometry::m_avgZ

mutableprivate

Definition at line 128 of file EcalEndcapGeometry.h.

Referenced by avgAbsZFrontFaceCenter().

◆ m_borderMgr

std::atomic<EZMgrFL<EBDetId>*> EcalEndcapGeometry::m_borderMgr

mutableprivate

Definition at line 124 of file EcalEndcapGeometry.h.

Referenced by getClosestBarrelCells(), and ~EcalEndcapGeometry().

◆ m_borderPtrVec

std::atomic<VecOrdListEBDetIdPtr*> EcalEndcapGeometry::m_borderPtrVec

mutableprivate

Definition at line 126 of file EcalEndcapGeometry.h.

Referenced by getClosestBarrelCells(), and ~EcalEndcapGeometry().

◆ m_cellVec

CellVec EcalEndcapGeometry::m_cellVec

private

Definition at line 131 of file EcalEndcapGeometry.h.

Referenced by avgAbsZFrontFaceCenter(), getCells(), getGeometryRawPtr(), initializeParms(), and newCell().

◆ m_check

std::atomic<bool> EcalEndcapGeometry::m_check

mutableprivate

Definition at line 129 of file EcalEndcapGeometry.h.

Referenced by avgAbsZFrontFaceCenter().

◆ m_del

CCGFloat EcalEndcapGeometry::m_del

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by initializeParms().

◆ m_nref

unsigned int EcalEndcapGeometry::m_nref

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by getCells(), initializeParms(), xindex(), and yindex().

◆ m_wref

CCGFloat EcalEndcapGeometry::m_wref

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by getCells(), initializeParms(), xindex(), and yindex().

◆ m_xhi

CCGFloat EcalEndcapGeometry::m_xhi[2]

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), and initializeParms().

◆ m_xlo

CCGFloat EcalEndcapGeometry::m_xlo[2]

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), initializeParms(), and xindex().

◆ m_xoff

CCGFloat EcalEndcapGeometry::m_xoff[2]

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), and initializeParms().

◆ m_yhi

CCGFloat EcalEndcapGeometry::m_yhi[2]

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), and initializeParms().

◆ m_ylo

CCGFloat EcalEndcapGeometry::m_ylo[2]

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), initializeParms(), and yindex().

◆ m_yoff

CCGFloat EcalEndcapGeometry::m_yoff[2]

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), and initializeParms().

◆ zeN

CCGFloat EcalEndcapGeometry::zeN

private

Definition at line 113 of file EcalEndcapGeometry.h.

Referenced by gId(), and initializeParms().

◆ zeP

CCGFloat EcalEndcapGeometry::zeP

private

Definition at line 113 of file EcalEndcapGeometry.h.

Referenced by getCells(), gId(), and initializeParms().

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ AlignedRecord

◆ AlignmentRecord

◆ CCGFloat

◆ CellVec

◆ DetIdType

◆ IdealRecord

◆ NumberingScheme

◆ OrderedListOfEBDetId

◆ PGeometryRecord

◆ Pt3D

◆ Pt3DVec

◆ VecOrdListEBDetIdPtr

Constructor & Destructor Documentation

◆ EcalEndcapGeometry()

◆ ~EcalEndcapGeometry()

Member Function Documentation

◆ alignmentTransformIndexGlobal()

◆ alignmentTransformIndexLocal()

◆ avgAbsZFrontFaceCenter()

◆ dbString()

◆ detIdFromLocalAlignmentIndex()

◆ getCells()

◆ getClosestBarrelCells()

◆ getClosestCell()

◆ getGeometryRawPtr()

◆ getNumberOfCrystalPerModule()

◆ getNumberOfModules()

◆ gId()

◆ hitString()

◆ initializeParms()

◆ localCorners()

◆ myPhi()

◆ newCell()

◆ numberOfAlignments()

◆ numberOfParametersPerShape()

◆ numberOfShapes()

◆ present()

◆ producerTag()

◆ setNumberOfCrystalPerModule()

◆ setNumberOfModules()

◆ xindex()

◆ yindex()

Member Data Documentation

◆ _nncrys

◆ _nnmods

◆ k_NumberOfCellsForCorners

◆ k_NumberOfParametersPerShape

◆ k_NumberOfShapes

◆ m_avgZ

◆ m_borderMgr

◆ m_borderPtrVec

◆ m_cellVec

◆ m_check

◆ m_del

◆ m_nref

◆ m_wref

◆ m_xhi

◆ m_xlo

◆ m_xoff

◆ m_yhi

◆ m_ylo

◆ m_yoff

◆ zeN

◆ zeP