| 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 // G4AdjointPosOnPhysVolGenerator 26 // G4AdjointPosOnPhysVolGenerator
27 // 27 //
28 // Class description: 28 // Class description:
29 // 29 //
30 // This class is responsible for the generatio 30 // This class is responsible for the generation of primary adjoint particles
31 // on the external surface of a user selected 31 // on the external surface of a user selected volume.
32 // The particles are generated uniformly on th 32 // The particles are generated uniformly on the surface with the angular
33 // distribution set to a cosine law relative t 33 // distribution set to a cosine law relative to normal of the surface.
34 // It is equivalent to the flux going in from 34 // It is equivalent to the flux going in from the surface if an isotropic flux
35 // is considered outside. 35 // is considered outside.
36 // It uses ray tracking technique and can be a 36 // It uses ray tracking technique and can be applied to all kind of convex
37 // volumes. Using the ray tracking technique t 37 // volumes. Using the ray tracking technique the area of the external surface
38 // is also computed. The area is needed to fix 38 // is also computed. The area is needed to fix the weight of the primary
39 // adjoint particle. 39 // adjoint particle.
40 // At the time of the development of this clas 40 // At the time of the development of this class, generation of points on
41 // volume surface and computation of surface w 41 // volume surface and computation of surface was limited in Geant4, therefore
42 // the general ray tracking technique was adop 42 // the general ray tracking technique was adopted. The direct method in
43 // G4VSolid could be now (2009) used instead. 43 // G4VSolid could be now (2009) used instead.
44 44
45 // Author: L. Desorgher, SpaceIT GmbH - 01.06. 45 // Author: L. Desorgher, SpaceIT GmbH - 01.06.2006
46 // Contract: ESA contract 21435/08/NL/AT 46 // Contract: ESA contract 21435/08/NL/AT
47 // Customer: ESA/ESTEC 47 // Customer: ESA/ESTEC
48 // ------------------------------------------- 48 // --------------------------------------------------------------------
49 #ifndef G4AdjointPosOnPhysVolGenerator_hh 49 #ifndef G4AdjointPosOnPhysVolGenerator_hh
50 #define G4AdjointPosOnPhysVolGenerator_hh 1 50 #define G4AdjointPosOnPhysVolGenerator_hh 1
51 51
52 #include "G4VPhysicalVolume.hh" 52 #include "G4VPhysicalVolume.hh"
53 #include "G4AffineTransform.hh" 53 #include "G4AffineTransform.hh"
54 #include "G4ThreeVector.hh" 54 #include "G4ThreeVector.hh"
55 55
56 class G4VSolid; 56 class G4VSolid;
57 57
58 class G4AdjointPosOnPhysVolGenerator 58 class G4AdjointPosOnPhysVolGenerator
59 { 59 {
60 //--------- 60 //---------
61 public: 61 public:
62 //--------- 62 //---------
63 63
64 static G4AdjointPosOnPhysVolGenerator* Ge 64 static G4AdjointPosOnPhysVolGenerator* GetInstance();
65 65
66 G4VPhysicalVolume* DefinePhysicalVolume(co 66 G4VPhysicalVolume* DefinePhysicalVolume(const G4String& aName);
67 void DefinePhysicalVolume1(const G4String& 67 void DefinePhysicalVolume1(const G4String& aName);
68 G4double ComputeAreaOfExtSurface(); 68 G4double ComputeAreaOfExtSurface();
69 G4double ComputeAreaOfExtSurface(G4int NSt 69 G4double ComputeAreaOfExtSurface(G4int NStat);
70 G4double ComputeAreaOfExtSurface(G4double 70 G4double ComputeAreaOfExtSurface(G4double epsilon);
71 G4double ComputeAreaOfExtSurface(G4VSolid* 71 G4double ComputeAreaOfExtSurface(G4VSolid* aSolid);
72 G4double ComputeAreaOfExtSurface(G4VSolid* 72 G4double ComputeAreaOfExtSurface(G4VSolid* aSolid,G4int NStat);
73 G4double ComputeAreaOfExtSurface(G4VSolid* 73 G4double ComputeAreaOfExtSurface(G4VSolid* aSolid,G4double epsilon);
74 74
75 void GenerateAPositionOnTheExtSurfaceOfASo 75 void GenerateAPositionOnTheExtSurfaceOfASolid(G4VSolid* aSolid,
76 76 G4ThreeVector& p,
77 77 G4ThreeVector& direction);
78 void GenerateAPositionOnTheExtSurfaceOfThe 78 void GenerateAPositionOnTheExtSurfaceOfTheSolid(G4ThreeVector& p,
79 79 G4ThreeVector& direction);
80 void GenerateAPositionOnTheExtSurfaceOfThe 80 void GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(G4ThreeVector& p,
81 81 G4ThreeVector& direction);
82 void GenerateAPositionOnTheExtSurfaceOfThe 82 void GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(G4ThreeVector& p,
83 83 G4ThreeVector& direction,
84 84 G4double& costh_to_normal);
85 85
86 inline void SetSolid(G4VSolid* aSolid) 86 inline void SetSolid(G4VSolid* aSolid)
87 { theSolid=aSolid; } 87 { theSolid=aSolid; }
88 inline G4double GetAreaOfExtSurfaceOfThePh 88 inline G4double GetAreaOfExtSurfaceOfThePhysicalVolume()
89 { return AreaOfExtSurfaceOfThePhysicalVo 89 { return AreaOfExtSurfaceOfThePhysicalVolume; }
90 inline G4double GetCosThDirComparedToNorma 90 inline G4double GetCosThDirComparedToNormal()
91 { return CosThDirComparedToNormal; } 91 { return CosThDirComparedToNormal; }
92 92
93 //--------- 93 //---------
94 private: // private methods 94 private: // private methods
95 //--------- 95 //---------
96 G4AdjointPosOnPhysVolGenerator() = default 96 G4AdjointPosOnPhysVolGenerator() = default;
97 ~G4AdjointPosOnPhysVolGenerator() = default 97 ~G4AdjointPosOnPhysVolGenerator() = default;
98 G4double ComputeAreaOfExtSurfaceStartingFr 98 G4double ComputeAreaOfExtSurfaceStartingFromSphere(G4VSolid* aSolid,
99 99 G4int NStat);
100 G4double ComputeAreaOfExtSurfaceStartingFr 100 G4double ComputeAreaOfExtSurfaceStartingFromBox(G4VSolid* aSolid,
101 101 G4int NStat);
102 void GenerateAPositionOnASolidBoundary(G4V 102 void GenerateAPositionOnASolidBoundary(G4VSolid* aSolid,
103 G4T 103 G4ThreeVector& p,
104 G4T 104 G4ThreeVector& direction);
105 G4double GenerateAPositionOnASphereBoundar 105 G4double GenerateAPositionOnASphereBoundary(G4VSolid* aSolid,
106 106 G4ThreeVector& p,
107 107 G4ThreeVector& direction);
108 G4double GenerateAPositionOnABoxBoundary(G 108 G4double GenerateAPositionOnABoxBoundary(G4VSolid* aSolid,
109 G 109 G4ThreeVector& p,
110 G 110 G4ThreeVector& direction);
111 void ComputeTransformationFromPhysVolToWor 111 void ComputeTransformationFromPhysVolToWorld();
112 112
113 //--------- 113 //---------
114 private: // attributes 114 private: // attributes
115 //--------- 115 //---------
116 116
117 static G4ThreadLocal G4AdjointPosOnPhysVolG 117 static G4ThreadLocal G4AdjointPosOnPhysVolGenerator* theInstance;
118 G4VSolid* theSolid = nullptr; 118 G4VSolid* theSolid = nullptr;
119 G4VPhysicalVolume* thePhysicalVolume = null 119 G4VPhysicalVolume* thePhysicalVolume = nullptr;
120 120
121 G4bool UseSphere{true}; 121 G4bool UseSphere{true};
122 G4String ModelOfSurfaceSource{"OnSolid"}; 122 G4String ModelOfSurfaceSource{"OnSolid"};
123 G4AffineTransform theTransformationFromPhys 123 G4AffineTransform theTransformationFromPhysVolToWorld;
124 G4double AreaOfExtSurfaceOfThePhysicalVolum 124 G4double AreaOfExtSurfaceOfThePhysicalVolume{0.};
125 G4double CosThDirComparedToNormal{0.}; 125 G4double CosThDirComparedToNormal{0.};
126 }; 126 };
127 127
128 #endif 128 #endif
129 129