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