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/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 <vector>

class EcalBarrelGeometry : public CaloSubdetectorGeometry 
{
   public:

      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 ;

      double 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 std::vector<HepGeom::Point3D<double> > localCorners( const double* pv, 
                           unsigned int  i,
                           HepGeom::Point3D<double> &   ref ) ;

      static CaloCellGeometry* newCell( const GlobalPoint& f1 ,
                    const GlobalPoint& f2 ,
                    const GlobalPoint& f3 ,
                    CaloCellGeometry::CornersMgr* mgr,
                    const double*      parm ,
                    const DetId&       detId ) ;
                    
   private:
      int _nnxtalEta;
  
      int _nnxtalPhi;
  
      std::vector<int> _EtaBaskets;
      
      int _PhiBaskets;

      mutable EZMgrFL<EEDetId>*     m_borderMgr ;

      mutable VecOrdListEEDetIdPtr* m_borderPtrVec ;

      mutable double m_radius ;
};


#endif