EcalPreshowerGeometry.h

Go to the documentation of this file.

#ifndef EcalPreshowerGeometry_h
#define EcalPreshowerGeometry_h

#include "Geometry/EcalCommonData/interface/EcalPreshowerNumberingScheme.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/PreshowerStrip.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/Records/interface/EcalPreshowerGeometryRecord.h"
#include "CondFormats/AlignmentRecord/interface/ESAlignmentRcd.h"
#include "DataFormats/EcalDetId/interface/ESDetId.h"
#include "Geometry/Records/interface/PEcalPreshowerRcd.h"
#include "Geometry/CaloGeometry/interface/CaloGenericDetId.h"
#include <vector>

class EcalPreshowerGeometry final : public CaloSubdetectorGeometry {
public:
  typedef std::vector<PreshowerStrip> CellVec;

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

  typedef IdealGeometryRecord IdealRecord;
  typedef EcalPreshowerGeometryRecord AlignedRecord;
  typedef ESAlignmentRcd AlignmentRecord;
  typedef PEcalPreshowerRcd PGeometryRecord;

  typedef EcalPreshowerNumberingScheme NumberingScheme;
  typedef CaloSubdetectorGeometry::ParVec ParVec;
  typedef CaloSubdetectorGeometry::ParVecVec ParVecVec;
  typedef ESDetId DetIdType;

  enum { k_NumberOfCellsForCorners = ESDetId::kSizeForDenseIndexing };

  enum { k_NumberOfShapes = 4 };

  enum { k_NumberOfParametersPerShape = 4 };

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

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

  EcalPreshowerGeometry();

  ~EcalPreshowerGeometry() override;

  void setzPlanes(CCGFloat z1minus, CCGFloat z2minus, CCGFloat z1plus, CCGFloat z2plus);

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

  // Get closest cell in arbitrary plane (1 or 2)
  virtual DetId getClosestCellInPlane(const GlobalPoint& r, int plane) const;

  void initializeParms() override;
  unsigned int numberOfTransformParms() const override { return 3; }

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

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

  static unsigned int numberOfAlignments() { return 8; }

  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 {
    if (id == DetId(0))
      return false;
    // not needed???
    auto index = CaloGenericDetId(id).denseIndex();
    return index < m_cellVec.size();
  }

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

private:
  const CCGFloat m_xWidWaf;
  const CCGFloat m_xInterLadGap;
  const CCGFloat m_xIntraLadGap;

  const CCGFloat m_yWidAct;
  const CCGFloat m_yCtrOff;

  CCGFloat m_zplane[4];

  CellVec m_cellVec;
};

#endif