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 G4CSGSolid 27 // 28 // 27.03.98 J.Apostolakis - First version. 29 // ------------------------------------------- 30 31 #include <cmath> 32 33 #include "G4CSGSolid.hh" 34 #include "G4QuickRand.hh" 35 #include "G4Polyhedron.hh" 36 37 #include "G4AutoLock.hh" 38 39 namespace 40 { 41 G4Mutex polyhedronMutex = G4MUTEX_INITIALIZE 42 } 43 44 ////////////////////////////////////////////// 45 // 46 // Constructor 47 // - Base class constructor 48 49 G4CSGSolid::G4CSGSolid(const G4String& name) : 50 G4VSolid(name) 51 { 52 } 53 54 ////////////////////////////////////////////// 55 // 56 // Fake default constructor - sets only member 57 // for usage restri 58 59 G4CSGSolid::G4CSGSolid( __void__& a ) 60 : G4VSolid(a) 61 { 62 } 63 64 ////////////////////////////////////////////// 65 // 66 // Destructor 67 // 68 69 G4CSGSolid::~G4CSGSolid() 70 { 71 delete fpPolyhedron; fpPolyhedron = nullptr; 72 } 73 74 ////////////////////////////////////////////// 75 // 76 // Copy constructor 77 // 78 79 G4CSGSolid::G4CSGSolid(const G4CSGSolid& rhs) 80 : G4VSolid(rhs), fCubicVolume(rhs.fCubicVolu 81 fSurfaceArea(rhs.fSurfaceArea) 82 { 83 } 84 85 ////////////////////////////////////////////// 86 // 87 // Assignment operator 88 89 G4CSGSolid& G4CSGSolid::operator = (const G4CS 90 { 91 // Check assignment to self 92 // 93 if (this == &rhs) { return *this; } 94 95 // Copy base class data 96 // 97 G4VSolid::operator=(rhs); 98 99 // Copy data 100 // 101 fCubicVolume = rhs.fCubicVolume; 102 fSurfaceArea = rhs.fSurfaceArea; 103 fRebuildPolyhedron = false; 104 delete fpPolyhedron; fpPolyhedron = nullptr 105 106 return *this; 107 } 108 109 G4double G4CSGSolid::GetRadiusInRing(G4double 110 { 111 G4double k = G4QuickRand(); 112 return (rmin <= 0) ? rmax*std::sqrt(k) 113 : std::sqrt(k*rmax*rmax + 114 } 115 116 std::ostream& G4CSGSolid::StreamInfo(std::ostr 117 { 118 os << "------------------------------------- 119 << " *** Dump for solid - " << GetName 120 << " ================================= 121 << " Solid type: " << GetEntityType() << 122 << " Parameters: \n" 123 << " NOT available !\n" 124 << "------------------------------------- 125 126 return os; 127 } 128 129 G4Polyhedron* G4CSGSolid::GetPolyhedron () con 130 { 131 if (fpPolyhedron == nullptr || 132 fRebuildPolyhedron || 133 fpPolyhedron->GetNumberOfRotationStepsAt 134 fpPolyhedron->GetNumberOfRotationSteps() 135 { 136 G4AutoLock l(&polyhedronMutex); 137 delete fpPolyhedron; 138 fpPolyhedron = CreatePolyhedron(); 139 fRebuildPolyhedron = false; 140 l.unlock(); 141 } 142 return fpPolyhedron; 143 } 144