---

EcalEndcapGeometry Class Reference

#include <Geometry/EcalAlgo/interface/EcalEndcapGeometry.h>

Inheritance diagram for EcalEndcapGeometry:

List of all members.

Public Types
typedef EcalEndcapGeometryRecord	AlignedRecord
typedef EEAlignmentRcd	AlignmentRecord
enum	CornersCount { k_NumberOfCellsForCorners = 14648 }
typedef IdealGeometryRecord	IdealRecord
typedef EcalEndcapNumberingScheme	NumberingScheme
typedef EZArrayFL< EBDetId >	OrderedListOfEBDetId
typedef std::vector < OrderedListOfEBDetId * >	VecOrdListEBDetIdPtr
Public Member Functions
	EcalEndcapGeometry ()
virtual CaloSubdetectorGeometry::DetIdSet	getCells (const GlobalPoint &r, double dR) const
	Get a list of all cells within a dR of the given cell.
const OrderedListOfEBDetId *	getClosestBarrelCells (EEDetId id) const
virtual DetId	getClosestCell (const GlobalPoint &r) const
int	getNumberOfCrystalPerModule () const
int	getNumberOfModules () const
void	initialize ()
void	setNumberOfCrystalPerModule (const int nncrys)
void	setNumberOfModules (const int nnmods)
virtual	~EcalEndcapGeometry ()
Static Public Member Functions
static std::string	hitString ()
static unsigned int	numberOfAlignments ()
static std::string	producerName ()
Private Member Functions
EEDetId	gId (float x, float y, float z) const
unsigned int	index (float x) const
Static Private Member Functions
static int	myPhi (int i)
Private Attributes
int	_nncrys
	number of crystals per module
int	_nnmods
	number of modules
EZMgrFL< EBDetId > *	m_borderMgr
VecOrdListEBDetIdPtr *	m_borderPtrVec
float	m_href
unsigned int	m_nref
float	m_wref
float	zeN
float	zeP

---

Detailed Description

Definition at line 17 of file EcalEndcapGeometry.h.

---

Member Typedef Documentation

typedef EcalEndcapGeometryRecord EcalEndcapGeometry::AlignedRecord

Definition at line 22 of file EcalEndcapGeometry.h.

typedef EEAlignmentRcd EcalEndcapGeometry::AlignmentRecord

Definition at line 23 of file EcalEndcapGeometry.h.

typedef IdealGeometryRecord EcalEndcapGeometry::IdealRecord

Definition at line 21 of file EcalEndcapGeometry.h.

typedef EcalEndcapNumberingScheme EcalEndcapGeometry::NumberingScheme

Definition at line 29 of file EcalEndcapGeometry.h.

typedef EZArrayFL<EBDetId> EcalEndcapGeometry::OrderedListOfEBDetId

Definition at line 25 of file EcalEndcapGeometry.h.

typedef std::vector<OrderedListOfEBDetId*> EcalEndcapGeometry::VecOrdListEBDetIdPtr

Definition at line 27 of file EcalEndcapGeometry.h.

---

Member Enumeration Documentation

enum EcalEndcapGeometry::CornersCount

Enumerator:

k_NumberOfCellsForCorners

Definition at line 31 of file EcalEndcapGeometry.h.

{ k_NumberOfCellsForCorners = 14648 } ;

---

Constructor & Destructor Documentation

EcalEndcapGeometry::EcalEndcapGeometry

(

)

Definition at line 8 of file EcalEndcapGeometry.cc.

                                       :
   _nnmods ( 0 ) ,
   _nncrys ( 0 ) ,
   m_borderMgr ( 0 ),
   m_borderPtrVec ( 0 ) 
{
}

EcalEndcapGeometry::~EcalEndcapGeometry

(

)

[virtual]

Definition at line 16 of file EcalEndcapGeometry.cc.

References m_borderMgr, and m_borderPtrVec.

{
   delete m_borderPtrVec ;
   delete m_borderMgr ;
}

---

Member Function Documentation

CaloSubdetectorGeometry::DetIdSet EcalEndcapGeometry::getCells	(	const GlobalPoint &	r,
		double	dR
	)			const [virtual]

Get a list of all cells within a dR of the given cell.

The default implementation makes a loop over all cell geometries. Cleverer implementations are suggested to use rough conversions between eta/phi and ieta/iphi and test on the boundaries.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 242 of file EcalEndcapGeometry.cc.

References reco::deltaR2(), eta, PV3DBase< T, PVType, FrameType >::eta(), CaloSubdetectorGeometry::getCells(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), id, m_href, m_nref, m_wref, p, PV3DBase< T, PVType, FrameType >::phi(), phi, funct::tan(), EEDetId::validDetId(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and zeP.

{
   CaloSubdetectorGeometry::DetIdSet dis ; // return object
   if( 0.000001 < dR )
   {
      if( dR > M_PI/2. ) // this code assumes small dR
      {
         dis = CaloSubdetectorGeometry::getCells( r, dR ) ; // base class version
      }
      else
      {
         const double dR2  ( dR*dR ) ;
         const double reta ( r.eta() ) ;
         const double rphi ( r.phi() ) ;
         const double rx   ( r.x() ) ;
         const double ry   ( r.y() ) ;
         const double rz   ( r.z() ) ;
         const double fac  ( fabs( zeP/rz ) ) ;
         const double xx   ( rx*fac ) ; // xyz at endcap z
         const double yy   ( ry*fac ) ; 
         const double zz   ( rz*fac ) ; 

         const double xang  ( atan( xx/zz ) ) ;
         const double lowX  ( zz>0 ? zz*tan( xang - dR ) : zz*tan( xang + dR ) ) ;
         const double highX ( zz>0 ? zz*tan( xang + dR ) : zz*tan( xang - dR ) ) ;
         const double yang  ( atan( yy/zz ) ) ;
         const double lowY  ( zz>0 ? zz*tan( yang - dR ) : zz*tan( yang + dR ) ) ;
         const double highY ( zz>0 ? zz*tan( yang + dR ) : zz*tan( yang - dR ) ) ;

         if( lowX  <  m_href &&   // proceed if any possible overlap with the endcap
             lowY  <  m_href &&
             highX > -m_href &&
             highY > -m_href    )
         {
            const int ix_ctr ( 1 + (int)floor( ( xx + m_href )/m_wref ) ) ;
            const int iy_ctr ( 1 + (int)floor( ( yy + m_href )/m_wref ) ) ;
            const int iz     ( rz>0 ? 1 : -1 ) ;
            
            const int ix_hi  ( ix_ctr + int( ( highX - xx )/m_wref ) + 2 ) ;
            const int ix_lo  ( ix_ctr - int( ( xx - lowX  )/m_wref ) - 2 ) ;
            
            const int iy_hi  ( iy_ctr + int( ( highY - yy )/m_wref ) + 2 ) ;
            const int iy_lo  ( iy_ctr - int( ( yy - lowY  )/m_wref ) - 2 ) ;
            
            for( int kx ( ix_lo ) ; kx <= ix_hi ; ++kx ) 
            {
               if( kx >  0      && 
                   kx <= m_nref    )
               {
                  for( int ky ( iy_lo ) ; ky <= iy_hi ; ++ky ) 
                  {
                     if( ky >  0      && 
                         ky <= m_nref    )
                     {
                        if( EEDetId::validDetId( kx, ky, iz ) ) // reject invalid ids
                        {
                           const EEDetId id ( kx, ky, iz ) ;
                           const CaloCellGeometry* cell ( getGeometry( id ) );
                           if( 0 != cell )
                           {
                              const GlobalPoint& p    ( cell->getPosition() ) ;
                              const double       eta  ( p.eta() ) ;
                              const double       phi  ( p.phi() ) ;
                              if( reco::deltaR2( eta, phi, reta, rphi ) < dR2 ) dis.insert( id ) ;
                           }
                        }
                     }
                  }
               }
            }
         }
      }
   }
   return dis;
}

const EcalEndcapGeometry::OrderedListOfEBDetId * EcalEndcapGeometry::getClosestBarrelCells

(

EEDetId

id

)

const

Definition at line 320 of file EcalEndcapGeometry.cc.

References degrees(), CaloSubdetectorGeometry::getGeometry(), i, m_borderMgr, m_borderPtrVec, myPhi(), phi, and ptr.

Referenced by PFRecHitProducerECAL::stdsimplemove().

{
   OrderedListOfEBDetId* ptr ( 0 ) ;
   if( 0 != id.rawId() )
   {
      const float phi ( 370. +
                        getGeometry( id )->getPosition().phi().degrees() );
      const int iPhi ( 1 + int(phi)%360 ) ;
      const int iz ( id.zside() ) ;
      if( 0 == m_borderMgr )
      {
         m_borderMgr = new EZMgrFL<EBDetId>( 720*9, 9 ) ;
      }
      if( 0 == m_borderPtrVec )
      {
         m_borderPtrVec = new VecOrdListEBDetIdPtr() ;
         m_borderPtrVec->reserve( 720 ) ;
         for( unsigned int i ( 0 ) ; i != 720 ; ++i )
         {
            const int kz     ( 360>i ? -1 : 1 ) ;
            const int iEta   ( kz*85 ) ;
            const int iEtam1 ( kz*84 ) ;
            const int iEtam2 ( kz*83 ) ;
            const int jPhi   ( i%360 + 1 ) ;
            OrderedListOfEBDetId& olist ( *new OrderedListOfEBDetId( m_borderMgr ) );
            olist[0]=EBDetId( iEta  ,        jPhi     ) ;
            olist[1]=EBDetId( iEta  , myPhi( jPhi+1 ) ) ;
            olist[2]=EBDetId( iEta  , myPhi( jPhi-1 ) ) ;
            olist[3]=EBDetId( iEtam1,        jPhi     ) ;
            olist[4]=EBDetId( iEtam1, myPhi( jPhi+1 ) ) ;
            olist[5]=EBDetId( iEtam1, myPhi( jPhi-1 ) ) ;
            olist[6]=EBDetId( iEta  , myPhi( jPhi+2 ) ) ;
            olist[7]=EBDetId( iEta  , myPhi( jPhi-2 ) ) ;
            olist[8]=EBDetId( iEtam2,        jPhi     ) ;
            m_borderPtrVec->push_back( &olist ) ;
         }
      }
      ptr = (*m_borderPtrVec)[ ( iPhi - 1 ) + ( 0>iz ? 0 : 360 ) ] ;
   }
   return ptr ;
}

DetId EcalEndcapGeometry::getClosestCell

(

const GlobalPoint &

r

)

const [virtual]

Reimplemented from CaloSubdetectorGeometry.

Definition at line 133 of file EcalEndcapGeometry.cc.

References funct::A, funct::C, e, CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), gId(), i, EEDetId::ix(), EEDetId::iy(), offset, EEDetId::validDetId(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), y, PV3DBase< T, PVType, FrameType >::z(), z, and EEDetId::zside().

Referenced by CaloGeometryHelper::getClosestCell().

{
   try
   {
      EEDetId mycellID ( gId( r.x(), r.y(), r.z() ) ) ; // educated guess

      if( EEDetId::validDetId( mycellID.ix(), 
                               mycellID.iy(),
                               mycellID.zside() ) )
      {
         // now get points in convenient ordering

         HepPoint3D  A;
         HepPoint3D  B;
         HepPoint3D  C;
         HepPoint3D  point(r.x(),r.y(),r.z());
         // D.K. : equation of plane : AA*x+BB*y+CC*z+DD=0;
         // finding equation for each edge
         
         // ================================================================
         double x,y,z;
         unsigned offset=0;
         int zsign=1;
         //================================================================
         std::vector<double> SS;
      
         // compute the distance of the point with respect of the 4 crystal lateral planes
         const GlobalPoint& myPosition=getGeometry(mycellID)->getPosition();
         
         x=myPosition.x();
         y=myPosition.y();
         z=myPosition.z();
         
         offset=0;
         // This will disappear when Andre has applied his fix
         zsign=1;
         
         if(z>0)
         {
            if(x>0&&y>0)
               offset=1;
            else  if(x<0&&y>0)
               offset=2;
            else if(x>0&&y<0)
               offset=0;
            else if (x<0&&y<0)
               offset=3;
            zsign=1;
         }
         else
         {
            if(x>0&&y>0)
               offset=3;
            else if(x<0&&y>0)
               offset=2;
            else if(x>0&&y<0)
               offset=0;
            else if(x<0&&y<0)
               offset=1;
            zsign=-1;
         }
         std::vector<GlobalPoint> corners;
         corners.clear();
         corners.resize(8);
         for(unsigned ic=0;ic<4;++ic)
         {
            corners[ic]=getGeometry(mycellID)->getCorners()[(unsigned)((zsign*ic+offset)%4)];
            corners[4+ic]=getGeometry(mycellID)->getCorners()[(unsigned)(4+(zsign*ic+offset)%4)];
         }
         
         for (short i=0; i < 4 ; ++i)
         {
            A = HepPoint3D(corners[i%4].x(),corners[i%4].y(),corners[i%4].z());
            B = HepPoint3D(corners[(i+1)%4].x(),corners[(i+1)%4].y(),corners[(i+1)%4].z());
            C = HepPoint3D(corners[4+(i+1)%4].x(),corners[4+(i+1)%4].y(),corners[4+(i+1)%4].z());
            HepPlane3D plane(A,B,C);
            plane.normalize();
            double distance = plane.distance(point);
            if (corners[0].z()<0.) distance=-distance;
            SS.push_back(distance);
         }
         
         // Only one move in necessary direction
         
         const bool yout ( 0 > SS[0]*SS[2] ) ;
         const bool xout ( 0 > SS[1]*SS[3] ) ;
         
         if( yout || xout )
         {
            const int ydel ( !yout ? 0 :  ( 0 < SS[0] ? -1 : 1 ) ) ;
            const int xdel ( !xout ? 0 :  ( 0 < SS[1] ? -1 : 1 ) ) ;
            const unsigned int ix ( mycellID.ix() + xdel ) ;
            const unsigned int iy ( mycellID.iy() + ydel ) ;
            const unsigned int iz ( mycellID.zside()     ) ;
            if( EEDetId::validDetId( ix, iy, iz ) ) 
               mycellID = EEDetId( ix, iy, iz ) ;
         }
  
         return mycellID;
      }
      return DetId(0);
   }
   catch ( cms::Exception &e ) 
   { 
      return DetId(0);
   }
}

int EcalEndcapGeometry::getNumberOfCrystalPerModule

(

)

const [inline]

Definition at line 39 of file EcalEndcapGeometry.h.

References _nncrys.

{ return _nncrys ; }

int EcalEndcapGeometry::getNumberOfModules

(

)

const [inline]

Definition at line 37 of file EcalEndcapGeometry.h.

References _nnmods.

{ return _nnmods ; }

EEDetId EcalEndcapGeometry::gId	(	float	x,
		float	y,
		float	z
	)			const [private]

Definition at line 97 of file EcalEndcapGeometry.cc.

References i, index(), k, EEDetId::validDetId(), and zeP.

Referenced by getClosestCell().

{
   const double       fac ( fabs(zeP/z) ) ;
   const unsigned int ix  ( index( x*fac ) ) ; 
   const unsigned int iy  ( index( y*fac ) ) ; 
   const unsigned int iz  ( z>0 ? 1 : -1 ) ;

   if( EEDetId::validDetId( ix, iy, iz ) ) 
   {
      return EEDetId( ix, iy, iz ) ; // first try is on target
   }
   else // try nearby coordinates, spiraling out from center
   {
      for( unsigned int i ( 1 ) ; i != 6 ; ++i )
      {
         for( unsigned int k ( 0 ) ; k != 8 ; ++k )
         {
            const int jx ( 0 == k || 4 == k || 5 == k ? +i :
                           ( 1 == k || 5 < k ? -i : 0 ) ) ;
            const int jy ( 2 == k || 4 == k || 6 == k ? +i :
                           ( 3 == k || 5 == k || 7 == k ? -i : 0 ) ) ;
            if( EEDetId::validDetId( ix + jx, iy + jy, iz ) ) 
            {
               return EEDetId( ix + jx, iy + jy, iz ) ;
            }
         }
      }
   }
   return EEDetId() ; // nowhere near any crystal
}

static std::string EcalEndcapGeometry::hitString

(

)

[inline, static]

Definition at line 55 of file EcalEndcapGeometry.h.

{ return "EcalHitsEE" ; }

unsigned int EcalEndcapGeometry::index

(

float

x

)

const [private]

Definition at line 82 of file EcalEndcapGeometry.cc.

References i, m_href, m_nref, and m_wref.

Referenced by gId().

{
   int i ( 1 + (int)floor( ( x + m_href )/m_wref ) ) ;
   if( 1      > i )
   {
      i = 1 ;
   }
   else
   {
      if( m_nref < (unsigned int)i ) i = m_nref ;
   }
   return i ;
}

void EcalEndcapGeometry::initialize

(

)

Definition at line 23 of file EcalEndcapGeometry.cc.

References funct::abs(), begin, CaloSubdetectorGeometry::cellGeometries(), i, EEDetId::ix(), EEDetId::iy(), m_href, m_nref, m_wref, p, x, z, zeN, and zeP.

{
  zeP=0.;
  zeN=0.;
  unsigned nP=0;
  unsigned nN=0;
  m_nref = 0 ;
  for( CaloSubdetectorGeometry::CellCont::const_iterator i ( cellGeometries().begin() ) ; 
       i != cellGeometries().end(); ++i )
  {
//      addCrystalToZGridmap(i->first,dynamic_cast<const TruncatedPyramid*>(i->second));
     float z=dynamic_cast<const TruncatedPyramid*>(i->second)->getPosition(0.).z();
     if(z>0.)
     {
        zeP+=z;
        ++nP;
     }
     else
     {
        zeN+=z;
        ++nN;
     }
     const EEDetId myId ( i->first ) ;
     const unsigned int ix ( myId.ix() ) ;
     const unsigned int iy ( myId.iy() ) ;
     if( abs( ix ) > m_nref ) m_nref = abs( ix ) ;
     if( abs( iy ) > m_nref ) m_nref = abs( iy ) ;
  }
  zeP/=(float)nP;
  zeN/=(float)nN;

  m_href = 0 ;
  for( CaloSubdetectorGeometry::CellCont::const_iterator i ( cellGeometries().begin() ) ; 
       i != cellGeometries().end(); ++i )
  {
     const EEDetId myId ( i->first ) ;
     const unsigned int ix ( myId.ix() ) ;
     const unsigned int iy ( myId.iy() ) ;
     if( ix == m_nref ) 
     {
        const GlobalPoint p ( dynamic_cast<const TruncatedPyramid*>
                              (i->second)->getPosition(0.)  ) ;
        const float x ( p.x()*fabs(zeP/p.z()) ) ;
        if( m_href < x ) m_href = x ;
     }
     if( iy == m_nref ) 
     {
        const GlobalPoint p ( dynamic_cast<const TruncatedPyramid*>
                              (i->second)->getPosition(0.)  ) ;
        const float x ( p.y()*fabs(zeP/p.z()) ) ;
        if( m_href < x ) m_href = x ;
     }
  }
  m_href = m_href*(1.*m_nref)/( 1.*m_nref - 1. ) ;
  m_wref = m_href/(0.5*m_nref) ;
}

static int EcalEndcapGeometry::myPhi

(

int

i

)

[inline, static, private]

Definition at line 63 of file EcalEndcapGeometry.h.

Referenced by getClosestBarrelCells().

{ i+=720; return ( 1 + (i-1)%360 ) ; }

static unsigned int EcalEndcapGeometry::numberOfAlignments

(

)

[inline, static]

Definition at line 59 of file EcalEndcapGeometry.h.

Referenced by FakeCaloAlignmentEP::produceEEAli().

{ return 4 ; }

static std::string EcalEndcapGeometry::producerName

(

)

[inline, static]

Definition at line 57 of file EcalEndcapGeometry.h.

{ return "EcalEndcap" ; }

void EcalEndcapGeometry::setNumberOfCrystalPerModule

(

const int

nncrys

)

[inline]

Definition at line 43 of file EcalEndcapGeometry.h.

References _nncrys.

{ _nncrys=nncrys ; }

void EcalEndcapGeometry::setNumberOfModules

(

const int

nnmods

)

[inline]

Definition at line 41 of file EcalEndcapGeometry.h.

References _nnmods.

{ _nnmods=nnmods ; }

---

Member Data Documentation

int EcalEndcapGeometry::_nncrys [private]

number of crystals per module

Definition at line 69 of file EcalEndcapGeometry.h.

Referenced by getNumberOfCrystalPerModule(), and setNumberOfCrystalPerModule().

int EcalEndcapGeometry::_nnmods [private]

number of modules

Definition at line 66 of file EcalEndcapGeometry.h.

Referenced by getNumberOfModules(), and setNumberOfModules().

EZMgrFL<EBDetId>* EcalEndcapGeometry::m_borderMgr [mutable, private]

Definition at line 80 of file EcalEndcapGeometry.h.

Referenced by getClosestBarrelCells(), and ~EcalEndcapGeometry().

VecOrdListEBDetIdPtr* EcalEndcapGeometry::m_borderPtrVec [mutable, private]

Definition at line 82 of file EcalEndcapGeometry.h.

Referenced by getClosestBarrelCells(), and ~EcalEndcapGeometry().

float EcalEndcapGeometry::m_href [private]

Definition at line 73 of file EcalEndcapGeometry.h.

Referenced by getCells(), index(), and initialize().

unsigned int EcalEndcapGeometry::m_nref [private]

Definition at line 75 of file EcalEndcapGeometry.h.

Referenced by getCells(), index(), and initialize().

float EcalEndcapGeometry::m_wref [private]

Definition at line 73 of file EcalEndcapGeometry.h.

Referenced by getCells(), index(), and initialize().

float EcalEndcapGeometry::zeN [private]

Definition at line 71 of file EcalEndcapGeometry.h.

Referenced by initialize().

float EcalEndcapGeometry::zeP [private]

Definition at line 71 of file EcalEndcapGeometry.h.

Referenced by getCells(), gId(), and initialize().

---

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

• Geometry/EcalAlgo/interface/EcalEndcapGeometry.h
• Geometry/EcalAlgo/src/EcalEndcapGeometry.cc

---