CaloCellGeometry Class Reference

#include <CaloCellGeometry.h>

Inheritance diagram for CaloCellGeometry:

Public Types
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
bool	emptyCorners () const
float	etaPos () const
float	etaSpan () const
const GlobalPoint &	getBackPoint () const
CornersVec const &	getCorners () const
	Returns the corner points of this cell's volume. More...
RepCorners const &	getCornersREP () const
const GlobalPoint &	getPosition () const
	Returns the position of reference for this cell. More...
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
virtual void	vocalCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const =0
virtual	~CaloCellGeometry ()

Static Public Member Functions
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 unsigned int	k_cornerSize = 8
static const CCGFloat	k_ScaleFromDDDtoGeant

Protected Member Functions
	CaloCellGeometry (CornersVec::const_reference gp, CornersMgr *mgr, const CCGFloat *par)
	CaloCellGeometry (const CornersVec &cv, const CCGFloat *par)
	CaloCellGeometry (void)
virtual void	initCorners (CornersVec &)=0
void	initSpan ()

Private Member Functions
void	initBack ()
void	initReps ()

Private Attributes
GlobalPoint	m_backPoint
CornersVec	m_corners
float	m_dEta
float	m_dPhi
const CCGFloat *	m_parms
GlobalPoint	m_refPoint
RhoEtaPhi	m_rep
std::array< RhoEtaPhi, k_cornerSize >	m_repCorners

Detailed Description

Abstract base class for an individual cell's geometry.

The base class declares a pure virtual function and also writes a definition (body) to force conscious acceptance of default behaviour.

If a derived class doesn't choose to override a normal virtual, it just inherits the base version's behaviour by default. If you want to provide a default behaviour but not let derived classes just inherit it "silently" like this, you can make it pure virtual but still provide a default that the derived class author has to call deliberately if he wants it:

class B {
  public:
      virtual bool f() = 0;
  };

  bool B::f() {
      return true;  // this is a good default, but
  }                 // shouldn't be used blindly

  class D : public B {
  public:
      bool f() {
          return B::f(); // if D wants the default
      }                  // behaviour, it has to say so
 };

Author

J. Mans, P. Meridiani

Definition at line 53 of file CaloCellGeometry.h.

Member Typedef Documentation

typedef float CaloCellGeometry::CCGFloat

Definition at line 57 of file CaloCellGeometry.h.

typedef EZMgrFL< GlobalPoint > CaloCellGeometry::CornersMgr

Definition at line 63 of file CaloCellGeometry.h.

typedef EZArrayFL< GlobalPoint > CaloCellGeometry::CornersVec

Definition at line 62 of file CaloCellGeometry.h.

typedef EZMgrFL< CCGFloat > CaloCellGeometry::ParMgr

Definition at line 67 of file CaloCellGeometry.h.

typedef EZArrayFL<CCGFloat> CaloCellGeometry::ParVec

Definition at line 65 of file CaloCellGeometry.h.

typedef std::vector<ParVec> CaloCellGeometry::ParVecVec

Definition at line 66 of file CaloCellGeometry.h.

typedef HepGeom::Point3D<CCGFloat> CaloCellGeometry::Pt3D

Definition at line 59 of file CaloCellGeometry.h.

typedef std::vector<Pt3D> CaloCellGeometry::Pt3DVec

Definition at line 60 of file CaloCellGeometry.h.

using CaloCellGeometry::RepCorners = std::array<RhoEtaPhi,k_cornerSize>

Definition at line 71 of file CaloCellGeometry.h.

typedef HepGeom::Transform3D CaloCellGeometry::Tr3D

Definition at line 58 of file CaloCellGeometry.h.

Constructor & Destructor Documentation

CaloCellGeometry::~CaloCellGeometry ( )

virtual

Definition at line 26 of file CaloCellGeometry.cc.

{}

CaloCellGeometry::CaloCellGeometry ( CornersVec::const_reference  gp, CornersMgr *  mgr, const CCGFloat *  par  )

protected

Definition at line 31 of file CaloCellGeometry.cc.

                                                      :
   m_refPoint ( gp  ),
   m_corners  ( mgr ),
   m_parms    ( par ),
   m_rep(gp.perp(),gp.eta(),gp.barePhi())
{}

CaloCellGeometry::CaloCellGeometry ( const CornersVec &  cv, const CCGFloat *  par  )

protected

Definition at line 40 of file CaloCellGeometry.cc.

                                            : 
   m_refPoint ( 0.25*( cv[0].x() + cv[1].x() + cv[2].x() + cv[3].x() ),
        0.25*( cv[0].y() + cv[1].y() + cv[2].y() + cv[3].y() ),
        0.25*( cv[0].z() + cv[1].z() + cv[2].z() + cv[3].z() )  ), 
   m_corners  ( cv ),
   m_parms    ( par ),
   m_rep( m_refPoint.perp(), m_refPoint.eta(), m_refPoint.barePhi())
{}

CaloCellGeometry::CaloCellGeometry ( void  )

protected

Definition at line 19 of file CaloCellGeometry.cc.

                                   :
   m_refPoint ( 0., 0., 0. ),
   m_corners  (  ) ,
   m_parms    ( (CCGFloat*) 0 )
{}

Member Function Documentation

const float * CaloCellGeometry::checkParmPtr ( const std::vector< CCGFloat > &  vd, CaloCellGeometry::ParVecVec &  pvv  )

static

Definition at line 119 of file CaloCellGeometry.cc.

References assert(), EZArrayFL< T >::begin(), cuy::ii, j, EZArrayFL< T >::size(), and findQualityFiles::v.

Referenced by getParmPtr().

{
   const float* pP ( 0 ) ;

   for( unsigned int ii ( 0 ) ; ii != pvv.size() ; ++ii )
   {
      const ParVec& v ( pvv[ii] ) ;
      assert( v.size() == vv.size() ) ;

      bool same ( true ) ;
      for( unsigned int j ( 0 ) ; j != vv.size() ; ++j )
      {
     same = same && ( fabs( vv[j] - v[j] )<1.e-6 ) ;
     if( !same ) break ;
      }
      if( same )
      {
     pP = &(*v.begin()) ;
     break ;
      }
   }
   return pP ;
}

bool CaloCellGeometry::emptyCorners ( ) const

inline

Definition at line 99 of file CaloCellGeometry.h.

References EZArrayFL< T >::empty(), and m_corners.

Referenced by operator<<().

{ return m_corners.empty() ;}

float CaloCellGeometry::etaPos ( ) const

inline

Definition at line 89 of file CaloCellGeometry.h.

References RhoEtaPhi::eta(), and m_rep.

{ return m_rep.eta();}

float CaloCellGeometry::etaSpan ( ) const

inline

Definition at line 92 of file CaloCellGeometry.h.

References m_dEta.

{ return m_dEta;}

const GlobalPoint& CaloCellGeometry::getBackPoint ( ) const

inline

Definition at line 85 of file CaloCellGeometry.h.

References m_backPoint.

{return m_backPoint;} 

CornersVec const& CaloCellGeometry::getCorners ( ) const

inline

Returns the corner points of this cell's volume.

Definition at line 79 of file CaloCellGeometry.h.

References assert(), m_corners, and EZArrayFL< T >::uninitialized().

Referenced by reco::PFRecHit::getCornersXYZ(), HCaloDetIdAssociator::getDetIdPoints(), getTransform(), initBack(), initReps(), initSpan(), inside(), operator<<(), and EcalGeomPhiSymHelper::setup().

{ assert(not m_corners.uninitialized()); return m_corners; }

RepCorners const& CaloCellGeometry::getCornersREP ( ) const

inline

Definition at line 80 of file CaloCellGeometry.h.

References m_repCorners.

Referenced by reco::PFRecHit::getCornersREP().

{  return m_repCorners;}

const float * CaloCellGeometry::getParmPtr ( const std::vector< CCGFloat > &  vd, CaloCellGeometry::ParMgr *  mgr, CaloCellGeometry::ParVecVec &  pvv  )

static

Definition at line 146 of file CaloCellGeometry.cc.

References checkParmPtr(), and i.

Referenced by CaloTowerHardcodeGeometryLoader::makeCell(), EcalTBHodoscopeGeometryLoaderFromDDD::makeGeometry(), and CaloGeometryDBEP< T, U >::produceAligned().

{
   const float* pP ( checkParmPtr( vv, pvv ) ) ;

   if( 0 == pP )
   {
      pvv.push_back( ParVec( mgr ) ) ;
      ParVec& back ( pvv.back() ) ;
      for( unsigned int i ( 0 ) ; i != vv.size() ; ++i )
      {
     back[i] = vv[i] ;
      }
      pP = &(*back.begin()) ;
   }
   return pP ;
}

const GlobalPoint& CaloCellGeometry::getPosition ( ) const

inline

Returns the position of reference for this cell.

Definition at line 84 of file CaloCellGeometry.h.

References m_refPoint.

Referenced by PreshowerHitMaker::addHit(), CaloTowerConstituentsMapBuilder::assignEEtoHE(), SuperClusterShapeAlgo::Calculate_Covariances(), DistanceToCell::DistanceToCell(), CaloGeometry::getPosition(), HCaloDetIdAssociator::getPosition(), CaloTowerGeometry::getSummary(), getTransform(), IdealZDCTrapezoid::initCorners(), PreshowerStrip::initCorners(), IdealZPrism::initCorners(), IdealCastorTrapezoid::initCorners(), PhotonFixCMS::initialise(), PreshowerPhiClusterAlgo::makeOneCluster(), DistanceToCell::operator()(), operator<<(), reco::PFRecHit::position(), and EcalGeomPhiSymHelper::setup().

{return m_refPoint;}

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

virtual

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

Reimplemented in TruncatedPyramid, FlatTrd, and PreshowerStrip.

Definition at line 70 of file CaloCellGeometry.cc.

References angle(), assert(), getCorners(), getPosition(), i, mag(), AlCaHLTBitMon_ParallelJobs::p, param(), csvLumiCalc::unit, vocalCorners(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by CaloTowerGeometry::getSummary().

{
   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( 0 != param() ) ;
   Pt3DVec lc ( 8, Pt3D(0,0,0) ) ;
   vocalCorners( lc, param(), lFront ) ;

   DPt3D dlFront ( lFront.x(), lFront.y(), lFront.z() ) ;

   const Pt3D lBack ( 0.25*(lc[4]+lc[5]+lc[6]+lc[7]) ) ;
   const DPt3D dlBack ( lBack.x(), lBack.y(), lBack.z() ) ;

   const Pt3D dlOne ( lc[0].x(), lc[0].y(), lc[0].z() ) ;

   const CornersVec& cor ( getCorners() ) ;
   DPt3DVec dkor ( 8, DPt3D(0,0,0) ) ;
   for( unsigned int i ( 0 ) ; i != 8 ; ++i )
   {
      dkor[i] = DPt3D ( cor[i].x(), cor[i].y(), cor[i].z() ) ;
   }

   DPt3D  dgBack ( 0.25*( dkor[4]+dkor[5]+dkor[6]+dkor[7] ) ) ;

   const DVec3D dgAxis ( (dgBack-dgFront).unit() ) ;

   dgBack = ( dgFront + (dlBack-dlFront).mag()*dgAxis ) ;
   const DPt3D  dgOneT ( dgFront + ( dlOne - dlFront ).mag()*( dkor[0] - dgFront ).unit() ) ;

   const double dlangle ( ( dlBack - dlFront).angle( dlOne - dlFront ) ) ;
   const double dgangle ( ( dgBack - dgFront).angle( dgOneT- dgFront ) ) ;
   const double ddangle ( dlangle - dgangle ) ;

   const DPlane3D dgPl ( dgFront, dgOneT, dgBack ) ;
   const DPt3D    dp2  ( dgFront + dgPl.normal().unit() ) ;

   const DPt3D  dgOne ( dgFront + HepGeom::Rotate3D( -ddangle, dgFront, dp2 )*
              DVec3D ( dgOneT - dgFront ) ) ;

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

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

void CaloCellGeometry::initBack ( )

inlineprivate

Definition at line 143 of file CaloCellGeometry.h.

References cuy::cv, getCorners(), m_backPoint, x, y, and z.

Referenced by initSpan().

                 {
    // from CaloTower code
    CornersVec const & cv = getCorners();
    m_backPoint = GlobalPoint(0.25 * (cv[4].x() + cv[5].x() + cv[6].x() + cv[7].x()),
                              0.25 * (cv[4].y() + cv[5].y() + cv[6].y() + cv[7].y()),
                              0.25 * (cv[4].z() + cv[5].z() + cv[6].z() + cv[7].z()));   
  }

virtual void CaloCellGeometry::initCorners ( CornersVec &  )

protectedpure virtual

Implemented in TruncatedPyramid, FlatTrd, IdealCastorTrapezoid, IdealZPrism, IdealObliquePrism, PreshowerStrip, and IdealZDCTrapezoid.

Referenced by initSpan().

void CaloCellGeometry::initReps ( )

inlineprivate

Definition at line 150 of file CaloCellGeometry.h.

References getCorners(), i, k_cornerSize, and m_repCorners.

Referenced by initSpan().

                  {
    for (auto i=0U;i<k_cornerSize; ++i) m_repCorners[i]= {getCorners()[i].perp(), getCorners()[i].eta(), getCorners()[i].barePhi()}; 
      
  }

void CaloCellGeometry::initSpan ( )

inlineprotected

Definition at line 131 of file CaloCellGeometry.h.

References funct::abs(), eta, getCorners(), initBack(), initCorners(), initReps(), m_corners, m_dEta, m_dPhi, and phi.

Referenced by IdealCastorTrapezoid::IdealCastorTrapezoid(), IdealZDCTrapezoid::IdealZDCTrapezoid(), IdealZPrism::IdealZPrism(), and PreshowerStrip::PreshowerStrip().

                  {
     initCorners(m_corners);
     m_dEta = std::abs(getCorners()[0].eta()-
                      getCorners()[2].eta());
     m_dPhi = std::abs(getCorners()[0].phi() -
                      getCorners()[2].phi());
     initBack();
     initReps();
  }

bool CaloCellGeometry::inside

(

const GlobalPoint &

point

)

const

Returns true if the specified point is inside this cell.

Definition at line 167 of file CaloCellGeometry.cc.

References EE, getCorners(), i, AlCaHLTBitMon_ParallelJobs::p, funct::true, x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by CastorGeometry::getClosestCell(), CaloDetIdAssociator::insideElement(), and HCaloDetIdAssociator::insideElement().

{
   bool ans ( false ) ;
   const Pt3D  p ( point.x(), point.y(), point.z() ) ;
   const CornersVec& cog ( getCorners() ) ;
   Pt3D  co[8] ;
   for( unsigned int i ( 0 ) ; i != 8 ; ++i )
   {
      co[i] = Pt3D ( cog[i].x(), cog[i].y(), cog[i].z() ) ;
   }

   const Plane3D  AA ( co[0], co[1], co[2] ) ; // z<0
   const Plane3D  BB ( co[6], co[5], co[4] ) ; // z>0

   if( AA.distance(p)*BB.distance(p) >= 0 )
   {
     const Plane3D  CC ( co[0], co[4], co[5] ) ; // x<0
     const Plane3D  DD ( co[2], co[6], co[7] ) ; // x>0
     if( CC.distance(p)*DD.distance(p) >= 0 )
     {
       const Plane3D  EE ( co[3], co[7], co[4] ) ; // y<0
       const Plane3D  FF ( co[1], co[5], co[6] ) ; // y>0
       if( EE.distance(p)*FF.distance(p) >= 0 )
       {
       ans = true ;
       }
      }
   }
   return ans ;
}

const CCGFloat* CaloCellGeometry::param ( ) const

inline

Definition at line 101 of file CaloCellGeometry.h.

References m_parms.

Referenced by IdealZDCTrapezoid::an(), IdealCastorTrapezoid::an(), ZdcGeometry::cellGeomPtr(), CastorGeometry::cellGeomPtr(), CaloTowerGeometry::cellGeomPtr(), IdealZPrism::dEta(), IdealCastorTrapezoid::dh(), IdealZPrism::dPhi(), IdealCastorTrapezoid::dR(), IdealZDCTrapezoid::dx(), PreshowerStrip::dx(), IdealCastorTrapezoid::dxh(), IdealCastorTrapezoid::dxl(), IdealZDCTrapezoid::dy(), PreshowerStrip::dy(), IdealZDCTrapezoid::dz(), PreshowerStrip::dz(), IdealZPrism::dz(), IdealCastorTrapezoid::dz(), IdealZPrism::eta(), CaloTowerGeometry::getSummary(), getTransform(), IdealZDCTrapezoid::initCorners(), IdealCastorTrapezoid::initCorners(), operator<<(), PreshowerStrip::tilt(), and IdealZPrism::z().

{ return m_parms ;}

float CaloCellGeometry::phiPos ( ) const

inline

Definition at line 90 of file CaloCellGeometry.h.

References m_rep, and RhoEtaPhi::phi().

{ return m_rep.phi();}

float CaloCellGeometry::phiSpan ( ) const

inline

Definition at line 93 of file CaloCellGeometry.h.

References m_dPhi.

{ return m_dPhi;}

RhoEtaPhi const& CaloCellGeometry::repPos ( ) const

inline

Definition at line 87 of file CaloCellGeometry.h.

References m_rep.

Referenced by reco::PFRecHit::positionREP().

{ return m_rep;}

float CaloCellGeometry::rhoPos ( ) const

inline

Definition at line 88 of file CaloCellGeometry.h.

References m_rep, and RhoEtaPhi::rho().

{ return m_rep.rho();}

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

pure virtual

Implemented in IdealCastorTrapezoid, IdealZPrism, TruncatedPyramid, FlatTrd, IdealObliquePrism, PreshowerStrip, and IdealZDCTrapezoid.

Referenced by getTransform().

Member Data Documentation

constexpr unsigned int CaloCellGeometry::k_cornerSize = 8

static

Definition at line 69 of file CaloCellGeometry.h.

Referenced by initReps().

const float CaloCellGeometry::k_ScaleFromDDDtoGeant

static

Definition at line 74 of file CaloCellGeometry.h.

Referenced by EcalTBHodoscopeGeometryLoaderFromDDD::makeGeometry().

GlobalPoint CaloCellGeometry::m_backPoint

private

Definition at line 157 of file CaloCellGeometry.h.

Referenced by getBackPoint(), and initBack().

CornersVec CaloCellGeometry::m_corners

private

Definition at line 158 of file CaloCellGeometry.h.

Referenced by emptyCorners(), getCorners(), and initSpan().

float CaloCellGeometry::m_dEta

private

Definition at line 161 of file CaloCellGeometry.h.

Referenced by etaSpan(), and initSpan().

float CaloCellGeometry::m_dPhi

private

Definition at line 162 of file CaloCellGeometry.h.

Referenced by initSpan(), and phiSpan().

const CCGFloat* CaloCellGeometry::m_parms

private

Definition at line 159 of file CaloCellGeometry.h.

Referenced by param().

GlobalPoint CaloCellGeometry::m_refPoint

private

Definition at line 156 of file CaloCellGeometry.h.

Referenced by getPosition().

RhoEtaPhi CaloCellGeometry::m_rep

private

Definition at line 160 of file CaloCellGeometry.h.

Referenced by etaPos(), phiPos(), repPos(), and rhoPos().

std::array<RhoEtaPhi,k_cornerSize> CaloCellGeometry::m_repCorners

private

Definition at line 163 of file CaloCellGeometry.h.

Referenced by getCornersREP(), and initReps().