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