CaloCellGeometry.h

Go to the documentation of this file.

#ifndef GEOMETRY_CALOGEOMETRY_CALOCELLGEOMETRY_H
#define GEOMETRY_CALOGEOMETRY_CALOCELLGEOMETRY_H 1

#include "Geometry/CaloGeometry/interface/EZArrayFL.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include <CLHEP/Geometry/Point3D.h>
#include <CLHEP/Geometry/Transform3D.h>
#include <vector>
#include <string>

#include "FWCore/Utilities/interface/GCC11Compatibility.h"



class CaloCellGeometry 
{
public:

  typedef float                      CCGFloat ;
  typedef HepGeom::Transform3D       Tr3D     ;
  typedef HepGeom::Point3D<CCGFloat> Pt3D     ;
  typedef std::vector<Pt3D>          Pt3DVec  ;

  typedef EZArrayFL< GlobalPoint > CornersVec ;
  typedef EZMgrFL< GlobalPoint >   CornersMgr ;

  typedef EZArrayFL<CCGFloat> ParVec ;
  typedef std::vector<ParVec> ParVecVec ;
  typedef EZMgrFL< CCGFloat > ParMgr ;

  enum CornersSize { k_cornerSize = 8 };

  static const CCGFloat k_ScaleFromDDDtoGeant ;

  virtual ~CaloCellGeometry() ;
      
  virtual const CornersVec& getCorners() const = 0 ;

  const GlobalPoint& getPosition() const {return m_refPoint ; }
  float etaPos() const { return m_eta;}
  float phiPos() const { return m_phi;}

  float etaSpan() const { return m_dEta;}
  float phiSpan() const { return m_dPhi;}


  bool inside( const GlobalPoint & point ) const ;  

  bool emptyCorners() const { return m_corners.empty() ;}

  const CCGFloat* param() const { return m_parms ;}

  static const CCGFloat* checkParmPtr( const std::vector<CCGFloat>& vd  ,
                       ParVecVec&                   pvv   ) ;

  static const CCGFloat* getParmPtr( const std::vector<CCGFloat>& vd  ,
                     ParMgr*                      mgr ,
                     ParVecVec&                   pvv   ) ;


  virtual void getTransform( Tr3D& tr, Pt3DVec* lptr ) const ;
  //------------------------------------------------------------------------

  virtual void vocalCorners( Pt3DVec&        vec ,
                 const CCGFloat* pv  ,
                 Pt3D&           ref  ) const = 0 ;

protected:

  CaloCellGeometry( CornersVec::const_reference gp ,
            const CornersMgr*           mgr,
            const CCGFloat*             par ) ;

  CaloCellGeometry( const CornersVec& cv,
            const CCGFloat*   par ) ;

  CornersVec& setCorners() const ;

  CaloCellGeometry( void );

  // MUST be called by children constructors
  void initSpan() const {
     m_dEta = std::abs(getCorners()[0].eta()-
                      getCorners()[2].eta());
     m_dPhi = std::abs(getCorners()[0].phi() -
                      getCorners()[2].phi());
  }

private:
  GlobalPoint         m_refPoint ;
  mutable CornersVec  m_corners  ;
  const CCGFloat*     m_parms    ;
  float m_eta, m_phi;
  mutable float m_dEta;
  mutable float m_dPhi;

};

std::ostream& operator<<( std::ostream& s, const CaloCellGeometry& cell ) ;

#endif