TruncatedPyramid.h

Go to the documentation of this file.

#ifndef TruncatedPyramid_h
#define TruncatedPyramid_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 <vector>
 
class TruncatedPyramid final : 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_Theta = 1;  //Polar angle of the line joining the
                                          //centres of the faces at -/+ Dz
  static constexpr uint32_t k_Phi = 2;    //Azimuthal angle of the line joing the
                                          //centres of the faces at -/+ Dz
  static constexpr uint32_t k_Dy1 = 3;    //Half-length along y of the face at -Dz
  static constexpr uint32_t k_Dx1 = 4;    //Half-length along x of the side at
                                          //y=-Dy1 of the face at -Dz
  static constexpr uint32_t k_Dx2 = 5;    //Half-length along x of the side at
                                          //y=+Dy1 of the face at -Dz
  static constexpr uint32_t k_Alp1 = 6;   //Angle w.r.t the y axis from the center
                                          //of the sides at y=-Dy1 to at y=+Dy1
  static constexpr uint32_t k_Dy2 = 7;    //Half-length along y of the face at +Dz
  static constexpr uint32_t k_Dx3 = 8;    //Half-length along x of the side at
                                          //y=-Dy2 of the face at +Dz
  static constexpr uint32_t k_Dx4 = 9;    //Half-length along x of the side at
                                          //y=+Dy2 of the face at +Dz
  static constexpr uint32_t k_Alp2 = 10;  //Angle w.r.t the y axis from the center
                                          //of the sides at y=-Dy2 to at y=+Dy2
 
  TruncatedPyramid(void);
 
  TruncatedPyramid(const TruncatedPyramid& tr);
 
  TruncatedPyramid& operator=(const TruncatedPyramid& tr);
 
  TruncatedPyramid(CornersMgr* cMgr,
                   const GlobalPoint& fCtr,
                   const GlobalPoint& bCtr,
                   const GlobalPoint& cor1,
                   const CCGFloat* parV);
 
  TruncatedPyramid(const CornersVec& corn, const CCGFloat* par);
 
  ~TruncatedPyramid() override;
 
  GlobalPoint getPosition(CCGFloat depth) const override;
 
  // Return thetaAxis polar angle of axis of the crystal
  CCGFloat getThetaAxis() const;
 
  // Return phiAxis azimuthal angle of axis of the crystal
  CCGFloat getPhiAxis() const;
 
  const GlobalVector& axis() const;
 
  // for geometry creation in other classes
  static void createCorners(const std::vector<CCGFloat>& pv, const Tr3D& tr, std::vector<GlobalPoint>& co);
 
  void vocalCorners(Pt3DVec& vec, const CCGFloat* pv, Pt3D& ref) const override;
 
  static void localCorners(Pt3DVec& vec, const CCGFloat* pv, Pt3D& ref);
 
  static void localCornersReflection(Pt3DVec& vec, const CCGFloat* pv, Pt3D& ref);
 
  static void localCornersSwap(Pt3DVec& vec, const CCGFloat* pv, Pt3D& ref);
 
  void getTransform(Tr3D& tr, Pt3DVec* lptr) const override;
 
private:
  void initCorners(CornersVec&) override;
 
  GlobalVector makeAxis(void);
 
  const GlobalPoint backCtr(void) const;
  GlobalVector m_axis;
  Pt3D m_corOne;
};
 
std::ostream& operator<<(std::ostream& s, const TruncatedPyramid& cell);
 
#endif