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