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 // G4AdjointCrossSurfChecker 26 // G4AdjointCrossSurfChecker 27 // 27 // 28 // Class Description: 28 // Class Description: 29 // 29 // 30 // This class is responsible for checking the 30 // This class is responsible for checking the crossing of a surface 31 // that could be the external boundary of a v 31 // that could be the external boundary of a volume or the external 32 // surface of a sphere. 32 // surface of a sphere. 33 // It is used to check if an adjoint particle 33 // It is used to check if an adjoint particle reaches the external 34 // surface or reenter the sensitive region del 34 // surface or reenter the sensitive region delimited by the adjoint 35 // source. 35 // source. 36 36 37 // Author: L. Desorgher, SpaceIT GmbH - 15 Jan 37 // Author: L. Desorgher, SpaceIT GmbH - 15 January 2007 38 // Contract: ESA contract 21435/08/NL/AT 38 // Contract: ESA contract 21435/08/NL/AT 39 // Customer: ESA/ESTEC 39 // Customer: ESA/ESTEC 40 //-------------------------------------------- 40 //--------------------------------------------------------------------- 41 #ifndef G4AdjointCrossSurfChecker_hh 41 #ifndef G4AdjointCrossSurfChecker_hh 42 #define G4AdjointCrossSurfChecker_hh 1 42 #define G4AdjointCrossSurfChecker_hh 1 43 43 44 #include "G4ThreeVector.hh" 44 #include "G4ThreeVector.hh" 45 #include "globals.hh" 45 #include "globals.hh" 46 46 47 #include <vector> 47 #include <vector> 48 48 49 class G4Step; 49 class G4Step; 50 50 51 class G4AdjointCrossSurfChecker 51 class G4AdjointCrossSurfChecker 52 { 52 { 53 public: 53 public: 54 static G4AdjointCrossSurfChecker* GetInstanc 54 static G4AdjointCrossSurfChecker* GetInstance(); 55 55 56 G4bool CrossingASphere(const G4Step* aStep, 56 G4bool CrossingASphere(const G4Step* aStep, G4double sphere_radius, G4ThreeVector sphere_center, 57 G4ThreeVector& cross_pos, G4double& cos_to 57 G4ThreeVector& cross_pos, G4double& cos_to_surface, G4bool& GoingIn); 58 58 59 G4bool GoingInOrOutOfaVolume( 59 G4bool GoingInOrOutOfaVolume( 60 const G4Step* aStep, const G4String& volum 60 const G4Step* aStep, const G4String& volume_name, G4double& cos_to_surface, G4bool& GoingIn); 61 61 62 G4bool GoingInOrOutOfaVolumeByExtSurface(con 62 G4bool GoingInOrOutOfaVolumeByExtSurface(const G4Step* aStep, const G4String& volume_name, 63 const G4String& mother_lvol_name, G4double 63 const G4String& mother_lvol_name, G4double& cos_to_surface, G4bool& GoingIn); 64 64 65 G4bool CrossingAGivenRegisteredSurface(const 65 G4bool CrossingAGivenRegisteredSurface(const G4Step* aStep, const G4String& surface_name, 66 G4ThreeVector& cross_pos, G4double& cos_to 66 G4ThreeVector& cross_pos, G4double& cos_to_surface, G4bool& GoingIn); 67 67 68 G4bool CrossingAGivenRegisteredSurface(const 68 G4bool CrossingAGivenRegisteredSurface(const G4Step* aStep, G4int ind, G4ThreeVector& cross_pos, 69 G4double& cos_to_surface, G4bool& GoingIn) 69 G4double& cos_to_surface, G4bool& GoingIn); 70 70 71 G4bool CrossingOneOfTheRegisteredSurface(con 71 G4bool CrossingOneOfTheRegisteredSurface(const G4Step* aStep, G4String& surface_name, 72 G4ThreeVector& cross_pos, G4double& cos_to 72 G4ThreeVector& cross_pos, G4double& cos_to_surface, G4bool& GoingIn); 73 73 74 G4bool CrossingAnInterfaceBetweenTwoVolumes( 74 G4bool CrossingAnInterfaceBetweenTwoVolumes(const G4Step* aStep, const G4String& vol1_name, 75 const G4String& vol2_name, G4ThreeVector& 75 const G4String& vol2_name, G4ThreeVector& cross_pos, G4double& cos_to_surface, G4bool& GoingIn); 76 76 77 G4bool AddaSphericalSurface( 77 G4bool AddaSphericalSurface( 78 const G4String& SurfaceName, G4double radi 78 const G4String& SurfaceName, G4double radius, G4ThreeVector pos, G4double& area); 79 G4bool AddaSphericalSurfaceWithCenterAtTheCe 79 G4bool AddaSphericalSurfaceWithCenterAtTheCenterOfAVolume(const G4String& SurfaceName, 80 G4double radius, const G4String& volume_na 80 G4double radius, const G4String& volume_name, G4ThreeVector& center, G4double& area); 81 G4bool AddanExtSurfaceOfAvolume( 81 G4bool AddanExtSurfaceOfAvolume( 82 const G4String& SurfaceName, const G4Strin 82 const G4String& SurfaceName, const G4String& volume_name, G4double& area); 83 G4bool AddanInterfaceBetweenTwoVolumes(const 83 G4bool AddanInterfaceBetweenTwoVolumes(const G4String& SurfaceName, const G4String& volume_name1, 84 const G4String& volume_name2, G4double& ar 84 const G4String& volume_name2, G4double& area); 85 void ClearListOfSelectedSurface(); 85 void ClearListOfSelectedSurface(); 86 86 87 private: 87 private: 88 G4AdjointCrossSurfChecker() = default; 88 G4AdjointCrossSurfChecker() = default; 89 ~G4AdjointCrossSurfChecker(); 89 ~G4AdjointCrossSurfChecker(); 90 90 91 G4int FindRegisteredSurface(const G4String& 91 G4int FindRegisteredSurface(const G4String& name); 92 92 93 private: 93 private: 94 static G4ThreadLocal G4AdjointCrossSurfCheck 94 static G4ThreadLocal G4AdjointCrossSurfChecker* instance; 95 95 96 std::vector<G4String> ListOfSurfaceName; 96 std::vector<G4String> ListOfSurfaceName; 97 std::vector<G4String> ListOfSurfaceType; 97 std::vector<G4String> ListOfSurfaceType; 98 std::vector<G4double> ListOfSphereRadius; 98 std::vector<G4double> ListOfSphereRadius; 99 std::vector<G4ThreeVector> ListOfSphereCente 99 std::vector<G4ThreeVector> ListOfSphereCenter; 100 std::vector<G4String> ListOfVol1Name; 100 std::vector<G4String> ListOfVol1Name; 101 std::vector<G4String> ListOfVol2Name; 101 std::vector<G4String> ListOfVol2Name; 102 std::vector<G4double> AreaOfSurface; 102 std::vector<G4double> AreaOfSurface; 103 }; 103 }; 104 104 105 #endif 105 #endif 106 106