| 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 // Gorad (Geant4 Open-source Radiation Analys
27 //
28 // Author : Makoto Asai (SLAC National Accele
29 //
30 // Development of Gorad is funded by NASA Joh
31 // under the contract NNJ15HK11B.
32 //
33 // *******************************************
34 //
35 // GRParallelWorldBiasingProcess.cc
36 // A process that takes care of the geometry
37 // This class assumes the existance of a par
38 // to the geometry importance biasing and bi
39 // associated to the region defined in that
40 //
41 // History
42 // September 8th, 2020 : first implementatio
43 //
44 // *******************************************
45
46 #include "G4ios.hh"
47 #include "GRParallelWorldBiasingProcess.hh"
48 #include "G4Step.hh"
49 #include "G4Track.hh"
50 #include "G4Navigator.hh"
51 #include "G4PathFinder.hh"
52 #include "G4VTouchable.hh"
53 #include "G4VPhysicalVolume.hh"
54 #include "G4ParticleChange.hh"
55 #include "G4ParticleChangeForNothing.hh"
56 #include "Randomize.hh"
57 #include "G4Region.hh"
58 #include "GRBiasingRegionInfo.hh"
59
60 GRParallelWorldBiasingProcess::
61 GRParallelWorldBiasingProcess(const G4String&
62 :G4ParallelWorldProcess(processName,theType)
63 {
64 particleChange = new G4ParticleChange();
65 emptyParticleChange = new G4ParticleChangeFo
66 }
67
68 GRParallelWorldBiasingProcess::~GRParallelWorl
69 {;}
70
71 G4VParticleChange* GRParallelWorldBiasingProce
72 const G4Track& track,
73 const G4Step& step)
74 {
75 fOldGhostTouchable = fGhostPostStepPoint->Ge
76 CopyStep(step);
77
78 if(fOnBoundary)
79 {
80 fNewGhostTouchable = fPathFinder->CreateTo
81 }
82 else
83 {
84 fNewGhostTouchable = fOldGhostTouchable;
85 }
86
87 fGhostPreStepPoint->SetTouchableHandle(fOldG
88 fGhostPostStepPoint->SetTouchableHandle(fNew
89
90 if(fNewGhostTouchable->GetVolume() == nullpt
91 {
92 // world volume boundary
93 emptyParticleChange->Initialize(track);
94 return emptyParticleChange;
95 }
96 if(fNewGhostTouchable->GetVolume() == fOldGh
97 {
98 // stayed in the same volume
99 emptyParticleChange->Initialize(track);
100 return emptyParticleChange;
101 }
102
103 G4int preStepCopyNo = fOldGhostTouchable->Ge
104 G4int postStepCopyNo = fNewGhostTouchable->G
105
106 GRBiasingRegionInfo* biasInfo = static_cast<
107 (fNewGhostTouchable->GetVolume()->GetLogic
108 auto probability = biasInfo->GetProbability(
109 if(probability < 1.)
110 {
111 // skip biasing to avoid over biasing
112 G4double ran = G4UniformRand();
113 if(ran > probability)
114 {
115 emptyParticleChange->Initialize(track);
116 return emptyParticleChange;
117 }
118 }
119
120 G4double trackWeight = track.GetWeight();
121 particleChange->Initialize(track);
122
123 auto biasFactor = biasInfo->GetBiasingFactor
124 if(biasFactor==1)
125 {
126 // not biasing
127 emptyParticleChange->Initialize(track);
128 return emptyParticleChange;
129 }
130
131 if(preStepCopyNo < postStepCopyNo)
132 {
133 // getting into more important volume, i.e
134 trackWeight /= biasFactor;
135 particleChange->ProposeParentWeight(trackW
136 for(G4int iClone=1;iClone<biasFactor;iClon
137 {
138 G4Track* clonedTrack = new G4Track(track
139 clonedTrack->SetWeight(trackWeight);
140 particleChange->AddSecondary(clonedTrack
141 }
142 particleChange->SetSecondaryWeightByProces
143 }
144 else if(preStepCopyNo > postStepCopyNo)
145 {
146 // getting into less important volume, i.e
147 G4double survivalRate = 1.0/biasFactor;
148 G4double ran = G4UniformRand();
149 if(ran > survivalRate)
150 {
151 // dead
152 particleChange->ProposeTrackStatus(fStop
153 }
154 else
155 {
156 // survive
157 trackWeight *= biasFactor;
158 particleChange->ProposeParentWeight(trac
159 }
160 }
161 else
162 {
163 // This should not happen!!
164 G4ExceptionDescription ed;
165 ed << "pre step point : "<<fOldGhostTouch
166 << "post step point : "<<fNewGhostTouch
167 G4Exception("GRParallelWorldBiasingProcess
168 }
169
170 return particleChange;
171 }
172
173