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

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 (const Pt3D &) 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 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)

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
Static Public Attributes inherited from CaloCellGeometry
static constexpr unsigned int	k_cornerSize = 8

static const CCGFloat	k_ScaleFromDDDtoGeant

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 22 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Pt3D TruncatedPyramid::Pt3D

Definition at line 23 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Pt3DVec TruncatedPyramid::Pt3DVec

Definition at line 24 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Tr3D TruncatedPyramid::Tr3D

Definition at line 25 of file TruncatedPyramid.h.

Constructor & Destructor Documentation

TruncatedPyramid::TruncatedPyramid ( void )

Definition at line 20 of file TruncatedPyramid.cc.

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

TruncatedPyramid::TruncatedPyramid ( const TruncatedPyramid & tr )

Definition at line 26 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 44 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 54 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 61 of file TruncatedPyramid.cc.

62 {}

Member Function Documentation

const GlobalVector & TruncatedPyramid::axis ( ) const

Definition at line 82 of file TruncatedPyramid.cc.

References m_axis.

Referenced by getTransform().

 { 
    return m_axis ; 
 }

const GlobalPoint TruncatedPyramid::backCtr ( void ) const

private

Definition at line 103 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 314 of file TruncatedPyramid.cc.

References popcon2dropbox::copy(), cross(), dot(), PVValHelper::dz, mps_fire::i, localCorners(), AlCaHLTBitMon_ParallelJobs::p, truncPyr::refl(), 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 76 of file TruncatedPyramid.cc.

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

Referenced by EcalClusterLocal::localCoordsEB(), EcalClusterLocal::localCoordsEE(), and operator<<().

 {
    return m_axis.phi() ; 
 } 

GlobalPoint TruncatedPyramid::getPosition ( CCGFloat depth ) const

overridevirtual

Reimplemented from CaloCellGeometry.

Definition at line 64 of file TruncatedPyramid.cc.

References CaloCellGeometry::getPosition(), and m_axis.

Referenced by EcalClusterLocal::localCoordsEB(), EcalClusterLocal::localCoordsEE(), and PFPhotonClusters::SetSeed().

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

CCGFloat TruncatedPyramid::getThetaAxis ( ) const

Definition at line 70 of file TruncatedPyramid.cc.

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

Referenced by EcalClusterLocal::localCoordsEB(), EcalClusterLocal::localCoordsEE(), 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 114 of file TruncatedPyramid.cc.

References angle(), axis(), PVValHelper::dz, MillePedeFileConverter_cfg::e, CaloCellGeometry::getPosition(), localCorners(), localCornersReflection(), localCornersSwap(), m_corOne, mag(), AlCaHLTBitMon_ParallelJobs::p, CaloCellGeometry::param(), csvLumiCalc::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 185 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 245 of file TruncatedPyramid.cc.

References funct::cos(), PVValHelper::dz, funct::sin(), reco::t2, 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[0] ) ;
    const CCGFloat th ( pv[1] ) ;
    const CCGFloat ph ( pv[2] ) ;
    const CCGFloat h1 ( pv[3] ) ;
    const CCGFloat b1 ( pv[4] ) ;
    const CCGFloat t1 ( pv[5] ) ;
    const CCGFloat a1 ( pv[6] ) ;
    const CCGFloat h2 ( pv[7] ) ;
    const CCGFloat b2 ( pv[8] ) ;
    const CCGFloat t2 ( pv[9] ) ;
    const CCGFloat a2 ( pv[10]) ;
   
    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 211 of file TruncatedPyramid.cc.

References localCorners(), truncPyr::refl(), and 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 286 of file TruncatedPyramid.cc.

References localCorners(), and 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 96 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 33 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 88 of file TruncatedPyramid.cc.

References localCorners().

 { 
    localCorners( vec, pv, ref ) ; 
 }

Member Data Documentation

GlobalVector TruncatedPyramid::m_axis

private

Definition at line 83 of file TruncatedPyramid.h.

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

Pt3D TruncatedPyramid::m_corOne

private

Definition at line 84 of file TruncatedPyramid.h.

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

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation