---

CaloSubdetectorGeometry.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOGEOMETRY_CALOSUBDETECTORGEOMETRY_H
#define GEOMETRY_CALOGEOMETRY_CALOSUBDETECTORGEOMETRY_H 1

#include <ext/hash_map>
#include <vector>
#include <set>
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "PhysicsTools/Utilities/interface/deltaR.h"

class CaloSubdetectorGeometry {

   public:

      typedef  __gnu_cxx::hash_map< unsigned int, const CaloCellGeometry *> CellCont;

      typedef std::set<DetId>       DetIdSet;

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

      CaloSubdetectorGeometry() : 
         m_parMgr ( 0 ) ,
         m_cmgr   ( 0 )   {}

      virtual ~CaloSubdetectorGeometry();

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

      void addCell( const DetId& id, 
                    const 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, 
                                                        int subdet  ) const ;

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

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

      //FIXME: Hcal implements its own  getValidDetId....

      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 ; }

      ParVecVec& parVecVec() { return m_parVecVec ; }

   protected:

      ParVecVec m_parVecVec ;

      mutable std::vector<DetId> m_validIds ;

      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 ;    

};


#endif

---