| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
>> 26 //
>> 27 // $Id: G4Para.hh 79491 2014-03-05 15:24:29Z gcosmo $
>> 28 //
>> 29 //
>> 30 // --------------------------------------------------------------------
>> 31 // GEANT 4 class header file
>> 32 //
26 // G4Para 33 // G4Para
27 // 34 //
28 // Class description: 35 // Class description:
29 // 36 //
30 // A G4Parallepiped, essentially a box with 37 // A G4Parallepiped, essentially a box with half lengths dx,dy,dz
31 // `skewed' so that there are angles theta & 38 // `skewed' so that there are angles theta & phi of the polar line
32 // joining the faces at +-dz in z, and alpha 39 // joining the faces at +-dz in z, and alpha formed by the y axis
33 // and the plane joinng the centre of the fa 40 // and the plane joinng the centre of the faces G4Parallel to the
34 // z-x plane at -dy and +dy. 41 // z-x plane at -dy and +dy.
35 // 42 //
36 // A G4Para is defined by: 43 // A G4Para is defined by:
37 // dx,dy,dz - Half-length in x,y,z 44 // dx,dy,dz - Half-length in x,y,z
38 // alpha - Angle formed by the y axis and 45 // alpha - Angle formed by the y axis and by the plane joining
39 // the centre of the faces G4Para 46 // the centre of the faces G4Parallel to the z-x plane
40 // at -dy and +dy 47 // at -dy and +dy
41 // theta - Polar angle of the line joinin 48 // theta - Polar angle of the line joining the centres of the
42 // faces at -dz and +dz in z 49 // faces at -dz and +dz in z
43 // phi - Azimuthal angle of the line jo 50 // phi - Azimuthal angle of the line joining the centres of the
44 // faces at -dz and +dz in z 51 // faces at -dz and +dz in z
45 // Member data: 52 // Member data:
46 // 53 //
47 // Note that the angles parameters are not s 54 // Note that the angles parameters are not stored - precomputed trig is
48 // stored instead. 55 // stored instead.
49 // 56 //
50 // fDx Half-length in x 57 // fDx Half-length in x
51 // fDy Half-length in y 58 // fDy Half-length in y
52 // fDz Half-length in z 59 // fDz Half-length in z
53 // 60 //
54 // fTalpha Tan of alpha 61 // fTalpha Tan of alpha
55 // fTthetaCphi Tan theta * Cos phi 62 // fTthetaCphi Tan theta * Cos phi
56 // fTthetaSphi Tan theta * Sin phi 63 // fTthetaSphi Tan theta * Sin phi
57 64
>> 65 // History:
58 // 21.3.94 P.Kent Old C++ code converted to to 66 // 21.3.94 P.Kent Old C++ code converted to tolerant geometry
59 // 31.10.96 V.Grichine Modifications according 67 // 31.10.96 V.Grichine Modifications according G4Box/Tubs before to commit
>> 68 // 18.11.99 V.Grichine , kUndefined was added to ESide
60 // ------------------------------------------- 69 // --------------------------------------------------------------------
61 #ifndef G4PARA_HH <<
62 #define G4PARA_HH <<
63 70
64 #include "G4GeomTypes.hh" << 71 #ifndef G4Para_HH
65 << 72 #define G4Para_HH
66 #if defined(G4GEOM_USE_USOLIDS) <<
67 #define G4GEOM_USE_UPARA 1 <<
68 #endif <<
69 <<
70 #if defined(G4GEOM_USE_UPARA) <<
71 #define G4UPara G4Para <<
72 #include "G4UPara.hh" <<
73 #else <<
74 73
75 #include "G4CSGSolid.hh" 74 #include "G4CSGSolid.hh"
76 #include "G4Polyhedron.hh" 75 #include "G4Polyhedron.hh"
77 76
78 class G4Para : public G4CSGSolid 77 class G4Para : public G4CSGSolid
79 { 78 {
80 public: // with description 79 public: // with description
81 80
82 G4Para(const G4String& pName, 81 G4Para(const G4String& pName,
83 G4double pDx, G4double pDy, G 82 G4double pDx, G4double pDy, G4double pDz,
84 G4double pAlpha, G4double pTh 83 G4double pAlpha, G4double pTheta, G4double pPhi);
85 << 84
86 G4Para(const G4String& pName, 85 G4Para(const G4String& pName,
87 const G4ThreeVector pt[8]); 86 const G4ThreeVector pt[8]);
88 87
89 ~G4Para() override; << 88 virtual ~G4Para();
90 << 89
91 // Accessors << 90 // Accessors
92 91
93 inline G4double GetZHalfLength() const; 92 inline G4double GetZHalfLength() const;
94 inline G4ThreeVector GetSymAxis() const; 93 inline G4ThreeVector GetSymAxis() const;
95 inline G4double GetYHalfLength() const; 94 inline G4double GetYHalfLength() const;
96 inline G4double GetXHalfLength() const; 95 inline G4double GetXHalfLength() const;
97 inline G4double GetTanAlpha() const; 96 inline G4double GetTanAlpha() const;
98 << 97
99 inline G4double GetAlpha() const; << 98 // Modifiers
100 inline G4double GetTheta() const; <<
101 inline G4double GetPhi() const; <<
102 // Obtain (re)computed values of original <<
103 <<
104 // Modifiers <<
105 99
106 inline void SetXHalfLength(G4double val); 100 inline void SetXHalfLength(G4double val);
107 inline void SetYHalfLength(G4double val); 101 inline void SetYHalfLength(G4double val);
108 inline void SetZHalfLength(G4double val); 102 inline void SetZHalfLength(G4double val);
109 inline void SetAlpha(G4double alpha); 103 inline void SetAlpha(G4double alpha);
110 inline void SetTanAlpha(G4double val); 104 inline void SetTanAlpha(G4double val);
111 inline void SetThetaAndPhi(G4double pTheta << 105 inline void SetThetaAndPhi(double pTheta, double pPhi);
112 <<
113 void SetAllParameters(G4double pDx, G4doub <<
114 G4double pAlpha, G4d <<
115 106
116 // Methods of solid << 107 void SetAllParameters(G4double pDx, G4double pDy, G4double pDz,
117 << 108 G4double pAlpha, G4double pTheta, G4double pPhi);
118 G4double GetCubicVolume() override; << 109
119 G4double GetSurfaceArea() override; << 110 // Other methods of solid
>> 111
>> 112 inline G4double GetCubicVolume();
>> 113 inline G4double GetSurfaceArea();
120 114
121 void ComputeDimensions(G4VPVParameterisati 115 void ComputeDimensions(G4VPVParameterisation* p,
122 const G4int n, 116 const G4int n,
123 const G4VPhysicalVo << 117 const G4VPhysicalVolume* pRep);
124 <<
125 void BoundingLimits(G4ThreeVector& pMin, G <<
126 118
127 G4bool CalculateExtent(const EAxis pAxis, 119 G4bool CalculateExtent(const EAxis pAxis,
128 const G4VoxelLimits 120 const G4VoxelLimits& pVoxelLimit,
129 const G4AffineTrans 121 const G4AffineTransform& pTransform,
130 G4double& pMi << 122 G4double& pMin, G4double& pMax) const;
>> 123
>> 124 EInside Inside(const G4ThreeVector& p) const;
131 125
132 EInside Inside(const G4ThreeVector& p) con << 126 G4ThreeVector SurfaceNormal( const G4ThreeVector& p) const;
133 <<
134 G4ThreeVector SurfaceNormal( const G4Three <<
135 127
136 G4double DistanceToIn(const G4ThreeVector& 128 G4double DistanceToIn(const G4ThreeVector& p,
137 const G4ThreeVector& << 129 const G4ThreeVector& v) const;
138 G4double DistanceToIn(const G4ThreeVector& << 130 G4double DistanceToIn(const G4ThreeVector& p) const;
139 << 131
140 G4double DistanceToOut(const G4ThreeVector 132 G4double DistanceToOut(const G4ThreeVector& p, const G4ThreeVector& v,
141 const G4bool calcNo << 133 const G4bool calcNorm=G4bool(false),
142 G4bool* valid << 134 G4bool *validNorm=0, G4ThreeVector *n=0) const;
143 G4ThreeVector << 135 G4double DistanceToOut(const G4ThreeVector& p) const;
144 G4double DistanceToOut(const G4ThreeVector <<
145 136
146 G4GeometryType GetEntityType() const overr << 137 G4GeometryType GetEntityType() const;
147 138
148 G4ThreeVector GetPointOnSurface() const ov << 139 G4ThreeVector GetPointOnSurface() const;
149 140
150 G4bool IsFaceted() const override; << 141 G4VSolid* Clone() const;
151 142
152 G4VSolid* Clone() const override; << 143 std::ostream& StreamInfo(std::ostream& os) const;
153 144
154 std::ostream& StreamInfo(std::ostream& os) << 145 // Visualisation functions
155 146
156 // Visualisation functions << 147 void DescribeYourselfTo (G4VGraphicsScene& scene) const;
>> 148 G4Polyhedron* CreatePolyhedron () const;
157 149
158 void DescribeYourselfTo (G4VGraph << 150 public: // without description
159 G4Polyhedron* CreatePolyhedron () const <<
160 151
161 G4Para(__void__&); 152 G4Para(__void__&);
162 // Fake default constructor for usage re 153 // Fake default constructor for usage restricted to direct object
163 // persistency for clients requiring pre 154 // persistency for clients requiring preallocation of memory for
164 // persistifiable objects << 155 // persistifiable objects.
165 156
166 G4Para(const G4Para& rhs); 157 G4Para(const G4Para& rhs);
167 G4Para& operator=(const G4Para& rhs); << 158 G4Para& operator=(const G4Para& rhs);
168 // Copy constructor and assignment opera << 159 // Copy constructor and assignment operator.
169 160
170 private: << 161 protected: // without description
171 162
172 void CheckParameters(); << 163 G4ThreeVectorList*
173 // Check parameters << 164 CreateRotatedVertices(const G4AffineTransform& pTransform) const;
174 165
175 void MakePlanes(); << 166 private:
176 // Set side planes <<
177 167
178 G4ThreeVector ApproxSurfaceNormal(const G4 << 168 G4ThreeVector ApproxSurfaceNormal( const G4ThreeVector& p) const;
179 // Algorithm for SurfaceNormal() followi 169 // Algorithm for SurfaceNormal() following the original
180 // specification for points not on the s 170 // specification for points not on the surface
181 171
182 << 172 G4ThreeVector GetPointOnPlane(G4ThreeVector p0, G4ThreeVector p1,
183 private: << 173 G4ThreeVector p2, G4ThreeVector p3,
>> 174 G4double& area) const;
>> 175 // Returns a random point on the surface of one of the faces.
>> 176
>> 177 private:
184 178
185 G4double halfCarTolerance; <<
186 G4double fDx,fDy,fDz; 179 G4double fDx,fDy,fDz;
187 G4double fTalpha,fTthetaCphi,fTthetaSphi; 180 G4double fTalpha,fTthetaCphi,fTthetaSphi;
188 struct { G4double a,b,c,d; } fPlanes[4]; <<
189 }; 181 };
190 182
191 #include "G4Para.icc" 183 #include "G4Para.icc"
192 <<
193 #endif <<
194 184
195 #endif 185 #endif
196 186