| 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 * File: G4FPYNormalFragmentDist.cc 27 * File: G4FPYNormalFragmentDist.cc
28 * Author: B. Wendt (wendbryc@isu.edu) 28 * Author: B. Wendt (wendbryc@isu.edu)
29 * 29 *
30 * Created on July 26, 2011, 12:26 PM 30 * Created on July 26, 2011, 12:26 PM
31 */ 31 */
32 << 32
33 #include "G4FPYNormalFragmentDist.hh" << 33 #include "G4Ions.hh"
>> 34 #include "Randomize.hh"
>> 35 #include "globals.hh"
34 36
35 #include "G4FFGDebuggingMacros.hh" 37 #include "G4FFGDebuggingMacros.hh"
36 #include "G4FFGEnumerations.hh" 38 #include "G4FFGEnumerations.hh"
>> 39 #include "G4FPYNormalFragmentDist.hh"
37 #include "G4FissionProductYieldDist.hh" 40 #include "G4FissionProductYieldDist.hh"
38 #include "G4Ions.hh" <<
39 #include "Randomize.hh" <<
40 #include "globals.hh" <<
41 41
42 G4FPYNormalFragmentDist::G4FPYNormalFragmentDi << 42 G4FPYNormalFragmentDist::
43 << 43 G4FPYNormalFragmentDist( G4int WhichIsotope,
44 << 44 G4FFGEnumerations::MetaState WhichMetaState,
45 << 45 G4FFGEnumerations::FissionCause WhichCause,
46 << 46 G4FFGEnumerations::YieldType WhichYieldType,
47 : G4FissionProductYieldDist(WhichIsotope, Wh << 47 std::istringstream& dataFile)
>> 48 : G4FissionProductYieldDist( WhichIsotope,
>> 49 WhichMetaState,
>> 50 WhichCause,
>> 51 WhichYieldType,
>> 52 dataFile)
48 { 53 {
49 // Initialize the class << 54 // Initialize the class
50 Initialize(); << 55 Initialize();
51 } 56 }
52 57
53 G4FPYNormalFragmentDist::G4FPYNormalFragmentDi << 58 G4FPYNormalFragmentDist::
54 << 59 G4FPYNormalFragmentDist( G4int WhichIsotope,
55 << 60 G4FFGEnumerations::MetaState WhichMetaState,
56 << 61 G4FFGEnumerations::FissionCause WhichCause,
57 << 62 G4FFGEnumerations::YieldType WhichYieldType,
58 : G4FissionProductYieldDist(WhichIsotope, Wh << 63 G4int Verbosity,
59 dataFile) << 64 std::istringstream& dataFile)
>> 65 : G4FissionProductYieldDist( WhichIsotope,
>> 66 WhichMetaState,
>> 67 WhichCause,
>> 68 WhichYieldType,
>> 69 Verbosity,
>> 70 dataFile)
60 { 71 {
61 // Initialize the class << 72 // Initialize the class
62 Initialize(); << 73 Initialize();
63 } 74 }
64 75
65 void G4FPYNormalFragmentDist::Initialize(){G4F << 76 void G4FPYNormalFragmentDist::
>> 77 Initialize( void )
>> 78 {
>> 79 G4FFG_FUNCTIONENTER__
66 80
67 / << 81 // Nothing here yet
68 82
69 G << 83 G4FFG_FUNCTIONLEAVE__
>> 84 }
70 85
71 G4Ions* G4FPYNormalFragmentDist::GetFissionPro << 86 G4Ions* G4FPYNormalFragmentDist::
>> 87 GetFissionProduct( void )
72 { 88 {
73 G4FFG_FUNCTIONENTER__ << 89 G4FFG_FUNCTIONENTER__
74 90
75 G4Ions* Particle = nullptr; << 91 G4Ions* Particle;
76 92
77 // Generate a (0, 1] random number and retur << 93 // Generate a (0, 1] random number and return the respective particle.
78 // The ENDF data tables lists 72172 as the l << 94 // The ENDF data tables lists 72172 as the largest fission fragment produced
79 // for any fission event. The maximum alpha << 95 // for any fission event. The maximum alpha production is 10 and the
80 // smallest fissile isotope is 90227. This m << 96 // smallest fissile isotope is 90227. This means that if isotope 72172 were
81 // selected as the first daughter product, t << 97 // selected as the first daughter product, then at 10 alpha particles only
82 // 15 nucleons and -2 protons would remain f << 98 // 15 nucleons and -2 protons would remain for the second daughter product.
83 // Although the actual probability of this o << 99 // Although the actual probability of this occurring is very small, or 0 in
84 // this case, a check should still be made t << 100 // this case, a check should still be made to ensure that the second
85 // daughter product can be physically realiz << 101 // daughter product can be physically realized. This would prevent a
86 // situation such as this extreme example wh << 102 // situation such as this extreme example which results in a nucleus of 13
87 // neutrons and 2 anti-protons. << 103 // neutrons and 2 anti-protons.
88 // This quick sanity check may become even m << 104 // This quick sanity check may become even more valid if the ENDF data
89 // tables are expanded in the future and inc << 105 // tables are expanded in the future and include larger fission products.
90 << 106
91 G4int icounter = 0; << 107 G4int icounter=0;
92 G4int icounter_max = 1024; << 108 G4int icounter_max=1024;
93 do { << 109 do
94 icounter++; << 110 {
95 if (icounter > icounter_max) { << 111 icounter++;
96 G4cout << "Loop-counter exceeded the thr << 112 if ( icounter > icounter_max ) {
97 << __FILE__ << "." << G4endl; << 113 G4cout << "Loop-counter exceeded the threshold value at " << __LINE__ << "th line of " << __FILE__ << "." << G4endl;
98 break; << 114 break;
99 } << 115 }
100 Particle = FindParticle(RandomEngine_->G4S << 116 Particle = FindParticle(RandomEngine_->G4SampleUniform());
101 } while (Particle->GetAtomicMass() > Remaini << 117 } while(Particle->GetAtomicMass() > RemainingA_ + 1
102 || Particle->GetAtomicNumber() > Re << 118 || Particle->GetAtomicNumber() > RemainingZ_ + 1);
103 // Loop checking, 11.05.2015, T. Koi << 119 // Loop checking, 11.05.2015, T. Koi
104 120
105 G4FFG_FUNCTIONLEAVE__ << 121 G4FFG_FUNCTIONLEAVE__
106 return Particle; << 122 return Particle;
107 } 123 }
108 124
109 G4FPYNormalFragmentDist::~G4FPYNormalFragmentD << 125 G4FPYNormalFragmentDist::~G4FPYNormalFragmentDist( void )
110 { 126 {
111 G4FFG_FUNCTIONENTER__ << 127 G4FFG_FUNCTIONENTER__
112 <<
113 // Empty - all the data elements to be decon <<
114 // ~G4FissionProductYieldDist() <<
115 128
116 G4FFG_FUNCTIONLEAVE__ << 129 // Empty - all the data elements to be deconstructed are removed by
>> 130 // ~G4FissionProductYieldDist()
>> 131
>> 132 G4FFG_FUNCTIONLEAVE__
117 } 133 }
118 134