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 // G4VoxelLimits << 26 // >> 27 // $Id: G4VoxelLimits.hh,v 1.9 2006/06/29 18:33:13 gunter Exp $ >> 28 // GEANT4 tag $Name: geant4-09-03-patch-02 $ >> 29 // >> 30 // class G4VoxelLimits 27 // 31 // 28 // Class description: 32 // Class description: 29 // 33 // 30 // Represents limitation/restrictions of space 34 // Represents limitation/restrictions of space, where restrictions 31 // are only made perpendicular to the cartesia 35 // are only made perpendicular to the cartesian axes. 32 // 36 // >> 37 // 33 // Member data: 38 // Member data: 34 // 39 // 35 // G4double fxAxisMin,fxAxisMax 40 // G4double fxAxisMin,fxAxisMax 36 // G4double fyAxisMin,fyAxisMax 41 // G4double fyAxisMin,fyAxisMax 37 // G4double fzAxisMin,fzAxisMax 42 // G4double fzAxisMin,fzAxisMax 38 // - The min and max values along each axis. + << 43 // - The min and max values along each axis. +-kInfinity if not restricted. >> 44 // >> 45 // >> 46 // Notes: >> 47 // >> 48 // Beware no break statements after returns in switch(pAxis)s. 39 49 40 // 13.07.95, P.Kent - Initial version. << 50 // History: >> 51 // 13.07.95 P.Kent Initial version. 41 // ------------------------------------------- 52 // -------------------------------------------------------------------- 42 #ifndef G4VOXELLIMITS_HH 53 #ifndef G4VOXELLIMITS_HH 43 #define G4VOXELLIMITS_HH 1 << 54 #define G4VOXELLIMITS_HH 44 55 45 #include "G4Types.hh" 56 #include "G4Types.hh" 46 #include "geomdefs.hh" 57 #include "geomdefs.hh" 47 58 48 #include "G4ThreeVector.hh" 59 #include "G4ThreeVector.hh" 49 60 50 #include <assert.h> 61 #include <assert.h> 51 62 52 class G4VoxelLimits 63 class G4VoxelLimits 53 { 64 { 54 public: << 65 public: // with description 55 66 56 G4VoxelLimits() = default; << 67 G4VoxelLimits(); 57 // Constructor - initialise to be unlimi 68 // Constructor - initialise to be unlimited. Volume unrestricted. 58 69 59 ~G4VoxelLimits() = default; << 70 ~G4VoxelLimits(); 60 // Destructor. No actions. 71 // Destructor. No actions. 61 72 62 void AddLimit(const EAxis pAxis, const G4d << 73 void AddLimit(const EAxis pAxis, const G4double pMin,const G4double pMax); 63 // Restrict the volume to between specif 74 // Restrict the volume to between specified min and max along the 64 // given axis. Cartesian axes only, pMin 75 // given axis. Cartesian axes only, pMin<=pMax. 65 76 66 G4double GetMaxXExtent() const; 77 G4double GetMaxXExtent() const; 67 // Return maximum x extent. 78 // Return maximum x extent. 68 G4double GetMaxYExtent() const; 79 G4double GetMaxYExtent() const; 69 // Return maximum y extent. 80 // Return maximum y extent. 70 G4double GetMaxZExtent() const; 81 G4double GetMaxZExtent() const; 71 // Return maximum z extent. 82 // Return maximum z extent. 72 83 73 G4double GetMinXExtent() const; 84 G4double GetMinXExtent() const; 74 // Return minimum x extent. 85 // Return minimum x extent. 75 G4double GetMinYExtent() const; 86 G4double GetMinYExtent() const; 76 // Return minimum y extent. 87 // Return minimum y extent. 77 G4double GetMinZExtent() const; 88 G4double GetMinZExtent() const; 78 // Return minimum z extent. 89 // Return minimum z extent. 79 90 80 G4double GetMaxExtent(const EAxis pAxis) c 91 G4double GetMaxExtent(const EAxis pAxis) const; 81 // Return maximum extent of volume along 92 // Return maximum extent of volume along specified axis. 82 G4double GetMinExtent(const EAxis pAxis) c 93 G4double GetMinExtent(const EAxis pAxis) const; 83 // Return minimum extent of volume along 94 // Return minimum extent of volume along specified axis. 84 95 85 G4bool IsXLimited() const; 96 G4bool IsXLimited() const; 86 // Return true if the x axis is limited. 97 // Return true if the x axis is limited. 87 G4bool IsYLimited() const; 98 G4bool IsYLimited() const; 88 // Return true if the y axis is limited. 99 // Return true if the y axis is limited. 89 G4bool IsZLimited() const; 100 G4bool IsZLimited() const; 90 // Return true if the z axis is limited. 101 // Return true if the z axis is limited. 91 102 92 G4bool IsLimited() const; 103 G4bool IsLimited() const; 93 // Return true if limited along any axis 104 // Return true if limited along any axis 94 G4bool IsLimited(const EAxis pAxis) const; 105 G4bool IsLimited(const EAxis pAxis) const; 95 // Return true if the specified axis is 106 // Return true if the specified axis is restricted/limited. 96 107 97 G4bool ClipToLimits(G4ThreeVector& pStart, << 108 G4bool ClipToLimits(G4ThreeVector& pStart,G4ThreeVector& pEnd) const; 98 // Clip the line segment pStart->pEnd to 109 // Clip the line segment pStart->pEnd to the volume described by the 99 // current limits. Return true if the li 110 // current limits. Return true if the line remains after clipping, 100 // else false, and leave the vectors in 111 // else false, and leave the vectors in an undefined state. 101 112 102 G4bool Inside(const G4ThreeVector& pVec) c 113 G4bool Inside(const G4ThreeVector& pVec) const; 103 // Return true if the specified vector i 114 // Return true if the specified vector is inside/on boundaries of limits. 104 115 105 G4int OutCode(const G4ThreeVector& pVec) c 116 G4int OutCode(const G4ThreeVector& pVec) const; 106 // Calculate the `outcode' for the speci 117 // Calculate the `outcode' for the specified vector. 107 // Intended for use during clipping agai 118 // Intended for use during clipping against the limits 108 // The bits are set given the following 119 // The bits are set given the following conditions: 109 // 0 pVec.x()<fxAxisMin && IsXLim 120 // 0 pVec.x()<fxAxisMin && IsXLimited() 110 // 1 pVec.x()>fxAxisMax && IsXLim 121 // 1 pVec.x()>fxAxisMax && IsXLimited() 111 // 2 pVec.y()<fyAxisMin && IsYLim 122 // 2 pVec.y()<fyAxisMin && IsYLimited() 112 // 3 pVec.y()>fyAxisMax && IsYLim 123 // 3 pVec.y()>fyAxisMax && IsYLimited() 113 // 4 pVec.z()<fzAxisMin && IsZLim 124 // 4 pVec.z()<fzAxisMin && IsZLimited() 114 // 5 pVec.z()>fzAxisMax && IsZLim 125 // 5 pVec.z()>fzAxisMax && IsZLimited() 115 126 116 private: 127 private: 117 128 118 G4double fxAxisMin = -kInfinity, fxAxisMax << 129 G4double fxAxisMin,fxAxisMax; 119 G4double fyAxisMin = -kInfinity, fyAxisMax << 130 G4double fyAxisMin,fyAxisMax; 120 G4double fzAxisMin = -kInfinity, fzAxisMax << 131 G4double fzAxisMin,fzAxisMax; 121 }; 132 }; 122 133 123 #include "G4VoxelLimits.icc" 134 #include "G4VoxelLimits.icc" 124 135 125 std::ostream& operator << (std::ostream& os, c 136 std::ostream& operator << (std::ostream& os, const G4VoxelLimits& pLim); 126 // Print the limits to the stream in the for 137 // Print the limits to the stream in the form: 127 // "{(xmin,xmax) (ymin,ymax) (zmin,zmax)}" 138 // "{(xmin,xmax) (ymin,ymax) (zmin,zmax)}" 128 // Replace (xmin,xmax) by (-,-) when not li 139 // Replace (xmin,xmax) by (-,-) when not limited. 129 140 130 #endif 141 #endif 131 142