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 // G4Para 26 // G4Para 27 // 27 // 28 // Class description: 28 // Class description: 29 // 29 // 30 // A G4Parallepiped, essentially a box with 30 // A G4Parallepiped, essentially a box with half lengths dx,dy,dz 31 // `skewed' so that there are angles theta & 31 // `skewed' so that there are angles theta & phi of the polar line 32 // joining the faces at +-dz in z, and alpha 32 // joining the faces at +-dz in z, and alpha formed by the y axis 33 // and the plane joinng the centre of the fa 33 // and the plane joinng the centre of the faces G4Parallel to the 34 // z-x plane at -dy and +dy. 34 // z-x plane at -dy and +dy. 35 // 35 // 36 // A G4Para is defined by: 36 // A G4Para is defined by: 37 // dx,dy,dz - Half-length in x,y,z 37 // dx,dy,dz - Half-length in x,y,z 38 // alpha - Angle formed by the y axis and 38 // alpha - Angle formed by the y axis and by the plane joining 39 // the centre of the faces G4Para 39 // the centre of the faces G4Parallel to the z-x plane 40 // at -dy and +dy 40 // at -dy and +dy 41 // theta - Polar angle of the line joinin 41 // theta - Polar angle of the line joining the centres of the 42 // faces at -dz and +dz in z 42 // faces at -dz and +dz in z 43 // phi - Azimuthal angle of the line jo 43 // phi - Azimuthal angle of the line joining the centres of the 44 // faces at -dz and +dz in z 44 // faces at -dz and +dz in z 45 // Member data: 45 // Member data: 46 // 46 // 47 // Note that the angles parameters are not s 47 // Note that the angles parameters are not stored - precomputed trig is 48 // stored instead. 48 // stored instead. 49 // 49 // 50 // fDx Half-length in x 50 // fDx Half-length in x 51 // fDy Half-length in y 51 // fDy Half-length in y 52 // fDz Half-length in z 52 // fDz Half-length in z 53 // 53 // 54 // fTalpha Tan of alpha 54 // fTalpha Tan of alpha 55 // fTthetaCphi Tan theta * Cos phi 55 // fTthetaCphi Tan theta * Cos phi 56 // fTthetaSphi Tan theta * Sin phi 56 // fTthetaSphi Tan theta * Sin phi 57 57 58 // 21.3.94 P.Kent Old C++ code converted to to 58 // 21.3.94 P.Kent Old C++ code converted to tolerant geometry 59 // 31.10.96 V.Grichine Modifications according 59 // 31.10.96 V.Grichine Modifications according G4Box/Tubs before to commit 60 // ------------------------------------------- 60 // -------------------------------------------------------------------- 61 #ifndef G4PARA_HH 61 #ifndef G4PARA_HH 62 #define G4PARA_HH 62 #define G4PARA_HH 63 63 64 #include "G4GeomTypes.hh" 64 #include "G4GeomTypes.hh" 65 65 66 #if defined(G4GEOM_USE_USOLIDS) 66 #if defined(G4GEOM_USE_USOLIDS) 67 #define G4GEOM_USE_UPARA 1 67 #define G4GEOM_USE_UPARA 1 68 #endif 68 #endif 69 69 70 #if defined(G4GEOM_USE_UPARA) 70 #if defined(G4GEOM_USE_UPARA) 71 #define G4UPara G4Para 71 #define G4UPara G4Para 72 #include "G4UPara.hh" 72 #include "G4UPara.hh" 73 #else 73 #else 74 74 75 #include "G4CSGSolid.hh" 75 #include "G4CSGSolid.hh" 76 #include "G4Polyhedron.hh" 76 #include "G4Polyhedron.hh" 77 77 78 class G4Para : public G4CSGSolid 78 class G4Para : public G4CSGSolid 79 { 79 { 80 public: // with description 80 public: // with description 81 81 82 G4Para(const G4String& pName, 82 G4Para(const G4String& pName, 83 G4double pDx, G4double pDy, G 83 G4double pDx, G4double pDy, G4double pDz, 84 G4double pAlpha, G4double pTh 84 G4double pAlpha, G4double pTheta, G4double pPhi); 85 85 86 G4Para(const G4String& pName, 86 G4Para(const G4String& pName, 87 const G4ThreeVector pt[8]); 87 const G4ThreeVector pt[8]); 88 88 89 ~G4Para() override; 89 ~G4Para() override; 90 90 91 // Accessors 91 // Accessors 92 92 93 inline G4double GetZHalfLength() const; 93 inline G4double GetZHalfLength() const; 94 inline G4ThreeVector GetSymAxis() const; 94 inline G4ThreeVector GetSymAxis() const; 95 inline G4double GetYHalfLength() const; 95 inline G4double GetYHalfLength() const; 96 inline G4double GetXHalfLength() const; 96 inline G4double GetXHalfLength() const; 97 inline G4double GetTanAlpha() const; 97 inline G4double GetTanAlpha() const; 98 98 99 inline G4double GetAlpha() const; 99 inline G4double GetAlpha() const; 100 inline G4double GetTheta() const; 100 inline G4double GetTheta() const; 101 inline G4double GetPhi() const; 101 inline G4double GetPhi() const; 102 // Obtain (re)computed values of original 102 // Obtain (re)computed values of original parameters 103 103 104 // Modifiers 104 // Modifiers 105 105 106 inline void SetXHalfLength(G4double val); 106 inline void SetXHalfLength(G4double val); 107 inline void SetYHalfLength(G4double val); 107 inline void SetYHalfLength(G4double val); 108 inline void SetZHalfLength(G4double val); 108 inline void SetZHalfLength(G4double val); 109 inline void SetAlpha(G4double alpha); 109 inline void SetAlpha(G4double alpha); 110 inline void SetTanAlpha(G4double val); 110 inline void SetTanAlpha(G4double val); 111 inline void SetThetaAndPhi(G4double pTheta 111 inline void SetThetaAndPhi(G4double pTheta, G4double pPhi); 112 112 113 void SetAllParameters(G4double pDx, G4doub 113 void SetAllParameters(G4double pDx, G4double pDy, G4double pDz, 114 G4double pAlpha, G4d 114 G4double pAlpha, G4double pTheta, G4double pPhi); 115 115 116 // Methods of solid 116 // Methods of solid 117 117 118 G4double GetCubicVolume() override; 118 G4double GetCubicVolume() override; 119 G4double GetSurfaceArea() override; 119 G4double GetSurfaceArea() override; 120 120 121 void ComputeDimensions(G4VPVParameterisati 121 void ComputeDimensions(G4VPVParameterisation* p, 122 const G4int n, 122 const G4int n, 123 const G4VPhysicalVo 123 const G4VPhysicalVolume* pRep) override; 124 124 125 void BoundingLimits(G4ThreeVector& pMin, G 125 void BoundingLimits(G4ThreeVector& pMin, G4ThreeVector& pMax) const override; 126 126 127 G4bool CalculateExtent(const EAxis pAxis, 127 G4bool CalculateExtent(const EAxis pAxis, 128 const G4VoxelLimits 128 const G4VoxelLimits& pVoxelLimit, 129 const G4AffineTrans 129 const G4AffineTransform& pTransform, 130 G4double& pMi 130 G4double& pMin, G4double& pMax) const override; 131 131 132 EInside Inside(const G4ThreeVector& p) con 132 EInside Inside(const G4ThreeVector& p) const override; 133 133 134 G4ThreeVector SurfaceNormal( const G4Three 134 G4ThreeVector SurfaceNormal( const G4ThreeVector& p) const override; 135 135 136 G4double DistanceToIn(const G4ThreeVector& 136 G4double DistanceToIn(const G4ThreeVector& p, 137 const G4ThreeVector& 137 const G4ThreeVector& v) const override; 138 G4double DistanceToIn(const G4ThreeVector& 138 G4double DistanceToIn(const G4ThreeVector& p) const override; 139 139 140 G4double DistanceToOut(const G4ThreeVector 140 G4double DistanceToOut(const G4ThreeVector& p, const G4ThreeVector& v, 141 const G4bool calcNo 141 const G4bool calcNorm = false, 142 G4bool* valid 142 G4bool* validNorm = nullptr, 143 G4ThreeVector 143 G4ThreeVector* n = nullptr) const override; 144 G4double DistanceToOut(const G4ThreeVector 144 G4double DistanceToOut(const G4ThreeVector& p) const override; 145 145 146 G4GeometryType GetEntityType() const overr 146 G4GeometryType GetEntityType() const override; 147 147 148 G4ThreeVector GetPointOnSurface() const ov 148 G4ThreeVector GetPointOnSurface() const override; 149 149 150 G4bool IsFaceted() const override; << 151 << 152 G4VSolid* Clone() const override; 150 G4VSolid* Clone() const override; 153 151 154 std::ostream& StreamInfo(std::ostream& os) 152 std::ostream& StreamInfo(std::ostream& os) const override; 155 153 156 // Visualisation functions 154 // Visualisation functions 157 155 158 void DescribeYourselfTo (G4VGraph 156 void DescribeYourselfTo (G4VGraphicsScene& scene) const override; 159 G4Polyhedron* CreatePolyhedron () const 157 G4Polyhedron* CreatePolyhedron () const override; 160 158 161 G4Para(__void__&); 159 G4Para(__void__&); 162 // Fake default constructor for usage re 160 // Fake default constructor for usage restricted to direct object 163 // persistency for clients requiring pre 161 // persistency for clients requiring preallocation of memory for 164 // persistifiable objects 162 // persistifiable objects 165 163 166 G4Para(const G4Para& rhs); 164 G4Para(const G4Para& rhs); 167 G4Para& operator=(const G4Para& rhs); 165 G4Para& operator=(const G4Para& rhs); 168 // Copy constructor and assignment opera 166 // Copy constructor and assignment operator 169 167 170 private: 168 private: 171 169 172 void CheckParameters(); 170 void CheckParameters(); 173 // Check parameters 171 // Check parameters 174 172 175 void MakePlanes(); 173 void MakePlanes(); 176 // Set side planes 174 // Set side planes 177 175 178 G4ThreeVector ApproxSurfaceNormal(const G4 176 G4ThreeVector ApproxSurfaceNormal(const G4ThreeVector& p) const; 179 // Algorithm for SurfaceNormal() followi 177 // Algorithm for SurfaceNormal() following the original 180 // specification for points not on the s 178 // specification for points not on the surface 181 179 182 180 183 private: 181 private: 184 182 185 G4double halfCarTolerance; 183 G4double halfCarTolerance; 186 G4double fDx,fDy,fDz; 184 G4double fDx,fDy,fDz; 187 G4double fTalpha,fTthetaCphi,fTthetaSphi; 185 G4double fTalpha,fTthetaCphi,fTthetaSphi; 188 struct { G4double a,b,c,d; } fPlanes[4]; 186 struct { G4double a,b,c,d; } fPlanes[4]; 189 }; 187 }; 190 188 191 #include "G4Para.icc" 189 #include "G4Para.icc" 192 190 193 #endif 191 #endif 194 192 195 #endif 193 #endif 196 194