| 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 // G4ILawTruncatedExp
27 // -------------------------------------------
28
29 #include "G4ILawTruncatedExp.hh"
30 #include "Randomize.hh"
31 #include "G4Track.hh"
32
33
34 G4ILawTruncatedExp::G4ILawTruncatedExp(const G
35 : G4VBiasingInteractionLaw(name)
36 {}
37
38 G4ILawTruncatedExp::~G4ILawTruncatedExp()
39 {}
40
41 void G4ILawTruncatedExp::SetForceCrossSection(
42 {
43 if (crossSection < 0.0)
44 {
45 G4Exception("G4ILawTruncatedExp::SetForceC
46 "BIAS.GEN.09", JustWarning,
47 "Cross-section value passed is
48 fIsSingular = true;
49 crossSection = 0.0;
50 }
51 fIsSingular = false;
52 fCrossSectionDefined = true;
53 fCrossSection = crossSection;
54 }
55
56 G4double G4ILawTruncatedExp::
57 ComputeEffectiveCrossSectionAt(G4double distan
58 {
59 if ( !fCrossSectionDefined )
60 {
61 G4Exception("G4ILawTruncatedExp::ComputeEf
62 "BIAS.GEN.10", JustWarning,
63 "Cross-section value requested
64 // -- zero cross-section, returns the limi
65 return 1.0 / ( fMaximumDistance - distance
66 }
67
68 G4double denum = 1.0 - std::exp(-fCrossSecti
69 return fCrossSection/denum;
70 }
71
72 G4double G4ILawTruncatedExp::
73 ComputeNonInteractionProbabilityAt(G4double di
74 {
75 if (!fCrossSectionDefined)
76 {
77 G4Exception("G4ILawTruncatedExp::ComputeNo
78 "BIAS.GEN.11", JustWarning,
79 "Non interaction probability v
80 // -- return limit case of null cross-sect
81 return 1.0 - distance / fMaximumDistance;
82 }
83 G4double num = 1.0 - std::exp( -fCrossSect
84 G4double denum = 1.0 - std::exp( -fCrossSect
85 return 1.0 - num/denum;
86 }
87
88 G4double G4ILawTruncatedExp::SampleInteraction
89 {
90 if ( !fCrossSectionDefined )
91 {
92 G4Exception("G4ILawTruncatedExp::Sample(..
93 "BIAS.GEN.12", JustWarning,
94 "Trying to sample while cross-
95 fInteractionDistance = G4UniformRand() * f
96 return fInteractionDistance;
97 }
98 fInteractionDistance = -std::log(1.0-G4Unifo
99 return fInteractionDistance;
100 }
101
102 G4double G4ILawTruncatedExp::
103 UpdateInteractionLengthForStep(G4double truePa
104 {
105 fInteractionDistance -= truePathLength;
106 fMaximumDistance -= truePathLength;
107
108 if ( fInteractionDistance < 0 )
109 {
110 G4ExceptionDescription ed;
111 ed << " Negative number of interaction len
112 << "' " << fInteractionDistance << ", s
113 G4Exception("G4ILawTruncatedExp::UpdateInt
114 "BIAS.GEN.13", JustWarning,
115 "Trying to sample while cross-
116 fInteractionDistance = 0.;
117 }
118
119 return fInteractionDistance;
120 }
121