TruncatedPyramid Class Reference

A base class to handle the particular shape of Ecal Xtals. Taken from ORCA Calorimetry Code. More...

#include <TruncatedPyramid.h>

Inheritance diagram for TruncatedPyramid:

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
CCGFloat	getPhiAxis () const
GlobalPoint	getPosition (CCGFloat depth) const override
CCGFloat	getThetaAxis () const
void	getTransform (Tr3D &tr, Pt3DVec *lptr) const override
	--------— only needed by specific utility; overloaded when needed -— More...
TruncatedPyramid &	operator= (const TruncatedPyramid &tr)
	TruncatedPyramid (void)
	TruncatedPyramid (const TruncatedPyramid &tr)
	TruncatedPyramid (CornersMgr *cMgr, const GlobalPoint &fCtr, const GlobalPoint &bCtr, const GlobalPoint &cor1, const CCGFloat *parV)
	TruncatedPyramid (const CornersVec &corn, const CCGFloat *par)
void	vocalCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const override
	~TruncatedPyramid () override
Public Member Functions inherited from CaloCellGeometry
bool	emptyCorners () const
virtual float	etaPos () const
virtual 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 const GlobalPoint &	getPosition () const
	Returns the position of reference for this cell. More...
virtual GlobalPoint	getPosition (const Pt3D &) const
bool	inside (const GlobalPoint &point) const
	Returns true if the specified point is inside this cell. More...
const CCGFloat *	param () const
virtual float	phiPos () const
virtual float	phiSpan () const
RhoEtaPhi const &	repPos () const
virtual 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 void	localCornersReflection (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref)
static void	localCornersSwap (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_Alp1 = 6
static constexpr uint32_t	k_Alp2 = 10
static constexpr uint32_t	k_Dx1 = 4
static constexpr uint32_t	k_Dx2 = 5
static constexpr uint32_t	k_Dx3 = 8
static constexpr uint32_t	k_Dx4 = 9
static constexpr uint32_t	k_Dy1 = 3
static constexpr uint32_t	k_Dy2 = 7
static constexpr uint32_t	k_Dz = 0
static constexpr uint32_t	k_Phi = 2
static constexpr uint32_t	k_Theta = 1
Static Public Attributes inherited from CaloCellGeometry
static constexpr unsigned int	k_cornerSize = 8
static const CCGFloat	k_ScaleFromDDDtoGeant

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

Private Attributes
GlobalVector	m_axis
Pt3D	m_corOne

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 particular shape of Ecal Xtals. Taken from ORCA Calorimetry Code.

Definition at line 19 of file TruncatedPyramid.h.

Member Typedef Documentation

typedef CaloCellGeometry::CCGFloat TruncatedPyramid::CCGFloat

Definition at line 21 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Pt3D TruncatedPyramid::Pt3D

Definition at line 22 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Pt3DVec TruncatedPyramid::Pt3DVec

Definition at line 23 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Tr3D TruncatedPyramid::Tr3D

Definition at line 24 of file TruncatedPyramid.h.

Constructor & Destructor Documentation

TruncatedPyramid::TruncatedPyramid

(

void

)

Definition at line 19 of file TruncatedPyramid.cc.

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

TruncatedPyramid::TruncatedPyramid

(

const TruncatedPyramid &

tr

)

Definition at line 21 of file TruncatedPyramid.cc.

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

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

Definition at line 32 of file TruncatedPyramid.cc.

References CaloCellGeometry::initSpan().

    : CaloCellGeometry(fCtr, cMgr, parV), m_axis((bCtr - fCtr).unit()), m_corOne(cor1.x(), cor1.y(), cor1.z()) {
  initSpan();
}

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

Definition at line 38 of file TruncatedPyramid.cc.

References CaloCellGeometry::initSpan().

    : CaloCellGeometry(corn, par), m_axis(makeAxis()), m_corOne(corn[0].x(), corn[0].y(), corn[0].z()) {
  initSpan();
}

TruncatedPyramid::~TruncatedPyramid ( )

override

Definition at line 43 of file TruncatedPyramid.cc.

{}

Member Function Documentation

const GlobalVector & TruncatedPyramid::axis

(

)

const

Definition at line 53 of file TruncatedPyramid.cc.

References m_axis.

Referenced by getTransform().

{ return m_axis; }

const GlobalPoint TruncatedPyramid::backCtr ( void  ) const

private

Definition at line 59 of file TruncatedPyramid.cc.

References CaloCellGeometry::getCorners(), x, y, and z.

Referenced by makeAxis().

                                                  {
  return GlobalPoint(0.25 * (getCorners()[4].x() + getCorners()[5].x() + getCorners()[6].x() + getCorners()[7].x()),
                     0.25 * (getCorners()[4].y() + getCorners()[5].y() + getCorners()[6].y() + getCorners()[7].y()),
                     0.25 * (getCorners()[4].z() + getCorners()[5].z() + getCorners()[6].z() + getCorners()[7].z()));
}

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

static

Definition at line 237 of file TruncatedPyramid.cc.

References cms::cuda::assert(), filterCSVwithJSON::copy, cross(), dot(), PVValHelper::dz, mps_fire::i, localCorners(), AlCaHLTBitMon_ParallelJobs::p, truncPyr::refl(), createJobs::tmp, x, y, and z.

                                                                                                              {
  assert(11 == pv.size());
  assert(8 == co.size());
  // to get the ordering right for fast sim, we have to use their convention
  // which were based on the old static geometry. Some gymnastics required here.

  const CCGFloat dz(pv[0]);
  const CCGFloat h1(pv[3]);
  const CCGFloat h2(pv[7]);
  Pt3DVec ko(8, Pt3D(0, 0, 0));

  // if reflection, different things for barrel and endcap
  static const FVec3D x(1, 0, 0);
  static const FVec3D y(0, 1, 0);
  static const FVec3D z(0, 0, 1);
  const bool refl(((tr * x).cross(tr * y)).dot(tr * z) < 0);  // has reflection!

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

  for (unsigned int i(0); i != 8; ++i) {
    ko[i] = tr * to[i];  // apply transformation
  }

  if (refl || h1 > h2) {
    if (11.2 < dz)  //barrel
    {
      if (!refl) {
        to[0] = ko[3];
        to[1] = ko[2];
        to[2] = ko[1];
        to[3] = ko[0];
        to[4] = ko[7];
        to[5] = ko[6];
        to[6] = ko[5];
        to[7] = ko[4];
      } else {
        to[0] = ko[0];
        to[1] = ko[1];
        to[2] = ko[2];
        to[3] = ko[3];
        to[4] = ko[4];
        to[5] = ko[5];
        to[6] = ko[6];
        to[7] = ko[7];
      }
    } else  //endcap
    {
      to[0] = ko[0];
      to[1] = ko[3];
      to[2] = ko[2];
      to[3] = ko[1];
      to[4] = ko[4];
      to[5] = ko[7];
      to[6] = ko[6];
      to[7] = ko[5];
    }
    copy(to.begin(), to.end(), ko.begin());  // faster than ko = to ? maybe.
  }
  for (unsigned int i(0); i != 8; ++i) {
    const Pt3D& p(ko[i]);
    co[i] = GlobalPoint(p.x(), p.y(), p.z());
  }
}

CCGFloat TruncatedPyramid::getPhiAxis

(

)

const

Definition at line 51 of file TruncatedPyramid.cc.

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

Referenced by egammaTools::localEcalClusterCoordsEB(), egammaTools::localEcalClusterCoordsEE(), and operator<<().

{ return m_axis.phi(); }

GlobalPoint TruncatedPyramid::getPosition ( CCGFloat  depth ) const

overridevirtual

Reimplemented from CaloCellGeometry.

Definition at line 45 of file TruncatedPyramid.cc.

References CaloCellGeometry::getPosition(), and m_axis.

Referenced by egammaTools::localEcalClusterCoordsEB(), and egammaTools::localEcalClusterCoordsEE().

                                                              {
  return CaloCellGeometry::getPosition() + depth * m_axis;
}

CCGFloat TruncatedPyramid::getThetaAxis

(

)

const

Definition at line 49 of file TruncatedPyramid.cc.

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

Referenced by egammaTools::localEcalClusterCoordsEB(), egammaTools::localEcalClusterCoordsEE(), and operator<<().

{ return m_axis.theta(); }

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

overridevirtual

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

Reimplemented from CaloCellGeometry.

Definition at line 65 of file TruncatedPyramid.cc.

References angle(), cms::cuda::assert(), axis(), PVValHelper::dz, alignCSCRings::e, CaloCellGeometry::getPosition(), localCorners(), localCornersReflection(), localCornersSwap(), m_corOne, mag(), 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 initCorners().

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

  const double dz(param()[0]);

  Pt3D lFront;
  assert(nullptr != param());
  std::vector<Pt3D> lc(8, Pt3D(0, 0, 0));
  if (11.2 > dz) {
    localCorners(lc, param(), lFront);
  } else {
    localCornersSwap(lc, param(), lFront);
  }

  // figure out if reflction volume or not

  Pt3D lBack(0.25 * (lc[4] + lc[5] + lc[6] + lc[7]));

  const double disl((lFront - lc[0]).mag());
  const double disr((lFront - lc[3]).mag());
  const double disg((gFront - m_corOne).mag());

  const double dell(fabs(disg - disl));
  const double delr(fabs(disg - disr));

  if (11.2 < dz && delr < dell)  // reflection volume if true
  {
    localCornersReflection(lc, param(), lFront);
    lBack = 0.25 * (lc[4] + lc[5] + lc[6] + lc[7]);
  }

  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 TruncatedPyramid::initCorners ( CaloCellGeometry::CornersVec &  corners )

overrideprivatevirtual

Implements CaloCellGeometry.

Definition at line 129 of file TruncatedPyramid.cc.

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

                                                                      {
  if (corners.uninitialized()) {
    Pt3DVec lc;

    Tr3D tr;
    getTransform(tr, &lc);

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

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

static

Definition at line 177 of file TruncatedPyramid.cc.

References isotrackTrainRegressor::a1, isotrackTrainRegressor::a2, cms::cuda::assert(), b1, b2, funct::cos(), PVValHelper::dz, k_Alp1, k_Alp2, k_Dx1, k_Dx2, k_Dx3, k_Dx4, k_Dy1, k_Dy2, k_Dz, k_Phi, k_Theta, funct::sin(), and funct::tan().

Referenced by createCorners(), getTransform(), EcalEndcapGeometry::localCorners(), localCornersReflection(), localCornersSwap(), makeEcalShape(), and vocalCorners().

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

  const CCGFloat dz(pv[TruncatedPyramid::k_Dz]);
  const CCGFloat th(pv[TruncatedPyramid::k_Theta]);
  const CCGFloat ph(pv[TruncatedPyramid::k_Phi]);
  const CCGFloat h1(pv[TruncatedPyramid::k_Dy1]);
  const CCGFloat b1(pv[TruncatedPyramid::k_Dx1]);
  const CCGFloat t1(pv[TruncatedPyramid::k_Dx2]);
  const CCGFloat a1(pv[TruncatedPyramid::k_Alp1]);
  const CCGFloat h2(pv[TruncatedPyramid::k_Dy2]);
  const CCGFloat b2(pv[TruncatedPyramid::k_Dx3]);
  const CCGFloat t2(pv[TruncatedPyramid::k_Dx4]);
  const CCGFloat a2(pv[TruncatedPyramid::k_Alp2]);

  const CCGFloat ta1(tan(a1));  // lower plane
  const CCGFloat ta2(tan(a2));  // upper plane

  const CCGFloat tth(tan(th));
  const CCGFloat tthcp(tth * cos(ph));
  const CCGFloat tthsp(tth * sin(ph));

  const unsigned int off(h1 < h2 ? 0 : 4);

  lc[0 + off] = Pt3D(-dz * tthcp - h1 * ta1 - b1, -dz * tthsp - h1, -dz);  // (-,-,-)
  lc[1 + off] = Pt3D(-dz * tthcp + h1 * ta1 - t1, -dz * tthsp + h1, -dz);  // (-,+,-)
  lc[2 + off] = Pt3D(-dz * tthcp + h1 * ta1 + t1, -dz * tthsp + h1, -dz);  // (+,+,-)
  lc[3 + off] = Pt3D(-dz * tthcp - h1 * ta1 + b1, -dz * tthsp - h1, -dz);  // (+,-,-)
  lc[4 - off] = Pt3D(dz * tthcp - h2 * ta2 - b2, dz * tthsp - h2, dz);     // (-,-,+)
  lc[5 - off] = Pt3D(dz * tthcp + h2 * ta2 - t2, dz * tthsp + h2, dz);     // (-,+,+)
  lc[6 - off] = Pt3D(dz * tthcp + h2 * ta2 + t2, dz * tthsp + h2, dz);     // (+,+,+)
  lc[7 - off] = Pt3D(dz * tthcp - h2 * ta2 + b2, dz * tthsp - h2, dz);     // (+,-,+)

  ref = 0.25 * (lc[0] + lc[1] + lc[2] + lc[3]);
}

void TruncatedPyramid::localCornersReflection ( Pt3DVec &  vec, const CCGFloat *  pv, Pt3D &  ref  )

static

Definition at line 147 of file TruncatedPyramid.cc.

References localCorners(), truncPyr::refl(), and createJobs::tmp.

Referenced by getTransform(), and EcalBarrelGeometry::localCorners().

                                                                                        {
  //   using namespace truncPyr ;
  localCorners(lc, pv, ref);
  Pt3D tmp;
  /*
   tmp   = lc[0] ;
   lc[0] = refl( lc[2] ) ;
   lc[2] = refl( tmp   ) ;
   tmp   = lc[1] ;
   lc[1] = refl( lc[3] ) ;
   lc[3] = refl( tmp   ) ;
   tmp   = lc[4] ;
   lc[4] = refl( lc[6] ) ;
   lc[6] = refl( tmp   ) ;
   tmp   = lc[5] ;
   lc[5] = refl( lc[7] ) ;
   lc[7] = refl( tmp   ) ;
*/
  lc[0] = truncPyr::refl(lc[0]);
  lc[1] = truncPyr::refl(lc[1]);
  lc[2] = truncPyr::refl(lc[2]);
  lc[3] = truncPyr::refl(lc[3]);
  lc[4] = truncPyr::refl(lc[4]);
  lc[5] = truncPyr::refl(lc[5]);
  lc[6] = truncPyr::refl(lc[6]);
  lc[7] = truncPyr::refl(lc[7]);

  ref = 0.25 * (lc[0] + lc[1] + lc[2] + lc[3]);
}

void TruncatedPyramid::localCornersSwap ( Pt3DVec &  vec, const CCGFloat *  pv, Pt3D &  ref  )

static

Definition at line 214 of file TruncatedPyramid.cc.

References localCorners(), and createJobs::tmp.

Referenced by getTransform(), and EcalBarrelGeometry::localCorners().

                                                                                  {
  localCorners(lc, pv, ref);

  Pt3D tmp;
  tmp = lc[0];
  lc[0] = lc[3];
  lc[3] = tmp;
  tmp = lc[1];
  lc[1] = lc[2];
  lc[2] = tmp;
  tmp = lc[4];
  lc[4] = lc[7];
  lc[7] = tmp;
  tmp = lc[5];
  lc[5] = lc[6];
  lc[6] = tmp;

  ref = 0.25 * (lc[0] + lc[1] + lc[2] + lc[3]);
}

GlobalVector TruncatedPyramid::makeAxis ( void  )

private

Definition at line 57 of file TruncatedPyramid.cc.

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

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

TruncatedPyramid & TruncatedPyramid::operator=

(

const TruncatedPyramid &

tr

)

Definition at line 23 of file TruncatedPyramid.cc.

References m_axis, and m_corOne.

                                                                        {
  CaloCellGeometry::operator=(tr);
  if (this != &tr) {
    m_axis = tr.m_axis;
    m_corOne = tr.m_corOne;
  }
  return *this;
}

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

overridevirtual

Implements CaloCellGeometry.

Definition at line 55 of file TruncatedPyramid.cc.

References localCorners().

{ localCorners(vec, pv, ref); }

Member Data Documentation

constexpr uint32_t TruncatedPyramid::k_Alp1 = 6

static

Definition at line 36 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Alp2 = 10

static

Definition at line 43 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dx1 = 4

static

Definition at line 32 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dx2 = 5

static

Definition at line 34 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dx3 = 8

static

Definition at line 39 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dx4 = 9

static

Definition at line 41 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dy1 = 3

static

Definition at line 31 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dy2 = 7

static

Definition at line 38 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Dz = 0

static

Definition at line 26 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Phi = 2

static

Definition at line 29 of file TruncatedPyramid.h.

Referenced by localCorners().

constexpr uint32_t TruncatedPyramid::k_Theta = 1

static

Definition at line 27 of file TruncatedPyramid.h.

Referenced by localCorners().

GlobalVector TruncatedPyramid::m_axis

private

Definition at line 91 of file TruncatedPyramid.h.

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

Pt3D TruncatedPyramid::m_corOne

private

Definition at line 92 of file TruncatedPyramid.h.

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