HcalGeometry.h

Go to the documentation of this file.

#ifndef HcalGeometry_h
#define HcalGeometry_h

#include "DataFormats/Common/interface/AtomicPtrCache.h"

#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/IdealObliquePrism.h"
#include "Geometry/CaloGeometry/interface/IdealZPrism.h"
#include "Geometry/CaloTopology/interface/HcalTopology.h"
#include "CondFormats/AlignmentRecord/interface/HcalAlignmentRcd.h"
#include "Geometry/Records/interface/HcalGeometryRecord.h"

class HcalFlexiHardcodeGeometryLoader;
class HcalHardcodeGeometryLoader;

class HcalGeometry : public CaloSubdetectorGeometry {

public:

  friend class HcalFlexiHardcodeGeometryLoader;
  friend class HcalHardcodeGeometryLoader;
  
  typedef std::vector<IdealObliquePrism> HBCellVec ;
  typedef std::vector<IdealObliquePrism> HECellVec ;
  typedef std::vector<IdealObliquePrism> HOCellVec ;
  typedef std::vector<IdealZPrism>       HFCellVec ;

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

  typedef HcalAlignmentRcd   AlignmentRecord ;
  typedef HcalGeometryRecord AlignedRecord   ;
  typedef PHcalRcd           PGeometryRecord ;
  typedef HcalDetId          DetIdType       ;

    
  enum { k_NumberOfParametersPerShape = 5 } ;

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

  virtual unsigned int numberOfShapes() const { return m_topology.getNumberOfShapes() ; }
  virtual unsigned int numberOfParametersPerShape() const { return k_NumberOfParametersPerShape ; }

  explicit HcalGeometry(const HcalTopology& topology);

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

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

  GlobalPoint                   getPosition(const DetId& id) const;
  GlobalPoint                   getBackPosition(const DetId& id) const;
  CaloCellGeometry::CornersVec  getCorners(const DetId& id) const;

  static std::string producerTag() { return "HCAL" ; }
  
  static unsigned int numberOfBarrelAlignments() { return 36 ; }

  static unsigned int numberOfEndcapAlignments() { return 36 ; }

  static unsigned int numberOfForwardAlignments() { return 36 ; }

  static unsigned int numberOfOuterAlignments() { return 60 ; }

  unsigned int getHxSize(const int type) const;

  static unsigned int numberOfAlignments() 
    { return ( numberOfBarrelAlignments() +
           numberOfEndcapAlignments() +
           numberOfOuterAlignments() +
           numberOfForwardAlignments() ) ; }

  static unsigned int alignmentBarrelIndexLocal(    const DetId& id ) ;
  static unsigned int alignmentEndcapIndexLocal(    const DetId& id ) ;
  static unsigned int alignmentForwardIndexLocal(   const DetId& id ) ;
  static unsigned int alignmentOuterIndexLocal(     const DetId& id ) ;
  static unsigned int alignmentTransformIndexLocal( const DetId& id ) ;

  static unsigned int alignmentBarEndForIndexLocal( const DetId& id,unsigned int nD ) ;

  static unsigned int alignmentTransformIndexGlobal( const DetId& id ) ;

  static DetId detIdFromLocalAlignmentIndex(   unsigned int i ) ;
  static DetId detIdFromBarrelAlignmentIndex(  unsigned int i ) ;
  static DetId detIdFromEndcapAlignmentIndex(  unsigned int i ) ;
  static DetId detIdFromForwardAlignmentIndex( unsigned int i ) ;
  static DetId detIdFromOuterAlignmentIndex(   unsigned int i ) ;

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

  virtual const CaloCellGeometry* getGeometry( const DetId& id ) const {
      return cellGeomPtr( m_topology.detId2denseId( id ) ) ;
  }

  virtual void getSummary( CaloSubdetectorGeometry::TrVec&  trVector,
               CaloSubdetectorGeometry::IVec&   iVector,
               CaloSubdetectorGeometry::DimVec& dimVector,
               CaloSubdetectorGeometry::IVec& dinsVector ) const ;

  const HcalTopology& topology() const { return m_topology; }

protected:

  virtual const CaloCellGeometry* cellGeomPtr( unsigned int index ) const ;

  virtual unsigned int indexFor(const DetId& id) const { return  m_topology.detId2denseId(id); }
  virtual unsigned int sizeForDenseIndex(const DetId& id) const { return m_topology.ncells(); }

private:

  //returns din
  unsigned int newCellImpl( const GlobalPoint& f1 ,
            const GlobalPoint& f2 ,
            const GlobalPoint& f3 ,
            const CCGFloat*    parm,
            const DetId&       detId     ) ;

  //can only be used by friend classes, to ensure sorting is done at the end
  void newCellFast( const GlobalPoint& f1 ,
            const GlobalPoint& f2 ,
            const GlobalPoint& f3 ,
            const CCGFloat*    parm,
            const DetId&       detId     ) ;

  void increaseReserve(unsigned int extra);
  void sortValidIds();

  void fillDetIds() const ;

  void init() ;

  int etaRing(HcalSubdetector bc, double abseta) const;
  int phiBin(HcalSubdetector bc, int etaring, double phi) const;
  DetId correctId(const DetId& id) const ;

  const HcalTopology& m_topology;
  bool                m_mergePosition;
  
  mutable edm::AtomicPtrCache<std::vector<DetId>> m_hbIds ;
  mutable edm::AtomicPtrCache<std::vector<DetId>> m_heIds ;
  mutable edm::AtomicPtrCache<std::vector<DetId>> m_hoIds ;
  mutable edm::AtomicPtrCache<std::vector<DetId>> m_hfIds ;
  mutable edm::AtomicPtrCache<std::vector<DetId>> m_emptyIds ;
  CaloSubdetectorGeometry::IVec m_dins;

  HBCellVec m_hbCellVec ;
  HECellVec m_heCellVec ;
  HOCellVec m_hoCellVec ;
  HFCellVec m_hfCellVec ;
};


#endif