| 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: G4CutTubs.icc 66356 2012-12-18 09:02:32Z gcosmo $
>> 28 //
>> 29 // --------------------------------------------------------------------
>> 30 // GEANT 4 inline definitions file
>> 31 //
>> 32 // G4CutTubs.icc
>> 33 //
26 // Implementation of inline methods of G4CutTu 34 // Implementation of inline methods of G4CutTubs
27 // ------------------------------------------- 35 // --------------------------------------------------------------------
28 36
29 inline << 37 inline G4double G4CutTubs::GetCubicVolume()
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 { 38 {
74 return sinEPhi; << 39 return G4VSolid::GetCubicVolume();
75 } 40 }
76 41
77 inline << 42 inline G4double G4CutTubs::GetSurfaceArea()
78 G4double G4CutTubs::GetCosEndPhi () const <<
79 { 43 {
80 return cosEPhi; << 44 return G4VSolid::GetSurfaceArea();
81 } 45 }
82 46
83 inline 47 inline
84 G4ThreeVector G4CutTubs::GetLowNorm () const 48 G4ThreeVector G4CutTubs::GetLowNorm () const
85 { 49 {
86 return fLowNorm; 50 return fLowNorm;
87 } 51 }
88 <<
89 inline 52 inline
90 G4ThreeVector G4CutTubs::GetHighNorm () const 53 G4ThreeVector G4CutTubs::GetHighNorm () const
91 { 54 {
92 return fHighNorm; 55 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 } 56 }
245 57