IdealObliquePrism Class Reference

#include <IdealObliquePrism.h>

Inheritance diagram for IdealObliquePrism:

Public Types
typedef CaloCellGeometry::CCGFloat	CCGFloat
typedef CaloCellGeometry::Pt3D	Pt3D
typedef CaloCellGeometry::Pt3DVec	Pt3DVec
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
CCGFloat	dEta () const
CCGFloat	dPhi () const
CCGFloat	dz () const
CCGFloat	eta () const
	IdealObliquePrism ()
	IdealObliquePrism (const IdealObliquePrism &idop)
	IdealObliquePrism (const GlobalPoint &faceCenter, CornersMgr *mgr, const CCGFloat *parm)
IdealObliquePrism &	operator= (const IdealObliquePrism &idop)
void	vocalCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const override
CCGFloat	z () const
	~IdealObliquePrism () 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 const GlobalPoint &	getPosition () const
	Returns the position of reference for this cell. More...
virtual GlobalPoint	getPosition (CCGFloat) const
virtual GlobalPoint	getPosition (const Pt3D &) const
virtual void	getTransform (Tr3D &tr, Pt3DVec *lptr) const
	--------— only needed by specific utility; overloaded when needed -— More...
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	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_dEta = 0
static constexpr uint32_t	k_dPhi = 1
static constexpr uint32_t	k_dZ = 2
static constexpr uint32_t	k_Eta = 3
static constexpr uint32_t	k_Z = 4
Static Public Attributes inherited from CaloCellGeometry
static constexpr unsigned int	k_cornerSize = 8
static const CCGFloat	k_ScaleFromDDDtoGeant

Private Member Functions
void	initCorners (CornersVec &) override

Static Private Member Functions
static GlobalPoint	etaPhiPerp (float eta, float phi, float perp)
static GlobalPoint	etaPhiZ (float eta, float phi, float z)

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

Oblique prism class used for HCAL (HB, HE, HO) volumes.

Required parameters for an ideal oblique prism:

• eta, phi of axis
• radial distance (along axis) to front and back faces
• single bit - faces parallel or perpendicular to z-axis
• eta width and phi width of faces (same for front/back)

Total: 6+1 parameters

Internally, the "point of reference" is the center (eta/phi) of the front face of the prism. Therefore, the only internally stored parameters are eta and phi widths, the axis tower thickness, and the parallel/perpendicular setting. The parallel/perpendicular setting is encoded in the sign of the thickness. (positive = parallel to z-axis, negative = perpendicular)

Author

J. Mans - Minnesota

Definition at line 28 of file IdealObliquePrism.h.

Member Typedef Documentation

typedef CaloCellGeometry::CCGFloat IdealObliquePrism::CCGFloat

Definition at line 30 of file IdealObliquePrism.h.

typedef CaloCellGeometry::Pt3D IdealObliquePrism::Pt3D

Definition at line 31 of file IdealObliquePrism.h.

typedef CaloCellGeometry::Pt3DVec IdealObliquePrism::Pt3DVec

Definition at line 32 of file IdealObliquePrism.h.

Constructor & Destructor Documentation

IdealObliquePrism::IdealObliquePrism

(

)

Definition at line 8 of file IdealObliquePrism.cc.

: CaloCellGeometry() {}

IdealObliquePrism::IdealObliquePrism

(

const IdealObliquePrism &

idop

)

Definition at line 10 of file IdealObliquePrism.cc.

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

IdealObliquePrism::IdealObliquePrism	(	const GlobalPoint &	faceCenter,
		CornersMgr *	mgr,
		const CCGFloat *	parm
	)

Definition at line 18 of file IdealObliquePrism.cc.

References CaloCellGeometry::initSpan().

    : CaloCellGeometry(faceCenter, mgr, parm) {
  initSpan();
}

IdealObliquePrism::~IdealObliquePrism ( )

override

Definition at line 23 of file IdealObliquePrism.cc.

{}

Member Function Documentation

CCGFloat IdealObliquePrism::dEta

(

)

const

Definition at line 25 of file IdealObliquePrism.cc.

References k_dEta, and CaloCellGeometry::param().

Referenced by initCorners(), localCorners(), and operator<<().

{ return param()[IdealObliquePrism::k_dEta]; }

CCGFloat IdealObliquePrism::dPhi

(

)

const

Definition at line 27 of file IdealObliquePrism.cc.

References k_dPhi, and CaloCellGeometry::param().

Referenced by initCorners(), localCorners(), and operator<<().

{ return param()[IdealObliquePrism::k_dPhi]; }

CCGFloat IdealObliquePrism::dz

(

)

const

Definition at line 29 of file IdealObliquePrism.cc.

References k_dZ, and CaloCellGeometry::param().

Referenced by Electron.Electron::cutBasedId(), ntupleDataFormat.Track::dzPull(), initCorners(), localCorners(), and operator<<().

{ return param()[IdealObliquePrism::k_dZ]; }

CCGFloat IdealObliquePrism::eta

(

void

)

const

Definition at line 31 of file IdealObliquePrism.cc.

References k_Eta, and CaloCellGeometry::param().

Referenced by Particle.Particle::__str__(), initCorners(), Jet.Jet::jetID(), localCorners(), and Jet.Jet::puJetId().

{ return param()[IdealObliquePrism::k_Eta]; }

GlobalPoint IdealObliquePrism::etaPhiPerp ( float  eta, float  phi, float  perp  )

staticprivate

Definition at line 37 of file IdealObliquePrism.cc.

Referenced by initCorners(), and localCorners().

                                                                          {
  return GlobalPoint(perp * cosf(phi), perp * sinf(phi), perp * sinhf(eta));
}

GlobalPoint IdealObliquePrism::etaPhiZ ( float  eta, float  phi, float  z  )

staticprivate

Definition at line 41 of file IdealObliquePrism.cc.

Referenced by initCorners(), and localCorners().

                                                                    {
  return GlobalPoint(z * cosf(phi) / sinhf(eta), z * sinf(phi) / sinhf(eta), z);
}

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

overrideprivatevirtual

Implements CaloCellGeometry.

Definition at line 89 of file IdealObliquePrism.cc.

References funct::cos(), dEta(), dPhi(), dz(), eta(), PV3DBase< T, PVType, FrameType >::eta(), etaPhiPerp(), etaPhiZ(), CaloCellGeometry::getPosition(), PV3DBase< T, PVType, FrameType >::mag(), mag(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), phi, PV3DBase< T, PVType, FrameType >::phi(), EZArrayFL< T >::uninitialized(), and PV3DBase< T, PVType, FrameType >::z().

                                                                  {
  if (co.uninitialized()) {
    CornersVec& corners(co);
    if (dz() > 0) {
      /* In this case, the faces are parallel to the zaxis.  
        This implies that all corners will have the same 
        cylindrical radius. 
     */
      const GlobalPoint p(getPosition());
      const CCGFloat r_near(p.perp() / cos(dPhi()));
      const CCGFloat r_far(r_near * ((p.mag() + 2 * dz()) / p.mag()));
      const CCGFloat eta(p.eta());
      const CCGFloat phi(p.phi());
      corners[0] = etaPhiPerp(eta + dEta(), phi + dPhi(), r_near);  // (+,+,near)
      corners[1] = etaPhiPerp(eta + dEta(), phi - dPhi(), r_near);  // (+,-,near)
      corners[2] = etaPhiPerp(eta - dEta(), phi - dPhi(), r_near);  // (-,-,near)
      corners[3] = etaPhiPerp(eta - dEta(), phi + dPhi(), r_near);  // (-,+,near)
      corners[4] = etaPhiPerp(eta + dEta(), phi + dPhi(), r_far);   // (+,+,far)
      corners[5] = etaPhiPerp(eta + dEta(), phi - dPhi(), r_far);   // (+,-,far)
      corners[6] = etaPhiPerp(eta - dEta(), phi - dPhi(), r_far);   // (-,-,far)
      corners[7] = etaPhiPerp(eta - dEta(), phi + dPhi(), r_far);   // (-,+,far)
    } else {
      /* In this case, the faces are perpendicular to the zaxis.  
        This implies that all corners will have the same 
        z-dimension. 
     */
      const GlobalPoint p(getPosition());
      const CCGFloat z_near(p.z());
      const CCGFloat mag(p.mag());
      const CCGFloat z_far(z_near * (1 - 2 * dz() / mag));  // negative to correct sign
      const CCGFloat eta(p.eta());
      const CCGFloat phi(p.phi());

      corners[0] = etaPhiZ(eta + dEta(), phi + dPhi(), z_near);  // (+,+,near)
      corners[1] = etaPhiZ(eta + dEta(), phi - dPhi(), z_near);  // (+,-,near)
      corners[2] = etaPhiZ(eta - dEta(), phi - dPhi(), z_near);  // (-,-,near)
      corners[3] = etaPhiZ(eta - dEta(), phi + dPhi(), z_near);  // (-,+,near)
      corners[4] = etaPhiZ(eta + dEta(), phi + dPhi(), z_far);   // (+,+,far)
      corners[5] = etaPhiZ(eta + dEta(), phi - dPhi(), z_far);   // (+,-,far)
      corners[6] = etaPhiZ(eta - dEta(), phi - dPhi(), z_far);   // (-,-,far)
      corners[7] = etaPhiZ(eta - dEta(), phi + dPhi(), z_far);   // (-,+,far)
    }
  }
}

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

static

Definition at line 45 of file IdealObliquePrism.cc.

References funct::cos(), dEta(), dPhi(), dz(), eta(), etaPhiPerp(), etaPhiZ(), mps_fire::i, k_dEta, k_dPhi, k_dZ, k_Eta, k_Z, PV3DBase< T, PVType, FrameType >::mag(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), x, y, and z().

Referenced by FWTGeoRecoGeometryESProducer::addCaloTowerGeometry(), FWTGeoRecoGeometryESProducer::addHcalCaloGeometryBarrel(), FWTGeoRecoGeometryESProducer::addHcalCaloGeometryEndcap(), FWTGeoRecoGeometryESProducer::addHcalCaloGeometryOuter(), CaloTowerGeometry::localCorners(), HcalGeometry::localCorners(), and vocalCorners().

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

  const CCGFloat dEta(pv[IdealObliquePrism::k_dEta]);
  const CCGFloat dPhi(pv[IdealObliquePrism::k_dPhi]);
  const CCGFloat dz(pv[IdealObliquePrism::k_dZ]);
  const CCGFloat eta(pv[IdealObliquePrism::k_Eta]);
  const CCGFloat z(pv[IdealObliquePrism::k_Z]);

  std::vector<GlobalPoint> gc(8, GlobalPoint(0, 0, 0));

  const GlobalPoint p(etaPhiZ(eta, 0, z));

  if (0 < dz) {
    const CCGFloat r_near(p.perp() / cos(dPhi));
    const CCGFloat r_far(r_near * ((p.mag() + 2 * dz) / p.mag()));
    gc[0] = etaPhiPerp(eta + dEta, +dPhi, r_near);  // (+,+,near)
    gc[1] = etaPhiPerp(eta + dEta, -dPhi, r_near);  // (+,-,near)
    gc[2] = etaPhiPerp(eta - dEta, -dPhi, r_near);  // (-,-,near)
    gc[3] = etaPhiPerp(eta - dEta, +dPhi, r_near);  // (-,+,near)
    gc[4] = etaPhiPerp(eta + dEta, +dPhi, r_far);   // (+,+,far)
    gc[5] = etaPhiPerp(eta + dEta, -dPhi, r_far);   // (+,-,far)
    gc[6] = etaPhiPerp(eta - dEta, -dPhi, r_far);   // (-,-,far)
    gc[7] = etaPhiPerp(eta - dEta, +dPhi, r_far);   // (-,+,far)
  } else {
    const CCGFloat z_near(z);
    const CCGFloat z_far(z * (1 - 2 * dz / p.mag()));
    gc[0] = etaPhiZ(eta + dEta, +dPhi, z_near);  // (+,+,near)
    gc[1] = etaPhiZ(eta + dEta, -dPhi, z_near);  // (+,-,near)
    gc[2] = etaPhiZ(eta - dEta, -dPhi, z_near);  // (-,-,near)
    gc[3] = etaPhiZ(eta - dEta, +dPhi, z_near);  // (-,+,near)
    gc[4] = etaPhiZ(eta + dEta, +dPhi, z_far);   // (+,+,far)
    gc[5] = etaPhiZ(eta + dEta, -dPhi, z_far);   // (+,-,far)
    gc[6] = etaPhiZ(eta - dEta, -dPhi, z_far);   // (-,-,far)
    gc[7] = etaPhiZ(eta - dEta, +dPhi, z_far);   // (-,+,far)
  }
  for (unsigned int i(0); i != 8; ++i) {
    lc[i] = Pt3D(gc[i].x(), gc[i].y(), gc[i].z());
  }

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

IdealObliquePrism & IdealObliquePrism::operator=

(

const IdealObliquePrism &

idop

)

Definition at line 12 of file IdealObliquePrism.cc.

                                                                             {
  if (&idop != this)
    CaloCellGeometry::operator=(idop);
  return *this;
}

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

overridevirtual

Implements CaloCellGeometry.

Definition at line 35 of file IdealObliquePrism.cc.

References localCorners().

{ localCorners(vec, pv, ref); }

CCGFloat IdealObliquePrism::z

(

)

const

Definition at line 33 of file IdealObliquePrism.cc.

References k_Z, and CaloCellGeometry::param().

Referenced by localCorners(), geometryXMLparser.Alignable::pos(), and ntupleDataFormat._HitObject::r3D().

{ return param()[IdealObliquePrism::k_Z]; }

Member Data Documentation

constexpr uint32_t IdealObliquePrism::k_dEta = 0

static

Definition at line 34 of file IdealObliquePrism.h.

Referenced by dEta(), and localCorners().

constexpr uint32_t IdealObliquePrism::k_dPhi = 1

static

Definition at line 35 of file IdealObliquePrism.h.

Referenced by dPhi(), and localCorners().

constexpr uint32_t IdealObliquePrism::k_dZ = 2

static

Definition at line 36 of file IdealObliquePrism.h.

Referenced by dz(), and localCorners().

constexpr uint32_t IdealObliquePrism::k_Eta = 3

static

Definition at line 37 of file IdealObliquePrism.h.

Referenced by eta(), and localCorners().

constexpr uint32_t IdealObliquePrism::k_Z = 4

static

Definition at line 38 of file IdealObliquePrism.h.

Referenced by localCorners(), and z().