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 for G4UEllipsoid wrapper cla 27 // 28 // 13-08-2019 Gabriele Cosmo, CERN 29 // ------------------------------------------- 30 31 #include "G4Ellipsoid.hh" 32 #include "G4UEllipsoid.hh" 33 34 #if ( defined(G4GEOM_USE_USOLIDS) || defined(G 35 36 #include "G4AffineTransform.hh" 37 #include "G4VPVParameterisation.hh" 38 #include "G4PhysicalConstants.hh" 39 #include "G4BoundingEnvelope.hh" 40 #include "G4Polyhedron.hh" 41 42 ////////////////////////////////////////////// 43 // 44 // Constructor - check & set half widths 45 46 47 G4UEllipsoid::G4UEllipsoid(const G4String& pNa 48 G4d 49 G4d 50 G4d 51 G4d 52 G4d 53 : Base_t(pName, dx, dy, dz, bcut, tcut) 54 { } 55 56 ////////////////////////////////////////////// 57 // 58 // Fake default constructor - sets only member 59 // for usage restri 60 61 G4UEllipsoid::G4UEllipsoid( __void__& a ) 62 : Base_t(a) 63 { } 64 65 ////////////////////////////////////////////// 66 // 67 // Destructor 68 69 G4UEllipsoid::~G4UEllipsoid() = default; 70 71 ////////////////////////////////////////////// 72 // 73 // Copy constructor 74 75 G4UEllipsoid::G4UEllipsoid(const G4UEllipsoid& 76 : Base_t(rhs) 77 { } 78 79 ////////////////////////////////////////////// 80 // 81 // Assignment operator 82 83 G4UEllipsoid& G4UEllipsoid::operator = (const 84 { 85 // Check assignment to self 86 // 87 if (this == &rhs) { return *this; } 88 89 // Copy base class data 90 // 91 Base_t::operator=(rhs); 92 93 return *this; 94 } 95 96 ////////////////////////////////////////////// 97 // 98 // Accessors 99 100 G4double G4UEllipsoid::GetDx() const 101 { 102 return Base_t::GetDx(); 103 } 104 105 G4double G4UEllipsoid::GetDy() const 106 { 107 return Base_t::GetDy(); 108 } 109 110 G4double G4UEllipsoid::GetDz() const 111 { 112 return Base_t::GetDz(); 113 } 114 115 G4double G4UEllipsoid::GetSemiAxisMax (G4int i 116 { 117 return (i==0) ? GetDx() 118 : (i==1) ? GetDy() 119 : GetDz(); 120 } 121 122 G4double G4UEllipsoid::GetZBottomCut() const 123 { 124 return Base_t::GetZBottomCut(); 125 } 126 127 G4double G4UEllipsoid::GetZTopCut() const 128 { 129 return Base_t::GetZTopCut(); 130 } 131 132 ////////////////////////////////////////////// 133 // 134 // Modifiers 135 136 void G4UEllipsoid::SetSemiAxis (G4double x, G4 137 { 138 Base_t::SetSemiAxes(x, y, z); 139 } 140 141 void G4UEllipsoid::SetZCuts (G4double newzBott 142 { 143 Base_t::SetZCuts(newzBottomCut, newzTopCut); 144 } 145 146 ////////////////////////////////////////////// 147 // 148 // Make a clone of the object 149 150 G4VSolid* G4UEllipsoid::Clone() const 151 { 152 return new G4UEllipsoid(*this); 153 } 154 155 ////////////////////////////////////////////// 156 // 157 // Get bounding box 158 159 void G4UEllipsoid::BoundingLimits(G4ThreeVecto 160 G4ThreeVecto 161 { 162 G4double dx = GetDx(); 163 G4double dy = GetDy(); 164 G4double dz = GetDz(); 165 G4double zmin = std::max(-dz,GetZBottomCut() 166 G4double zmax = std::min( dz,GetZTopCut()); 167 pMin.set(-dx,-dy,zmin); 168 pMax.set( dx, dy,zmax); 169 170 // Check correctness of the bounding box 171 // 172 if (pMin.x() >= pMax.x() || pMin.y() >= pMax 173 { 174 std::ostringstream message; 175 message << "Bad bounding box (min >= max) 176 << GetName() << " !" 177 << "\npMin = " << pMin 178 << "\npMax = " << pMax; 179 G4Exception("G4UEllipsoid::BoundingLimits( 180 JustWarning, message); 181 StreamInfo(G4cout); 182 } 183 } 184 185 ////////////////////////////////////////////// 186 // 187 // Calculate extent under transform and specif 188 189 G4bool 190 G4UEllipsoid::CalculateExtent(const EAxis pAxi 191 const G4VoxelLim 192 const G4AffineTr 193 G4double& 194 { 195 G4ThreeVector bmin, bmax; 196 197 // Get bounding box 198 BoundingLimits(bmin,bmax); 199 200 // Find extent 201 G4BoundingEnvelope bbox(bmin,bmax); 202 return bbox.CalculateExtent(pAxis,pVoxelLimi 203 } 204 205 ////////////////////////////////////////////// 206 // 207 // CreatePolyhedron 208 209 G4Polyhedron* G4UEllipsoid::CreatePolyhedron() 210 { 211 return new G4PolyhedronEllipsoid(GetDx(), Ge 212 GetZBottomC 213 } 214 215 #endif // G4GEOM_USE_USOLIDS 216