|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
#ifndef GeometryCaloGeometryFlatTrd_h
#define GeometryCaloGeometryFlatTrd_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 FlatTrd
\brief A base class to handle the particular shape of HGCal volumes.
*/
class FlatTrd : 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
static constexpr uint32_t k_Cell = 11; //Cell size
//Assumes dY2=dY1; dX3=dX1; dX4=dX2; Alp2=Alp1=Theta=Phi=0
FlatTrd(void);
FlatTrd(const FlatTrd& tr);
FlatTrd& operator=(const FlatTrd& tr);
FlatTrd(CornersMgr* cMgr,
const GlobalPoint& fCtr,
const GlobalPoint& bCtr,
const GlobalPoint& cor1,
const CCGFloat* parV);
FlatTrd(const CornersVec& corn, const CCGFloat* par);
FlatTrd(const FlatTrd& tr, const Pt3D& local);
~FlatTrd() override;
GlobalPoint const& getPosition() const override { return m_global; }
GlobalPoint getPosition(const Pt3D& local) const override;
float etaPos() const override { return m_global.eta(); }
float phiPos() const override { return m_global.phi(); }
float etaSpan() const override;
float phiSpan() const override;
float zPos() const { return m_global.z(); }
float dz() const { return param()[k_dZ]; }
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 unsigned int ncorner_ = 8;
static constexpr unsigned int ncornerBy2_ = 4;
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 FlatTrd& cell);
#endif
|