EcalEndcapGeometry.h

Go to the documentation of this file.

#ifndef EcalEndcapGeometry_h
#define EcalEndcapGeometry_h

#include "Geometry/CaloGeometry/interface/EZArrayFL.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/Records/interface/EcalEndcapGeometryRecord.h"
#include "CondFormats/AlignmentRecord/interface/EEAlignmentRcd.h"
#include "Geometry/EcalCommonData/interface/EcalEndcapNumberingScheme.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/Records/interface/PEcalEndcapRcd.h"
#include "FWCore/Utilities/interface/thread_safety_macros.h"
#include <vector>
#include <map>
#include <atomic>

class TruncatedPyramid;

class EcalEndcapGeometry final : public CaloSubdetectorGeometry {
public:
  typedef std::vector<TruncatedPyramid> CellVec;

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

  typedef IdealGeometryRecord IdealRecord;
  typedef EcalEndcapGeometryRecord AlignedRecord;
  typedef EEAlignmentRcd AlignmentRecord;
  typedef PEcalEndcapRcd PGeometryRecord;

  typedef EZArrayFL<EBDetId> OrderedListOfEBDetId;  // like an stl vector: begin(), end(), [i]

  typedef std::vector<OrderedListOfEBDetId*> VecOrdListEBDetIdPtr;

  typedef EcalEndcapNumberingScheme NumberingScheme;

  typedef EEDetId DetIdType;

  enum { k_NumberOfCellsForCorners = EEDetId::kSizeForDenseIndexing };

  enum { k_NumberOfShapes = 1 };

  enum { k_NumberOfParametersPerShape = 11 };

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

  unsigned int numberOfShapes() const override { return k_NumberOfShapes; }
  unsigned int numberOfParametersPerShape() const override { return k_NumberOfParametersPerShape; }

  EcalEndcapGeometry();

  ~EcalEndcapGeometry() override;

  int getNumberOfModules() const { return _nnmods; }

  int getNumberOfCrystalPerModule() const { return _nncrys; }

  void setNumberOfModules(const int nnmods) { _nnmods = nnmods; }

  void setNumberOfCrystalPerModule(const int nncrys) { _nncrys = nncrys; }

  const OrderedListOfEBDetId* getClosestBarrelCells(EEDetId id) const;

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

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

  void initializeParms() override;

  CCGFloat avgAbsZFrontFaceCenter() const;  // average over both endcaps. Positive!

  static std::string hitString() { return "EcalHitsEE"; }

  static std::string producerTag() { return "EcalEndcap"; }

  static unsigned int numberOfAlignments() { return 4; }

  static unsigned int alignmentTransformIndexLocal(const DetId& id);

  static unsigned int alignmentTransformIndexGlobal(const DetId& id);

  static DetId detIdFromLocalAlignmentIndex(unsigned int iLoc);

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

  bool present(const DetId& id) const override;

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

private:
  static int myPhi(int i) {
    i += 720;
    return (1 + (i - 1) % 360);
  }

  int _nnmods;

  int _nncrys;

  CCGFloat zeP, zeN;

  CCGFloat m_wref, m_xlo[2], m_xhi[2], m_ylo[2], m_yhi[2], m_xoff[2], m_yoff[2], m_del;

  unsigned int m_nref;

  unsigned int xindex(CCGFloat x, CCGFloat z) const;
  unsigned int yindex(CCGFloat y, CCGFloat z) const;

  EEDetId gId(float x, float y, float z) const;

  mutable std::atomic<EZMgrFL<EBDetId>*> m_borderMgr;

  mutable std::atomic<VecOrdListEBDetIdPtr*> m_borderPtrVec;

  CMS_THREAD_GUARD(m_check) mutable CCGFloat m_avgZ;
  mutable std::atomic<bool> m_check;

  CellVec m_cellVec;
};

#endif