FlatTrd Class Reference

A base class to handle the particular shape of HGCal volumes. More...

#include <FlatTrd.h>

Inheritance diagram for FlatTrd:

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 float	etaPos () const
	FlatTrd (void)
	FlatTrd (const FlatTrd &tr)
	FlatTrd (CornersMgr *cMgr, const GlobalPoint &fCtr, const GlobalPoint &bCtr, const GlobalPoint &cor1, const CCGFloat *parV)
	FlatTrd (const CornersVec &corn, const CCGFloat *par)
	FlatTrd (const FlatTrd &tr, const Pt3D &local)
Pt3D	getLocal (const GlobalPoint &global) const
CCGFloat	getPhiAxis () const
virtual const GlobalPoint &	getPosition () const
GlobalPoint	getPosition (const Pt3D &local) const
CCGFloat	getThetaAxis () const
virtual void	getTransform (Tr3D &tr, Pt3DVec *lptr) const
	--------— only needed by specific utility; overloaded when needed -— More...
FlatTrd &	operator= (const FlatTrd &tr)
virtual float	phiPos () const
virtual void	vocalCorners (Pt3DVec &vec, const CCGFloat *pv, Pt3D &ref) const
virtual	~FlatTrd ()
Public Member Functions inherited from CaloCellGeometry
bool	emptyCorners () const
float	etaPos () const
float	etaSpan () const
const GlobalPoint &	getBackPoint () const
const CornersVec &	getCorners () const
	Returns the corner points of this cell's volume. More...
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
float	phiPos () const
float	phiSpan () 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 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
virtual void	initCorners (CornersVec &) override
GlobalVector	makeAxis (void)

Private Attributes
GlobalVector	m_axis
Pt3D	m_corOne
GlobalPoint	m_global
Pt3D	m_local
Tr3D	m_tr

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, CornersMgr *mgr, const CCGFloat *par)
	CaloCellGeometry (const CornersVec &cv, const CCGFloat *par)
	CaloCellGeometry (void)
void	initSpan ()

Detailed Description

A base class to handle the particular shape of HGCal volumes.

Definition at line 19 of file FlatTrd.h.

Member Typedef Documentation

typedef CaloCellGeometry::CCGFloat FlatTrd::CCGFloat

Definition at line 22 of file FlatTrd.h.

typedef CaloCellGeometry::Pt3D FlatTrd::Pt3D

Definition at line 23 of file FlatTrd.h.

typedef CaloCellGeometry::Pt3DVec FlatTrd::Pt3DVec

Definition at line 24 of file FlatTrd.h.

typedef CaloCellGeometry::Tr3D FlatTrd::Tr3D

Definition at line 25 of file FlatTrd.h.

Constructor & Destructor Documentation

FlatTrd::FlatTrd

(

void

)

Definition at line 22 of file FlatTrd.cc.

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

FlatTrd::FlatTrd

(

const FlatTrd &

tr

)

Definition at line 27 of file FlatTrd.cc.

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

FlatTrd::FlatTrd	(	CornersMgr *	cMgr,
		const GlobalPoint &	fCtr,
		const GlobalPoint &	bCtr,
		const GlobalPoint &	cor1,
		const CCGFloat *	parV
	)

Definition at line 48 of file FlatTrd.cc.

References gather_cfg::cout, etaPos(), getTransform(), glb, m_global, m_local, m_tr, and phiPos().

                                    :
  CaloCellGeometry ( fCtr, cMgr, parV ) ,
  m_axis           ( ( bCtr - fCtr ).unit() ) ,
  m_corOne         ( cor1.x(), cor1.y(), cor1.z() ),
  m_local          (0., 0., 0.) {
  getTransform(m_tr,0);
  Pt3D glb = m_tr*m_local;
  m_global = GlobalPoint(glb.x(),glb.y(),glb.z());
#ifdef DebugLog
  std::cout << "FlatTrd: Local " << m_local << " Global " << glb << " eta "
        << etaPos() << " phi " << phiPos() << " Translation "
        << m_tr.getTranslation() << " and rotation " << m_tr.getRotation();
#endif
} 

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

Definition at line 68 of file FlatTrd.cc.

References gather_cfg::cout, etaPos(), getTransform(), glb, m_axis, m_global, m_local, m_tr, makeAxis(), and phiPos().

                                     :
  CaloCellGeometry ( corn, par  ) , 
  m_corOne         ( corn[0].x(), corn[0].y(), corn[0].z() ),
  m_local          (0., 0., 0.) {
  getTransform(m_tr,0);
  m_axis   = makeAxis();
  Pt3D glb = m_tr*m_local;
  m_global = GlobalPoint(glb.x(),glb.y(),glb.z());
#ifdef DebugLog
  std::cout << "FlatTrd: Local " << m_local << " Global " << glb << " eta "
        << etaPos() << " phi " << phiPos() << " Translation "
        << m_tr.getTranslation() << " and rotation " << m_tr.getRotation();
#endif
} 

FlatTrd::FlatTrd	(	const FlatTrd &	tr,
		const Pt3D &	local
	)

Definition at line 84 of file FlatTrd.cc.

References gather_cfg::cout, etaPos(), glb, m_global, m_local, m_tr, and phiPos().

                                                        : 
  CaloCellGeometry( tr ),  m_local(local) {
  *this = tr;
  Pt3D glb = m_tr*m_local;
  m_global = GlobalPoint(glb.x(),glb.y(),glb.z());
#ifdef DebugLog
  std::cout << "FlatTrd: Local " << m_local << " Global " << glb << " eta "
        << etaPos() << " phi " << phiPos() << " Translation "
        << m_tr.getTranslation() << " and rotation " << m_tr.getRotation();
#endif
}

FlatTrd::~FlatTrd ( )

virtual

Definition at line 96 of file FlatTrd.cc.

{}

Member Function Documentation

const GlobalVector & FlatTrd::axis

(

)

const

Definition at line 116 of file FlatTrd.cc.

References m_axis.

Referenced by getTransform().

                                        { 
  return m_axis ; 
}

const GlobalPoint FlatTrd::backCtr ( void  ) const

private

Definition at line 244 of file FlatTrd.cc.

References CaloCellGeometry::getCorners(), m_local, m_tr, and CaloCellGeometry::param().

Referenced by makeAxis().

                                         {
  float dz = (getCorners()[4].z() > getCorners()[0].z()) ? 
    param()[0] : -param()[0];
  Pt3D local_b(m_local.x(),m_local.y(),m_local.z()+dz);
  Pt3D global_b = m_tr*local_b;
  GlobalPoint global(global_b.x(),global_b.y(),global_b.z());
  return global;
}

void FlatTrd::createCorners ( const std::vector< CCGFloat > &  pv, const Tr3D &  tr, std::vector< GlobalPoint > &  co  )

static

Definition at line 126 of file FlatTrd.cc.

References assert(), gather_cfg::cout, i, localCorners(), AlCaHLTBitMon_ParallelJobs::p, and tmp.

Referenced by HGCalGeometryLoader::build().

                                                     {

  assert( 11 <= pv.size() ) ;
  assert( 8 == co.size() ) ;

  Pt3DVec        ko ( 8, Pt3D(0,0,0) ) ;

  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
    const Pt3D & p ( ko[i] ) ;
    co[ i ] = GlobalPoint( p.x(), p.y(), p.z() ) ;
#ifdef DebugLog
    std::cout << "Corner[" << i << "] = " << co[i] << std::endl;
#endif
  }
}

virtual float FlatTrd::etaPos ( ) const

inlinevirtual

Definition at line 48 of file FlatTrd.h.

References PV3DBase< T, PVType, FrameType >::eta(), and m_global.

Referenced by FlatTrd(), operator<<(), and operator=().

{ return m_global.eta(); }

Pt3D FlatTrd::getLocal

(

const GlobalPoint &

global

)

const

Definition at line 103 of file FlatTrd.cc.

References m_tr, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                      {
  Pt3D local = m_tr.inverse()*Pt3D(global.x(),global.y(),global.z());
  return local;
}

CCGFloat FlatTrd::getPhiAxis

(

)

const

Definition at line 112 of file FlatTrd.cc.

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

Referenced by operator<<().

                                   {
  return m_axis.phi() ; 
} 

virtual const GlobalPoint& FlatTrd::getPosition ( ) const

inlinevirtual

Definition at line 46 of file FlatTrd.h.

References m_global.

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

{ return m_global; }

GlobalPoint FlatTrd::getPosition

(

const Pt3D &

local

)

const

Definition at line 98 of file FlatTrd.cc.

References glb, and m_tr.

                                                         {
  Pt3D glb = m_tr*local;
  return GlobalPoint(glb.x(),glb.y(),glb.z());
}

CCGFloat FlatTrd::getThetaAxis

(

)

const

Definition at line 108 of file FlatTrd.cc.

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

Referenced by operator<<().

                                     { 
  return m_axis.theta() ; 
} 

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

virtual

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

Reimplemented from CaloCellGeometry.

Definition at line 180 of file FlatTrd.cc.

References angle(), assert(), axis(), alignCSCRings::e, CaloCellGeometry::getPosition(), localCorners(), 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 FlatTrd(), and initCorners().

                                                          {
  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() ) ;
  std::vector<Pt3D > lc( 8, Pt3D(0,0,0) ) ;
  localCorners( lc, param(), lFront ) ;

  // figure out if reflction volume or not

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

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

overrideprivatevirtual

Implements CaloCellGeometry.

Definition at line 225 of file FlatTrd.cc.

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

                                                        {

  if( co.uninitialized() ) {
    CornersVec& corners ( co ) ;
    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 FlatTrd::localCorners ( Pt3DVec &  vec, const CCGFloat *  pv, Pt3D &  ref  )

static

Definition at line 149 of file FlatTrd.cc.

References assert(), gather_cfg::cout, h, and funct::tan().

Referenced by createCorners(), getTransform(), and vocalCorners().

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

   const CCGFloat dz ( pv[0] ) ;
   const CCGFloat h  ( pv[3] ) ;
   const CCGFloat bl ( pv[4] ) ;
   const CCGFloat tl ( pv[5] ) ;
   const CCGFloat a1 ( pv[6] ) ;
  
   const CCGFloat ta1 ( tan( a1 ) ) ;

   lc[0] = Pt3D ( - h*ta1 - bl, - h , -dz ); // (-,-,-)
   lc[1] = Pt3D ( + h*ta1 - tl, + h , -dz ); // (-,+,-)
   lc[2] = Pt3D ( + h*ta1 + tl, + h , -dz ); // (+,+,-)
   lc[3] = Pt3D ( - h*ta1 + bl, - h , -dz ); // (+,-,-)
   lc[4] = Pt3D ( - h*ta1 - bl, - h ,  dz ); // (-,-,+)
   lc[5] = Pt3D ( + h*ta1 - tl, + h ,  dz ); // (-,+,+)
   lc[6] = Pt3D ( + h*ta1 + tl, + h ,  dz ); // (+,+,+)
   lc[7] = Pt3D ( - h*ta1 + bl, - h ,  dz ); // (+,-,+)

   ref   = 0.25*( lc[0] + lc[1] + lc[2] + lc[3] ) ;
#ifdef DebugLog
   std::cout << "Ref " << ref << " Local Corners " << lc[0] << "|" << lc[1] 
         << "|" << lc[2] << "|" << lc[3] << "|" << lc[4] << "|" << lc[5]
         << "|" << lc[6] << "|" << lc[7] << std::endl;
#endif
}

GlobalVector FlatTrd::makeAxis ( void  )

private

Definition at line 240 of file FlatTrd.cc.

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

Referenced by FlatTrd().

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

FlatTrd & FlatTrd::operator=

(

const FlatTrd &

tr

)

Definition at line 31 of file FlatTrd.cc.

References gather_cfg::cout, etaPos(), m_axis, m_corOne, m_global, m_local, m_tr, and phiPos().

                                               {
  CaloCellGeometry::operator=( tr ) ;
  if ( this != &tr ) {
    m_axis   = tr.m_axis ;
    m_corOne = tr.m_corOne ; 
    m_local  = tr.m_local;
    m_global = tr.m_global;
    m_tr     = tr.m_tr;
  }
#ifdef DebugLog
  std::cout << "FlatTrd(Copy): Local " << m_local << " Global " << m_global
        << " eta " << etaPos() << " phi " << phiPos() << " Translation "
        << m_tr.getTranslation() << " and rotation " << m_tr.getRotation();
#endif
  return *this ; 
}

virtual float FlatTrd::phiPos ( ) const

inlinevirtual

Definition at line 49 of file FlatTrd.h.

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

Referenced by FlatTrd(), operator<<(), and operator=().

{ return m_global.phi(); }

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

virtual

Implements CaloCellGeometry.

Definition at line 120 of file FlatTrd.cc.

References localCorners().

                                             { 
  localCorners( vec, pv, ref ) ; 
}

Member Data Documentation

GlobalVector FlatTrd::m_axis

private

Definition at line 81 of file FlatTrd.h.

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

Pt3D FlatTrd::m_corOne

private

Definition at line 82 of file FlatTrd.h.

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

GlobalPoint FlatTrd::m_global

private

Definition at line 83 of file FlatTrd.h.

Referenced by etaPos(), FlatTrd(), getPosition(), operator=(), and phiPos().

Pt3D FlatTrd::m_local

private

Definition at line 82 of file FlatTrd.h.

Referenced by backCtr(), FlatTrd(), and operator=().

Tr3D FlatTrd::m_tr

private

Definition at line 84 of file FlatTrd.h.

Referenced by backCtr(), FlatTrd(), getLocal(), getPosition(), and operator=().