---

TruncatedPyramid Class Reference

A base class to handle the particular shape of Ecal Xtals. More...

#include <Geometry/CaloGeometry/interface/TruncatedPyramid.h>

Inheritance diagram for TruncatedPyramid:

List of all members.

Public Types
typedef std::vector< HepPlane3D >	BoundaryVec
Public Member Functions
const GlobalVector &	axis () const
virtual const CornersVec &	getCorners () const
float	getPhiAxis () const
const GlobalPoint	getPosition (float depth) const
	Position corresponding to the center of the front face at a certain depth (default is zero) along the crystal axis.
float	getThetaAxis () const
virtual bool	inside (const GlobalPoint &point) const
	TruncatedPyramid (const CornersVec &corn)
virtual	~TruncatedPyramid ()
Static Public Member Functions
static void	createCorners (const std::vector< double > &pv, const HepTransform3D &tr, CornersVec &co)
Private Member Functions
GlobalVector	axis ()
	print out the element, with an optional string prefix, maybe OVAL identifier
const GlobalPoint	backCtr () const
Private Attributes
GlobalVector	m_axis

---

Detailed Description

A base class to handle the particular shape of Ecal Xtals.

Taken from ORCA Calorimetry Code

Definition at line 20 of file TruncatedPyramid.h.

---

Member Typedef Documentation

typedef std::vector< HepPlane3D > TruncatedPyramid::BoundaryVec

Definition at line 24 of file TruncatedPyramid.h.

---

Constructor & Destructor Documentation

TruncatedPyramid::TruncatedPyramid

(

const CornersVec &

corn

)

[inline]

Definition at line 26 of file TruncatedPyramid.h.

00026                                                   :
00027          CaloCellGeometry ( corn ) , 
00028          m_axis           ( axis() ) {} 
/*,

virtual TruncatedPyramid::~TruncatedPyramid

(

)

[inline, virtual]

Definition at line 32 of file TruncatedPyramid.h.

{}// delete m_bou ; }

---

Member Function Documentation

GlobalVector TruncatedPyramid::axis

(

)

[inline, private]

print out the element, with an optional string prefix, maybe OVAL identifier

Definition at line 64 of file TruncatedPyramid.h.

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

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

const GlobalVector& TruncatedPyramid::axis

(

)

const [inline]

Definition at line 52 of file TruncatedPyramid.h.

References m_axis.

{ return m_axis ; }

const GlobalPoint TruncatedPyramid::backCtr

(

)

const [inline, private]

Definition at line 70 of file TruncatedPyramid.h.

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

Referenced by axis().

      {
         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< double > &	pv,
		const HepTransform3D &	tr,
		CaloCellGeometry::CornersVec &	co
	)			[static]

Definition at line 12 of file TruncatedPyramid.cc.

References a1, a2, b1, b2, edmNew::copy(), funct::cos(), dot(), h1, h2, i, CaloCellGeometry::k_ScaleFromDDDtoGeant, off, p, funct::sin(), funct::tan(), x, y, and z.

{
   assert( 11 == pv.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.

   static const float kf ( CaloCellGeometry::k_ScaleFromDDDtoGeant ) ;

   const double dz ( kf*pv[0] ) ;
   const double th (    pv[1] ) ;
   const double ph (    pv[2] ) ;
   const double h1 ( kf*pv[3] ) ;
   const double b1 ( kf*pv[4] ) ;
   const double t1 ( kf*pv[5] ) ;
   const double a1 (    pv[6] ) ;
   const double h2 ( kf*pv[7] ) ;
   const double b2 ( kf*pv[8] ) ;
   const double t2 ( kf*pv[9] ) ;
   const double a2 (    pv[10]) ;

   std::vector<HepPoint3D> to ( 8, HepPoint3D(0,0,0) ) ;
   std::vector<HepPoint3D> ko ( 8, HepPoint3D(0,0,0) ) ;
  
   const double ta1 ( tan( a1 ) ) ; // lower plane
   const double ta2 ( tan( a2 ) ) ; // upper plane

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

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

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

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

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

   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 HepPoint3D& p ( ko[i] ) ;
      co[ i ] = GlobalPoint( p.x(), p.y(), p.z() ) ;
   }
}

virtual const CornersVec& TruncatedPyramid::getCorners

(

)

const [inline, virtual]

Implements CaloCellGeometry.

Definition at line 44 of file TruncatedPyramid.h.

References CaloCellGeometry::getCorners().

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

{ return CaloCellGeometry::getCorners() ; }

float TruncatedPyramid::getPhiAxis

(

)

const [inline]

Definition at line 50 of file TruncatedPyramid.h.

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

Referenced by operator<<().

{ return m_axis.phi() ; } 

const GlobalPoint TruncatedPyramid::getPosition

(

float

depth

)

const [inline]

Position corresponding to the center of the front face at a certain depth (default is zero) along the crystal axis.

If "depth" is <=0, the nomial position of the cell is returned (center of the front face).

Definition at line 41 of file TruncatedPyramid.h.

References CaloCellGeometry::getPosition(), and m_axis.

Referenced by PFRecHitProducerECAL::createEcalRecHit(), and PFRecHitProducerECAL::findEcalRecHitGeometry().

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

float TruncatedPyramid::getThetaAxis

(

)

const [inline]

Definition at line 47 of file TruncatedPyramid.h.

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

Referenced by operator<<().

{ return m_axis.theta() ; } 

bool TruncatedPyramid::inside

(

const GlobalPoint &

point

)

const [virtual]

Implements CaloCellGeometry.

Definition at line 117 of file TruncatedPyramid.cc.

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

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

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

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

---

Member Data Documentation

GlobalVector TruncatedPyramid::m_axis [private]

Definition at line 80 of file TruncatedPyramid.h.

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

---

The documentation for this class was generated from the following files:

• Geometry/CaloGeometry/interface/TruncatedPyramid.h
• Geometry/CaloGeometry/src/TruncatedPyramid.cc

---