HcalDDDGeometry.h

Go to the documentation of this file.

#ifndef Geometry_HcalTowerAlgo_HcalDDDGeometry_h
#define Geometry_HcalTowerAlgo_HcalDDDGeometry_h

#include "Geometry/HcalCommonData/interface/HcalCellType.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"
#include "FWCore/Utilities/interface/thread_safety_macros.h"
#include <atomic>

class HcalDDDGeometryLoader;

class HcalDDDGeometry : public CaloSubdetectorGeometry {

public:

  friend class HcalDDDGeometryLoader;

  typedef std::vector<IdealObliquePrism> HBCellVec ;
  typedef std::vector<IdealObliquePrism> HECellVec ;
  typedef std::vector<IdealObliquePrism> HOCellVec ;
  typedef std::vector<IdealZPrism>       HFCellVec ;

  explicit HcalDDDGeometry(const HcalTopology& theTopo);
  ~HcalDDDGeometry() override;
  
  const std::vector<DetId>& getValidDetIds( DetId::Detector det    = DetId::Detector ( 0 ) , 
                            int             subdet = 0   ) const override;

  DetId getClosestCell(const GlobalPoint& r) const override ;

  int insertCell (std::vector<HcalCellType> const & );

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

protected:

  // Modify the RawPtr class
  const CaloCellGeometry* getGeometryRawPtr (uint32_t index) const override;

private:

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

  std::vector<HcalCellType> hcalCells_;
  CMS_THREAD_GUARD(m_filledDetIds) mutable std::vector<DetId> m_hbIds ;
  CMS_THREAD_GUARD(m_filledDetIds) mutable std::vector<DetId> m_heIds ;
  CMS_THREAD_GUARD(m_filledDetIds) mutable std::vector<DetId> m_hoIds ;
  CMS_THREAD_GUARD(m_filledDetIds) mutable std::vector<DetId> m_hfIds ;
  CMS_THREAD_GUARD(m_filledDetIds) mutable std::vector<DetId> m_emptyIds ;

  const HcalTopology& topo_;
  double              etaMax_;

  HBCellVec m_hbCellVec ;
  HECellVec m_heCellVec ;
  HOCellVec m_hoCellVec ;
  HFCellVec m_hfCellVec ;
  mutable std::atomic<bool> m_filledDetIds;
};

#endif