| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // Implementation of inline methods of G4CutTu
27 // -------------------------------------------
28
29 inline
30 G4double G4CutTubs::GetInnerRadius () const
31 {
32 return fRMin;
33 }
34
35 inline
36 G4double G4CutTubs::GetOuterRadius () const
37 {
38 return fRMax;
39 }
40
41 inline
42 G4double G4CutTubs::GetZHalfLength () const
43 {
44 return fDz;
45 }
46
47 inline
48 G4double G4CutTubs::GetStartPhiAngle () const
49 {
50 return fSPhi;
51 }
52
53 inline
54 G4double G4CutTubs::GetDeltaPhiAngle () const
55 {
56 return fDPhi;
57 }
58
59 inline
60 G4double G4CutTubs::GetSinStartPhi () const
61 {
62 return sinSPhi;
63 }
64
65 inline
66 G4double G4CutTubs::GetCosStartPhi () const
67 {
68 return cosSPhi;
69 }
70
71 inline
72 G4double G4CutTubs::GetSinEndPhi () const
73 {
74 return sinEPhi;
75 }
76
77 inline
78 G4double G4CutTubs::GetCosEndPhi () const
79 {
80 return cosEPhi;
81 }
82
83 inline
84 G4ThreeVector G4CutTubs::GetLowNorm () const
85 {
86 return fLowNorm;
87 }
88
89 inline
90 G4ThreeVector G4CutTubs::GetHighNorm () const
91 {
92 return fHighNorm;
93 }
94
95 inline
96 void G4CutTubs::Initialize()
97 {
98 fZMin = 0.;
99 fZMax = 0.;
100 fCubicVolume = 0.;
101 fSurfaceArea = 0.;
102 fRebuildPolyhedron = true;
103 }
104
105 inline
106 void G4CutTubs::InitializeTrigonometry()
107 {
108 G4double hDPhi = 0.5*fDPhi;
109 G4double cPhi = fSPhi + hDPhi;
110 G4double ePhi = fSPhi + fDPhi;
111
112 sinCPhi = std::sin(cPhi);
113 cosCPhi = std::cos(cPhi);
114 cosHDPhi = std::cos(hDPhi);
115 cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolera
116 cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolera
117 sinSPhi = std::sin(fSPhi);
118 cosSPhi = std::cos(fSPhi);
119 sinEPhi = std::sin(ePhi);
120 cosEPhi = std::cos(ePhi);
121 }
122
123 inline void G4CutTubs::CheckSPhiAngle(G4double
124 {
125 // Ensure fSphi in 0-2PI or -2PI-0 range if
126
127 if ( sPhi < 0 )
128 {
129 fSPhi = CLHEP::twopi - std::fmod(std::fabs
130 }
131 else
132 {
133 fSPhi = std::fmod(sPhi,CLHEP::twopi) ;
134 }
135 if ( fSPhi+fDPhi > CLHEP::twopi )
136 {
137 fSPhi -= CLHEP::twopi ;
138 }
139 }
140
141 inline void G4CutTubs::CheckDPhiAngle(G4double
142 {
143 fPhiFullCutTube = true;
144 if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5
145 {
146 fDPhi=CLHEP::twopi;
147 fSPhi=0;
148 }
149 else
150 {
151 fPhiFullCutTube = false;
152 if ( dPhi > 0 )
153 {
154 fDPhi = dPhi;
155 }
156 else
157 {
158 std::ostringstream message;
159 message << "Invalid dphi." << G4endl
160 << "Negative or zero delta-Phi (
161 << GetName();
162 G4Exception("G4CutTubs::CheckDPhiAngle()
163 FatalException, message);
164 }
165 }
166 }
167
168 inline void G4CutTubs::CheckPhiAngles(G4double
169 {
170 CheckDPhiAngle(dPhi);
171 if ( (fDPhi<CLHEP::twopi) && ((sPhi) != 0.0)
172 InitializeTrigonometry();
173 }
174
175 inline
176 void G4CutTubs::SetInnerRadius (G4double newRM
177 {
178 if ( newRMin < 0 ) // Check radii
179 {
180 std::ostringstream message;
181 message << "Invalid radii." << G4endl
182 << "Invalid values for radii in so
183 << " newRMin = " << newRMin
184 << ", fRMax = " << fRMax << G4endl
185 << " Negative inner radius!
186 G4Exception("G4CutTubs::SetInnerRadius()",
187 FatalException, message);
188 }
189 fRMin= newRMin;
190 Initialize();
191 }
192
193 inline
194 void G4CutTubs::SetOuterRadius (G4double newRM
195 {
196 if ( newRMax <= 0 ) // Check radii
197 {
198 std::ostringstream message;
199 message << "Invalid radii." << G4endl
200 << "Invalid values for radii in so
201 << " fRMin = " << fRMin
202 << ", newRMax = " << newRMax << G4
203 << " Invalid outer radius!"
204 G4Exception("G4CutTubs::SetOuterRadius()",
205 FatalException, message);
206 }
207 fRMax= newRMax;
208 Initialize();
209 }
210
211 inline
212 void G4CutTubs::SetZHalfLength (G4double newDz
213 {
214 if (newDz<=0) // Check z-len
215 {
216 std::ostringstream message;
217 message << "Invalid Z half-length." << G4e
218 << "Negative Z half-length (" << n
219 << GetName();
220 G4Exception("G4CutTubs::SetZHalfLength()",
221 FatalException, message);
222 }
223 fDz= newDz;
224 Initialize();
225 }
226
227 inline
228 void G4CutTubs::SetStartPhiAngle (G4double new
229 {
230 // Flag 'compute' can be used to explicitely
231 // trigonometry in case SetDeltaPhiAngle() i
232
233 CheckSPhiAngle(newSPhi);
234 fPhiFullCutTube = false;
235 if (compute) { InitializeTrigonometry(); }
236 Initialize();
237 }
238
239 inline
240 void G4CutTubs::SetDeltaPhiAngle (G4double new
241 {
242 CheckPhiAngles(fSPhi, newDPhi);
243 Initialize();
244 }
245