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