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 G4CutTu 27 // ------------------------------------------- 28 29 inline 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 { 74 return sinEPhi; 75 } 76 77 inline 78 G4double G4CutTubs::GetCosEndPhi () const 79 { 80 return cosEPhi; 81 } 82 83 inline 84 G4ThreeVector G4CutTubs::GetLowNorm () const 85 { 86 return fLowNorm; 87 } 88 89 inline 90 G4ThreeVector G4CutTubs::GetHighNorm () const 91 { 92 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 } 245