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 27 // ------------------------------------------- 28 // GEANT 4 class header file 29 // 30 // 31 // FastAerosolSolid 32 // 33 // Class description: 34 // 35 // A FastAerosolSolid is a collection of fDr 36 // with positions set randomly by FastAeroso 37 // volume given by a bulk shape given by Fas 38 // 39 // The FastAerosol member fCloud handles the 40 // (finding nearest droplet, populating the 41 // efficient simulations. 42 // 43 // This class is heavily based on system sol 44 // 45 // Author: A.Knaian (ara@nklabs.com), N.MacFad 46 // ------------------------------------------- 47 48 #ifndef FastAerosolSolid_HH 49 #define FastAerosolSolid_HH 50 51 #include "FastAerosol.hh" 52 53 #include "G4Polyhedron.hh" 54 55 // rotations 56 #include <functional> 57 #include "G4RotationMatrix.hh" 58 59 class FastAerosolSolid : public G4VSolid 60 { 61 public: 62 FastAerosolSolid(const G4String& pName, FastA 63 G4VSolid* pDroplet, 64 std::function<G4RotationMatrix ( 65 66 FastAerosolSolid(const G4String& pName, FastA 67 G4VSolid* pDroplet); 68 ~FastAerosolSolid()=default; 69 70 // Access functions 71 inline G4double GetCubicVolume(); 72 inline G4double GetSurfaceArea(); 73 74 // Solid standard methods 75 G4bool CalculateExtent(const EAxis pAxis, 76 const G4VoxelLimits& pVoxelLimit, 77 const G4AffineTransform& pTransfor 78 G4double& pmin, G4double& pmax) co 79 80 EInside Inside(const G4ThreeVector& p) const; 81 82 G4ThreeVector SurfaceNormal(const G4ThreeVect 83 84 G4double DistanceToIn(const G4ThreeVector& p, 85 const G4ThreeVector& v) const; 86 87 G4double DistanceToIn(const G4ThreeVector& p) 88 G4double DistanceToOut(const G4ThreeVector& p 89 const G4ThreeVector& v, 90 const G4bool calcNorm=G 91 G4bool *validNorm=0, G4 92 93 G4double DistanceToOut(const G4ThreeVector& p 94 95 G4GeometryType GetEntityType() const; 96 97 G4VSolid* Clone() const; 98 99 std::ostream& StreamInfo(std::ostream& os) co 100 101 G4ThreeVector GetPointOnSurface() const; 102 103 G4Polyhedron* GetPolyhedron () const; 104 void DescribeYourselfTo(G4VGraphicsScene& sce 105 G4VisExtent GetExtent() const; 106 G4Polyhedron* CreatePolyhedron() const; 107 108 public: // without description 109 FastAerosolSolid(__void__&); 110 // 111 // Fake default constructor for usage restrict 112 // persistency for clients requiring prealloca 113 // persistifiable objects. 114 115 FastAerosolSolid(const FastAerosolSolid& rhs) 116 FastAerosolSolid& operator=(const FastAerosol 117 // Copy constructor and assignment operator. 118 119 inline void SetStepLim(G4double newLim); 120 121 private: 122 123 inline void Initialize(); 124 // 125 // Reset relevant values to zero 126 127 G4double fStepLim = DBL_MAX; // Maximum step 128 129 FastAerosol* fCloud; // FastAerosol which han 130 G4VSolid* fDroplet; // Droplet shape 131 G4VSolid* fBulk; // Aerosol bulk 132 133 G4double fR = 0.0; // Droplet bounding radius 134 135 G4double fVisDx, fVisDy, fVisDz; // Visual ex 136 137 G4double fCubicVolume = 0.0; // Cubic volume 138 G4double fSurfaceArea = 0.0; // Surface area 139 140 G4double farFromCloudDist; 141 142 std::function<G4RotationMatrix (G4ThreeVector 143 // rotation function 144 145 protected: // without description 146 147 mutable G4bool fRebuildPolyhedron; 148 mutable G4Polyhedron* fpPolyhedron; 149 }; 150 151 #include "FastAerosolSolid.icc" 152 153 #endif 154