#include <EcalBarrelGeometry.h>

Inheritance diagram for EcalBarrelGeometry:

Public Types
enum	{ k_NumberOfCellsForCorners = EBDetId::kSizeForDenseIndexing }

enum	{ k_NumberOfShapes = 17 }

enum	{ k_NumberOfParametersPerShape = 11 }

typedef EcalBarrelGeometryRecord	AlignedRecord

typedef EBAlignmentRcd	AlignmentRecord

typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< TruncatedPyramid >	CellVec

typedef EBDetId	DetIdType

typedef IdealGeometryRecord	IdealRecord

typedef EcalBarrelNumberingScheme	NumberingScheme

typedef EZArrayFL< EEDetId >	OrderedListOfEEDetId

typedef PEcalBarrelRcd	PGeometryRecord

typedef CaloCellGeometry::Pt3D	Pt3D

typedef CaloCellGeometry::Pt3DVec	Pt3DVec

typedef std::vector< OrderedListOfEEDetId * >	VecOrdListEEDetIdPtr

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	avgRadiusXYFrontFaceCenter () const

	EcalBarrelGeometry ()

int	getBasketSizeInPhi () const

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

DetId	getClosestCell (const GlobalPoint &r) const override

const OrderedListOfEEDetId *	getClosestEndcapCells (EBDetId id) const

const std::vector< int > &	getEtaBaskets () const

int	getNumXtalsEtaDirection () const

int	getNumXtalsPhiDirection () const

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	setBasketSizeInPhi (const int &PhiBaskets)

void	setEtaBaskets (const std::vector< int > &EtaBaskets)

void	setNumXtalsEtaDirection (const int &nnxtalEta)

void	setNumXtalsPhiDirection (const int &nnxtalPhi)

	~EcalBarrelGeometry () 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 void	initializeParms ()

virtual unsigned int	numberOfTransformParms () const

CaloSubdetectorGeometry &	operator= (const CaloSubdetectorGeometry &)=delete

ParMgr *	parMgr ()

const ParMgr *	parMgrConst () const

ParVecVec &	parVecVec ()

const ParVecVec &	parVecVec () 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 ()

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 Attributes
std::vector< int >	_EtaBaskets

int	_nnxtalEta

int	_nnxtalPhi

int	_PhiBaskets

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

std::atomic< VecOrdListEEDetIdPtr * >	m_borderPtrVec

CellVec	m_cellVec

std::atomic< bool >	m_check

CCGFloat	m_radius

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 19 of file EcalBarrelGeometry.h.

Member Typedef Documentation

◆ AlignedRecord

typedef EcalBarrelGeometryRecord EcalBarrelGeometry::AlignedRecord

Definition at line 28 of file EcalBarrelGeometry.h.

◆ AlignmentRecord

typedef EBAlignmentRcd EcalBarrelGeometry::AlignmentRecord

Definition at line 29 of file EcalBarrelGeometry.h.

◆ CCGFloat

typedef CaloCellGeometry::CCGFloat EcalBarrelGeometry::CCGFloat

Definition at line 23 of file EcalBarrelGeometry.h.

◆ CellVec

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

Definition at line 21 of file EcalBarrelGeometry.h.

◆ DetIdType

typedef EBDetId EcalBarrelGeometry::DetIdType

Definition at line 38 of file EcalBarrelGeometry.h.

◆ IdealRecord

typedef IdealGeometryRecord EcalBarrelGeometry::IdealRecord

Definition at line 27 of file EcalBarrelGeometry.h.

◆ NumberingScheme

typedef EcalBarrelNumberingScheme EcalBarrelGeometry::NumberingScheme

Definition at line 36 of file EcalBarrelGeometry.h.

◆ OrderedListOfEEDetId

typedef EZArrayFL<EEDetId> EcalBarrelGeometry::OrderedListOfEEDetId

Definition at line 32 of file EcalBarrelGeometry.h.

◆ PGeometryRecord

typedef PEcalBarrelRcd EcalBarrelGeometry::PGeometryRecord

Definition at line 30 of file EcalBarrelGeometry.h.

◆ Pt3D

typedef CaloCellGeometry::Pt3D EcalBarrelGeometry::Pt3D

Definition at line 24 of file EcalBarrelGeometry.h.

◆ Pt3DVec

typedef CaloCellGeometry::Pt3DVec EcalBarrelGeometry::Pt3DVec

Definition at line 25 of file EcalBarrelGeometry.h.

◆ VecOrdListEEDetIdPtr

typedef std::vector<OrderedListOfEEDetId*> EcalBarrelGeometry::VecOrdListEEDetIdPtr

Definition at line 34 of file EcalBarrelGeometry.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

Enumerator
k_NumberOfCellsForCorners

Definition at line 40 of file EcalBarrelGeometry.h.

40 { k_NumberOfCellsForCorners = EBDetId::kSizeForDenseIndexing };

◆ anonymous enum

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 42 of file EcalBarrelGeometry.h.

42 { k_NumberOfShapes = 17 };

◆ anonymous enum

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 44 of file EcalBarrelGeometry.h.

44 { k_NumberOfParametersPerShape = 11 };

Constructor & Destructor Documentation

◆ EcalBarrelGeometry()

EcalBarrelGeometry::EcalBarrelGeometry ( )

Definition at line 20 of file EcalBarrelGeometry.cc.

     : _nnxtalEta(85),
       _nnxtalPhi(360),
       _PhiBaskets(18),
       m_borderMgr(nullptr),
       m_borderPtrVec(nullptr),
       m_radius(-1.),
       m_check(false),
       m_cellVec(k_NumberOfCellsForCorners) {
   const int neba[] = {25, 45, 65, 85};
   _EtaBaskets = std::vector<int>(neba, neba + 4);
 }

References _EtaBaskets.

◆ ~EcalBarrelGeometry()

EcalBarrelGeometry::~EcalBarrelGeometry ( )

override

Definition at line 33 of file EcalBarrelGeometry.cc.

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

References m_borderMgr, m_borderPtrVec, and findQualityFiles::v.

Member Function Documentation

◆ alignmentTransformIndexGlobal()

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

static

Definition at line 59 of file EcalBarrelGeometry.cc.

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

References DetId::Ecal.

◆ alignmentTransformIndexLocal()

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

static

Definition at line 45 of file EcalBarrelGeometry.cc.

                                                                              {
   const CaloGenericDetId gid(id);
  
   assert(gid.isEB());
  
   unsigned int index(EBDetId(id).ism() - 1);
  
   return index;
 }

References cms::cuda::assert(), CaloGenericDetId::isEB(), and ecalpyutils::ism().

◆ avgRadiusXYFrontFaceCenter()

CCGFloat EcalBarrelGeometry::avgRadiusXYFrontFaceCenter ( ) const

Definition at line 419 of file EcalBarrelGeometry.cc.

                                                               {
   if (!m_check.load(std::memory_order_acquire)) {
     CCGFloat sum(0);
     for (uint32_t i(0); i != m_cellVec.size(); ++i) {
       auto cell(cellGeomPtr(i));
       if (nullptr != cell) {
         const GlobalPoint& pos(cell->getPosition());
         sum += pos.perp();
       }
     }
     m_radius = sum / m_cellVec.size();
     m_check.store(true, std::memory_order_release);
   }
   return m_radius;
 }

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

◆ dbString()

static std::string EcalBarrelGeometry::dbString ( )

inlinestatic

Definition at line 46 of file EcalBarrelGeometry.h.

46 { return "PEcalBarrelRcd"; }

◆ detIdFromLocalAlignmentIndex()

DetId EcalBarrelGeometry::detIdFromLocalAlignmentIndex ( unsigned int iLoc )

static

Definition at line 55 of file EcalBarrelGeometry.cc.

                                                                         {
   return EBDetId(iLoc + 1, 1, EBDetId::SMCRYSTALMODE);
 }

References EBDetId::SMCRYSTALMODE.

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

◆ getBasketSizeInPhi()

int EcalBarrelGeometry::getBasketSizeInPhi ( ) const

inline

Definition at line 61 of file EcalBarrelGeometry.h.

61 { return _PhiBaskets; }

References _PhiBaskets.

Referenced by ClusterShapeAlgo::Calculate_BarrelBasketEnergyFraction().

◆ getCells()

CaloSubdetectorGeometry::DetIdSet EcalBarrelGeometry::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 266 of file EcalBarrelGeometry.cc.

                                                                                                   {
   constexpr int maxphi(EBDetId::MAX_IPHI);
   constexpr int maxeta(EBDetId::MAX_IETA);
   constexpr float scale(maxphi / (2 * M_PI));  // angle to index
  
   CaloSubdetectorGeometry::DetIdSet dis;  // this is the return object
  
   if (0.000001 < dR) {
     if (dR > M_PI / 2.)  // this version needs "small" dR
     {
       dis = CaloSubdetectorGeometry::getCells(r, dR);  // base class version
     } else {
       const float dR2(dR * dR);
       const float reta(r.eta());
       const float rz(r.z());
       const float rphi(r.phi());
       const float lowEta(reta - dR);
       const float highEta(reta + dR);
  
       if (highEta > -1.5 && lowEta < 1.5)  // in barrel
       {
         const int ieta_center(int(reta * scale + ((rz < 0) ? (-1) : (1))));
         const float phi(rphi < 0 ? rphi + float(2 * M_PI) : rphi);
         const int iphi_center(int(phi * scale + 11.f));  // phi=-9.4deg is iphi=1
  
         const float fr(dR * scale);         // # crystal widths in dR
         const float frp(1.08f * fr + 1.f);  // conservatively above fr
         const float frm(0.92f * fr - 1.f);  // conservatively below fr
         const int idr((int)frp);            // integerize
         const int idr2p((int)(frp * frp));
         const int idr2m(frm > 0 ? int(frm * frm) : 0);
  
         for (int de(-idr); de <= idr; ++de)  // over eta limits
         {
           int ieta(de + ieta_center);
  
           if (std::abs(ieta) <= maxeta && ieta != 0)  // eta is in EB
           {
             const int de2(de * de);
             for (int dp(-idr); dp <= idr; ++dp)  // over phi limits
             {
               const int irange2(dp * dp + de2);
  
               if (irange2 <= idr2p)  // cut off corners that must be too far away
               {
                 const int iphi((iphi_center + dp + maxphi - 1) % maxphi + 1);
  
                 if (iphi != 0) {
                   const EBDetId id(ieta, iphi);
  
                   bool ok(irange2 < idr2m);  // no more calculation necessary if inside this radius
  
                   if (!ok)  // if not ok, then we have to test this cell for being inside cone
                   {
                     const CaloCellGeometry* cell(&m_cellVec[id.denseIndex()]);
                     const float eta(cell->etaPos());
                     const float phi(cell->phiPos());
                     ok = (reco::deltaR2(eta, phi, reta, rphi) < dR2);
                   }
                   if (ok)
                     dis.insert(id);
                 }
               }
             }
           }
         }
       }
     }
   }
   return dis;
 }

References funct::abs(), reco::deltaR2(), Calorimetry_cff::dp, HGC3DClusterGenMatchSelector_cfi::dR, PVValHelper::eta, CaloCellGeometry::etaPos(), f, CaloSubdetectorGeometry::getCells(), triggerObjects_cff::id, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, m_cellVec, M_PI, EBDetId::MAX_IETA, EBDetId::MAX_IPHI, convertSQLiteXML::ok, phi, CaloCellGeometry::phiPos(), alignCSCRings::r, and Scenarios_cff::scale.

◆ getClosestCell()

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 63 of file EcalBarrelGeometry.cc.

                                                                    {
   // z is the easy one
   int leverx = 1;
   int levery = 1;
   CCGFloat pointz = r.z();
   int zbin = 1;
   if (pointz < 0)
     zbin = -1;
  
   // Now find the closest eta
   CCGFloat pointeta = r.eta();
   //  double eta;
   CCGFloat deta = 999.;
   int etabin = 1;
  
   int guessed_eta = (int)(fabs(pointeta) / 0.0174) + 1;
   int guessed_eta_begin = guessed_eta - 1;
   int guessed_eta_end = guessed_eta + 1;
   if (guessed_eta_begin < 1)
     guessed_eta_begin = 1;
   if (guessed_eta_end > 85)
     guessed_eta_end = 85;
  
   for (int bin = guessed_eta_begin; bin <= guessed_eta_end; bin++) {
     try {
       if (!present(EBDetId(zbin * bin, 1, EBDetId::ETAPHIMODE)))
         continue;
  
       CCGFloat eta = getGeometry(EBDetId(zbin * bin, 1, EBDetId::ETAPHIMODE))->etaPos();
  
       if (fabs(pointeta - eta) < deta) {
         deta = fabs(pointeta - eta);
         etabin = bin;
       } else
         break;
     } catch (cms::Exception& e) {
     }
   }
  
   // Now the closest phi. always same number of phi bins(!?)
   constexpr CCGFloat twopi = M_PI + M_PI;
   // 10 degree tilt
   constexpr CCGFloat tilt = twopi / 36.;
  
   CCGFloat pointphi = r.phi() + tilt;
  
   // put phi in correct range (0->2pi)
   if (pointphi > twopi)
     pointphi -= twopi;
   if (pointphi < 0)
     pointphi += twopi;
  
   //calculate phi bin, distinguish + and - eta
   int phibin = static_cast<int>(pointphi / (twopi / _nnxtalPhi)) + 1;
   //   if(point.z()<0.0)
   //     {
   //       phibin = nxtalPhi/2 - 1 - phibin;
   //       if(phibin<0)
   //         phibin += nxtalPhi;
   //     }
   try {
     EBDetId myCell(zbin * etabin, phibin, EBDetId::ETAPHIMODE);
  
     if (!present(myCell))
       return DetId(0);
  
     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
  
     // Since the point can lie between crystals, it is necessary to keep track of the movements
     // to avoid infinite loops
     CCGFloat history[4]{0.f};
  
     //
     // stop movement in eta direction when closest cell was found (point between crystals)
     int start = 1;
     int counter = 0;
     // Moving until find closest crystal in eta and phi directions (leverx and levery)
     while (leverx == 1 || levery == 1) {
       leverx = 0;
       levery = 0;
       const CaloCellGeometry::CornersVec& corners(getGeometry(myCell)->getCorners());
       CCGFloat SS[4];
  
       // compute the distance of the point with respect of the 4 crystal lateral planes
       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 (plane.d() > 0.)
           distance = -distance;
         if (corners[0].z() < 0.)
           distance = -distance;
         SS[i] = distance;
       }
  
       // SS's - normals
       // check position of the point with respect to opposite side of crystal
       // if SS's have opposite sign, the  point lies inside that crystal
  
       if ((SS[0] > 0. && SS[2] > 0.) || (SS[0] < 0. && SS[2] < 0.)) {
         levery = 1;
         if (history[0] > 0. && history[2] > 0. && SS[0] < 0 && SS[2] < 0 &&
             (fabs(SS[0]) + fabs(SS[2])) > (fabs(history[0]) + fabs(history[2])))
           levery = 0;
         if (history[0] < 0. && history[2] < 0. && SS[0] > 0 && SS[2] > 0 &&
             (fabs(SS[0]) + fabs(SS[2])) > (fabs(history[0]) + fabs(history[2])))
           levery = 0;
  
         if (SS[0] > 0.) {
           EBDetId nextPoint;
           if (myCell.iphi() == EBDetId::MIN_IPHI)
             nextPoint = EBDetId(myCell.ieta(), EBDetId::MAX_IPHI);
           else
             nextPoint = EBDetId(myCell.ieta(), myCell.iphi() - 1);
           if (present(nextPoint))
             myCell = nextPoint;
           else
             levery = 0;
         } else {
           EBDetId nextPoint;
           if (myCell.iphi() == EBDetId::MAX_IPHI)
             nextPoint = EBDetId(myCell.ieta(), EBDetId::MIN_IPHI);
           else
             nextPoint = EBDetId(myCell.ieta(), myCell.iphi() + 1);
           if (present(nextPoint))
             myCell = nextPoint;
           else
             levery = 0;
         }
       }
  
       if (((SS[1] > 0. && SS[3] > 0.) || (SS[1] < 0. && SS[3] < 0.)) && start == 1) {
         leverx = 1;
  
         if (history[1] > 0. && history[3] > 0. && SS[1] < 0 && SS[3] < 0 &&
             (fabs(SS[1]) + fabs(SS[3])) > (fabs(history[1]) + fabs(history[3]))) {
           leverx = 0;
           start = 0;
         }
  
         if (history[1] < 0. && history[3] < 0. && SS[1] > 0 && SS[3] > 0 &&
             (fabs(SS[1]) + fabs(SS[3])) > (fabs(history[1]) + fabs(history[3]))) {
           leverx = 0;
           start = 0;
         }
  
         if (SS[1] > 0.) {
           EBDetId nextPoint;
           if (myCell.ieta() == -1)
             nextPoint = EBDetId(1, myCell.iphi());
           else {
             int nieta = myCell.ieta() + 1;
             if (nieta == 86)
               nieta = 85;
             nextPoint = EBDetId(nieta, myCell.iphi());
           }
           if (present(nextPoint))
             myCell = nextPoint;
           else
             leverx = 0;
         } else {
           EBDetId nextPoint;
           if (myCell.ieta() == 1)
             nextPoint = EBDetId(-1, myCell.iphi());
           else {
             int nieta = myCell.ieta() - 1;
             if (nieta == -86)
               nieta = -85;
             nextPoint = EBDetId(nieta, myCell.iphi());
           }
           if (present(nextPoint))
             myCell = nextPoint;
           else
             leverx = 0;
         }
       }
  
       // Update the history. If the point lies between crystals, the closest one
       // is returned
       std::copy(SS, SS + 4, history);
  
       counter++;
       if (counter == 10) {
         leverx = 0;
         levery = 0;
       }
     }
     // D.K. if point lies netween cells, take a closest cell.
     return DetId(myCell);
   } catch (cms::Exception& e) {
     return DetId(0);
   }
 }

References _nnxtalPhi, MaterialEffects_cfi::A, TtFullHadDaughter::B, newFWLiteAna::bin, gen::C, filterCSVwithJSON::copy, HLT_2018_cff::distance, MillePedeFileConverter_cfg::e, PVValHelper::eta, EBDetId::ETAPHIMODE, CaloSubdetectorGeometry::getGeometry(), mps_fire::i, EBDetId::ieta(), createfilelist::int, EBDetId::iphi(), M_PI, EBDetId::MAX_IPHI, EBDetId::MIN_IPHI, point, present(), alignCSCRings::r, x, y, and z.

Referenced by CaloGeometryHelper::getClosestCell().

◆ getClosestEndcapCells()

const EcalBarrelGeometry::OrderedListOfEEDetId * EcalBarrelGeometry::getClosestEndcapCells ( EBDetId id ) const

Definition at line 338 of file EcalBarrelGeometry.cc.

                                                                                                         {
   OrderedListOfEEDetId* ptr(nullptr);
   auto ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
   if (!ptrVec) {
     if (0 != id.rawId()) {
       const int iPhi(id.iphi());
       const int iz(id.ieta() > 0 ? 1 : -1);
       const EEDetId eeid(EEDetId::idOuterRing(iPhi, iz));
       const int iq(eeid.iquadrant());
       const int xout(1 == iq || 4 == iq ? 1 : -1);
       const int yout(1 == iq || 2 == iq ? 1 : -1);
       if (!m_borderMgr.load(std::memory_order_acquire)) {
         EZMgrFL<EEDetId>* expect = nullptr;
         auto ptrMgr = new EZMgrFL<EEDetId>(720 * 9, 9);
         bool exchanged = m_borderMgr.compare_exchange_strong(expect, ptrMgr, std::memory_order_acq_rel);
         if (!exchanged)
           delete ptrMgr;
       }
       VecOrdListEEDetIdPtr* expect = nullptr;
       auto ptrVec = new VecOrdListEEDetIdPtr();
       ptrVec->reserve(720);
       for (unsigned int i(0); i != 720; ++i) {
         const int kz(360 > i ? -1 : 1);
         const EEDetId eeid(EEDetId::idOuterRing(i % 360 + 1, kz));
  
         const int jx(eeid.ix());
         const int jy(eeid.iy());
  
         OrderedListOfEEDetId& olist(*new OrderedListOfEEDetId(m_borderMgr.load(std::memory_order_acquire)));
         int il(0);
  
         for (unsigned int k(1); k <= 25; ++k) {
           const int kx(1 == k || 2 == k || 3 == k || 12 == k || 13 == k
                            ? 0
                            : (4 == k || 6 == k || 8 == k || 15 == k || 20 == k
                                   ? 1
                                   : (5 == k || 7 == k || 9 == k || 16 == k || 19 == k
                                          ? -1
                                          : (10 == k || 14 == k || 21 == k || 22 == k || 25 == k ? 2 : -2))));
           const int ky(1 == k || 4 == k || 5 == k || 10 == k || 11 == k
                            ? 0
                            : (2 == k || 6 == k || 7 == k || 14 == k || 17 == k
                                   ? 1
                                   : (3 == k || 8 == k || 9 == k || 18 == k || 21 == k
                                          ? -1
                                          : (12 == k || 15 == k || 16 == k || 22 == k || 23 == k ? 2 : -2))));
  
           if (8 >= il && EEDetId::validDetId(jx + kx * xout, jy + ky * yout, kz)) {
             olist[il++] = EEDetId(jx + kx * xout, jy + ky * yout, kz);
           }
         }
         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;
 }

References mps_fire::i, EEDetId::idOuterRing(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, EEDetId::iquadrant(), EEDetId::ix(), EEDetId::iy(), dqmdumpme::k, m_borderMgr, m_borderPtrVec, and EEDetId::validDetId().

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

◆ getEtaBaskets()

const std::vector<int>& EcalBarrelGeometry::getEtaBaskets ( ) const

inline

Definition at line 59 of file EcalBarrelGeometry.h.

59 { return _EtaBaskets; }

References _EtaBaskets.

Referenced by ClusterShapeAlgo::Calculate_BarrelBasketEnergyFraction().

◆ getGeometryRawPtr()

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

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 435 of file EcalBarrelGeometry.cc.

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

References m_cellVec, and CaloCellGeometry::param().

◆ getNumXtalsEtaDirection()

int EcalBarrelGeometry::getNumXtalsEtaDirection ( ) const

inline

Definition at line 57 of file EcalBarrelGeometry.h.

57 { return _nnxtalEta; }

References _nnxtalEta.

◆ getNumXtalsPhiDirection()

int EcalBarrelGeometry::getNumXtalsPhiDirection ( ) const

inline

Definition at line 55 of file EcalBarrelGeometry.h.

55 { return _nnxtalPhi; }

References _nnxtalPhi.

◆ hitString()

static std::string EcalBarrelGeometry::hitString ( )

inlinestatic

Definition at line 81 of file EcalBarrelGeometry.h.

81 { return "EcalHitsEB"; }

◆ localCorners()

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

static

Definition at line 401 of file EcalBarrelGeometry.cc.

                                                                                                 {
   const bool negz(EBDetId::kSizeForDenseIndexing / 2 > i);
   const bool odd(1 == i % 2);
  
   if (((negz && !odd) || (!negz && odd))) {
     TruncatedPyramid::localCornersReflection(lc, pv, ref);
   } else {
     TruncatedPyramid::localCornersSwap(lc, pv, ref);
   }
 }

References mps_fire::i, EBDetId::kSizeForDenseIndexing, TruncatedPyramid::localCornersReflection(), TruncatedPyramid::localCornersSwap(), and MetAnalyzer::pv().

◆ newCell()

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

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 412 of file EcalBarrelGeometry.cc.

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

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

◆ numberOfAlignments()

static unsigned int EcalBarrelGeometry::numberOfAlignments ( )

inlinestatic

Definition at line 85 of file EcalBarrelGeometry.h.

85 { return 36; }

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

◆ numberOfParametersPerShape()

unsigned int EcalBarrelGeometry::numberOfParametersPerShape ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 49 of file EcalBarrelGeometry.h.

49 { return k_NumberOfParametersPerShape; }

References k_NumberOfParametersPerShape.

◆ numberOfShapes()

unsigned int EcalBarrelGeometry::numberOfShapes ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 48 of file EcalBarrelGeometry.h.

48 { return k_NumberOfShapes; }

References k_NumberOfShapes.

◆ present()

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

overridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 441 of file EcalBarrelGeometry.cc.

                                                       {
   if (id.det() == DetId::Ecal && id.subdetId() == EcalBarrel) {
     EBDetId ebId(id);
     if (EBDetId::validDetId(ebId.ieta(), ebId.iphi()))
       return true;
   }
   return false;
 }

References DetId::Ecal, EcalBarrel, EBDetId::ieta(), EBDetId::iphi(), and EBDetId::validDetId().

Referenced by spr::extraIds(), and getClosestCell().

◆ producerTag()

static std::string EcalBarrelGeometry::producerTag ( )

inlinestatic

Definition at line 83 of file EcalBarrelGeometry.h.

83 { return "EcalBarrel"; }

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

◆ setBasketSizeInPhi()

void EcalBarrelGeometry::setBasketSizeInPhi ( const int & PhiBaskets )

inline

Definition at line 69 of file EcalBarrelGeometry.h.

69 { _PhiBaskets = PhiBaskets; }

References _PhiBaskets.

◆ setEtaBaskets()

void EcalBarrelGeometry::setEtaBaskets ( const std::vector< int > & EtaBaskets )

inline

Definition at line 67 of file EcalBarrelGeometry.h.

67 { _EtaBaskets = EtaBaskets; }

References _EtaBaskets.

◆ setNumXtalsEtaDirection()

void EcalBarrelGeometry::setNumXtalsEtaDirection ( const int & nnxtalEta )

inline

Definition at line 65 of file EcalBarrelGeometry.h.

65 { _nnxtalEta = nnxtalEta; }

References _nnxtalEta.

◆ setNumXtalsPhiDirection()

void EcalBarrelGeometry::setNumXtalsPhiDirection ( const int & nnxtalPhi )

inline

Definition at line 63 of file EcalBarrelGeometry.h.

63 { _nnxtalPhi = nnxtalPhi; }

References _nnxtalPhi.

Member Data Documentation

◆ _EtaBaskets

std::vector<int> EcalBarrelGeometry::_EtaBaskets

private

size of the baskets in the eta direction. This is needed to find out whether two adjacent crystals lie in the same basked ('module') or not (e.g. this can be used for correcting cluster energies etc.)

Definition at line 118 of file EcalBarrelGeometry.h.

Referenced by EcalBarrelGeometry(), getEtaBaskets(), and setEtaBaskets().

◆ _nnxtalEta

int EcalBarrelGeometry::_nnxtalEta

private

number of crystals in eta direction

Definition at line 109 of file EcalBarrelGeometry.h.

Referenced by getNumXtalsEtaDirection(), and setNumXtalsEtaDirection().

◆ _nnxtalPhi

int EcalBarrelGeometry::_nnxtalPhi

private

number of crystals in phi direction

Definition at line 112 of file EcalBarrelGeometry.h.

Referenced by getClosestCell(), getNumXtalsPhiDirection(), and setNumXtalsPhiDirection().

◆ _PhiBaskets

int EcalBarrelGeometry::_PhiBaskets

private

size of one basket in phi

Definition at line 121 of file EcalBarrelGeometry.h.

Referenced by getBasketSizeInPhi(), and setBasketSizeInPhi().

◆ m_borderMgr

std::atomic<EZMgrFL<EEDetId>*> EcalBarrelGeometry::m_borderMgr

mutableprivate

Definition at line 123 of file EcalBarrelGeometry.h.

Referenced by getClosestEndcapCells(), and ~EcalBarrelGeometry().

◆ m_borderPtrVec

std::atomic<VecOrdListEEDetIdPtr*> EcalBarrelGeometry::m_borderPtrVec

mutableprivate

Definition at line 125 of file EcalBarrelGeometry.h.

Referenced by getClosestEndcapCells(), and ~EcalBarrelGeometry().

◆ m_cellVec

CellVec EcalBarrelGeometry::m_cellVec

private

Definition at line 129 of file EcalBarrelGeometry.h.

Referenced by avgRadiusXYFrontFaceCenter(), getCells(), getGeometryRawPtr(), and newCell().

◆ m_check

std::atomic<bool> EcalBarrelGeometry::m_check

mutableprivate

Definition at line 127 of file EcalBarrelGeometry.h.

Referenced by avgRadiusXYFrontFaceCenter().

◆ m_radius

CCGFloat EcalBarrelGeometry::m_radius

mutableprivate

Definition at line 126 of file EcalBarrelGeometry.h.

Referenced by avgRadiusXYFrontFaceCenter().

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ AlignedRecord

◆ AlignmentRecord

◆ CCGFloat

◆ CellVec

◆ DetIdType

◆ IdealRecord

◆ NumberingScheme

◆ OrderedListOfEEDetId

◆ PGeometryRecord

◆ Pt3D

◆ Pt3DVec

◆ VecOrdListEEDetIdPtr

Member Enumeration Documentation

◆ anonymous enum

◆ anonymous enum

◆ anonymous enum

Constructor & Destructor Documentation

◆ EcalBarrelGeometry()

◆ ~EcalBarrelGeometry()

Member Function Documentation

◆ alignmentTransformIndexGlobal()

◆ alignmentTransformIndexLocal()

◆ avgRadiusXYFrontFaceCenter()

◆ dbString()

◆ detIdFromLocalAlignmentIndex()

◆ getBasketSizeInPhi()

◆ getCells()

◆ getClosestCell()

◆ getClosestEndcapCells()

◆ getEtaBaskets()

◆ getGeometryRawPtr()

◆ getNumXtalsEtaDirection()

◆ getNumXtalsPhiDirection()

◆ hitString()

◆ localCorners()

◆ newCell()

◆ numberOfAlignments()

◆ numberOfParametersPerShape()

◆ numberOfShapes()

◆ present()

◆ producerTag()

◆ setBasketSizeInPhi()

◆ setEtaBaskets()

◆ setNumXtalsEtaDirection()

◆ setNumXtalsPhiDirection()

Member Data Documentation

◆ _EtaBaskets

◆ _nnxtalEta

◆ _nnxtalPhi

◆ _PhiBaskets

◆ m_borderMgr

◆ m_borderPtrVec

◆ m_cellVec

◆ m_check

◆ m_radius