Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // Implementation of inline methods of G4Spher 27 // ------------------------------------------- 28 29 inline 30 G4double G4Sphere::GetInnerRadius() const 31 { 32 return fRmin; 33 } 34 35 inline 36 G4double G4Sphere::GetOuterRadius() const 37 { 38 return fRmax; 39 } 40 41 inline 42 G4double G4Sphere::GetStartPhiAngle() const 43 { 44 return fSPhi; 45 } 46 47 inline 48 G4double G4Sphere::GetDeltaPhiAngle() const 49 { 50 return fDPhi; 51 } 52 53 inline 54 G4double G4Sphere::GetStartThetaAngle() const 55 { 56 return fSTheta; 57 } 58 59 G4double G4Sphere::GetDeltaThetaAngle() const 60 { 61 return fDTheta; 62 } 63 64 inline 65 G4double G4Sphere::GetSinStartPhi () const 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 { 103 return sinETheta; 104 } 105 106 inline 107 G4double G4Sphere::GetCosEndTheta () const 108 { 109 return cosETheta; 110 } 111 112 inline 113 void G4Sphere::Initialize() 114 { 115 fCubicVolume = 0.; 116 fSurfaceArea = 0.; 117 fRebuildPolyhedron = true; 118 } 119 120 inline 121 void G4Sphere::InitializePhiTrigonometry() 122 { 123 hDPhi = 0.5*fDPhi; // 124 cPhi = fSPhi + hDPhi; 125 ePhi = fSPhi + fDPhi; 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 } 137 138 inline 139 void G4Sphere::InitializeThetaTrigonometry() 140 { 141 eTheta = fSTheta + fDTheta; 142 143 sinSTheta = std::sin(fSTheta); 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 } 153 154 inline 155 void G4Sphere::CheckThetaAngles(G4double sThet 156 { 157 if ( (sTheta<0) || (sTheta>CLHEP::pi) ) 158 { 159 std::ostringstream message; 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 } 191 192 inline 193 void G4Sphere::CheckSPhiAngle(G4double sPhi) 194 { 195 // Ensure fSphi in 0-2PI or -2PI-0 range if 196 197 if ( sPhi < 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 } 210 211 inline 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 238 inline 239 void G4Sphere::CheckPhiAngles(G4double sPhi, G 240 { 241 CheckDPhiAngle(dPhi); 242 if (!fFullPhiSphere && (sPhi != 0.0)) { Chec 243 fFullSphere = fFullPhiSphere && fFullThetaSp 244 245 InitializePhiTrigonometry(); 246 } 247 248 inline 249 void G4Sphere::SetInnerRadius(G4double newRmin 250 { 251 fRmin= newRmin; 252 fRminTolerance = (fRmin) != 0.0 ? std::max( 253 Initialize(); 254 } 255 256 inline 257 void G4Sphere::SetOuterRadius(G4double newRmax 258 { 259 fRmax= newRmax; 260 fRmaxTolerance = std::max( kRadTolerance, fE 261 Initialize(); 262 } 263 264 inline 265 void G4Sphere::SetStartPhiAngle(G4double newSP 266 { 267 // Flag 'compute' can be used to explicitely 268 // trigonometry in case SetDeltaPhiAngle() i 269 270 CheckSPhiAngle(newSPhi); 271 fFullPhiSphere = false; 272 if (compute) { InitializePhiTrigonometry(); 273 Initialize(); 274 } 275 276 inline 277 void G4Sphere::SetDeltaPhiAngle(G4double newDP 278 { 279 CheckPhiAngles(fSPhi, newDPhi); 280 Initialize(); 281 } 282 283 inline 284 void G4Sphere::SetStartThetaAngle(G4double new 285 { 286 CheckThetaAngles(newSTheta, fDTheta); 287 Initialize(); 288 } 289 290 inline 291 void G4Sphere::SetDeltaThetaAngle(G4double new 292 { 293 CheckThetaAngles(fSTheta, newDTheta); 294 Initialize(); 295 } 296