FlatHexagon Class Reference

A base class to handle the hexagonal shape of HGCal silicon volumes. More...

#include <FlatHexagon.h>

Inheritance diagram for FlatHexagon:

Public Types
typedef CaloCellGeometry::CCGFloat	CCGFloat
typedef CaloCellGeometry::Pt3D	Pt3D
typedef CaloCellGeometry::Pt3DVec	Pt3DVec
typedef CaloCellGeometry::Tr3D	Tr3D
Public Types inherited from CaloCellGeometry
typedef float	CCGFloat
typedef EZMgrFL< GlobalPoint >	CornersMgr
typedef EZArrayFL< GlobalPoint >	CornersVec
typedef EZMgrFL< CCGFloat >	ParMgr
typedef EZArrayFL< CCGFloat >	ParVec
typedef std::vector< ParVec >	ParVecVec
typedef HepGeom::Point3D< CCGFloat >	Pt3D
typedef std::vector< Pt3D >	Pt3DVec
using	RepCorners = std::array< RhoEtaPhi, k_cornerSize >
typedef HepGeom::Transform3D	Tr3D

Public Member Functions
const GlobalVector &	axis () const
virtual float	etaPos () const
	FlatHexagon (void)
	FlatHexagon (const FlatHexagon &tr)
	FlatHexagon (CornersMgr *cMgr, const GlobalPoint &fCtr, const GlobalPoint &bCtr, const GlobalPoint &cor1, const CCGFloat *parV)
	FlatHexagon (const CornersVec &corn, const CCGFloat *par)
	FlatHexagon (const FlatHexagon &tr, const Pt3D &local)
Pt3D	getLocal (const GlobalPoint &global) const
CCGFloat	getPhiAxis () const
GlobalPoint const &	getPosition () const override
	Returns the position of reference for this cell. More...
GlobalPoint	getPosition (const Pt3D &local) const override
CCGFloat	getThetaAxis () const
void	getTransform (Tr3D &tr, Pt3DVec *lptr) const override
	--------— only needed by specific utility; overloaded when needed -— More...
FlatHexagon &	operator= (const FlatHexagon &tr)
virtual float	phiPos () const
void	setPosition (const GlobalPoint &p)
void	vocalCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const override
	~FlatHexagon () override
Public Member Functions inherited from CaloCellGeometry
bool	emptyCorners () const
float	etaPos () const
float	etaSpan () const
GlobalPoint const &	getBackPoint () const
CornersVec const &	getCorners () const
	Returns the corner points of this cell's volume. More...
RepCorners const &	getCornersREP () const
virtual GlobalPoint	getPosition (CCGFloat) const
bool	inside (const GlobalPoint &point) const
	Returns true if the specified point is inside this cell. More...
const CCGFloat *	param () const
float	phiPos () const
float	phiSpan () const
RhoEtaPhi const &	repPos () const
float	rhoPos () const
void	setBackPoint (const GlobalPoint &pos)
virtual	~CaloCellGeometry ()

Static Public Member Functions
static void	createCorners (const std::vector< CCGFloat > &pv, const Tr3D &tr, std::vector< GlobalPoint > &co)
static void	localCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref)
Static Public Member Functions inherited from CaloCellGeometry
static const CCGFloat *	checkParmPtr (const std::vector< CCGFloat > &vd, ParVecVec &pvv)
static const CCGFloat *	getParmPtr (const std::vector< CCGFloat > &vd, ParMgr *mgr, ParVecVec &pvv)

Static Public Attributes
static constexpr uint32_t	k_dZ = 0
static constexpr uint32_t	k_r = 1
static constexpr uint32_t	k_R = 2
static constexpr unsigned int	ncorner_ = 12
static constexpr unsigned int	ncornerBy2_ = 6
static constexpr double	oneBySix_ = 1.0 / 6.0
Static Public Attributes inherited from CaloCellGeometry
static constexpr unsigned int	k_cornerSize = 8
static const CCGFloat	k_ScaleFromDDDtoGeant

Private Member Functions
GlobalPoint	backCtr (void) const
void	initCorners (CornersVec &) override
GlobalVector	makeAxis (void)

Private Attributes
GlobalVector	m_axis
Pt3D	m_corOne
GlobalPoint	m_global
Pt3D	m_local
Tr3D	m_tr

Additional Inherited Members
Protected Member Functions inherited from CaloCellGeometry
	CaloCellGeometry (CornersVec::const_reference gp, CornersMgr *mgr, const CCGFloat *par)
	CaloCellGeometry (const CornersVec &cv, const CCGFloat *par)
	CaloCellGeometry (void)
void	initSpan ()
void	setCornerVec (const std::vector< GlobalPoint > &cor)
void	setRefPoint (const GlobalPoint &pos)

Detailed Description

A base class to handle the hexagonal shape of HGCal silicon volumes.

Definition at line 20 of file FlatHexagon.h.

Member Typedef Documentation

typedef CaloCellGeometry::CCGFloat FlatHexagon::CCGFloat

Definition at line 22 of file FlatHexagon.h.

typedef CaloCellGeometry::Pt3D FlatHexagon::Pt3D

Definition at line 23 of file FlatHexagon.h.

typedef CaloCellGeometry::Pt3DVec FlatHexagon::Pt3DVec

Definition at line 24 of file FlatHexagon.h.

typedef CaloCellGeometry::Tr3D FlatHexagon::Tr3D

Definition at line 25 of file FlatHexagon.h.

Constructor & Destructor Documentation

FlatHexagon::FlatHexagon

(

void

)

Definition at line 23 of file FlatHexagon.cc.

    : CaloCellGeometry(), m_axis(0., 0., 0.), m_corOne(0., 0., 0.), m_local(0., 0., 0.), m_global(0., 0., 0.) {}

FlatHexagon::FlatHexagon

(

const FlatHexagon &

tr

)

Definition at line 26 of file FlatHexagon.cc.

: CaloCellGeometry(tr) { *this = tr; }

FlatHexagon::FlatHexagon	(	CornersMgr *	cMgr,
		const GlobalPoint &	fCtr,
		const GlobalPoint &	bCtr,
		const GlobalPoint &	cor1,
		const CCGFloat *	parV
	)

Definition at line 45 of file FlatHexagon.cc.

References etaPos(), getTransform(), glb, m_global, m_local, m_tr, and phiPos().

    : CaloCellGeometry(fCtr, cMgr, parV),
      m_axis((bCtr - fCtr).unit()),
      m_corOne(cor1.x(), cor1.y(), cor1.z()),
      m_local(0., 0., 0.) {
  getTransform(m_tr, nullptr);
  Pt3D glb = m_tr * m_local;
  m_global = GlobalPoint(glb.x(), glb.y(), glb.z());
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "FlatHexagon: Local " << m_local << " Global " << glb << " eta " << etaPos()
                                   << " phi " << phiPos() << " Translation " << m_tr.getTranslation()
                                   << " and rotation " << m_tr.getRotation();
#endif
}

FlatHexagon::FlatHexagon	(	const CornersVec &	corn,
		const CCGFloat *	par
	)

Definition at line 61 of file FlatHexagon.cc.

References etaPos(), getTransform(), glb, m_axis, m_global, m_local, m_tr, makeAxis(), and phiPos().

    : CaloCellGeometry(corn, par), m_corOne(corn[0].x(), corn[0].y(), corn[0].z()), m_local(0., 0., 0.) {
  getTransform(m_tr, nullptr);
  m_axis = makeAxis();
  Pt3D glb = m_tr * m_local;
  m_global = GlobalPoint(glb.x(), glb.y(), glb.z());
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "FlatHexagon: Local " << m_local << " Global " << glb << " eta " << etaPos()
                                   << " phi " << phiPos() << " Translation " << m_tr.getTranslation()
                                   << " and rotation " << m_tr.getRotation();
#endif
}

FlatHexagon::FlatHexagon	(	const FlatHexagon &	tr,
		const Pt3D &	local
	)

Definition at line 74 of file FlatHexagon.cc.

References etaPos(), glb, m_global, m_local, m_tr, and phiPos().

                                                                 : CaloCellGeometry(tr), m_local(local) {
  *this = tr;
  Pt3D glb = m_tr * m_local;
  m_global = GlobalPoint(glb.x(), glb.y(), glb.z());
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "FlatHexagon: Local " << m_local << " Global " << glb << " eta " << etaPos()
                                   << " phi " << phiPos() << " Translation " << m_tr.getTranslation()
                                   << " and rotation " << m_tr.getRotation();
#endif
}

FlatHexagon::~FlatHexagon ( )

override

Definition at line 85 of file FlatHexagon.cc.

{}

Member Function Documentation

const GlobalVector & FlatHexagon::axis

(

)

const

Definition at line 112 of file FlatHexagon.cc.

References m_axis.

Referenced by getTransform(), and phiPos().

{ return m_axis; }

GlobalPoint FlatHexagon::backCtr ( void  ) const

private

Definition at line 225 of file FlatHexagon.cc.

References PVValHelper::dz, CaloCellGeometry::getCorners(), k_dZ, m_local, m_tr, ncornerBy2_, and CaloCellGeometry::param().

Referenced by makeAxis().

                                       {
  float dz =
      (getCorners()[ncornerBy2_].z() > getCorners()[0].z()) ? param()[FlatHexagon::k_dZ] : -param()[FlatHexagon::k_dZ];
  Pt3D local_b(m_local.x(), m_local.y(), m_local.z() + dz);
  Pt3D global_b = m_tr * local_b;
  GlobalPoint global(global_b.x(), global_b.y(), global_b.z());
  return global;
}

void FlatHexagon::createCorners ( const std::vector< CCGFloat > &  pv, const Tr3D &  tr, std::vector< GlobalPoint > &  co  )

static

Definition at line 116 of file FlatHexagon.cc.

References mps_fire::i, localCorners(), ncorner_, AlCaHLTBitMon_ParallelJobs::p, and createJobs::tmp.

Referenced by HGCalGeometryLoader::buildGeom(), phiPos(), and CaloGeometryDBEP< T, U >::produceAligned().

                                                                                                         {
  assert(3 == pv.size());
  assert(ncorner_ == co.size());

  Pt3DVec ko(ncorner_, Pt3D(0, 0, 0));

  Pt3D tmp;
  Pt3DVec to(ncorner_, Pt3D(0, 0, 0));
  localCorners(to, &pv.front(), tmp);

  for (unsigned int i(0); i != ncorner_; ++i) {
    ko[i] = tr * to[i];  // apply transformation
    const Pt3D& p(ko[i]);
    co[i] = GlobalPoint(p.x(), p.y(), p.z());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CaloGeometry") << "Corner[" << i << "] = " << co[i];
#endif
  }
}

virtual float FlatHexagon::etaPos ( ) const

inlinevirtual

Definition at line 51 of file FlatHexagon.h.

References PV3DBase< T, PVType, FrameType >::eta(), and m_global.

Referenced by FlatHexagon(), operator<<(), and operator=().

{ return m_global.eta(); }

Pt3D FlatHexagon::getLocal

(

const GlobalPoint &

global

)

const

Definition at line 99 of file FlatHexagon.cc.

References DTRecHitClients_cfi::local, m_tr, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by phiPos().

                                                          {
  Pt3D local = m_tr.inverse() * Pt3D(global.x(), global.y(), global.z());
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "FlatHexagon::Global " << global.x() << ":" << global.y() << ":" << global.z()
                                   << " Local " << local.x() << ":" << local.y() << ":" << local.z();
#endif
  return local;
}

CCGFloat FlatHexagon::getPhiAxis

(

)

const

Definition at line 110 of file FlatHexagon.cc.

References m_axis, and PV3DBase< T, PVType, FrameType >::phi().

Referenced by operator<<(), and phiPos().

{ return m_axis.phi(); }

GlobalPoint const& FlatHexagon::getPosition ( ) const

inlineoverridevirtual

Returns the position of reference for this cell.

Reimplemented from CaloCellGeometry.

Definition at line 49 of file FlatHexagon.h.

References DTRecHitClients_cfi::local, and m_global.

Referenced by makeAxis(), and operator<<().

{ return m_global; }

GlobalPoint FlatHexagon::getPosition ( const Pt3D &  local ) const

overridevirtual

Reimplemented from CaloCellGeometry.

Definition at line 87 of file FlatHexagon.cc.

References glb, DTRecHitClients_cfi::local, and m_tr.

                                                            {
  Pt3D glb = m_tr * local;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "FlatHexagon::Local " << local.x() << ":" << local.y() << ":" << local.z()
                                   << " Global " << glb.x() << ":" << glb.y() << ":" << glb.z() << " TR " << m_tr.xx()
                                   << ":" << m_tr.xy() << ":" << m_tr.xz() << ":" << m_tr.yx() << ":" << m_tr.yy()
                                   << ":" << m_tr.yz() << ":" << m_tr.zx() << ":" << m_tr.zy() << ":" << m_tr.zz()
                                   << ":" << m_tr.dx() << ":" << m_tr.dy() << ":" << m_tr.dz();
#endif
  return GlobalPoint(glb.x(), glb.y(), glb.z());
}

CCGFloat FlatHexagon::getThetaAxis

(

)

const

Definition at line 108 of file FlatHexagon.cc.

References m_axis, and PV3DBase< T, PVType, FrameType >::theta().

Referenced by operator<<(), and phiPos().

{ return m_axis.theta(); }

void FlatHexagon::getTransform ( Tr3D &  tr, Pt3DVec *  lptr  ) const

overridevirtual

--------— only needed by specific utility; overloaded when needed -—

Reimplemented from CaloCellGeometry.

Definition at line 165 of file FlatHexagon.cc.

References angle(), axis(), MillePedeFileConverter_cfg::e, CaloCellGeometry::getPosition(), localCorners(), m_corOne, mag(), ncorner_, AlCaHLTBitMon_ParallelJobs::p, CaloCellGeometry::param(), unit(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by FlatHexagon(), initCorners(), and phiPos().

                                                            {
  const GlobalPoint& p(CaloCellGeometry::getPosition());
  const Pt3D gFront(p.x(), p.y(), p.z());
  const DPt3D dgFront(p.x(), p.y(), p.z());

  Pt3D lFront;
  assert(nullptr != param());
  std::vector<Pt3D> lc(ncorner_, Pt3D(0, 0, 0));
  localCorners(lc, param(), lFront);

  // figure out if reflction volume or not

  Pt3D lBack((lc[6] + lc[7] + lc[8] + lc[9] + lc[10] + lc[11]) / 6.0);

  const DPt3D dlFront(lFront.x(), lFront.y(), lFront.z());
  const DPt3D dlBack(lBack.x(), lBack.y(), lBack.z());
  const DPt3D dlOne(lc[0].x(), lc[0].y(), lc[0].z());

  const FVec3D dgAxis(axis().x(), axis().y(), axis().z());

  const DPt3D dmOne(m_corOne.x(), m_corOne.y(), m_corOne.z());

  const DPt3D dgBack(dgFront + (dlBack - dlFront).mag() * dgAxis);
  DPt3D dgOne(dgFront + (dlOne - dlFront).mag() * (dmOne - dgFront).unit());

  const double dlangle((dlBack - dlFront).angle(dlOne - dlFront));
  const double dgangle((dgBack - dgFront).angle(dgOne - dgFront));
  const double dangle(dlangle - dgangle);

  if (1.e-6 < fabs(dangle)) {  //guard against precision problems
    const DPlane3D dgPl(dgFront, dgOne, dgBack);
    const DPt3D dp2(dgFront + dgPl.normal().unit());

    DPt3D dgOld(dgOne);

    dgOne = (dgFront + HepGeom::Rotate3D(-dangle, dgFront, dp2) * DVec3D(dgOld - dgFront));
  }

  tr = Tr3D(dlFront, dlBack, dlOne, dgFront, dgBack, dgOne);

  if (nullptr != lptr)
    (*lptr) = lc;
}

void FlatHexagon::initCorners ( CaloCellGeometry::CornersVec &  co )

overrideprivatevirtual

Implements CaloCellGeometry.

Definition at line 209 of file FlatHexagon.cc.

References getTransform(), mps_fire::i, ncorner_, and EZArrayFL< T >::uninitialized().

                                                            {
  if (co.uninitialized()) {
    CornersVec& corners(co);
    Pt3DVec lc;
    Tr3D tr;
    getTransform(tr, &lc);

    for (unsigned int i(0); i != ncorner_; ++i) {
      const Pt3D corn(tr * lc[i]);
      corners[i] = GlobalPoint(corn.x(), corn.y(), corn.z());
    }
  }
}

void FlatHexagon::localCorners ( Pt3DVec &  vec, const CCGFloat *  pv, Pt3D &  ref  )

static

Definition at line 136 of file FlatHexagon.cc.

References PVValHelper::dz, k_dZ, k_r, k_R, ncorner_, oneBySix_, dttmaxenums::R, and alignCSCRings::r.

Referenced by createCorners(), getTransform(), HGCalGeometry::localCorners(), phiPos(), and vocalCorners().

                                                                         {
  assert(nullptr != pv);
  assert(ncorner_ == lc.size());

  const CCGFloat dz(pv[FlatHexagon::k_dZ]);
  const CCGFloat r(pv[FlatHexagon::k_r]);
  const CCGFloat R(pv[FlatHexagon::k_R]);

  lc[0] = Pt3D(0, -2 * R, -dz);  // (0,-,-)
  lc[1] = Pt3D(-r, -R, -dz);     // (-,-,-)
  lc[2] = Pt3D(-r, R, -dz);      // (-,+,-)
  lc[3] = Pt3D(0, 2 * R, -dz);   // (0,+,-)
  lc[4] = Pt3D(r, R, -dz);       // (+,+,-)
  lc[5] = Pt3D(r, -R, -dz);      // (+,-,-)
  lc[6] = Pt3D(0, -2 * R, dz);   // (0,-,+)
  lc[7] = Pt3D(-r, -R, dz);      // (-,-,+)
  lc[8] = Pt3D(-r, R, dz);       // (-,+,+)
  lc[9] = Pt3D(0, 2 * R, dz);    // (0,+,+)
  lc[10] = Pt3D(r, R, dz);       // (+,+,+)
  lc[11] = Pt3D(r, -R, dz);      // (+,-,+)

  ref = oneBySix_ * (lc[0] + lc[1] + lc[2] + lc[3] + lc[4] + lc[5]);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "Ref " << ref << " Local Corners " << lc[0] << "|" << lc[1] << "|" << lc[2] << "|"
                                   << lc[3] << "|" << lc[4] << "|" << lc[5] << "|" << lc[6] << "|" << lc[7] << "|"
                                   << lc[8] << "|" << lc[9] << "|" << lc[10] << "|" << lc[11];
#endif
}

GlobalVector FlatHexagon::makeAxis ( void  )

private

Definition at line 223 of file FlatHexagon.cc.

References backCtr(), getPosition(), and Vector3DBase< T, FrameTag >::unit().

Referenced by FlatHexagon().

{ return GlobalVector(backCtr() - getPosition()).unit(); }

FlatHexagon & FlatHexagon::operator=

(

const FlatHexagon &

tr

)

Definition at line 28 of file FlatHexagon.cc.

References etaPos(), m_axis, m_corOne, m_global, m_local, m_tr, and phiPos().

                                                         {
  CaloCellGeometry::operator=(tr);
  if (this != &tr) {
    m_axis = tr.m_axis;
    m_corOne = tr.m_corOne;
    m_local = tr.m_local;
    m_global = tr.m_global;
    m_tr = tr.m_tr;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CaloGeometry") << "FlatHexagon(Copy): Local " << m_local << " Global " << m_global << " eta "
                                   << etaPos() << " phi " << phiPos() << " Translation " << m_tr.getTranslation()
                                   << " and rotation " << m_tr.getRotation();
#endif
  return *this;
}

virtual float FlatHexagon::phiPos ( ) const

inlinevirtual

Definition at line 52 of file FlatHexagon.h.

References axis(), createCorners(), getLocal(), getPhiAxis(), getThetaAxis(), getTransform(), localCorners(), m_global, PV3DBase< T, PVType, FrameType >::phi(), MetAnalyzer::pv(), and vocalCorners().

Referenced by FlatHexagon(), operator<<(), and operator=().

{ return m_global.phi(); }

void FlatHexagon::setPosition ( const GlobalPoint &  p )

inline

Definition at line 71 of file FlatHexagon.h.

References m_global, AlCaHLTBitMon_ParallelJobs::p, and CaloCellGeometry::setRefPoint().

                                         {
    m_global = p;
    setRefPoint(p);
  }

void FlatHexagon::vocalCorners ( Pt3DVec &  vec, const CCGFloat *  pv, Pt3D &  ref  ) const

overridevirtual

Implements CaloCellGeometry.

Definition at line 114 of file FlatHexagon.cc.

References localCorners().

Referenced by phiPos().

{ localCorners(vec, pv, ref); }

Member Data Documentation

constexpr uint32_t FlatHexagon::k_dZ = 0

static

Definition at line 27 of file FlatHexagon.h.

Referenced by backCtr(), HGCalGeometryLoader::build(), HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), HGCalGeometry::getNewCorners(), HGCalGeometry::getSummary(), and localCorners().

constexpr uint32_t FlatHexagon::k_r = 1

static

Definition at line 28 of file FlatHexagon.h.

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::get8Corners(), HGCalGeometry::getCorners(), HGCalGeometry::getNewCorners(), HGCalGeometry::getSummary(), and localCorners().

constexpr uint32_t FlatHexagon::k_R = 2

static

Definition at line 29 of file FlatHexagon.h.

Referenced by HGCalGeometryLoader::build(), HGCalGeometry::getCorners(), HGCalGeometry::getNewCorners(), HGCalGeometry::getSummary(), and localCorners().

GlobalVector FlatHexagon::m_axis

private

Definition at line 86 of file FlatHexagon.h.

Referenced by axis(), FlatHexagon(), getPhiAxis(), getThetaAxis(), and operator=().

Pt3D FlatHexagon::m_corOne

private

Definition at line 87 of file FlatHexagon.h.

Referenced by getTransform(), and operator=().

GlobalPoint FlatHexagon::m_global

private

Definition at line 88 of file FlatHexagon.h.

Referenced by etaPos(), FlatHexagon(), getPosition(), operator=(), phiPos(), and setPosition().

Pt3D FlatHexagon::m_local

private

Definition at line 87 of file FlatHexagon.h.

Referenced by backCtr(), FlatHexagon(), and operator=().

Tr3D FlatHexagon::m_tr

private

Definition at line 89 of file FlatHexagon.h.

Referenced by backCtr(), FlatHexagon(), getLocal(), getPosition(), and operator=().

constexpr unsigned int FlatHexagon::ncorner_ = 12

static

Definition at line 77 of file FlatHexagon.h.

Referenced by HGCalGeometryLoader::buildGeom(), createCorners(), HGCalGeometry::getCorners(), getTransform(), initCorners(), localCorners(), and CaloGeometryDBEP< T, U >::produceAligned().

constexpr unsigned int FlatHexagon::ncornerBy2_ = 6

static

Definition at line 78 of file FlatHexagon.h.

Referenced by backCtr(), HGCalGeometryLoader::buildGeom(), and CaloGeometryDBEP< T, U >::produceAligned().

constexpr double FlatHexagon::oneBySix_ = 1.0 / 6.0

static

Definition at line 76 of file FlatHexagon.h.

Referenced by HGCalGeometryLoader::buildGeom(), localCorners(), and CaloGeometryDBEP< T, U >::produceAligned().