| 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 //
27 /// \file GB05BOptnSplitAndKillByCrossSection.
28 /// \brief Implementation of the GB05BOptnSpli
29
30 #include "GB05BOptnSplitAndKillByCrossSection.
31
32 #include "Randomize.hh"
33
34 //....oooOO0OOooo........oooOO0OOooo........oo
35
36 GB05BOptnSplitAndKillByCrossSection::GB05BOptn
37 : G4VBiasingOperation(name), fParticleChange
38 {}
39
40 //....oooOO0OOooo........oooOO0OOooo........oo
41
42 GB05BOptnSplitAndKillByCrossSection::~GB05BOpt
43
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 G4double GB05BOptnSplitAndKillByCrossSection::
47
48 {
49 *condition = NotForced;
50
51 // -- Sample the exponential law using the t
52 // -- to couterbalance for:
53 G4double proposedStepLength = -std::log(G4Un
54
55 return proposedStepLength;
56 }
57
58 //....oooOO0OOooo........oooOO0OOooo........oo
59
60 G4VParticleChange*
61 GB05BOptnSplitAndKillByCrossSection::GenerateB
62 {
63 // -- This method is called if we have limit
64 // -- We hence make the splitting or killing
65
66 // Get track weight:
67 G4double initialWeight = track->GetWeight();
68
69 // The "particle change" is the object to be
70 // the tracking the update of the primary st
71 // secondary tracks.
72 fParticleChange.Initialize(*track);
73
74 // -- Splitting and killing factors.
75 // -- They are taken the same, but the killi
76 G4double splittingFactor = 2.0;
77 G4double killingFactor = 2.0;
78
79 if (track->GetMomentumDirection().z() > 0) {
80 // -- We split if the track is moving forw
81
82 // Define the tracks weight:
83 G4double weightOfTrack = initialWeight / s
84
85 // Ask currect track weight to be changed
86 fParticleChange.ProposeParentWeight(weight
87 // Now make clones of this track (this is
88 // we will then have the primary and clone
89 // splitting by a factor 2:
90 G4Track* clone = new G4Track(*track);
91 clone->SetWeight(weightOfTrack);
92 fParticleChange.AddSecondary(clone);
93 // -- Below's call added for safety & illu
94 // -- modify the clone (ie : daughter) wei
95 // -- primary. Here call not mandatory and
96 fParticleChange.SetSecondaryWeightByProces
97 }
98 else {
99 // -- We apply Russian roulette if the tra
100
101 // Shoot a random number (in ]0,1[ segment
102 G4double random = G4UniformRand();
103 G4double killingProbability = 1.0 - 1.0 /
104 if (random < killingProbability) {
105 // We ask for the the track to be killed
106 fParticleChange.ProposeTrackStatus(fStop
107 }
108 else {
109 // In this case, the track survives. We
110 // to conserve weight among killed and s
111 fParticleChange.ProposeParentWeight(init
112 }
113 }
114
115 return &fParticleChange;
116 }
117
118 //....oooOO0OOooo........oooOO0OOooo........oo
119