FlatHexagon.h

Go to the documentation of this file.

#ifndef GeometryCaloGeometryFlatHexagon_h
#define GeometryCaloGeometryFlatHexagon_h

#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include <CLHEP/Geometry/Point3D.h>
#include <CLHEP/Geometry/Plane3D.h>
#include <CLHEP/Geometry/Vector3D.h>
#include <CLHEP/Geometry/Transform3D.h>
#include <cmath>
#include <vector>

class FlatHexagon : public CaloCellGeometry {
public:

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

  static constexpr uint32_t k_dZ   = 0;//Half-length along the z-axis
  static constexpr uint32_t k_r    = 1;//Half length along x-axis
  static constexpr uint32_t k_R    = 2;//Quarter of maximum dimension along y
  
  FlatHexagon( void );
  
  FlatHexagon( const FlatHexagon& tr ) ;
  
  FlatHexagon& operator=( const FlatHexagon& tr ) ;
  
  FlatHexagon( CornersMgr*  cMgr ,
           const GlobalPoint& fCtr ,
           const GlobalPoint& bCtr ,
           const GlobalPoint& cor1 ,
           const CCGFloat*    parV   ) ;
  
  FlatHexagon( const CornersVec& corn ,
           const CCGFloat*   par    ) ;
  
  FlatHexagon( const FlatHexagon& tr, const Pt3D& local) ;

  ~FlatHexagon() override ;
  
  GlobalPoint const & getPosition() const override { return m_global; }
  GlobalPoint getPosition( const Pt3D& local ) const override;
  virtual float etaPos() const { return m_global.eta(); }
  virtual float phiPos() const { return m_global.phi(); }
  Pt3D getLocal( const GlobalPoint& global ) const;

  // Return thetaAxis polar angle of axis of the crystal
  CCGFloat getThetaAxis() const ;
  
  // Return phiAxis azimuthal angle of axis of the crystal
  CCGFloat getPhiAxis() const ;
  
  void vocalCorners( Pt3DVec&        vec ,
             const CCGFloat* pv  ,
             Pt3D&           ref  ) const override;
  
  const GlobalVector& axis() const ;
  
  static void createCorners( const std::vector<CCGFloat>& pv ,
                             const Tr3D&                  tr ,
                             std::vector<GlobalPoint>&    co   ) ;
  
  static void localCorners( Pt3DVec&        vec ,
                            const CCGFloat* pv  ,
                            Pt3D&           ref  ) ;
  
  void getTransform( Tr3D& tr, Pt3DVec* lptr ) const override;
  
  void setPosition ( const GlobalPoint& p ) { m_global = p;  setRefPoint(p); }

  static constexpr double       oneBySix_   = 1.0/6.0;
  static constexpr unsigned int ncorner_    = 12;
  static constexpr unsigned int ncornerBy2_ = 6;

private:
  
  void initCorners(CornersVec& ) override;
  
  GlobalVector makeAxis( void );
  
  GlobalPoint backCtr( void ) const;    
  GlobalVector m_axis;
  Pt3D         m_corOne, m_local;
  GlobalPoint  m_global;
  Tr3D         m_tr;
};

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

#endif