| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // G4AdjointPrimaryGenerator class implementat
27 //
28 // Author: L. Desorgher, SpaceIT GmbH - Novemb
29 // Contract: ESA contract 21435/08/NL/AT
30 // Customer: ESA/ESTEC
31 // -------------------------------------------
32
33 #include "G4AdjointPrimaryGenerator.hh"
34 #include "G4PhysicalConstants.hh"
35 #include "G4Event.hh"
36 #include "G4SingleParticleSource.hh"
37 #include "G4ParticleDefinition.hh"
38 #include "G4AdjointPosOnPhysVolGenerator.hh"
39 #include "G4Navigator.hh"
40 #include "G4TransportationManager.hh"
41 #include "G4VPhysicalVolume.hh"
42 #include "G4Material.hh"
43 #include "Randomize.hh"
44
45 // -------------------------------------------
46 //
47 G4AdjointPrimaryGenerator::G4AdjointPrimaryGen
48 {
49 center_spherical_source = G4ThreeVector(0.,0
50 type_of_adjoint_source="Spherical";
51 theSingleParticleSource = new G4SingleParti
52
53 theSingleParticleSource->GetEneDist()->SetEn
54 theSingleParticleSource->GetEneDist()->SetAl
55 theSingleParticleSource->GetPosDist()->SetPo
56 theSingleParticleSource->GetAngDist()->SetAn
57
58 theG4AdjointPosOnPhysVolGenerator = G4Adjoin
59 }
60
61 // -------------------------------------------
62 //
63 G4AdjointPrimaryGenerator::~G4AdjointPrimaryGe
64 {
65 delete theSingleParticleSource;
66 }
67
68 // -------------------------------------------
69 //
70 void G4AdjointPrimaryGenerator::
71 GenerateAdjointPrimaryVertex(G4Event* anEvent,
72 G4double E1, G4do
73 {
74 if (type_of_adjoint_source == "ExternalSurfa
75 {
76 // Generate position and direction relativ
77 // of sensitive volume
78
79 G4double costh_to_normal=1.;
80 G4ThreeVector pos =G4ThreeVector(0.,0.,0.)
81 G4ThreeVector direction = G4ThreeVector(0.
82 theG4AdjointPosOnPhysVolGenerator
83 ->GenerateAPositionOnTheExtSurfaceOfTheP
84
85 if (costh_to_normal <1.e-4) { costh_to_no
86
87 // compute now the position along the ray
88 //
89 theSingleParticleSource->GetAngDist()
90 ->SetParticleMomentumDirection(-directio
91 theSingleParticleSource->GetPosDist()->Set
92 }
93
94 theSingleParticleSource->GetEneDist()->SetEm
95 theSingleParticleSource->GetEneDist()->SetEm
96
97 theSingleParticleSource->SetParticleDefiniti
98 theSingleParticleSource->GeneratePrimaryVert
99 }
100
101 // -------------------------------------------
102 //
103 void G4AdjointPrimaryGenerator::
104 GenerateFwdPrimaryVertex(G4Event* anEvent,G4Pa
105 G4double E1, G4double
106 {
107 if (type_of_adjoint_source == "ExternalSurfa
108 {
109 // Generate position and direction relativ
110 // of sensitive volume
111
112 G4double costh_to_normal=1.;
113 G4ThreeVector pos =G4ThreeVector(0.,0.,0.)
114 G4ThreeVector direction = G4ThreeVector(0.
115 theG4AdjointPosOnPhysVolGenerator
116 ->GenerateAPositionOnTheExtSurfaceOfTheP
117
118 if (costh_to_normal <1.e-4) { costh_to_no
119 theSingleParticleSource->GetAngDist()
120 ->SetParticleMomentumDirection(direction
121 theSingleParticleSource->GetPosDist()->Set
122 }
123
124 theSingleParticleSource->GetEneDist()->SetEm
125 theSingleParticleSource->GetEneDist()->SetEm
126
127 theSingleParticleSource->SetParticleDefiniti
128 theSingleParticleSource->GeneratePrimaryVert
129 }
130
131 // -------------------------------------------
132 //
133 void G4AdjointPrimaryGenerator::
134 SetSphericalAdjointPrimarySource(G4double radi
135 {
136 radius_spherical_source = radius;
137 center_spherical_source = center_pos;
138 type_of_adjoint_source = "Spherical";
139 theSingleParticleSource->GetPosDist()->SetPo
140 theSingleParticleSource->GetPosDist()->SetPo
141 theSingleParticleSource->GetPosDist()->SetCe
142 theSingleParticleSource->GetPosDist()->SetRa
143 theSingleParticleSource->GetAngDist()->SetAn
144 theSingleParticleSource->GetAngDist()->SetMa
145 theSingleParticleSource->GetAngDist()->SetMi
146 }
147
148 // -------------------------------------------
149 //
150 void G4AdjointPrimaryGenerator::
151 SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume
152 {
153 theG4AdjointPosOnPhysVolGenerator->DefinePhy
154 type_of_adjoint_source ="ExternalSurfaceOfAV
155 theSingleParticleSource->GetPosDist()->SetPo
156 theSingleParticleSource->GetAngDist()->SetAn
157 }
158
159 // -------------------------------------------
160 //
161 void G4AdjointPrimaryGenerator::
162 ComputeAccumulatedDepthVectorAlongBackRay(G4Th
163 G4Th
164 G4do
165 {
166 if (fLinearNavigator == nullptr)
167 {
168 fLinearNavigator = G4TransportationManager
169 ->GetNavigatorForTracking
170 }
171 G4ThreeVector position = glob_pos;
172 G4double safety=1.;
173 G4VPhysicalVolume* thePhysVolume =
174 fLinearNavigator->LocateGlobalPointAnd
175 G4double newStep = fLinearNavigator->Compute
176
177 delete theAccumulatedDepthVector;
178 theAccumulatedDepthVector = new G4PhysicsFre
179
180 G4double acc_depth=0.;
181 G4double acc_length=0.;
182 theAccumulatedDepthVector->InsertValues(acc_
183
184 while (newStep > 0. && thePhysVolume != null
185 {
186 acc_length+=newStep;
187 acc_depth+=newStep*thePhysVolume->GetLogic
188 ->GetMater
189 theAccumulatedDepthVector->InsertValues(ac
190 position=position+newStep*direction;
191 thePhysVolume = fLinearNavigator
192 ->LocateGlobalPointAndSetup(
193 newStep = fLinearNavigator->ComputeStep(po
194 }
195 }
196
197 // -------------------------------------------
198 //
199 G4double G4AdjointPrimaryGenerator::
200 SampleDistanceAlongBackRayAndComputeWeightCorr
201 {
202 G4double rand = G4UniformRand();
203 G4double distance = theAccumulatedDepthVecto
204 weight_corr=1.;
205 return distance;
206 }
207