EcalBarrelGeometry.h

Go to the documentation of this file.

#ifndef EcalBarrelGeometry_h
#define EcalBarrelGeometry_h

#include "Geometry/CaloGeometry/interface/EZArrayFL.h"
#include "Geometry/EcalCommonData/interface/EcalBarrelNumberingScheme.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/TruncatedPyramid.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/Records/interface/EcalBarrelGeometryRecord.h"
#include "CondFormats/AlignmentRecord/interface/EBAlignmentRcd.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "Geometry/Records/interface/PEcalBarrelRcd.h"
#include "FWCore/Utilities/interface/thread_safety_macros.h"
#include <vector>
#include <atomic>

class EcalBarrelGeometry final : public CaloSubdetectorGeometry 
{
   public:

      typedef std::vector<TruncatedPyramid> CellVec ;

      typedef CaloCellGeometry::CCGFloat CCGFloat ;
      typedef CaloCellGeometry::Pt3D     Pt3D     ;
      typedef CaloCellGeometry::Pt3DVec  Pt3DVec  ;

      typedef IdealGeometryRecord      IdealRecord   ;
      typedef EcalBarrelGeometryRecord AlignedRecord ;
      typedef EBAlignmentRcd           AlignmentRecord ;
      typedef PEcalBarrelRcd           PGeometryRecord ;

      typedef EZArrayFL<EEDetId> OrderedListOfEEDetId ; // like an stl vector: begin(), end(), [i]

      typedef std::vector<OrderedListOfEEDetId*>  VecOrdListEEDetIdPtr ;

      typedef EcalBarrelNumberingScheme NumberingScheme ;

      typedef EBDetId DetIdType ;

      enum { k_NumberOfCellsForCorners = EBDetId::kSizeForDenseIndexing } ;

      enum { k_NumberOfShapes = 17 } ;

      enum { k_NumberOfParametersPerShape = 11 } ;

      static std::string dbString() { return "PEcalBarrelRcd" ; }

      virtual unsigned int numberOfShapes() const { return k_NumberOfShapes ; }
      virtual unsigned int numberOfParametersPerShape() const { return k_NumberOfParametersPerShape ; }

      EcalBarrelGeometry() ;
  
      virtual ~EcalBarrelGeometry();

      int getNumXtalsPhiDirection()           const { return _nnxtalPhi ; }

      int getNumXtalsEtaDirection()           const { return _nnxtalEta ; }

      const std::vector<int>& getEtaBaskets() const { return _EtaBaskets ; }

      int getBasketSizeInPhi()                const { return _PhiBaskets ; }  

      void setNumXtalsPhiDirection( const int& nnxtalPhi )     { _nnxtalPhi=nnxtalPhi ; }

      void setNumXtalsEtaDirection( const int& nnxtalEta )     { _nnxtalEta=nnxtalEta ; }

      void setEtaBaskets( const std::vector<int>& EtaBaskets ) { _EtaBaskets=EtaBaskets ; }

      void setBasketSizeInPhi( const int& PhiBaskets )         { _PhiBaskets=PhiBaskets ; }

      const OrderedListOfEEDetId* getClosestEndcapCells( EBDetId id ) const ;

      // Get closest cell, etc...
      virtual DetId getClosestCell( const GlobalPoint& r ) const ;

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

      CCGFloat avgRadiusXYFrontFaceCenter() const ;

      static std::string hitString() { return "EcalHitsEB" ; }

      static std::string producerTag() { return "EcalBarrel" ; }

      static unsigned int numberOfAlignments() { return 36 ; }

      static unsigned int alignmentTransformIndexLocal( const DetId& id ) ;

      static unsigned int alignmentTransformIndexGlobal( const DetId& id ) ;

      static DetId detIdFromLocalAlignmentIndex( unsigned int iLoc ) ;

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

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

      virtual const CaloCellGeometry* cellGeomPtr( uint32_t index ) const ;

   private:

      int _nnxtalEta;
  
      int _nnxtalPhi;
  
      std::vector<int> _EtaBaskets;
      
      int _PhiBaskets;

      mutable std::atomic<EZMgrFL<EEDetId>*>     m_borderMgr ;

      mutable std::atomic<VecOrdListEEDetIdPtr*> m_borderPtrVec ;
      CMS_THREAD_GUARD(m_check) mutable CCGFloat m_radius ; 
      mutable std::atomic<bool> m_check;

      CellVec m_cellVec ;
};


#endif