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: G4Sphere.icc,v 1.6 2006/06/29 18:44:29 gunter Exp $ >> 28 // GEANT4 tag $Name: geant4-08-01 $ >> 29 // >> 30 // -------------------------------------------------------------------- >> 31 // GEANT 4 inline definitions file >> 32 // >> 33 // G4Sphere.icc >> 34 // 26 // Implementation of inline methods of G4Spher 35 // Implementation of inline methods of G4Sphere 27 // ------------------------------------------- 36 // -------------------------------------------------------------------- 28 37 29 inline 38 inline 30 G4double G4Sphere::GetInnerRadius() const << 39 G4double G4Sphere::GetInsideRadius() const 31 { 40 { 32 return fRmin; 41 return fRmin; 33 } 42 } 34 43 35 inline 44 inline 36 G4double G4Sphere::GetOuterRadius() const 45 G4double G4Sphere::GetOuterRadius() const 37 { 46 { 38 return fRmax; 47 return fRmax; 39 } 48 } 40 49 41 inline 50 inline 42 G4double G4Sphere::GetStartPhiAngle() const 51 G4double G4Sphere::GetStartPhiAngle() const 43 { 52 { 44 return fSPhi; 53 return fSPhi; 45 } 54 } 46 55 47 inline 56 inline 48 G4double G4Sphere::GetDeltaPhiAngle() const 57 G4double G4Sphere::GetDeltaPhiAngle() const 49 { 58 { 50 return fDPhi; 59 return fDPhi; 51 } 60 } 52 61 53 inline 62 inline 54 G4double G4Sphere::GetStartThetaAngle() const 63 G4double G4Sphere::GetStartThetaAngle() const 55 { 64 { 56 return fSTheta; 65 return fSTheta; 57 } 66 } 58 67 59 G4double G4Sphere::GetDeltaThetaAngle() const 68 G4double G4Sphere::GetDeltaThetaAngle() const 60 { 69 { 61 return fDTheta; 70 return fDTheta; 62 } 71 } 63 72 64 inline 73 inline 65 G4double G4Sphere::GetSinStartPhi () const << 74 void G4Sphere::SetInsideRadius(G4double newRmin) 66 { << 67 return sinSPhi; << 68 } << 69 << 70 inline << 71 G4double G4Sphere::GetCosStartPhi () const << 72 { << 73 return cosSPhi; << 74 } << 75 << 76 inline << 77 G4double G4Sphere::GetSinEndPhi () const << 78 { << 79 return sinEPhi; << 80 } << 81 << 82 inline << 83 G4double G4Sphere::GetCosEndPhi () const << 84 { << 85 return cosEPhi; << 86 } << 87 << 88 inline << 89 G4double G4Sphere::GetSinStartTheta () const << 90 { << 91 return sinSTheta; << 92 } << 93 << 94 inline << 95 G4double G4Sphere::GetCosStartTheta () const << 96 { << 97 return cosSTheta; << 98 } << 99 << 100 inline << 101 G4double G4Sphere::GetSinEndTheta () const << 102 { 75 { 103 return sinETheta; << 76 fRmin= newRmin; 104 } << 77 fCubicVolume= 0.; 105 << 78 fpPolyhedron = 0; 106 inline << 107 G4double G4Sphere::GetCosEndTheta () const << 108 { << 109 return cosETheta; << 110 } 79 } 111 80 112 inline 81 inline 113 void G4Sphere::Initialize() << 82 void G4Sphere::SetOuterRadius(G4double newRmax) 114 { 83 { 115 fCubicVolume = 0.; << 84 fRmax= newRmax; 116 fSurfaceArea = 0.; << 85 fCubicVolume= 0.; 117 fRebuildPolyhedron = true; << 86 fpPolyhedron = 0; 118 } 87 } 119 88 120 inline 89 inline 121 void G4Sphere::InitializePhiTrigonometry() << 90 void G4Sphere::SetStartPhiAngle(G4double newSphi) 122 { 91 { 123 hDPhi = 0.5*fDPhi; // << 92 fSPhi= newSphi; 124 cPhi = fSPhi + hDPhi; << 93 fCubicVolume= 0.; 125 ePhi = fSPhi + fDPhi; << 94 fpPolyhedron = 0; 126 << 127 sinCPhi = std::sin(cPhi); << 128 cosCPhi = std::cos(cPhi); << 129 cosHDPhi = std::cos(hDPhi); << 130 cosHDPhiIT = std::cos(hDPhi - 0.5*kAngTolera << 131 cosHDPhiOT = std::cos(hDPhi + 0.5*kAngTolera << 132 sinSPhi = std::sin(fSPhi); << 133 cosSPhi = std::cos(fSPhi); << 134 sinEPhi = std::sin(ePhi); << 135 cosEPhi = std::cos(ePhi); << 136 } 95 } 137 96 138 inline 97 inline 139 void G4Sphere::InitializeThetaTrigonometry() << 98 void G4Sphere::SetDeltaPhiAngle(G4double newDphi) 140 { 99 { 141 eTheta = fSTheta + fDTheta; << 100 fDPhi= newDphi; 142 << 101 fCubicVolume= 0.; 143 sinSTheta = std::sin(fSTheta); << 102 fpPolyhedron = 0; 144 cosSTheta = std::cos(fSTheta); << 145 sinETheta = std::sin(eTheta); << 146 cosETheta = std::cos(eTheta); << 147 << 148 tanSTheta = sinSTheta/cosSTheta; << 149 tanSTheta2 = tanSTheta*tanSTheta; << 150 tanETheta = sinETheta/cosETheta; << 151 tanETheta2 = tanETheta*tanETheta; << 152 } 103 } 153 104 154 inline 105 inline 155 void G4Sphere::CheckThetaAngles(G4double sThet << 106 void G4Sphere::SetStartThetaAngle(G4double newSTheta) 156 { 107 { 157 if ( (sTheta<0) || (sTheta>CLHEP::pi) ) << 108 fSTheta=newSTheta; 158 { << 109 fCubicVolume= 0.; 159 std::ostringstream message; << 110 fpPolyhedron = 0; 160 message << "sTheta outside 0-PI range." << << 161 << "Invalid starting Theta angle f << 162 G4Exception("G4Sphere::CheckThetaAngles()" << 163 FatalException, message); << 164 } << 165 else << 166 { << 167 fSTheta=sTheta; << 168 } << 169 if ( dTheta+sTheta >= CLHEP::pi ) << 170 { << 171 fDTheta=CLHEP::pi-sTheta; << 172 } << 173 else if ( dTheta > 0 ) << 174 { << 175 fDTheta=dTheta; << 176 } << 177 else << 178 { << 179 std::ostringstream message; << 180 message << "Invalid dTheta." << G4endl << 181 << "Negative delta-Theta (" << dTh << 182 << GetName(); << 183 G4Exception("G4Sphere::CheckThetaAngles()" << 184 FatalException, message); << 185 } << 186 fFullThetaSphere = fDTheta-fSTheta >= CLHEP: << 187 fFullSphere = fFullPhiSphere && fFullThetaSp << 188 << 189 InitializeThetaTrigonometry(); << 190 } 111 } 191 112 192 inline 113 inline 193 void G4Sphere::CheckSPhiAngle(G4double sPhi) << 114 void G4Sphere::SetDeltaThetaAngle(G4double newDTheta) 194 { 115 { 195 // Ensure fSphi in 0-2PI or -2PI-0 range if << 116 fDTheta=newDTheta; 196 << 117 fCubicVolume= 0.; 197 if ( sPhi < 0 ) << 118 fpPolyhedron = 0; 198 { << 199 fSPhi = CLHEP::twopi - std::fmod(std::fabs << 200 } << 201 else << 202 { << 203 fSPhi = std::fmod(sPhi,CLHEP::twopi) ; << 204 } << 205 if ( fSPhi+fDPhi > CLHEP::twopi ) << 206 { << 207 fSPhi -= CLHEP::twopi ; << 208 } << 209 } 119 } 210 120 211 inline << 121 // Old access functions 212 void G4Sphere::CheckDPhiAngle(G4double dPhi) << 213 { << 214 fFullPhiSphere = true; << 215 if ( dPhi >= CLHEP::twopi-kAngTolerance*0.5 << 216 { << 217 fDPhi=CLHEP::twopi; << 218 } << 219 else << 220 { << 221 fFullPhiSphere = false; << 222 if ( dPhi > 0 ) << 223 { << 224 fDPhi = dPhi; << 225 } << 226 else << 227 { << 228 std::ostringstream message; << 229 message << "Invalid dphi." << G4endl << 230 << "Negative delta-Phi (" << dPh << 231 << GetName(); << 232 G4Exception("G4Sphere::CheckDPhiAngle()" << 233 FatalException, message); << 234 } << 235 } << 236 } << 237 122 238 inline 123 inline 239 void G4Sphere::CheckPhiAngles(G4double sPhi, G << 124 G4double G4Sphere::GetRmin() const 240 { 125 { 241 CheckDPhiAngle(dPhi); << 126 return GetInsideRadius(); 242 if (!fFullPhiSphere && (sPhi != 0.0)) { Chec << 243 fFullSphere = fFullPhiSphere && fFullThetaSp << 244 << 245 InitializePhiTrigonometry(); << 246 } 127 } 247 128 248 inline 129 inline 249 void G4Sphere::SetInnerRadius(G4double newRmin << 130 G4double G4Sphere::GetRmax() const 250 { 131 { 251 fRmin= newRmin; << 132 return GetOuterRadius(); 252 fRminTolerance = (fRmin) != 0.0 ? std::max( << 253 Initialize(); << 254 } 133 } 255 134 256 inline 135 inline 257 void G4Sphere::SetOuterRadius(G4double newRmax << 136 G4double G4Sphere::GetSPhi() const 258 { 137 { 259 fRmax= newRmax; << 138 return GetStartPhiAngle(); 260 fRmaxTolerance = std::max( kRadTolerance, fE << 261 Initialize(); << 262 } 139 } 263 140 264 inline 141 inline 265 void G4Sphere::SetStartPhiAngle(G4double newSP << 142 G4double G4Sphere::GetDPhi() const 266 { 143 { 267 // Flag 'compute' can be used to explicitely << 144 return GetDeltaPhiAngle(); 268 // trigonometry in case SetDeltaPhiAngle() i << 269 << 270 CheckSPhiAngle(newSPhi); << 271 fFullPhiSphere = false; << 272 if (compute) { InitializePhiTrigonometry(); << 273 Initialize(); << 274 } 145 } 275 146 276 inline 147 inline 277 void G4Sphere::SetDeltaPhiAngle(G4double newDP << 148 G4double G4Sphere::GetSTheta() const 278 { 149 { 279 CheckPhiAngles(fSPhi, newDPhi); << 150 return GetStartThetaAngle(); 280 Initialize(); << 281 } 151 } 282 152 283 inline 153 inline 284 void G4Sphere::SetStartThetaAngle(G4double new << 154 G4double G4Sphere::GetDTheta() const 285 { 155 { 286 CheckThetaAngles(newSTheta, fDTheta); << 156 return GetDeltaThetaAngle(); 287 Initialize(); << 288 } 157 } >> 158 /////////////////////////////////////////////////// 289 159 290 inline 160 inline 291 void G4Sphere::SetDeltaThetaAngle(G4double new << 161 G4double G4Sphere::GetCubicVolume() 292 { 162 { 293 CheckThetaAngles(fSTheta, newDTheta); << 163 if(fCubicVolume != 0.) ; 294 Initialize(); << 164 else fCubicVolume = fDPhi*(std::cos(fSTheta)-std::cos(fSTheta+fDTheta))* >> 165 (fRmax*fRmax*fRmax-fRmin*fRmin*fRmin)/3.; >> 166 return fCubicVolume; 295 } 167 } 296 168