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
enum	CornersSize { k_cornerSize = 8 }
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
typedef HepGeom::Transform3D	Tr3D

Public Member Functions
const GlobalVector &	axis () const
virtual const CornersVec &	getCorners () const
	Returns the corner points of this cell's volume. More...
CCGFloat	getPhiAxis () const
const GlobalPoint	getPosition (CCGFloat depth) const
CCGFloat	getThetaAxis () const
virtual void	getTransform (Tr3D &tr, Pt3DVec *lptr) const
	--------— only needed by specific utility; overloaded when needed -— More...
TruncatedPyramid &	operator= (const TruncatedPyramid &tr)
	TruncatedPyramid (void)
	TruncatedPyramid (const TruncatedPyramid &tr)
	TruncatedPyramid (const CornersMgr *cMgr, const GlobalPoint &fCtr, const GlobalPoint &bCtr, const GlobalPoint &cor1, const CCGFloat *parV)
	TruncatedPyramid (const CornersVec &corn, const CCGFloat *par)
virtual void	vocalCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const
virtual	~TruncatedPyramid ()
Public Member Functions inherited from CaloCellGeometry
bool	emptyCorners () const
const GlobalPoint &	getPosition () const
	Returns the position of reference for this cell. More...
bool	inside (const GlobalPoint &point) const
	Returns true if the specified point is inside this cell. More...
const CCGFloat *	param () const
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
GlobalVector	makeAxis (void)

Private Attributes
GlobalVector	m_axis
Pt3D	m_corOne

Additional Inherited Members
Static Public Attributes inherited from CaloCellGeometry
static const CCGFloat	k_ScaleFromDDDtoGeant
Protected Member Functions inherited from CaloCellGeometry
	CaloCellGeometry (CornersVec::const_reference gp, const CornersMgr *mgr, const CCGFloat *par)
	CaloCellGeometry (const CornersVec &cv, const CCGFloat *par)
	CaloCellGeometry (void)
	CaloCellGeometry (const CaloCellGeometry &cell)
CaloCellGeometry &	operator= (const CaloCellGeometry &cell)
CornersVec &	setCorners () const

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

typedef CaloCellGeometry::Pt3D TruncatedPyramid::Pt3D

Definition at line 24 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Pt3DVec TruncatedPyramid::Pt3DVec

Definition at line 25 of file TruncatedPyramid.h.

typedef CaloCellGeometry::Tr3D TruncatedPyramid::Tr3D

Definition at line 26 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	(	const CornersMgr *	cMgr,
		const GlobalPoint &	fCtr,
		const GlobalPoint &	bCtr,
		const GlobalPoint &	cor1,
		const CCGFloat *	parV
	)

Definition at line 44 of file TruncatedPyramid.cc.

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

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

Definition at line 54 of file TruncatedPyramid.cc.

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

TruncatedPyramid::~TruncatedPyramid ( )

virtual

Definition at line 61 of file TruncatedPyramid.cc.

{}

Member Function Documentation

const GlobalVector & TruncatedPyramid::axis

(

)

const

Definition at line 83 of file TruncatedPyramid.cc.

References m_axis.

Referenced by getTransform().

{ 
   return m_axis ; 
}

const GlobalPoint TruncatedPyramid::backCtr ( void  ) const

private

Definition at line 104 of file TruncatedPyramid.cc.

References getCorners(), x, detailsBasic3DVector::y, and detailsBasic3DVector::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 319 of file TruncatedPyramid.cc.

References filterCSVwithJSON::copy, cross(), dot(), i, localCorners(), AlCaHLTBitMon_ParallelJobs::p, truncPyr::refl(), tmp, x, detailsBasic3DVector::y, and detailsBasic3DVector::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 Vec3D  x ( 1, 0, 0 ) ;
   static const Vec3D  y ( 0, 1, 0 ) ;
   static const Vec3D  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() ) ;
   }
}

const CaloCellGeometry::CornersVec & TruncatedPyramid::getCorners ( ) const

virtual

Returns the corner points of this cell's volume.

Implements CaloCellGeometry.

Definition at line 186 of file TruncatedPyramid.cc.

References CaloCellGeometry::getCorners(), getTransform(), i, CaloCellGeometry::setCorners(), and EZArrayFL< T >::uninitialized().

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

{
   const CornersVec& co ( CaloCellGeometry::getCorners() ) ;
   if( co.uninitialized() ) 
   {
      CornersVec& corners ( setCorners() ) ;

      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() ) ;
      }
   }
   return co ;
}

CCGFloat TruncatedPyramid::getPhiAxis

(

)

const

Definition at line 77 of file TruncatedPyramid.cc.

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

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

{
   return m_axis.phi() ; 
} 

const GlobalPoint TruncatedPyramid::getPosition

(

CCGFloat

depth

)

const

Definition at line 65 of file TruncatedPyramid.cc.

References CaloCellGeometry::getPosition(), and m_axis.

Referenced by EcalSimple2007H4TBAnalyzer::analyze(), PFRecHitProducerECAL::createEcalRecHit(), PFRecHitProducerECAL::findEcalRecHitGeometry(), EcalClusterLocal::localCoordsEB(), ggPFClusters::localCoordsEB(), EcalClusterLocal::localCoordsEE(), and ggPFClusters::localCoordsEE().

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

CCGFloat TruncatedPyramid::getThetaAxis

(

)

const

Definition at line 71 of file TruncatedPyramid.cc.

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

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

{ 
   return m_axis.theta() ; 
} 

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

virtual

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

Reimplemented from CaloCellGeometry.

Definition at line 115 of file TruncatedPyramid.cc.

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

Referenced by getCorners().

{
   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( 0 != 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 Vec3D 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( 0 != lptr ) (*lptr) = lc ;
}

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

static

Definition at line 250 of file TruncatedPyramid.cc.

References funct::cos(), funct::sin(), and funct::tan().

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

{
   assert( 0 != 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 216 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 291 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 97 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, m_corOne, and CaloCellGeometry::operator=().

{
   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

virtual

Implements CaloCellGeometry.

Definition at line 89 of file TruncatedPyramid.cc.

References localCorners().

{ 
   localCorners( vec, pv, ref ) ; 
}

Member Data Documentation

GlobalVector TruncatedPyramid::m_axis

private

Definition at line 84 of file TruncatedPyramid.h.

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

Pt3D TruncatedPyramid::m_corOne

private

Definition at line 85 of file TruncatedPyramid.h.

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