| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
>> 23 //
>> 24 // $Id: G4Tubs.icc,v 1.4 2004/10/10 10:44:19 johna Exp $
>> 25 // GEANT4 tag $Name: geant4-07-01 $
>> 26 //
>> 27 // --------------------------------------------------------------------
>> 28 // GEANT 4 inline definitions file
>> 29 //
>> 30 // G4Tubs.icc
>> 31 //
26 // Implementation of inline methods of G4Tubs 32 // Implementation of inline methods of G4Tubs
27 // ------------------------------------------- 33 // --------------------------------------------------------------------
28 34
29 inline 35 inline
30 G4double G4Tubs::GetInnerRadius () const 36 G4double G4Tubs::GetInnerRadius () const
31 { 37 {
32 return fRMin; 38 return fRMin;
33 } 39 }
34 40
35 inline 41 inline
36 G4double G4Tubs::GetOuterRadius () const 42 G4double G4Tubs::GetOuterRadius () const
37 { 43 {
38 return fRMax; 44 return fRMax;
39 } 45 }
40 46
41 inline 47 inline
42 G4double G4Tubs::GetZHalfLength () const 48 G4double G4Tubs::GetZHalfLength () const
43 { 49 {
44 return fDz; 50 return fDz;
45 } 51 }
46 52
47 inline 53 inline
48 G4double G4Tubs::GetStartPhiAngle () const 54 G4double G4Tubs::GetStartPhiAngle () const
49 { 55 {
50 return fSPhi; 56 return fSPhi;
51 } 57 }
52 58
53 inline 59 inline
54 G4double G4Tubs::GetDeltaPhiAngle () const 60 G4double G4Tubs::GetDeltaPhiAngle () const
55 { 61 {
56 return fDPhi; 62 return fDPhi;
57 } 63 }
58 64
59 inline 65 inline
60 G4double G4Tubs::GetSinStartPhi () const << 66 void G4Tubs::SetInnerRadius (G4double newRMin)
61 { <<
62 return sinSPhi; <<
63 } <<
64 <<
65 inline <<
66 G4double G4Tubs::GetCosStartPhi () const <<
67 { 67 {
68 return cosSPhi; << 68 fRMin= newRMin;
>> 69 fCubicVolume= 0.;
>> 70 fpPolyhedron = 0;
69 } 71 }
70 72
71 inline 73 inline
72 G4double G4Tubs::GetSinEndPhi () const << 74 void G4Tubs::SetOuterRadius (G4double newRMax)
73 { 75 {
74 return sinEPhi; << 76 fRMax= newRMax;
>> 77 fCubicVolume= 0.;
>> 78 fpPolyhedron = 0;
75 } 79 }
76 80
77 inline 81 inline
78 G4double G4Tubs::GetCosEndPhi () const << 82 void G4Tubs::SetZHalfLength (G4double newDz)
79 { 83 {
80 return cosEPhi; << 84 fDz= newDz;
>> 85 fCubicVolume= 0.;
>> 86 fpPolyhedron = 0;
81 } 87 }
82 88
83 inline 89 inline
84 void G4Tubs::Initialize() << 90 void G4Tubs::SetStartPhiAngle (G4double newSPhi)
85 { 91 {
86 fCubicVolume = 0.; << 92 fSPhi= newSPhi;
87 fSurfaceArea = 0.; << 93 fCubicVolume= 0.;
88 fInvRmax= 1.0 / fRMax; << 94 fpPolyhedron = 0;
89 fInvRmin= fRMin > 0. ? 1.0 / fRMin : 0.0; <<
90 fRebuildPolyhedron = true; <<
91 } 95 }
92 96
93 inline 97 inline
94 void G4Tubs::InitializeTrigonometry() << 98 void G4Tubs::SetDeltaPhiAngle (G4double newDPhi)
95 { <<
96 G4double hDPhi = 0.5*fDPhi; <<
97 G4double cPhi = fSPhi + hDPhi; <<
98 G4double ePhi = fSPhi + fDPhi; <<
99 <<
100 sinCPhi = std::sin(cPhi); <<
101 cosCPhi = std::cos(cPhi); <<
102 cosHDPhi = std::cos(hDPhi); <<
103 cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolera <<
104 cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolera <<
105 sinSPhi = std::sin(fSPhi); <<
106 cosSPhi = std::cos(fSPhi); <<
107 sinEPhi = std::sin(ePhi); <<
108 cosEPhi = std::cos(ePhi); <<
109 } <<
110 <<
111 inline void G4Tubs::CheckSPhiAngle(G4double sP <<
112 { <<
113 // Ensure fSphi in 0-2PI or -2PI-0 range if <<
114 <<
115 if ( sPhi < 0 ) <<
116 { <<
117 fSPhi = CLHEP::twopi - std::fmod(std::fabs <<
118 } <<
119 else <<
120 { <<
121 fSPhi = std::fmod(sPhi,CLHEP::twopi) ; <<
122 } <<
123 if ( fSPhi+fDPhi > CLHEP::twopi ) <<
124 { <<
125 fSPhi -= CLHEP::twopi ; <<
126 } <<
127 } <<
128 <<
129 inline void G4Tubs::CheckDPhiAngle(G4double dP <<
130 { 99 {
131 fPhiFullTube = true; << 100 fDPhi= newDPhi;
132 if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5 << 101 fCubicVolume= 0.;
133 { << 102 fpPolyhedron = 0;
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 } 103 }
155 104
156 inline void G4Tubs::CheckPhiAngles(G4double sP << 105 // Older names for access functions
157 { <<
158 CheckDPhiAngle(dPhi); <<
159 if ( (fDPhi<CLHEP::twopi) && ((sPhi) != 0.0) <<
160 InitializeTrigonometry(); <<
161 } <<
162 106
163 inline 107 inline
164 void G4Tubs::SetInnerRadius (G4double newRMin) << 108 G4double G4Tubs::GetRMin () const
165 { 109 {
166 if ( newRMin < 0 ) // Check radii << 110 return GetInnerRadius();
167 { <<
168 std::ostringstream message; <<
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 } <<
177 fRMin= newRMin; <<
178 Initialize(); <<
179 } 111 }
180 112
181 inline 113 inline
182 void G4Tubs::SetOuterRadius (G4double newRMax) << 114 G4double G4Tubs::GetRMax () const
183 { 115 {
184 if ( newRMax <= 0 ) // Check radii << 116 return GetOuterRadius();
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 } 117 }
198 118
199 inline 119 inline
200 void G4Tubs::SetZHalfLength (G4double newDz) << 120 G4double G4Tubs::GetDz () const
201 { 121 {
202 if (newDz<=0) // Check z-len << 122 return GetZHalfLength() ;
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 } 123 }
214 124
215 inline 125 inline
216 void G4Tubs::SetStartPhiAngle (G4double newSPh << 126 G4double G4Tubs::GetSPhi () const
217 { 127 {
218 // Flag 'compute' can be used to explicitely << 128 return GetStartPhiAngle();
219 // trigonometry in case SetDeltaPhiAngle() i <<
220 <<
221 CheckSPhiAngle(newSPhi); <<
222 fPhiFullTube = false; <<
223 if (compute) { InitializeTrigonometry(); } <<
224 Initialize(); <<
225 } 129 }
226 130
227 inline 131 inline
228 void G4Tubs::SetDeltaPhiAngle (G4double newDPh << 132 G4double G4Tubs::GetDPhi () const
229 { 133 {
230 CheckPhiAngles(fSPhi, newDPhi); << 134 return GetDeltaPhiAngle();
231 Initialize(); <<
232 } 135 }
233 136
>> 137 ///////////////////////////////////////////////////
>> 138
234 inline 139 inline
235 G4double G4Tubs::GetCubicVolume() 140 G4double G4Tubs::GetCubicVolume()
236 { 141 {
237 if(fCubicVolume != 0.) {;} << 142 if(fCubicVolume != 0.) ;
238 else { fCubicVolume = fDPhi*fDz*(fRMax*fRM << 143 else fCubicVolume = fDPhi*fDz*(fRMax*fRMax-fRMin*fRMin);
239 return fCubicVolume; 144 return fCubicVolume;
240 } <<
241 <<
242 inline <<
243 G4double G4Tubs::GetSurfaceArea() <<
244 { <<
245 if(fSurfaceArea != 0.) {;} <<
246 else <<
247 { <<
248 fSurfaceArea = fDPhi*(fRMin+fRMax)*(2*fDz+ <<
249 if (!fPhiFullTube) <<
250 { <<
251 fSurfaceArea = fSurfaceArea + 4*fDz*(fRM <<
252 } <<
253 } <<
254 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 } 145 }
267 146