| 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 // 26 //
27 // 27 //
28 // Hadronic Process: Nuclear De-excitations 28 // Hadronic Process: Nuclear De-excitations
29 // by V. Lara (Oct 1998) 29 // by V. Lara (Oct 1998)
30 30
31 #ifndef G4CompetitiveFission_h 31 #ifndef G4CompetitiveFission_h
32 #define G4CompetitiveFission_h 1 32 #define G4CompetitiveFission_h 1
33 33
34 #include "G4VEvaporationChannel.hh" 34 #include "G4VEvaporationChannel.hh"
35 #include "G4Fragment.hh" 35 #include "G4Fragment.hh"
36 #include "G4VEmissionProbability.hh" 36 #include "G4VEmissionProbability.hh"
37 #include "G4FissionParameters.hh" 37 #include "G4FissionParameters.hh"
38 #include <CLHEP/Units/SystemOfUnits.h> 38 #include <CLHEP/Units/SystemOfUnits.h>
39 #include "G4Exp.hh" 39 #include "G4Exp.hh"
40 40
41 class G4VFissionBarrier; 41 class G4VFissionBarrier;
42 class G4VEmissionProbability; 42 class G4VEmissionProbability;
43 class G4VLevelDensityParameter; 43 class G4VLevelDensityParameter;
44 class G4PairingCorrection; 44 class G4PairingCorrection;
45 45
46 class G4CompetitiveFission : public G4VEvapora 46 class G4CompetitiveFission : public G4VEvaporationChannel
47 { 47 {
48 public: 48 public:
49 49
50 G4CompetitiveFission(); << 50 explicit G4CompetitiveFission();
51 ~G4CompetitiveFission() override; 51 ~G4CompetitiveFission() override;
52 52
53 void Initialise() override; <<
54 <<
55 G4Fragment* EmittedFragment(G4Fragment* theN 53 G4Fragment* EmittedFragment(G4Fragment* theNucleus) override;
56 54
57 G4double GetEmissionProbability(G4Fragment* 55 G4double GetEmissionProbability(G4Fragment* theNucleus) override;
58 56
59 void SetFissionBarrier(G4VFissionBarrier * a 57 void SetFissionBarrier(G4VFissionBarrier * aBarrier);
60 58
61 void SetEmissionStrategy(G4VEmissionProbabil 59 void SetEmissionStrategy(G4VEmissionProbability * aFissionProb);
62 60
63 void SetLevelDensityParameter(G4VLevelDensit 61 void SetLevelDensityParameter(G4VLevelDensityParameter * aLevelDensity);
64 62
65 inline G4double GetFissionBarrier(void) cons 63 inline G4double GetFissionBarrier(void) const;
66 64
67 inline G4double GetLevelDensityParameter(voi 65 inline G4double GetLevelDensityParameter(void) const;
68 66
69 inline G4double GetMaximalKineticEnergy(void 67 inline G4double GetMaximalKineticEnergy(void) const;
70 68
71 G4CompetitiveFission(const G4CompetitiveFiss <<
72 const G4CompetitiveFission & operator=(const <<
73 G4bool operator==(const G4CompetitiveFission <<
74 G4bool operator!=(const G4CompetitiveFission <<
75 <<
76 private: 69 private:
77 70
78 // Sample AtomicNumber of Fission products 71 // Sample AtomicNumber of Fission products
79 G4int FissionAtomicNumber(G4int A); 72 G4int FissionAtomicNumber(G4int A);
80 73
81 G4double MassDistribution(G4double x, G4int 74 G4double MassDistribution(G4double x, G4int A);
82 75
83 // Sample Charge of fission products 76 // Sample Charge of fission products
84 G4int FissionCharge(G4int A, G4int Z, G4doub 77 G4int FissionCharge(G4int A, G4int Z, G4double Af);
85 78
86 // Sample Kinetic energy of fission products 79 // Sample Kinetic energy of fission products
87 G4double FissionKineticEnergy(G4int A, G4int 80 G4double FissionKineticEnergy(G4int A, G4int Z,
88 G4int Af1, G4int Zf1, 81 G4int Af1, G4int Zf1,
89 G4int Af2, G4int Zf2, 82 G4int Af2, G4int Zf2,
90 G4double U, G4double Tmax); 83 G4double U, G4double Tmax);
91 84
92 inline G4double Ratio(G4double A, G4double A 85 inline G4double Ratio(G4double A, G4double A11,
93 G4double B1, G4double 86 G4double B1, G4double A00) const;
94 87
95 inline G4double SymmetricRatio(G4int A, G4do 88 inline G4double SymmetricRatio(G4int A, G4double A11) const;
96 89
97 inline G4double AsymmetricRatio(G4int A, G4d 90 inline G4double AsymmetricRatio(G4int A, G4double A11) const;
98 91
99 inline G4double LocalExp(G4double x) const; 92 inline G4double LocalExp(G4double x) const;
100 93
>> 94 G4CompetitiveFission(const G4CompetitiveFission &right);
>> 95 const G4CompetitiveFission & operator=(const G4CompetitiveFission &right);
>> 96 G4bool operator==(const G4CompetitiveFission &right) const;
>> 97 G4bool operator!=(const G4CompetitiveFission &right) const;
>> 98
101 // Maximal Kinetic Energy that can be carrie 99 // Maximal Kinetic Energy that can be carried by fragment
102 G4double maxKineticEnergy{0.0}; << 100 G4double maxKineticEnergy;
103 G4double fissionBarrier{0.0}; << 101 G4double fissionBarrier;
104 G4double fissionProbability{0.0}; << 102 G4double fissionProbability;
105 G4double fFactor{1.0}; <<
106 103
107 // For Fission barrier 104 // For Fission barrier
108 G4VFissionBarrier* theFissionBarrierPtr; 105 G4VFissionBarrier* theFissionBarrierPtr;
109 106
110 // For Fission probability emission 107 // For Fission probability emission
111 G4VEmissionProbability* theFissionProbabilit 108 G4VEmissionProbability* theFissionProbabilityPtr;
112 109
113 // For Level Density calculation 110 // For Level Density calculation
114 G4VLevelDensityParameter* theLevelDensityPtr 111 G4VLevelDensityParameter* theLevelDensityPtr;
115 G4PairingCorrection* pairingCorrection; 112 G4PairingCorrection* pairingCorrection;
116 113
117 G4bool myOwnFissionProbability{true}; << 114 G4bool myOwnFissionProbability;
118 G4bool myOwnFissionBarrier{true}; << 115 G4bool myOwnFissionBarrier;
119 G4bool myOwnLevelDensity{true}; << 116 G4bool myOwnLevelDensity;
120 117
121 G4FissionParameters theParam; 118 G4FissionParameters theParam;
122 119
123 G4int theSecID; // Creator model ID for the 120 G4int theSecID; // Creator model ID for the secondaries created by this model
124 G4bool isInitialised{false}; <<
125 }; 121 };
126 122
127 inline G4double G4CompetitiveFission::GetFissi 123 inline G4double G4CompetitiveFission::GetFissionBarrier(void) const
128 { 124 {
129 return fissionBarrier; 125 return fissionBarrier;
130 } 126 }
131 127
132 inline G4double G4CompetitiveFission::GetMaxim 128 inline G4double G4CompetitiveFission::GetMaximalKineticEnergy(void) const
133 { 129 {
134 return maxKineticEnergy; 130 return maxKineticEnergy;
135 } 131 }
136 132
137 inline 133 inline
138 G4double G4CompetitiveFission::Ratio(G4double 134 G4double G4CompetitiveFission::Ratio(G4double A, G4double A11,
139 G4double B1, G4double A00) const 135 G4double B1, G4double A00) const
140 { 136 {
141 G4double res; 137 G4double res;
142 if (A11 >= A*0.5 && A11 <= (A00+10.0)) { 138 if (A11 >= A*0.5 && A11 <= (A00+10.0)) {
143 G4double x = (A11-A00)/A; 139 G4double x = (A11-A00)/A;
144 res = 1.0 - B1*x*x; 140 res = 1.0 - B1*x*x;
145 } else { 141 } else {
146 G4double x = 10.0/A; 142 G4double x = 10.0/A;
147 res = 1.0 - B1*x*x - 2.0*x*B1*(A11-A00-10. 143 res = 1.0 - B1*x*x - 2.0*x*B1*(A11-A00-10.0)/A;
148 } 144 }
149 return res; 145 return res;
150 } 146 }
151 147
152 inline 148 inline
153 G4double G4CompetitiveFission::AsymmetricRatio 149 G4double G4CompetitiveFission::AsymmetricRatio(G4int A, G4double A11) const
154 { 150 {
155 return Ratio(G4double(A),A11,23.5,134.0); 151 return Ratio(G4double(A),A11,23.5,134.0);
156 } 152 }
157 153
158 inline 154 inline
159 G4double G4CompetitiveFission::SymmetricRatio( 155 G4double G4CompetitiveFission::SymmetricRatio(G4int A, G4double A11) const
160 { 156 {
161 G4double A0 = G4double(A); 157 G4double A0 = G4double(A);
162 return Ratio(A0,A11,5.32,A0*0.5); 158 return Ratio(A0,A11,5.32,A0*0.5);
163 } 159 }
164 160
165 inline G4double G4CompetitiveFission::LocalExp 161 inline G4double G4CompetitiveFission::LocalExp(G4double x) const
166 { 162 {
167 return (std::abs(x) < 8.) ? G4Exp(-0.5*x*x) 163 return (std::abs(x) < 8.) ? G4Exp(-0.5*x*x) : 0.0;
168 } 164 }
169 165
170 #endif 166 #endif
171 167
172 168
173 169