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 G4Para 26 // Implementation of inline methods of G4Para 27 // ------------------------------------------- 27 // -------------------------------------------------------------------- 28 28 29 inline 29 inline 30 G4double G4Para::GetZHalfLength() const 30 G4double G4Para::GetZHalfLength() const 31 { 31 { 32 return fDz; 32 return fDz; 33 } 33 } 34 34 35 inline 35 inline 36 G4ThreeVector G4Para::GetSymAxis() const 36 G4ThreeVector G4Para::GetSymAxis() const 37 { 37 { 38 return G4ThreeVector(fTthetaCphi,fTthetaSphi 38 return G4ThreeVector(fTthetaCphi,fTthetaSphi,1.).unit(); 39 } 39 } 40 40 41 inline 41 inline 42 G4double G4Para::GetYHalfLength() const 42 G4double G4Para::GetYHalfLength() const 43 { 43 { 44 return fDy; 44 return fDy; 45 } 45 } 46 46 47 inline 47 inline 48 G4double G4Para::GetXHalfLength() const 48 G4double G4Para::GetXHalfLength() const 49 { 49 { 50 return fDx; 50 return fDx; 51 } 51 } 52 52 53 inline 53 inline 54 G4double G4Para::GetTanAlpha() const 54 G4double G4Para::GetTanAlpha() const 55 { 55 { 56 return fTalpha; 56 return fTalpha; 57 } 57 } 58 58 59 inline 59 inline 60 void G4Para::SetXHalfLength(G4double val) 60 void G4Para::SetXHalfLength(G4double val) 61 { 61 { 62 fDx = val; 62 fDx = val; 63 63 64 fCubicVolume = 0.; 64 fCubicVolume = 0.; 65 fSurfaceArea = 0.; 65 fSurfaceArea = 0.; 66 fRebuildPolyhedron = true; 66 fRebuildPolyhedron = true; 67 67 68 CheckParameters(); 68 CheckParameters(); 69 MakePlanes(); 69 MakePlanes(); 70 } 70 } 71 71 72 inline 72 inline 73 void G4Para::SetYHalfLength(G4double val) 73 void G4Para::SetYHalfLength(G4double val) 74 { 74 { 75 fDy = val; 75 fDy = val; 76 76 77 fCubicVolume = 0.; 77 fCubicVolume = 0.; 78 fSurfaceArea = 0.; 78 fSurfaceArea = 0.; 79 fRebuildPolyhedron = true; 79 fRebuildPolyhedron = true; 80 80 81 CheckParameters(); 81 CheckParameters(); 82 MakePlanes(); 82 MakePlanes(); 83 } 83 } 84 84 85 inline 85 inline 86 void G4Para::SetZHalfLength(G4double val) 86 void G4Para::SetZHalfLength(G4double val) 87 { 87 { 88 fDz = val; 88 fDz = val; 89 89 90 fCubicVolume = 0.; 90 fCubicVolume = 0.; 91 fSurfaceArea = 0.; 91 fSurfaceArea = 0.; 92 fRebuildPolyhedron = true; 92 fRebuildPolyhedron = true; 93 93 94 CheckParameters(); 94 CheckParameters(); 95 MakePlanes(); 95 MakePlanes(); 96 } 96 } 97 97 98 inline 98 inline 99 void G4Para::SetAlpha(G4double alpha) 99 void G4Para::SetAlpha(G4double alpha) 100 { 100 { 101 fTalpha = std::tan(alpha); 101 fTalpha = std::tan(alpha); 102 102 103 fSurfaceArea = 0.; 103 fSurfaceArea = 0.; 104 fCubicVolume = 0.; 104 fCubicVolume = 0.; 105 fRebuildPolyhedron = true; 105 fRebuildPolyhedron = true; 106 106 107 MakePlanes(); 107 MakePlanes(); 108 } 108 } 109 109 110 inline 110 inline 111 void G4Para::SetTanAlpha(G4double val) 111 void G4Para::SetTanAlpha(G4double val) 112 { 112 { 113 fTalpha = val; 113 fTalpha = val; 114 114 115 fCubicVolume = 0.; 115 fCubicVolume = 0.; 116 fSurfaceArea = 0.; 116 fSurfaceArea = 0.; 117 fRebuildPolyhedron = true; 117 fRebuildPolyhedron = true; 118 118 119 MakePlanes(); 119 MakePlanes(); 120 } 120 } 121 121 122 inline 122 inline 123 void G4Para::SetThetaAndPhi(G4double pTheta, G << 123 void G4Para::SetThetaAndPhi(double pTheta, double pPhi) 124 { 124 { 125 G4double tanTheta = std::tan(pTheta); 125 G4double tanTheta = std::tan(pTheta); 126 fTthetaCphi = tanTheta*std::cos(pPhi); 126 fTthetaCphi = tanTheta*std::cos(pPhi); 127 fTthetaSphi = tanTheta*std::sin(pPhi); 127 fTthetaSphi = tanTheta*std::sin(pPhi); 128 128 129 fCubicVolume = 0.; 129 fCubicVolume = 0.; 130 fSurfaceArea = 0.; 130 fSurfaceArea = 0.; 131 fRebuildPolyhedron = true; 131 fRebuildPolyhedron = true; 132 132 133 MakePlanes(); 133 MakePlanes(); 134 } 134 } 135 135 136 inline G4double G4Para::GetPhi() const 136 inline G4double G4Para::GetPhi() const 137 { 137 { 138 return std::atan2(fTthetaSphi,fTthetaCphi); 138 return std::atan2(fTthetaSphi,fTthetaCphi); 139 } 139 } 140 140 141 inline G4double G4Para::GetTheta() const 141 inline G4double G4Para::GetTheta() const 142 { 142 { 143 return std::atan(std::sqrt(fTthetaCphi*fTth 143 return std::atan(std::sqrt(fTthetaCphi*fTthetaCphi 144 +fTthetaSphi*fTt 144 +fTthetaSphi*fTthetaSphi)); 145 } 145 } 146 146 147 inline G4double G4Para::GetAlpha() const 147 inline G4double G4Para::GetAlpha() const 148 { 148 { 149 return std::atan(fTalpha); 149 return std::atan(fTalpha); 150 } 150 } 151 151