| 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 // This software was developed by Lawrence Liv
28 //
29 // Redistribution and use in source and binary
30 // modification, are permitted provided that t
31 //
32 // 1. Redistributions of source code must reta
33 // this list of conditions and the following
34 // 2. Redistributions in binary form must repr
35 // this list of conditions and the following
36 // and/or other materials provided with the
37 // 3. The name of the author may not be used t
38 // derived from this software without specif
39 //
40 // THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``A
41 // WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
42 // MERCHANTABILITY AND FITNESS FOR A PARTICULA
43 // EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DI
44 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGE
45 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES
46 // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AN
47 // WHETHER IN CONTRACT, STRICT LIABILITY, OR T
48 // OTHERWISE) ARISING IN ANY WAY OUT OF THE US
49 // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
50 //
51 // Copyright (c) 2006 The Regents of the Unive
52 // All rights reserved.
53 // UCRL-CODE-224807
54 //
55 //
56 //
57 // neutron_hp -- source file
58 // J.M. Verbeke, Jan-2007
59 // A prototype of the low energy neutron trans
60 //
61 #include "G4FissLib.hh"
62 #include "G4SystemOfUnits.hh"
63
64 G4FissLib::G4FissLib()
65 :xSec(0)
66 {
67 SetMinEnergy(0.0);
68 SetMaxEnergy(20.*MeV);
69 if(!G4FindDataDir("G4NEUTRONHPDATA")) {
70 G4cout << "Please setenv G4NEUTRONHPDATA
71 throw G4HadronicException(__FILE__, __LIN
72 }
73 dirName = G4FindDataDir("G4NEUTRONHPDATA");
74 G4String tString = "/Fission/";
75 dirName = dirName + tString;
76 numEle = G4Element::GetNumberOfElements();
77 theFission = new G4ParticleHPChannel[numEle]
78
79 for (std::size_t i=0; i<numEle; ++i)
80 {
81 if((*(G4Element::GetElementTable()))[i]->G
82 {
83 theFission[i].Init((*(G4Element::GetElem
84 theFission[i].Register(&theLibrary);
85 }
86 }
87 }
88
89 G4FissLib::~G4FissLib()
90 {
91 delete [] theFission;
92 }
93
94 G4HadFinalState*
95 G4FissLib::ApplyYourself(const G4HadProjectile
96 {
97 G4ParticleHPManager::GetInstance()->OpenReac
98
99 const G4Material* theMaterial = aTrack.GetMa
100 std::size_t n = theMaterial->GetNumberOfElem
101 std::size_t index = theMaterial->GetElement(
102
103 if (n != 1) {
104 xSec = new G4double[n];
105 G4double sum = 0;
106 G4int i;
107 std::size_t imat;
108 const G4double * NumAtomsPerVolume = theMa
109 G4double rWeight;
110 G4ParticleHPThermalBoost aThermalE;
111 for (i = 0; i < (G4int)n; ++i) {
112 imat = theMaterial->GetElement(i)->GetIn
113 rWeight = NumAtomsPerVolume[i];
114 xSec[i] = theFission[imat].GetXsec(aTher
115
116 theMaterial->GetTemperatur
117 xSec[i] *= rWeight;
118 sum+=xSec[i];
119 }
120
121 G4double random = G4UniformRand();
122 G4double running = 0;
123 for (i = 0; i < (G4int)n; ++i) {
124 running += xSec[i];
125 index = theMaterial->GetElement(i)->GetI
126 if(random<=running/sum) break;
127 }
128 delete [] xSec;
129 }
130
131 //return theFission[index].ApplyYourself(aTr
132 //Overwrite target parameters
133 G4HadFinalState* result = theFission[index]
134 aNucleus.SetParameters(G4ParticleHPManager:
135 const G4Element* target_element = (*G4Eleme
136 const G4Isotope* target_isotope=NULL;
137 G4int iele = (G4int)target_element->GetNumb
138 for ( G4int j = 0 ; j != iele ; ++j ) {
139 target_isotope=target_element->GetIsotop
140 if ( target_isotope->GetN() == G4Particl
141 }
142 aNucleus.SetIsotope( target_isotope );
143 G4ParticleHPManager::GetInstance()->CloseRe
144 return result;
145 }
146
147 const std::pair<G4double, G4double> G4FissLib:
148 {
149 // max energy non-conservation is mass of he
150 return std::pair<G4double, G4double>(10.0*pe
151 }
152