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 {
  friend class HcalGeometryPlan1Tester;

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

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

  explicit HcalGeometry(const HcalTopology& topology);

  ~HcalGeometry() override;

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

  std::shared_ptr<const CaloCellGeometry> getGeometry(const DetId& id) const override;

  DetId getClosestCell(const GlobalPoint& r) const override;
  DetId getClosestCell(const GlobalPoint& r, bool ignoreCorrect) const;

  CaloSubdetectorGeometry::DetIdSet getCells(const GlobalPoint& r, double dR) const override;

  GlobalPoint getPosition(const DetId& id) const;
  GlobalPoint getPosition(uint32_t id, bool type) const;
  GlobalPoint getBackPosition(const DetId& id) const;
  GlobalPoint getBackPosition(uint32_t id, bool type) 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);

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

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

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

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

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

private:
  // Base clas for getting geometry
  std::shared_ptr<const CaloCellGeometry> getGeometryBase(const DetId& id) const {
    return cellGeomPtr(m_topology.detId2denseId(id));
  }

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