CaloSubdetectorGeometry.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOGEOMETRY_CALOSUBDETECTORGEOMETRY_H
#define GEOMETRY_CALOGEOMETRY_CALOSUBDETECTORGEOMETRY_H 1

#include <vector>
#include <set>
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "DataFormats/Math/interface/deltaR.h"

class CaloSubdetectorGeometry {

   public:

      typedef  std::vector< const CaloCellGeometry * > CellCont;

      typedef std::set<DetId>       DetIdSet;

      typedef CaloCellGeometry::ParMgr    ParMgr ;
      typedef CaloCellGeometry::ParVec    ParVec ;
      typedef CaloCellGeometry::ParVecVec ParVecVec ;

      typedef std::vector<double>         TrVec  ;
      typedef std::vector<unsigned int>   IVec   ;
      typedef std::vector<double>         DimVec ;

      CaloSubdetectorGeometry() ;

      virtual ~CaloSubdetectorGeometry();

      const CellCont& cellGeometries() const { return m_cellG ; }  

      void addCell( const DetId&      id, 
            CaloCellGeometry* ccg ) ;

      virtual bool present( const DetId& id ) const;

      virtual const CaloCellGeometry* getGeometry( const DetId& id ) const ;

      virtual const std::vector<DetId>& getValidDetIds( DetId::Detector det    = DetId::Detector(0) , 
                            int             subdet = 0                   ) const ;

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

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


      double deltaPhi( const DetId& detId ) const ;

      double deltaEta( const DetId& detId ) const ;



      void allocateCorners( CaloCellGeometry::CornersVec::size_type n ) ;

      CaloCellGeometry::CornersMgr* cornersMgr() { return m_cmgr ; }

      void allocatePar( ParVec::size_type n, unsigned int m ) ;

      ParMgr* parMgr() { return m_parMgr ; }
      const ParMgr* parMgrConst() const { return m_parMgr ; }

      ParVecVec&       parVecVec()       { return m_parVecVec ; }
      const ParVecVec& parVecVec() const { return m_parVecVec ; }
      
      virtual unsigned int numberOfShapes() const { return 1 ; }
      virtual unsigned int numberOfParametersPerShape() const { return 1 ; }
      virtual unsigned int numberOfTransformParms() const { return 6 ; }

      virtual void fillDefaultNamedParameters() const { return ; }
      
      void getSummary( TrVec&  trVector,
               IVec&   iVector,
               DimVec& dimVector   ) const ;

      virtual void initializeParms() { return ; } 

   protected:

      ParVecVec m_parVecVec ;

      static double deltaR( const GlobalPoint& p1,
                const GlobalPoint& p2  ) 
      { return reco::deltaR( p1, p2 ) ; }

   private:

      ParMgr*   m_parMgr ;

      CaloCellGeometry::CornersMgr* m_cmgr ;

      CaloSubdetectorGeometry(            const CaloSubdetectorGeometry& ) ;
      CaloSubdetectorGeometry& operator=( const CaloSubdetectorGeometry& ) ;

      CellCont m_cellG ;    

      mutable bool m_sortedIds ;

      mutable std::vector<DetId> m_validIds ;

      mutable std::vector<double>*  m_deltaPhi ;
      mutable std::vector<double>*  m_deltaEta ;
};


#endif