EcalPreshowerGeometry Class Reference

#include <EcalPreshowerGeometry.h>

Inheritance diagram for EcalPreshowerGeometry:

Public Types
enum	{ k_NumberOfCellsForCorners = ESDetId::kSizeForDenseIndexing }
enum	{ k_NumberOfShapes = 4 }
enum	{ k_NumberOfParametersPerShape = 4 }
typedef EcalPreshowerGeometryRecord	AlignedRecord
typedef ESAlignmentRcd	AlignmentRecord
typedef CaloCellGeometry::CCGFloat	CCGFloat
typedef std::vector < PreshowerStrip >	CellVec
typedef ESDetId	DetIdType
typedef IdealGeometryRecord	IdealRecord
typedef EcalPreshowerNumberingScheme	NumberingScheme
typedef CaloSubdetectorGeometry::ParVec	ParVec
typedef CaloSubdetectorGeometry::ParVecVec	ParVecVec
typedef PEcalPreshowerRcd	PGeometryRecord
typedef CaloCellGeometry::Pt3D	Pt3D
typedef CaloCellGeometry::Pt3DVec	Pt3DVec
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
	EcalPreshowerGeometry ()
DetId	getClosestCell (const GlobalPoint &r) const override
virtual DetId	getClosestCellInPlane (const GlobalPoint &r, int plane) 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
unsigned int	numberOfTransformParms () const override
bool	present (const DetId &id) const override
	is this detid present in the geometry? More...
void	setzPlanes (CCGFloat z1minus, CCGFloat z2minus, CCGFloat z1plus, CCGFloat z2plus)
	~EcalPreshowerGeometry () override
	The EcalPreshowerGeometry will delete all its cell geometries at destruction time. More...
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 DetIdSet	getCells (const GlobalPoint &r, double dR) const
	Get a list of all cells within a dR of the given cell. More...
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...
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
CellVec	m_cellVec
const CCGFloat	m_xInterLadGap
const CCGFloat	m_xIntraLadGap
const CCGFloat	m_xWidWaf
const CCGFloat	m_yCtrOff
const CCGFloat	m_yWidAct
CCGFloat	m_zplane [4]

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 15 of file EcalPreshowerGeometry.h.

Member Typedef Documentation

typedef EcalPreshowerGeometryRecord EcalPreshowerGeometry::AlignedRecord

Definition at line 24 of file EcalPreshowerGeometry.h.

typedef ESAlignmentRcd EcalPreshowerGeometry::AlignmentRecord

Definition at line 25 of file EcalPreshowerGeometry.h.

typedef CaloCellGeometry::CCGFloat EcalPreshowerGeometry::CCGFloat

Definition at line 19 of file EcalPreshowerGeometry.h.

typedef std::vector<PreshowerStrip> EcalPreshowerGeometry::CellVec

Definition at line 17 of file EcalPreshowerGeometry.h.

typedef ESDetId EcalPreshowerGeometry::DetIdType

Definition at line 31 of file EcalPreshowerGeometry.h.

typedef IdealGeometryRecord EcalPreshowerGeometry::IdealRecord

Definition at line 23 of file EcalPreshowerGeometry.h.

typedef EcalPreshowerNumberingScheme EcalPreshowerGeometry::NumberingScheme

Definition at line 28 of file EcalPreshowerGeometry.h.

typedef CaloSubdetectorGeometry::ParVec EcalPreshowerGeometry::ParVec

Definition at line 29 of file EcalPreshowerGeometry.h.

typedef CaloSubdetectorGeometry::ParVecVec EcalPreshowerGeometry::ParVecVec

Definition at line 30 of file EcalPreshowerGeometry.h.

typedef PEcalPreshowerRcd EcalPreshowerGeometry::PGeometryRecord

Definition at line 26 of file EcalPreshowerGeometry.h.

typedef CaloCellGeometry::Pt3D EcalPreshowerGeometry::Pt3D

Definition at line 20 of file EcalPreshowerGeometry.h.

typedef CaloCellGeometry::Pt3DVec EcalPreshowerGeometry::Pt3DVec

Definition at line 21 of file EcalPreshowerGeometry.h.

Member Enumeration Documentation

anonymous enum

Enumerator
k_NumberOfCellsForCorners

Definition at line 33 of file EcalPreshowerGeometry.h.

{ k_NumberOfCellsForCorners = ESDetId::kSizeForDenseIndexing };

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 35 of file EcalPreshowerGeometry.h.

{ k_NumberOfShapes = 4 };

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 37 of file EcalPreshowerGeometry.h.

{ k_NumberOfParametersPerShape = 4 };

Constructor & Destructor Documentation

EcalPreshowerGeometry::EcalPreshowerGeometry

(

)

Definition at line 16 of file EcalPreshowerGeometry.cc.

References m_zplane.

    : m_xWidWaf(6.3),
      m_xInterLadGap(0.05),  // additional gap between wafers in adj ladders
      m_xIntraLadGap(0.04),  // gap between wafers in same ladder
      m_yWidAct(6.1),
      m_yCtrOff(0.05),  // gap at center
      m_cellVec(k_NumberOfCellsForCorners) {
  m_zplane[0] = 0.;
  m_zplane[1] = 0.;
  m_zplane[2] = 0.;
  m_zplane[3] = 0.;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("EcalGeom") << "EcalPreshowerGeometry::Creating an instance";
#endif
}

EcalPreshowerGeometry::~EcalPreshowerGeometry ( )

override

The EcalPreshowerGeometry will delete all its cell geometries at destruction time.

Definition at line 32 of file EcalPreshowerGeometry.cc.

{}

Member Function Documentation

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

static

Definition at line 75 of file EcalPreshowerGeometry.cc.

References DetId::Ecal.

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

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

static

Definition at line 34 of file EcalPreshowerGeometry.cc.

References cms::cuda::assert(), CaloGenericDetId::isES(), and edm::second().

                                                                                {
  const CaloGenericDetId gid(id);

  assert(gid.isES());

  // plane 2 is split into 2 dees along x=0 for both signs of z

  // plane 1 at zsign=-1 is split into 2 dees between six=19 and six=20 for siy<=20,
  //                                             and six=21 and 22 for siy>=21

  // plane 1 at zsign=+1 is split into 2 dees between six=20 and six=21 for siy<=20,
  //                                             and six=19 and 20 for siy>=21

  // Desired numbering
  //                LEFT    RIGHT (as one faces the Dee from the IP)
  //  ES-  pl=2     0       1
  //       pl=1     2       3    the reversal of pl=2 and pl=1 is intentional here (CM Kuo)
  //  ES+  pl=1     4       5
  //       pl=2     6       7

  const ESDetId esid(id);
  const int jx(esid.six() - 1);
  const int jy(esid.siy() - 1);
  const int jz(esid.zside() + 1);
  const int pl(esid.plane() - 1);
  const bool second(1 == pl);
  const bool top(19 < jy);
  const bool negz(0 == jz);
  const int lrl(19 > jx ? 0 : 1);
  const int lrr(21 > jx ? 0 : 1);

  return (second ? jx / 20 + 3 * jz :        // 2nd plane split along middle
              (negz && !top ? lrl + 2 :      // 1st plane at neg z and bottom half split at six=19&20
                   (negz && top ? lrr + 2 :  // 1st plane at neg z and top half split at six=21&22
                        (!negz && !top ? lrr + 4 : lrl + 4))));  // opposite at positive z
}

static std::string EcalPreshowerGeometry::dbString ( )

inlinestatic

Definition at line 39 of file EcalPreshowerGeometry.h.

{ return "PEcalPreshowerRcd"; }

DetId EcalPreshowerGeometry::detIdFromLocalAlignmentIndex ( unsigned int  iLoc )

static

Definition at line 71 of file EcalPreshowerGeometry.cc.

Referenced by WriteESAlignments::convert(), FakeCaloAlignmentEP::produceESAli(), and TestCaloAlignmentEP::produceESAli().

                                                                           {
  return ESDetId(1, 10 + 20 * (iLoc % 2), 10, 2 > iLoc || 5 < iLoc ? 2 : 1, 2 * (iLoc / 4) - 1);
}

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 137 of file EcalPreshowerGeometry.cc.

References getClosestCellInPlane().

{ return getClosestCellInPlane(point, 2); }

DetId EcalPreshowerGeometry::getClosestCellInPlane ( const GlobalPoint &  r, int  plane  ) const

virtual

Definition at line 139 of file EcalPreshowerGeometry.cc.

References cuy::col, CaloSubdetectorGeometry::getGeometry(), ESDetId::ISTRIP_MAX, ESDetId::ISTRIP_MIN, m_xInterLadGap, m_xIntraLadGap, m_xWidWaf, m_yCtrOff, m_yWidAct, m_zplane, AlCaHLTBitMon_ParallelJobs::p, present(), ESDetId::validDetId(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by PreshowerHitMaker::addHit(), getClosestCell(), and ReducedESRecHitCollectionProducer::produce().

                                                                                            {
  const CCGFloat x(point.x());
  const CCGFloat y(point.y());
  const CCGFloat z(point.z());

  if (0 == z || 1 > plane || 2 < plane)
    return DetId(0);

  const unsigned int iz((0 > z ? 0 : 2) + plane - 1);

  const CCGFloat ze(m_zplane[iz]);
  const CCGFloat xe(x * ze / z);
  const CCGFloat ye(y * ze / z);

  const CCGFloat x0(1 == plane ? xe : ye);
  const CCGFloat y0(1 == plane ? ye : xe);

  static const CCGFloat xWid(m_xWidWaf + m_xIntraLadGap + m_xInterLadGap);

  const int row(1 + int(y0 + 20. * m_yWidAct - m_yCtrOff) / m_yWidAct);
  const int col(1 + int((x0 + 20. * xWid) / xWid));

  CCGFloat closest(1e9);

  DetId detId(0);

  const int jz(0 > ze ? -1 : 1);

  //   std::cout<<"** p="<<point<<", ("<<xe<<", "<<ye<<", "<<ze<<"), row="<<row<<", col="<<col<<std::endl;

  for (int ix(-1); ix != 2; ++ix)  // search within +-1 in row and col
  {
    for (int iy(-1); iy != 2; ++iy) {
      for (int jstrip(ESDetId::ISTRIP_MIN); jstrip <= ESDetId::ISTRIP_MAX; ++jstrip) {
        const int jx(1 == plane ? col + ix : row + iy);
        const int jy(1 == plane ? row + iy : col + ix);
        if (ESDetId::validDetId(jstrip, jx, jy, plane, jz)) {
          const ESDetId esId(jstrip, jx, jy, plane, jz);
          auto cell = getGeometry(esId);
          if (nullptr != cell) {
            const GlobalPoint& p(cell->getPosition());
            const CCGFloat dist2((p.x() - xe) * (p.x() - xe) + (p.y() - ye) * (p.y() - ye));
            if (dist2 < closest && present(esId)) {
              closest = dist2;
              detId = esId;
            }
          }
        }
      }
    }
  }
  return detId;
}

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

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 210 of file EcalPreshowerGeometry.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);
}

static std::string EcalPreshowerGeometry::hitString ( )

inlinestatic

Definition at line 60 of file EcalPreshowerGeometry.h.

{ return "EcalHitsES"; }

void EcalPreshowerGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 79 of file EcalPreshowerGeometry.cc.

References cms::cuda::assert(), CaloSubdetectorGeometry::getGeometry(), CaloSubdetectorGeometry::getValidDetIds(), mps_fire::i, ESDetId::plane(), setzPlanes(), and ESDetId::zside().

                                            {
  unsigned int n1minus(0);
  unsigned int n2minus(0);
  unsigned int n1plus(0);
  unsigned int n2plus(0);
  CCGFloat z1minus(0);
  CCGFloat z2minus(0);
  CCGFloat z1plus(0);
  CCGFloat z2plus(0);
  const std::vector<DetId>& esDetIds(getValidDetIds());
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("EcalGeom") << "EcalPreshowerGeometry:: Get " << esDetIds.size() << " valid DetIds";
#endif

  for (unsigned int i(0); i != esDetIds.size(); ++i) {
    const ESDetId esid(esDetIds[i]);
    auto cell = getGeometry(esid);
    if (nullptr != cell) {
      const CCGFloat zz(cell->getPosition().z());
      if (1 == esid.plane()) {
        if (0 > esid.zside()) {
          z1minus += zz;
          ++n1minus;
        } else {
          z1plus += zz;
          ++n1plus;
        }
      }
      if (2 == esid.plane()) {
        if (0 > esid.zside()) {
          z2minus += zz;
          ++n2minus;
        } else {
          z2plus += zz;
          ++n2plus;
        }
      }
    }
  }
  assert(0 != n1minus && 0 != n2minus && 0 != n1plus && 0 != n2plus);
  z1minus /= (1. * n1minus);
  z2minus /= (1. * n2minus);
  z1plus /= (1. * n1plus);
  z2plus /= (1. * n2plus);
  assert(0 != z1minus && 0 != z2minus && 0 != z1plus && 0 != z2plus);
  setzPlanes(z1minus, z2minus, z1plus, z2plus);
}

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

static

Definition at line 193 of file EcalPreshowerGeometry.cc.

References PreshowerStrip::localCorners().

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

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

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 197 of file EcalPreshowerGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), CaloSubdetectorGeometry::cornersMgr(), and m_cellVec.

                                                        {
  const unsigned int cellIndex(ESDetId(detId).denseIndex());
  m_cellVec[cellIndex] = PreshowerStrip(f1, cornersMgr(), parm);
  addValidID(detId);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("EcalGeom") << "EcalPreshowerGeometry::Add a new DetId " << ESDetId(detId);
#endif
}

static unsigned int EcalPreshowerGeometry::numberOfAlignments ( )

inlinestatic

Definition at line 64 of file EcalPreshowerGeometry.h.

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

{ return 8; }

unsigned int EcalPreshowerGeometry::numberOfParametersPerShape ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 42 of file EcalPreshowerGeometry.h.

References k_NumberOfParametersPerShape.

{ return k_NumberOfParametersPerShape; }

unsigned int EcalPreshowerGeometry::numberOfShapes ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 41 of file EcalPreshowerGeometry.h.

References k_NumberOfShapes.

{ return k_NumberOfShapes; }

unsigned int EcalPreshowerGeometry::numberOfTransformParms ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 58 of file EcalPreshowerGeometry.h.

{ return 3; }

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

inlineoverridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 81 of file EcalPreshowerGeometry.h.

References CaloGenericDetId::denseIndex(), and m_cellVec.

Referenced by getClosestCellInPlane().

                                               {
    if (id == DetId(0))
      return false;
    // not needed???
    auto index = CaloGenericDetId(id).denseIndex();
    return index < m_cellVec.size();
  }

static std::string EcalPreshowerGeometry::producerTag ( )

inlinestatic

Definition at line 62 of file EcalPreshowerGeometry.h.

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

{ return "EcalPreshower"; }

void EcalPreshowerGeometry::setzPlanes	(	CCGFloat	z1minus,
		CCGFloat	z2minus,
		CCGFloat	z1plus,
		CCGFloat	z2plus
	)

Definition at line 127 of file EcalPreshowerGeometry.cc.

References cms::cuda::assert(), and m_zplane.

Referenced by initializeParms().

                                                                                                           {
  assert(0 > z1minus && 0 > z2minus && 0 < z1plus && 0 < z2plus);

  m_zplane[0] = z1minus;
  m_zplane[1] = z2minus;
  m_zplane[2] = z1plus;
  m_zplane[3] = z2plus;
}

Member Data Documentation

CellVec EcalPreshowerGeometry::m_cellVec

private

Definition at line 103 of file EcalPreshowerGeometry.h.

Referenced by getGeometryRawPtr(), newCell(), and present().

const CCGFloat EcalPreshowerGeometry::m_xInterLadGap

private

Definition at line 95 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_xIntraLadGap

private

Definition at line 96 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_xWidWaf

private

Definition at line 94 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_yCtrOff

private

Definition at line 99 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

const CCGFloat EcalPreshowerGeometry::m_yWidAct

private

Definition at line 98 of file EcalPreshowerGeometry.h.

Referenced by getClosestCellInPlane().

CCGFloat EcalPreshowerGeometry::m_zplane[4]

private

Definition at line 101 of file EcalPreshowerGeometry.h.

Referenced by EcalPreshowerGeometry(), getClosestCellInPlane(), and setzPlanes().