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" << 45 #include "globals.hh" 44 #include "globals.hh" 46 << 45 #include "G4ThreeVector.hh" 47 #include <vector> 46 #include <vector> 48 47 49 class G4Step; 48 class G4Step; 50 49 51 class G4AdjointCrossSurfChecker << 50 class G4AdjointCrossSurfChecker 52 { << 51 { 53 public: << 52 public: 54 static G4AdjointCrossSurfChecker* GetInstanc << 53 55 << 54 static G4AdjointCrossSurfChecker* GetInstance(); 56 G4bool CrossingASphere(const G4Step* aStep, << 55 57 G4ThreeVector& cross_pos, G4double& cos_to << 56 G4bool CrossingASphere(const G4Step* aStep, 58 << 57 G4double sphere_radius, 59 G4bool GoingInOrOutOfaVolume( << 58 G4ThreeVector sphere_center, 60 const G4Step* aStep, const G4String& volum << 59 G4ThreeVector& cross_pos, 61 << 60 G4double& cos_to_surface, 62 G4bool GoingInOrOutOfaVolumeByExtSurface(con << 61 G4bool& GoingIn); 63 const G4String& mother_lvol_name, G4double << 62 64 << 63 G4bool GoingInOrOutOfaVolume(const G4Step* aStep, 65 G4bool CrossingAGivenRegisteredSurface(const << 64 const G4String& volume_name, 66 G4ThreeVector& cross_pos, G4double& cos_to << 65 G4double& cos_to_surface, 67 << 66 G4bool& GoingIn); 68 G4bool CrossingAGivenRegisteredSurface(const << 67 69 G4double& cos_to_surface, G4bool& GoingIn) << 68 G4bool GoingInOrOutOfaVolumeByExtSurface(const G4Step* aStep, 70 << 69 const G4String& volume_name, 71 G4bool CrossingOneOfTheRegisteredSurface(con << 70 const G4String& mother_lvol_name, 72 G4ThreeVector& cross_pos, G4double& cos_to << 71 G4double& cos_to_surface, 73 << 72 G4bool& GoingIn); 74 G4bool CrossingAnInterfaceBetweenTwoVolumes( << 73 75 const G4String& vol2_name, G4ThreeVector& << 74 G4bool CrossingAGivenRegisteredSurface(const G4Step* aStep, 76 << 75 const G4String& surface_name, 77 G4bool AddaSphericalSurface( << 76 G4ThreeVector& cross_pos, 78 const G4String& SurfaceName, G4double radi << 77 G4double& cos_to_surface, 79 G4bool AddaSphericalSurfaceWithCenterAtTheCe << 78 G4bool& GoingIn); 80 G4double radius, const G4String& volume_na << 79 81 G4bool AddanExtSurfaceOfAvolume( << 80 G4bool CrossingAGivenRegisteredSurface(const G4Step* aStep, 82 const G4String& SurfaceName, const G4Strin << 81 G4int ind, 83 G4bool AddanInterfaceBetweenTwoVolumes(const << 82 G4ThreeVector& cross_pos, 84 const G4String& volume_name2, G4double& ar << 83 G4double& cos_to_surface, 85 void ClearListOfSelectedSurface(); << 84 G4bool& GoingIn); 86 << 85 87 private: << 86 G4bool CrossingOneOfTheRegisteredSurface(const G4Step* aStep, 88 G4AdjointCrossSurfChecker() = default; << 87 G4String& surface_name, 89 ~G4AdjointCrossSurfChecker(); << 88 G4ThreeVector& cross_pos, 90 << 89 G4double& cos_to_surface, 91 G4int FindRegisteredSurface(const G4String& << 90 G4bool& GoingIn); 92 << 91 93 private: << 92 G4bool CrossingAnInterfaceBetweenTwoVolumes(const G4Step* aStep, 94 static G4ThreadLocal G4AdjointCrossSurfCheck << 93 const G4String& vol1_name, 95 << 94 const G4String& vol2_name, 96 std::vector<G4String> ListOfSurfaceName; << 95 G4ThreeVector& cross_pos, 97 std::vector<G4String> ListOfSurfaceType; << 96 G4double& cos_to_surface, 98 std::vector<G4double> ListOfSphereRadius; << 97 G4bool& GoingIn); 99 std::vector<G4ThreeVector> ListOfSphereCente << 98 100 std::vector<G4String> ListOfVol1Name; << 99 G4bool AddaSphericalSurface(const G4String& SurfaceName, 101 std::vector<G4String> ListOfVol2Name; << 100 G4double radius, 102 std::vector<G4double> AreaOfSurface; << 101 G4ThreeVector pos, >> 102 G4double& area); >> 103 G4bool AddaSphericalSurfaceWithCenterAtTheCenterOfAVolume( >> 104 const G4String& SurfaceName, >> 105 G4double radius, >> 106 const G4String& volume_name, >> 107 G4ThreeVector& center, >> 108 G4double& area); >> 109 G4bool AddanExtSurfaceOfAvolume(const G4String& SurfaceName, >> 110 const G4String& volume_name, >> 111 G4double& area); >> 112 G4bool AddanInterfaceBetweenTwoVolumes(const G4String& SurfaceName, >> 113 const G4String& volume_name1, >> 114 const G4String& volume_name2, >> 115 G4double& area); >> 116 void ClearListOfSelectedSurface(); >> 117 >> 118 private: >> 119 >> 120 G4AdjointCrossSurfChecker(); >> 121 ~G4AdjointCrossSurfChecker(); >> 122 >> 123 G4int FindRegisteredSurface(const G4String& name); >> 124 >> 125 private: >> 126 >> 127 static G4ThreadLocal G4AdjointCrossSurfChecker* instance; >> 128 >> 129 std::vector<G4String> ListOfSurfaceName; >> 130 std::vector<G4String> ListOfSurfaceType; >> 131 std::vector<G4double> ListOfSphereRadius; >> 132 std::vector<G4ThreeVector> ListOfSphereCenter; >> 133 std::vector<G4String> ListOfVol1Name; >> 134 std::vector<G4String> ListOfVol2Name; >> 135 std::vector<G4double> AreaOfSurface; 103 }; 136 }; 104 137 105 #endif 138 #endif 106 139