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