EcalEndcapGeometry Class Reference

#include <EcalEndcapGeometry.h>

Inheritance diagram for EcalEndcapGeometry:

Public Types
enum	{ k_NumberOfCellsForCorners = EEDetId::kSizeForDenseIndexing }
enum	{ k_NumberOfShapes = 1 }
enum	{ k_NumberOfParametersPerShape = 11 }
typedef EcalEndcapGeometryRecord	AlignedRecord
typedef EEAlignmentRcd	AlignmentRecord
typedef CaloCellGeometry::CCGFloat	CCGFloat
typedef std::vector< TruncatedPyramid >	CellVec
typedef EEDetId	DetIdType
typedef IdealGeometryRecord	IdealRecord
typedef EcalEndcapNumberingScheme	NumberingScheme
typedef EZArrayFL< EBDetId >	OrderedListOfEBDetId
typedef PEcalEndcapRcd	PGeometryRecord
typedef CaloCellGeometry::Pt3D	Pt3D
typedef CaloCellGeometry::Pt3DVec	Pt3DVec
typedef std::vector< OrderedListOfEBDetId * >	VecOrdListEBDetIdPtr
Public Types inherited from CaloSubdetectorGeometry
typedef CaloCellGeometry::CCGFloat	CCGFloat
typedef std::vector< CaloCellGeometry const * >	CellSet
typedef std::set< DetId >	DetIdSet
typedef std::vector< CCGFloat >	DimVec
typedef std::vector< unsigned int >	IVec
typedef CaloCellGeometry::ParMgr	ParMgr
typedef CaloCellGeometry::ParVec	ParVec
typedef CaloCellGeometry::ParVecVec	ParVecVec
typedef std::vector< CCGFloat >	TrVec

Public Member Functions
CCGFloat	avgAbsZFrontFaceCenter () const
	EcalEndcapGeometry ()
CaloSubdetectorGeometry::DetIdSet	getCells (const GlobalPoint &r, double dR) const override
	Get a list of all cells within a dR of the given cell. More...
const OrderedListOfEBDetId *	getClosestBarrelCells (EEDetId id) const
DetId	getClosestCell (const GlobalPoint &r) const override
int	getNumberOfCrystalPerModule () const
int	getNumberOfModules () const
void	initializeParms () override
void	newCell (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId) override
unsigned int	numberOfParametersPerShape () const override
unsigned int	numberOfShapes () const override
void	setNumberOfCrystalPerModule (const int nncrys)
void	setNumberOfModules (const int nnmods)
	~EcalEndcapGeometry () override
Public Member Functions inherited from CaloSubdetectorGeometry
void	allocateCorners (CaloCellGeometry::CornersVec::size_type n)
void	allocatePar (ParVec::size_type n, unsigned int m)
	CaloSubdetectorGeometry ()
	CaloSubdetectorGeometry (const CaloSubdetectorGeometry &)=delete
	avoid copies More...
CaloCellGeometry::CornersMgr *	cornersMgr ()
CCGFloat	deltaEta (const DetId &detId) const
CCGFloat	deltaPhi (const DetId &detId) const
virtual void	fillDefaultNamedParameters () const
virtual CellSet	getCellSet (const GlobalPoint &r, double dR) const
virtual const CaloCellGeometry *	getGeometry (const DetId &id) const
	Get the cell geometry of a given detector id. Should return false if not found. More...
virtual void	getSummary (TrVec &trVector, IVec &iVector, DimVec &dimVector, IVec &dinsVector) const
virtual const std::vector< DetId > &	getValidDetIds (DetId::Detector det=DetId::Detector(0), int subdet=0) const
	Get a list of valid detector ids (for the given subdetector) More...
virtual unsigned int	numberOfTransformParms () const
CaloSubdetectorGeometry &	operator= (const CaloSubdetectorGeometry &)=delete
ParMgr *	parMgr ()
const ParMgr *	parMgrConst () const
ParVecVec &	parVecVec ()
const ParVecVec &	parVecVec () const
virtual bool	present (const DetId &id) const
	is this detid present in the geometry? More...
virtual	~CaloSubdetectorGeometry ()
	The base class DOES assume that it owns the CaloCellGeometry objects. More...

Static Public Member Functions
static unsigned int	alignmentTransformIndexGlobal (const DetId &id)
static unsigned int	alignmentTransformIndexLocal (const DetId &id)
static std::string	dbString ()
static DetId	detIdFromLocalAlignmentIndex (unsigned int iLoc)
static std::string	hitString ()
static void	localCorners (Pt3DVec &lc, const CCGFloat *pv, unsigned int i, Pt3D &ref)
static unsigned int	numberOfAlignments ()
static std::string	producerTag ()

Protected Member Functions
const CaloCellGeometry *	cellGeomPtr (uint32_t index) const override
Protected Member Functions inherited from CaloSubdetectorGeometry
void	addValidID (const DetId &id)
virtual unsigned int	indexFor (const DetId &id) const
virtual unsigned int	sizeForDenseIndex (const DetId &id) const

Private Member Functions
EEDetId	gId (float x, float y, float z) const
unsigned int	xindex (CCGFloat x, CCGFloat z) const
unsigned int	yindex (CCGFloat y, CCGFloat z) const

Static Private Member Functions
static int	myPhi (int i)

Private Attributes
int	_nncrys
	number of crystals per module More...
int	_nnmods
	number of modules More...
CCGFloat	m_avgZ
std::atomic< EZMgrFL< EBDetId > * >	m_borderMgr
std::atomic< VecOrdListEBDetIdPtr * >	m_borderPtrVec
CellVec	m_cellVec
std::atomic< bool >	m_check
CCGFloat	m_del
unsigned int	m_nref
CCGFloat	m_wref
CCGFloat	m_xhi [2]
CCGFloat	m_xlo [2]
CCGFloat	m_xoff [2]
CCGFloat	m_yhi [2]
CCGFloat	m_ylo [2]
CCGFloat	m_yoff [2]
CCGFloat	zeN
CCGFloat	zeP

Additional Inherited Members
Static Protected Member Functions inherited from CaloSubdetectorGeometry
static CCGFloat	deltaR (const GlobalPoint &p1, const GlobalPoint &p2)
Protected Attributes inherited from CaloSubdetectorGeometry
ParVecVec	m_parVecVec
std::vector< DetId >	m_validIds

Detailed Description

Definition at line 20 of file EcalEndcapGeometry.h.

Member Typedef Documentation

typedef EcalEndcapGeometryRecord EcalEndcapGeometry::AlignedRecord

Definition at line 31 of file EcalEndcapGeometry.h.

typedef EEAlignmentRcd EcalEndcapGeometry::AlignmentRecord

Definition at line 32 of file EcalEndcapGeometry.h.

typedef CaloCellGeometry::CCGFloat EcalEndcapGeometry::CCGFloat

Definition at line 26 of file EcalEndcapGeometry.h.

typedef std::vector<TruncatedPyramid> EcalEndcapGeometry::CellVec

Definition at line 24 of file EcalEndcapGeometry.h.

typedef EEDetId EcalEndcapGeometry::DetIdType

Definition at line 41 of file EcalEndcapGeometry.h.

typedef IdealGeometryRecord EcalEndcapGeometry::IdealRecord

Definition at line 30 of file EcalEndcapGeometry.h.

typedef EcalEndcapNumberingScheme EcalEndcapGeometry::NumberingScheme

Definition at line 39 of file EcalEndcapGeometry.h.

typedef EZArrayFL<EBDetId> EcalEndcapGeometry::OrderedListOfEBDetId

Definition at line 35 of file EcalEndcapGeometry.h.

typedef PEcalEndcapRcd EcalEndcapGeometry::PGeometryRecord

Definition at line 33 of file EcalEndcapGeometry.h.

typedef CaloCellGeometry::Pt3D EcalEndcapGeometry::Pt3D

Definition at line 27 of file EcalEndcapGeometry.h.

typedef CaloCellGeometry::Pt3DVec EcalEndcapGeometry::Pt3DVec

Definition at line 28 of file EcalEndcapGeometry.h.

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

Definition at line 37 of file EcalEndcapGeometry.h.

Member Enumeration Documentation

anonymous enum

Enumerator
k_NumberOfCellsForCorners

Definition at line 43 of file EcalEndcapGeometry.h.

{ k_NumberOfCellsForCorners = EEDetId::kSizeForDenseIndexing } ;

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 45 of file EcalEndcapGeometry.h.

{ k_NumberOfShapes = 1 } ;

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 47 of file EcalEndcapGeometry.h.

{ k_NumberOfParametersPerShape = 11 } ;

Constructor & Destructor Documentation

EcalEndcapGeometry::EcalEndcapGeometry

(

void

)

Definition at line 14 of file EcalEndcapGeometry.cc.

References m_xhi, m_xlo, m_xoff, m_yhi, m_ylo, and m_yoff.

Referenced by numberOfParametersPerShape().

  : _nnmods( 316 ),
    _nncrys( 25 ),
    zeP( 0. ),
    zeN( 0. ),
    m_wref( 0. ),
    m_del( 0. ),
    m_nref( 0 ),
    m_borderMgr( nullptr ),
    m_borderPtrVec( nullptr ),
    m_avgZ( -1 ),
    m_check( false ),
    m_cellVec( k_NumberOfCellsForCorners )
{
  m_xlo[0] = 999.;
  m_xlo[1] = 999.;
  m_xhi[0] = -999.;
  m_xhi[1] = -999.;
  m_ylo[0] = 999.;
  m_ylo[1] = 999.;
  m_yhi[0] = -999.;
  m_yhi[1] = -999.;
  m_xoff[0] = 0.;
  m_xoff[1] = 0.;
  m_yoff[0] = 0.;
  m_yoff[1] = 0.;
}

EcalEndcapGeometry::~EcalEndcapGeometry ( )

override

Definition at line 42 of file EcalEndcapGeometry.cc.

References m_borderMgr, m_borderPtrVec, and findQualityFiles::v.

Referenced by numberOfParametersPerShape().

{
  if(m_borderPtrVec)
  {
    auto ptr = m_borderPtrVec.load(std::memory_order_acquire);
    for(auto& v: (*ptr)) {
        delete v;
        v = nullptr;
    }
    delete m_borderPtrVec.load() ;
  }
  delete m_borderMgr.load() ;
}

Member Function Documentation

unsigned int EcalEndcapGeometry::alignmentTransformIndexGlobal ( const DetId &  id )

static

Definition at line 74 of file EcalEndcapGeometry.cc.

References DetId::Ecal.

Referenced by numberOfAlignments().

{
   return (unsigned int)DetId::Ecal - 1 ;
}

unsigned int EcalEndcapGeometry::alignmentTransformIndexLocal ( const DetId &  id )

static

Definition at line 57 of file EcalEndcapGeometry.cc.

References diffTreeTool::index, CaloGenericDetId::isEE(), and HcalDigiParam_cfi::zside.

Referenced by numberOfAlignments().

{
   const CaloGenericDetId gid ( id ) ;

   assert( gid.isEE() ) ;
   unsigned int index ( EEDetId(id).ix()/51 + ( EEDetId(id).zside()<0 ? 0 : 2 ) ) ;

   return index ;
}

CCGFloat EcalEndcapGeometry::avgAbsZFrontFaceCenter

(

)

const

Definition at line 500 of file EcalEndcapGeometry.cc.

References cellGeomPtr(), CaloCellGeometry::getPosition(), mps_fire::i, m_avgZ, m_cellVec, m_check, and PV3DBase< T, PVType, FrameType >::z().

Referenced by setNumberOfCrystalPerModule().

{
   if(!m_check.load(std::memory_order_acquire))
   {
      CCGFloat sum ( 0 ) ;
      for( unsigned int i ( 0 ) ; i != m_cellVec.size() ; ++i )
      {
     const CaloCellGeometry* cell ( cellGeomPtr(i) ) ;
     if( nullptr != cell )
     {
        sum += fabs( cell->getPosition().z() ) ;
     }
      }
      m_avgZ = sum/m_cellVec.size();
      m_check.store(true, std::memory_order_release);
   }
   return m_avgZ;
}

const CaloCellGeometry * EcalEndcapGeometry::cellGeomPtr ( uint32_t  index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 520 of file EcalEndcapGeometry.cc.

References m_cellVec, and CaloCellGeometry::param().

Referenced by avgAbsZFrontFaceCenter(), initializeParms(), and numberOfAlignments().

{
   const CaloCellGeometry* cell ( &m_cellVec[ index ] ) ;
   return ( m_cellVec.size() < index ||
        nullptr == cell->param() ? nullptr : cell ) ;
}

static std::string EcalEndcapGeometry::dbString ( )

inlinestatic

Definition at line 50 of file EcalEndcapGeometry.h.

{ return "PEcalEndcapRcd" ; }

DetId EcalEndcapGeometry::detIdFromLocalAlignmentIndex ( unsigned int  iLoc )

static

Definition at line 68 of file EcalEndcapGeometry.cc.

Referenced by numberOfAlignments(), FakeCaloAlignmentEP::produceEEAli(), and TestCaloAlignmentEP::produceEEAli().

{
   return EEDetId( 20 + 50*( iLoc%2 ), 50, 2*( iLoc/2 ) - 1 ) ;
}

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

overridevirtual

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 351 of file EcalEndcapGeometry.cc.

References funct::abs(), reco::deltaR2(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), CaloCellGeometry::etaPos(), CaloSubdetectorGeometry::getCells(), hcalTTPDigis_cfi::id, m_cellVec, m_nref, M_PI, m_wref, m_xhi, m_xlo, m_yhi, m_ylo, phi, PV3DBase< T, PVType, FrameType >::phi(), CaloCellGeometry::phiPos(), corrVsCorr::rx, corrVsCorr::ry, funct::tan(), EEDetId::validDetId(), PV3DBase< T, PVType, FrameType >::x(), xindex(), geometryCSVtoXML::xx, PV3DBase< T, PVType, FrameType >::y(), yindex(), geometryCSVtoXML::yy, PV3DBase< T, PVType, FrameType >::z(), zeP, and geometryCSVtoXML::zz.

Referenced by setNumberOfCrystalPerModule().

{
   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 float dR2  ( dR*dR ) ;
     const float reta ( r.eta() ) ;
     const float rphi ( r.phi() ) ;
     const float rx   ( r.x() ) ;
     const float ry   ( r.y() ) ;
     const float rz   ( r.z() ) ;
     const float fac  ( std::abs( zeP/rz ) ) ;
     const float xx   ( rx*fac ) ; // xyz at endcap z
     const float yy   ( ry*fac ) ; 
     const float zz   ( rz*fac ) ; 

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

     const float refxlo ( 0 > rz ? m_xlo[0] : m_xlo[1] ) ;
     const float refxhi ( 0 > rz ? m_xhi[0] : m_xhi[1] ) ;
     const float refylo ( 0 > rz ? m_ylo[0] : m_ylo[1] ) ;
     const float refyhi ( 0 > rz ? m_yhi[0] : m_yhi[1] ) ;

     if( lowX  <  refxhi &&   // proceed if any possible overlap with the endcap
         lowY  <  refyhi &&
         highX >  refxlo &&
         highY >  refylo    )
     {
        const int ix_ctr ( xindex( xx, rz ) ) ;
        const int iy_ctr ( yindex( yy, rz ) ) ;
        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 <= (int) m_nref    )
           {
          for( int ky ( iy_lo ) ; ky <= iy_hi ; ++ky ) 
          {
             if( ky >  0      && 
             ky <= (int) m_nref    )
             {
            if( EEDetId::validDetId( kx, ky, iz ) ) // reject invalid ids
            {
              const EEDetId id ( kx, ky, iz ) ;
              const CaloCellGeometry* cell  = &m_cellVec[ id.denseIndex()];
              const float       eta  (cell->etaPos() ) ;
              const float       phi  (cell->phiPos() ) ;
              if( reco::deltaR2( eta, phi, reta, rphi ) < dR2 ) dis.insert( id ) ;
            }
             }
          }
           }
        }
     }
      }
   }
   return dis;
}

const EcalEndcapGeometry::OrderedListOfEBDetId * EcalEndcapGeometry::getClosestBarrelCells

(

EEDetId

id

)

const

Definition at line 430 of file EcalEndcapGeometry.cc.

References degrees(), CaloSubdetectorGeometry::getGeometry(), mps_fire::i, m_borderMgr, m_borderPtrVec, myPhi(), phi, and HcalDigiParam_cfi::zside.

Referenced by setNumberOfCrystalPerModule(), spr::simpleMove(), and PFECALHashNavigator::stdsimplemove().

{
   OrderedListOfEBDetId* ptr ( nullptr ) ;
   auto ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
   if(!ptrVec) {
       if(0 != id.rawId() && nullptr != getGeometry(id)) {
           const float phi(370.+getGeometry(id)->getPosition().phi().degrees());
           const int iPhi(1+int(phi)%360) ;
           const int iz(id.zside()) ;
           if (!m_borderMgr.load(std::memory_order_acquire)) {
               EZMgrFL<EBDetId>* expect = nullptr;
               auto ptrMgr = new EZMgrFL<EBDetId>( 720*9, 9 );
               bool exchanged = m_borderMgr.compare_exchange_strong(expect, ptrMgr, std::memory_order_acq_rel);
               if(!exchanged) delete ptrMgr;
           }
           VecOrdListEBDetIdPtr* expect = nullptr;
           auto ptrVec = new VecOrdListEBDetIdPtr();
           ptrVec->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.load(std::memory_order_acquire) ) );
               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     ) ;
               ptrVec->emplace_back( &olist ) ;
           }
           bool exchanged = m_borderPtrVec.compare_exchange_strong(expect, ptrVec, std::memory_order_acq_rel);
           if(!exchanged) delete ptrVec;
           ptrVec = m_borderPtrVec.load(std::memory_order_acquire);
           ptr = (*ptrVec)[ ( iPhi - 1 ) + ( 0>iz ? 0 : 360 ) ] ;
       }
   }
   return ptr;
}

DetId EcalEndcapGeometry::getClosestCell ( const GlobalPoint &  r ) const

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 239 of file EcalEndcapGeometry.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, patCaloMETCorrections_cff::C, SoftLeptonByDistance_cfi::distance, MillePedeFileConverter_cfg::e, CaloCellGeometry::getCorners(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), gId(), mps_fire::i, EEDetId::ix(), EEDetId::iy(), PFRecoTauDiscriminationByIsolation_cfi::offset, point, EEDetId::validDetId(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, PV3DBase< T, PVType, FrameType >::z(), and EEDetId::zside().

Referenced by spr::chargeIsolationEcal(), CaloGeometryHelper::getClosestCell(), spr::propagateCALO(), spr::propagateCosmicCALO(), and setNumberOfCrystalPerModule().

{
   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

     Pt3D A;
     Pt3D B;
     Pt3D C;
     Pt3D 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
     
     // ================================================================
     CCGFloat x,y,z;
     unsigned offset=0;
     int zsign=1;
     //================================================================
      
     // compute the distance of the point with respect of the 4 crystal lateral planes

     
     if( nullptr != getGeometry(mycellID) )
     {
        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
     
        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;
        }
        GlobalPoint corners[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)];
        }

        CCGFloat SS[4];
        for (short i=0; i < 4 ; ++i)
        {
           A = Pt3D(corners[i%4].x(),corners[i%4].y(),corners[i%4].z());
           B = Pt3D(corners[(i+1)%4].x(),corners[(i+1)%4].y(),corners[(i+1)%4].z());
           C = Pt3D(corners[4+(i+1)%4].x(),corners[4+(i+1)%4].y(),corners[4+(i+1)%4].z());
           Pl3D plane(A,B,C);
           plane.normalize();
           CCGFloat distance = plane.distance(point);
           if (corners[0].z()<0.) distance=-distance;
           SS[i] = 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);
   }
   return DetId(0);
}

int EcalEndcapGeometry::getNumberOfCrystalPerModule ( ) const

inline

Definition at line 61 of file EcalEndcapGeometry.h.

References _nncrys.

{ return _nncrys ; }

int EcalEndcapGeometry::getNumberOfModules ( ) const

inline

Definition at line 59 of file EcalEndcapGeometry.h.

References _nnmods.

{ return _nnmods ; }

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

private

Definition at line 203 of file EcalEndcapGeometry.cc.

References mps_fire::i, gen::k, EEDetId::validDetId(), xindex(), yindex(), zeN, and zeP.

Referenced by getClosestCell().

{
   const CCGFloat     fac ( fabs( ( 0 > z ? zeN : zeP )/z ) ) ;
   const unsigned int ix  ( xindex( x*fac, z ) ) ; 
   const unsigned int iy  ( yindex( y*fac, z ) ) ; 
   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 ( )

inlinestatic

Definition at line 78 of file EcalEndcapGeometry.h.

{ return "EcalHitsEE" ; }

void EcalEndcapGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 80 of file EcalEndcapGeometry.cc.

References cellGeomPtr(), EEDetId::detIdFromDenseIndex(), CaloCellGeometry::getPosition(), mps_fire::i, EEDetId::ix(), EEDetId::iy(), m_cellVec, m_del, m_nref, m_wref, m_xhi, m_xlo, m_xoff, m_yhi, m_ylo, m_yoff, AlCaHLTBitMon_ParallelJobs::p, x, y, z, PV3DBase< T, PVType, FrameType >::z(), zeN, zeP, and geometryCSVtoXML::zz.

Referenced by setNumberOfCrystalPerModule().

{
  zeP=0.;
  zeN=0.;
  unsigned nP=0;
  unsigned nN=0;
  m_nref = 0 ;

  for( uint32_t i ( 0 ) ; i != m_cellVec.size() ; ++i )
  {
     const CaloCellGeometry* cell ( cellGeomPtr(i) ) ;
     if( nullptr != cell )
     {
    const CCGFloat z ( cell->getPosition().z() ) ;
    if(z>0.)
    {
       zeP+=z;
       ++nP;
    }
    else
    {
       zeN+=z;
       ++nN;
    }
    const EEDetId myId ( EEDetId::detIdFromDenseIndex(i) ) ;
    const unsigned int ix ( myId.ix() ) ;
    const unsigned int iy ( myId.iy() ) ;
    if( ix > m_nref ) m_nref = ix ;
    if( iy > m_nref ) m_nref = iy ;
     }
  }
  if( 0 < nP ) zeP/=(CCGFloat)nP;
  if( 0 < nN ) zeN/=(CCGFloat)nN;

  m_xlo[0] =  999 ;
  m_xhi[0] = -999 ;
  m_ylo[0] =  999 ;
  m_yhi[0] = -999 ;
  m_xlo[1] =  999 ;
  m_xhi[1] = -999 ;
  m_ylo[1] =  999 ;
  m_yhi[1] = -999 ;
  for( uint32_t i ( 0 ) ; i != m_cellVec.size() ; ++i )
  {
     const CaloCellGeometry* cell ( cellGeomPtr(i) ) ;
     if( nullptr != cell )
     {
    const GlobalPoint& p ( cell->getPosition()  ) ;
    const CCGFloat z ( p.z() ) ;
    const CCGFloat zz ( 0 > z ? zeN : zeP ) ;
    const CCGFloat x ( p.x()*zz/z ) ;
    const CCGFloat y ( p.y()*zz/z ) ;

    if( 0 > z && x < m_xlo[0] ) m_xlo[0] = x ;
    if( 0 < z && x < m_xlo[1] ) m_xlo[1] = x ;
    if( 0 > z && y < m_ylo[0] ) m_ylo[0] = y ;
    if( 0 < z && y < m_ylo[1] ) m_ylo[1] = y ;
     
    if( 0 > z && x > m_xhi[0] ) m_xhi[0] = x ;
    if( 0 < z && x > m_xhi[1] ) m_xhi[1] = x ;
    if( 0 > z && y > m_yhi[0] ) m_yhi[0] = y ;
    if( 0 < z && y > m_yhi[1] ) m_yhi[1] = y ;
     }
  }

  m_xoff[0] = ( m_xhi[0] + m_xlo[0] )/2. ;
  m_xoff[1] = ( m_xhi[1] + m_xlo[1] )/2. ;
  m_yoff[0] = ( m_yhi[0] + m_ylo[0] )/2. ;
  m_yoff[1] = ( m_yhi[1] + m_ylo[1] )/2. ;

  m_del = ( m_xhi[0] - m_xlo[0] + m_xhi[1] - m_xlo[1] +
        m_yhi[0] - m_ylo[0] + m_yhi[1] - m_ylo[1]   ) ;

  if( 1 != m_nref ) m_wref = m_del/(4.*(m_nref-1)) ;

  m_xlo[0] -= m_wref/2 ;
  m_xlo[1] -= m_wref/2 ;
  m_xhi[0] += m_wref/2 ;
  m_xhi[1] += m_wref/2 ;

  m_ylo[0] -= m_wref/2 ;
  m_ylo[1] -= m_wref/2 ;
  m_yhi[0] += m_wref/2 ;
  m_yhi[1] += m_wref/2 ;

  m_del += m_wref ;
/*
  std::cout<<"zeP="<<zeP<<", zeN="<<zeN<<", nP="<<nP<<", nN="<<nN<<std::endl ;

  std::cout<<"xlo[0]="<<m_xlo[0]<<", xlo[1]="<<m_xlo[1]<<", xhi[0]="<<m_xhi[0]<<", xhi[1]="<<m_xhi[1]
       <<"\nylo[0]="<<m_ylo[0]<<", ylo[1]="<<m_ylo[1]<<", yhi[0]="<<m_yhi[0]<<", yhi[1]="<<m_yhi[1]<<std::endl ;

  std::cout<<"xoff[0]="<<m_xoff[0]<<", xoff[1]"<<m_xoff[1]<<", yoff[0]="<<m_yoff[0]<<", yoff[1]"<<m_yoff[1]<<std::endl ;

  std::cout<<"nref="<<m_nref<<", m_wref="<<m_wref<<std::endl ;
*/  
}

void EcalEndcapGeometry::localCorners ( Pt3DVec &  lc, const CCGFloat *  pv, unsigned int  i, Pt3D &  ref  )

static

Definition at line 477 of file EcalEndcapGeometry.cc.

References TruncatedPyramid::localCorners().

Referenced by numberOfAlignments().

{
   TruncatedPyramid::localCorners( lc, pv, ref ) ;
}

static int EcalEndcapGeometry::myPhi ( int  i )

inlinestaticprivate

Definition at line 107 of file EcalEndcapGeometry.h.

Referenced by getClosestBarrelCells().

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

void EcalEndcapGeometry::newCell ( const GlobalPoint &  f1, const GlobalPoint &  f2, const GlobalPoint &  f3, const CCGFloat *  parm, const DetId &  detId  )

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 486 of file EcalEndcapGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), CaloSubdetectorGeometry::cornersMgr(), and m_cellVec.

Referenced by numberOfAlignments().

{
   const unsigned int cellIndex ( EEDetId( detId ).denseIndex() ) ;
   m_cellVec[ cellIndex ] =
      TruncatedPyramid( cornersMgr(), f1, f2, f3, parm ) ;
   addValidID( detId ) ;
}

static unsigned int EcalEndcapGeometry::numberOfAlignments ( )

inlinestatic

Definition at line 82 of file EcalEndcapGeometry.h.

References alignmentTransformIndexGlobal(), alignmentTransformIndexLocal(), cellGeomPtr(), detIdFromLocalAlignmentIndex(), connectstrParser::f1, connectstrParser::f2, connectstrParser::f3, mps_fire::i, diffTreeTool::index, localCorners(), newCell(), and MetAnalyzer::pv().

Referenced by FakeCaloAlignmentEP::produceEEAli(), and TestCaloAlignmentEP::produceEEAli().

{ return 4 ; }

unsigned int EcalEndcapGeometry::numberOfParametersPerShape ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 53 of file EcalEndcapGeometry.h.

References EcalEndcapGeometry(), k_NumberOfParametersPerShape, and ~EcalEndcapGeometry().

{ return k_NumberOfParametersPerShape ; }

unsigned int EcalEndcapGeometry::numberOfShapes ( ) const

inlineoverridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 52 of file EcalEndcapGeometry.h.

References k_NumberOfShapes.

{ return k_NumberOfShapes ; }

static std::string EcalEndcapGeometry::producerTag ( )

inlinestatic

Definition at line 80 of file EcalEndcapGeometry.h.

Referenced by PCaloGeometryBuilder::beginRun(), and CaloGeometryBuilder::produceAligned().

{ return "EcalEndcap" ; }

void EcalEndcapGeometry::setNumberOfCrystalPerModule ( const int  nncrys )

inline

Definition at line 65 of file EcalEndcapGeometry.h.

References _nncrys, avgAbsZFrontFaceCenter(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, getCells(), getClosestBarrelCells(), getClosestCell(), initializeParms(), and alignCSCRings::r.

{ _nncrys=nncrys ; }

void EcalEndcapGeometry::setNumberOfModules ( const int  nnmods )

inline

Definition at line 63 of file EcalEndcapGeometry.h.

References _nnmods.

{ _nnmods=nnmods ; }

unsigned int EcalEndcapGeometry::xindex ( CCGFloat  x, CCGFloat  z  ) const

private

Definition at line 180 of file EcalEndcapGeometry.cc.

References mps_fire::i, m_nref, m_wref, and m_xlo.

Referenced by getCells(), and gId().

{
   const CCGFloat xlo ( 0 > z ? m_xlo[0]  : m_xlo[1]  ) ;
   const int i ( 1 + int( ( x - xlo )/m_wref ) ) ;

   return ( 1 > i ? 1 :
        ( m_nref < (unsigned int) i ? m_nref : (unsigned int) i ) ) ;

}

unsigned int EcalEndcapGeometry::yindex ( CCGFloat  y, CCGFloat  z  ) const

private

Definition at line 192 of file EcalEndcapGeometry.cc.

References mps_fire::i, m_nref, m_wref, and m_ylo.

Referenced by getCells(), and gId().

{
   const CCGFloat ylo ( 0 > z ? m_ylo[0]  : m_ylo[1]  ) ;
   const int i ( 1 + int( ( y - ylo )/m_wref ) ) ;

   return ( 1 > i ? 1 :
        ( m_nref < (unsigned int) i ? m_nref : (unsigned int) i ) ) ;
}

Member Data Documentation

int EcalEndcapGeometry::_nncrys

private

number of crystals per module

Definition at line 113 of file EcalEndcapGeometry.h.

Referenced by getNumberOfCrystalPerModule(), and setNumberOfCrystalPerModule().

int EcalEndcapGeometry::_nnmods

private

number of modules

Definition at line 110 of file EcalEndcapGeometry.h.

Referenced by getNumberOfModules(), and setNumberOfModules().

CCGFloat EcalEndcapGeometry::m_avgZ

mutableprivate

Definition at line 130 of file EcalEndcapGeometry.h.

Referenced by avgAbsZFrontFaceCenter().

std::atomic<EZMgrFL<EBDetId>*> EcalEndcapGeometry::m_borderMgr

mutableprivate

Definition at line 126 of file EcalEndcapGeometry.h.

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

std::atomic<VecOrdListEBDetIdPtr*> EcalEndcapGeometry::m_borderPtrVec

mutableprivate

Definition at line 128 of file EcalEndcapGeometry.h.

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

CellVec EcalEndcapGeometry::m_cellVec

private

Definition at line 133 of file EcalEndcapGeometry.h.

Referenced by avgAbsZFrontFaceCenter(), cellGeomPtr(), getCells(), initializeParms(), and newCell().

std::atomic<bool> EcalEndcapGeometry::m_check

mutableprivate

Definition at line 131 of file EcalEndcapGeometry.h.

Referenced by avgAbsZFrontFaceCenter().

CCGFloat EcalEndcapGeometry::m_del

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by initializeParms().

unsigned int EcalEndcapGeometry::m_nref

private

Definition at line 119 of file EcalEndcapGeometry.h.

Referenced by getCells(), initializeParms(), xindex(), and yindex().

CCGFloat EcalEndcapGeometry::m_wref

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by getCells(), initializeParms(), xindex(), and yindex().

CCGFloat EcalEndcapGeometry::m_xhi[2]

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), and initializeParms().

CCGFloat EcalEndcapGeometry::m_xlo[2]

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), initializeParms(), and xindex().

CCGFloat EcalEndcapGeometry::m_xoff[2]

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), and initializeParms().

CCGFloat EcalEndcapGeometry::m_yhi[2]

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), and initializeParms().

CCGFloat EcalEndcapGeometry::m_ylo[2]

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), getCells(), initializeParms(), and yindex().

CCGFloat EcalEndcapGeometry::m_yoff[2]

private

Definition at line 117 of file EcalEndcapGeometry.h.

Referenced by EcalEndcapGeometry(), and initializeParms().

CCGFloat EcalEndcapGeometry::zeN

private

Definition at line 115 of file EcalEndcapGeometry.h.

Referenced by gId(), and initializeParms().

CCGFloat EcalEndcapGeometry::zeP

private

Definition at line 115 of file EcalEndcapGeometry.h.

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