| 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 //
>> 28 // --------------------------------------------------------------------
>> 29 // GEANT 4 inline definitions file
>> 30 //
>> 31 // G4Tubs.icc
>> 32 //
26 // Implementation of inline methods of G4Tubs 33 // Implementation of inline methods of G4Tubs
27 // ------------------------------------------- 34 // --------------------------------------------------------------------
28 35
29 inline 36 inline
30 G4double G4Tubs::GetInnerRadius () const 37 G4double G4Tubs::GetInnerRadius () const
31 { 38 {
32 return fRMin; 39 return fRMin;
33 } 40 }
34 41
35 inline 42 inline
36 G4double G4Tubs::GetOuterRadius () const 43 G4double G4Tubs::GetOuterRadius () const
37 { 44 {
38 return fRMax; 45 return fRMax;
39 } 46 }
40 47
41 inline 48 inline
42 G4double G4Tubs::GetZHalfLength () const 49 G4double G4Tubs::GetZHalfLength () const
43 { 50 {
44 return fDz; 51 return fDz;
45 } 52 }
46 53
47 inline 54 inline
48 G4double G4Tubs::GetStartPhiAngle () const 55 G4double G4Tubs::GetStartPhiAngle () const
49 { 56 {
50 return fSPhi; 57 return fSPhi;
51 } 58 }
52 59
53 inline 60 inline
54 G4double G4Tubs::GetDeltaPhiAngle () const 61 G4double G4Tubs::GetDeltaPhiAngle () const
55 { 62 {
56 return fDPhi; 63 return fDPhi;
57 } 64 }
58 65
59 inline 66 inline
60 G4double G4Tubs::GetSinStartPhi () const 67 G4double G4Tubs::GetSinStartPhi () const
61 { 68 {
62 return sinSPhi; 69 return sinSPhi;
63 } 70 }
64 71
65 inline 72 inline
66 G4double G4Tubs::GetCosStartPhi () const 73 G4double G4Tubs::GetCosStartPhi () const
67 { 74 {
68 return cosSPhi; 75 return cosSPhi;
69 } 76 }
70 77
71 inline 78 inline
72 G4double G4Tubs::GetSinEndPhi () const 79 G4double G4Tubs::GetSinEndPhi () const
73 { 80 {
74 return sinEPhi; 81 return sinEPhi;
75 } 82 }
76 83
77 inline 84 inline
78 G4double G4Tubs::GetCosEndPhi () const 85 G4double G4Tubs::GetCosEndPhi () const
79 { 86 {
80 return cosEPhi; 87 return cosEPhi;
81 } 88 }
82 89
83 inline << 90 inline
84 void G4Tubs::Initialize() 91 void G4Tubs::Initialize()
85 { 92 {
86 fCubicVolume = 0.; 93 fCubicVolume = 0.;
87 fSurfaceArea = 0.; 94 fSurfaceArea = 0.;
88 fInvRmax= 1.0 / fRMax; <<
89 fInvRmin= fRMin > 0. ? 1.0 / fRMin : 0.0; <<
90 fRebuildPolyhedron = true; 95 fRebuildPolyhedron = true;
91 } 96 }
92 97
93 inline << 98 inline
94 void G4Tubs::InitializeTrigonometry() 99 void G4Tubs::InitializeTrigonometry()
95 { 100 {
96 G4double hDPhi = 0.5*fDPhi; 101 G4double hDPhi = 0.5*fDPhi; // half delta phi
97 G4double cPhi = fSPhi + hDPhi; << 102 G4double cPhi = fSPhi + hDPhi;
98 G4double ePhi = fSPhi + fDPhi; 103 G4double ePhi = fSPhi + fDPhi;
99 104
100 sinCPhi = std::sin(cPhi); 105 sinCPhi = std::sin(cPhi);
101 cosCPhi = std::cos(cPhi); 106 cosCPhi = std::cos(cPhi);
102 cosHDPhi = std::cos(hDPhi); <<
103 cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolera 107 cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolerance); // inner/outer tol half dphi
104 cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolera 108 cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolerance);
105 sinSPhi = std::sin(fSPhi); 109 sinSPhi = std::sin(fSPhi);
106 cosSPhi = std::cos(fSPhi); 110 cosSPhi = std::cos(fSPhi);
107 sinEPhi = std::sin(ePhi); 111 sinEPhi = std::sin(ePhi);
108 cosEPhi = std::cos(ePhi); 112 cosEPhi = std::cos(ePhi);
109 } 113 }
110 114
111 inline void G4Tubs::CheckSPhiAngle(G4double sP 115 inline void G4Tubs::CheckSPhiAngle(G4double sPhi)
112 { 116 {
113 // Ensure fSphi in 0-2PI or -2PI-0 range if 117 // Ensure fSphi in 0-2PI or -2PI-0 range if shape crosses 0
114 118
115 if ( sPhi < 0 ) 119 if ( sPhi < 0 )
116 { 120 {
117 fSPhi = CLHEP::twopi - std::fmod(std::fabs 121 fSPhi = CLHEP::twopi - std::fmod(std::fabs(sPhi),CLHEP::twopi);
118 } 122 }
119 else 123 else
120 { 124 {
121 fSPhi = std::fmod(sPhi,CLHEP::twopi) ; 125 fSPhi = std::fmod(sPhi,CLHEP::twopi) ;
122 } 126 }
123 if ( fSPhi+fDPhi > CLHEP::twopi ) 127 if ( fSPhi+fDPhi > CLHEP::twopi )
124 { 128 {
125 fSPhi -= CLHEP::twopi ; 129 fSPhi -= CLHEP::twopi ;
126 } 130 }
127 } 131 }
128 132
129 inline void G4Tubs::CheckDPhiAngle(G4double dP 133 inline void G4Tubs::CheckDPhiAngle(G4double dPhi)
130 { 134 {
131 fPhiFullTube = true; 135 fPhiFullTube = true;
132 if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5 136 if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5 )
133 { 137 {
134 fDPhi=CLHEP::twopi; 138 fDPhi=CLHEP::twopi;
135 fSPhi=0; 139 fSPhi=0;
136 } 140 }
137 else 141 else
138 { 142 {
139 fPhiFullTube = false; 143 fPhiFullTube = false;
140 if ( dPhi > 0 ) 144 if ( dPhi > 0 )
141 { 145 {
142 fDPhi = dPhi; 146 fDPhi = dPhi;
143 } 147 }
144 else 148 else
145 { 149 {
146 std::ostringstream message; 150 std::ostringstream message;
147 message << "Invalid dphi." << G4endl 151 message << "Invalid dphi." << G4endl
148 << "Negative or zero delta-Phi ( 152 << "Negative or zero delta-Phi (" << dPhi << "), for solid: "
149 << GetName(); 153 << GetName();
150 G4Exception("G4Tubs::CheckDPhiAngle()", 154 G4Exception("G4Tubs::CheckDPhiAngle()", "GeomSolids0002",
151 FatalException, message); 155 FatalException, message);
152 } 156 }
153 } 157 }
154 } 158 }
155 159
156 inline void G4Tubs::CheckPhiAngles(G4double sP 160 inline void G4Tubs::CheckPhiAngles(G4double sPhi, G4double dPhi)
157 { 161 {
158 CheckDPhiAngle(dPhi); 162 CheckDPhiAngle(dPhi);
159 if ( (fDPhi<CLHEP::twopi) && ((sPhi) != 0.0) << 163 if ( (fDPhi<CLHEP::twopi) && (sPhi) ) { CheckSPhiAngle(sPhi); }
160 InitializeTrigonometry(); 164 InitializeTrigonometry();
161 } 165 }
162 166
163 inline 167 inline
164 void G4Tubs::SetInnerRadius (G4double newRMin) 168 void G4Tubs::SetInnerRadius (G4double newRMin)
165 { 169 {
166 if ( newRMin < 0 ) // Check radii 170 if ( newRMin < 0 ) // Check radii
167 { 171 {
168 std::ostringstream message; 172 std::ostringstream message;
169 message << "Invalid radii." << G4endl 173 message << "Invalid radii." << G4endl
170 << "Invalid values for radii in so 174 << "Invalid values for radii in solid " << GetName() << G4endl
171 << " newRMin = " << newRMin 175 << " newRMin = " << newRMin
172 << ", fRMax = " << fRMax << G4endl 176 << ", fRMax = " << fRMax << G4endl
173 << " Negative inner radius! 177 << " Negative inner radius!";
174 G4Exception("G4Tubs::SetInnerRadius()", "G 178 G4Exception("G4Tubs::SetInnerRadius()", "GeomSolids0002",
175 FatalException, message); 179 FatalException, message);
176 } 180 }
177 fRMin= newRMin; 181 fRMin= newRMin;
178 Initialize(); 182 Initialize();
179 } 183 }
180 184
181 inline 185 inline
182 void G4Tubs::SetOuterRadius (G4double newRMax) 186 void G4Tubs::SetOuterRadius (G4double newRMax)
183 { 187 {
184 if ( newRMax <= 0 ) // Check radii 188 if ( newRMax <= 0 ) // Check radii
185 { 189 {
186 std::ostringstream message; 190 std::ostringstream message;
187 message << "Invalid radii." << G4endl 191 message << "Invalid radii." << G4endl
188 << "Invalid values for radii in so 192 << "Invalid values for radii in solid " << GetName() << G4endl
189 << " fRMin = " << fRMin 193 << " fRMin = " << fRMin
190 << ", newRMax = " << newRMax << G4 194 << ", newRMax = " << newRMax << G4endl
191 << " Invalid outer radius!" 195 << " Invalid outer radius!";
192 G4Exception("G4Tubs::SetOuterRadius()", "G 196 G4Exception("G4Tubs::SetOuterRadius()", "GeomSolids0002",
193 FatalException, message); 197 FatalException, message);
194 } 198 }
195 fRMax= newRMax; 199 fRMax= newRMax;
196 Initialize(); 200 Initialize();
197 } 201 }
198 202
199 inline 203 inline
200 void G4Tubs::SetZHalfLength (G4double newDz) 204 void G4Tubs::SetZHalfLength (G4double newDz)
201 { 205 {
202 if (newDz<=0) // Check z-len 206 if (newDz<=0) // Check z-len
203 { 207 {
204 std::ostringstream message; 208 std::ostringstream message;
205 message << "Invalid Z half-length." << G4e 209 message << "Invalid Z half-length." << G4endl
206 << "Negative Z half-length (" << n 210 << "Negative Z half-length (" << newDz << "), for solid: "
207 << GetName(); 211 << GetName();
208 G4Exception("G4Tubs::SetZHalfLength()", "G 212 G4Exception("G4Tubs::SetZHalfLength()", "GeomSolids0002",
209 FatalException, message); 213 FatalException, message);
210 } 214 }
211 fDz= newDz; 215 fDz= newDz;
212 Initialize(); 216 Initialize();
213 } 217 }
214 218
215 inline 219 inline
216 void G4Tubs::SetStartPhiAngle (G4double newSPh 220 void G4Tubs::SetStartPhiAngle (G4double newSPhi, G4bool compute)
217 { 221 {
218 // Flag 'compute' can be used to explicitely 222 // Flag 'compute' can be used to explicitely avoid recomputation of
219 // trigonometry in case SetDeltaPhiAngle() i 223 // trigonometry in case SetDeltaPhiAngle() is invoked afterwards
220 224
221 CheckSPhiAngle(newSPhi); 225 CheckSPhiAngle(newSPhi);
222 fPhiFullTube = false; 226 fPhiFullTube = false;
223 if (compute) { InitializeTrigonometry(); } 227 if (compute) { InitializeTrigonometry(); }
224 Initialize(); 228 Initialize();
225 } 229 }
226 230
227 inline 231 inline
228 void G4Tubs::SetDeltaPhiAngle (G4double newDPh 232 void G4Tubs::SetDeltaPhiAngle (G4double newDPhi)
229 { 233 {
230 CheckPhiAngles(fSPhi, newDPhi); 234 CheckPhiAngles(fSPhi, newDPhi);
231 Initialize(); 235 Initialize();
232 } 236 }
233 237
>> 238 // Older names for access functions
>> 239
>> 240 inline
>> 241 G4double G4Tubs::GetRMin () const
>> 242 {
>> 243 return GetInnerRadius();
>> 244 }
>> 245
>> 246 inline
>> 247 G4double G4Tubs::GetRMax () const
>> 248 {
>> 249 return GetOuterRadius();
>> 250 }
>> 251
>> 252 inline
>> 253 G4double G4Tubs::GetDz () const
>> 254 {
>> 255 return GetZHalfLength() ;
>> 256 }
>> 257
>> 258 inline
>> 259 G4double G4Tubs::GetSPhi () const
>> 260 {
>> 261 return GetStartPhiAngle();
>> 262 }
>> 263
>> 264 inline
>> 265 G4double G4Tubs::GetDPhi () const
>> 266 {
>> 267 return GetDeltaPhiAngle();
>> 268 }
>> 269
234 inline 270 inline
235 G4double G4Tubs::GetCubicVolume() 271 G4double G4Tubs::GetCubicVolume()
236 { 272 {
237 if(fCubicVolume != 0.) {;} 273 if(fCubicVolume != 0.) {;}
238 else { fCubicVolume = fDPhi*fDz*(fRMax*fRM 274 else { fCubicVolume = fDPhi*fDz*(fRMax*fRMax-fRMin*fRMin); }
239 return fCubicVolume; 275 return fCubicVolume;
240 } 276 }
241 277
242 inline 278 inline
243 G4double G4Tubs::GetSurfaceArea() 279 G4double G4Tubs::GetSurfaceArea()
244 { 280 {
245 if(fSurfaceArea != 0.) {;} 281 if(fSurfaceArea != 0.) {;}
246 else 282 else
247 { 283 {
248 fSurfaceArea = fDPhi*(fRMin+fRMax)*(2*fDz+ 284 fSurfaceArea = fDPhi*(fRMin+fRMax)*(2*fDz+fRMax-fRMin);
249 if (!fPhiFullTube) 285 if (!fPhiFullTube)
250 { 286 {
251 fSurfaceArea = fSurfaceArea + 4*fDz*(fRM 287 fSurfaceArea = fSurfaceArea + 4*fDz*(fRMax-fRMin);
252 } 288 }
253 } 289 }
254 return fSurfaceArea; 290 return fSurfaceArea;
255 } <<
256 <<
257 inline <<
258 G4double G4Tubs::FastInverseRxy( const G4Three <<
259 G4doubl <<
260 G4doubl <<
261 { <<
262 G4double rho2 = sqr( pos.x() ) + sqr( pos.y( <<
263 G4bool onSurface = std::fabs( rho2 * invRad <<
264 G4double invRho = onSurface ? invRad : 1.0 / <<
265 return invRho; <<
266 } 291 }
267 292