1 | // |
---|
2 | // ******************************************************************** |
---|
3 | // * License and Disclaimer * |
---|
4 | // * * |
---|
5 | // * The Geant4 software is copyright of the Copyright Holders of * |
---|
6 | // * the Geant4 Collaboration. It is provided under the terms and * |
---|
7 | // * conditions of the Geant4 Software License, included in the file * |
---|
8 | // * LICENSE and available at http://cern.ch/geant4/license . These * |
---|
9 | // * include a list of copyright holders. * |
---|
10 | // * * |
---|
11 | // * Neither the authors of this software system, nor their employing * |
---|
12 | // * institutes,nor the agencies providing financial support for this * |
---|
13 | // * work make any representation or warranty, express or implied, * |
---|
14 | // * regarding this software system or assume any liability for its * |
---|
15 | // * use. Please see the license in the file LICENSE and URL above * |
---|
16 | // * for the full disclaimer and the limitation of liability. * |
---|
17 | // * * |
---|
18 | // * This code implementation is the result of the scientific and * |
---|
19 | // * technical work of the GEANT4 collaboration. * |
---|
20 | // * By using, copying, modifying or distributing the software (or * |
---|
21 | // * any work based on the software) you agree to acknowledge its * |
---|
22 | // * use in resulting scientific publications, and indicate your * |
---|
23 | // * acceptance of all terms of the Geant4 Software license. * |
---|
24 | // ******************************************************************** |
---|
25 | // |
---|
26 | // |
---|
27 | // This software was developed by Lawrence Livermore National Laboratory. |
---|
28 | // |
---|
29 | // Redistribution and use in source and binary forms, with or without |
---|
30 | // modification, are permitted provided that the following conditions are met: |
---|
31 | // |
---|
32 | // 1. Redistributions of source code must retain the above copyright notice, |
---|
33 | // this list of conditions and the following disclaimer. |
---|
34 | // 2. Redistributions in binary form must reproduce the above copyright notice, |
---|
35 | // this list of conditions and the following disclaimer in the documentation |
---|
36 | // and/or other materials provided with the distribution. |
---|
37 | // 3. The name of the author may not be used to endorse or promote products |
---|
38 | // derived from this software without specific prior written permission. |
---|
39 | // |
---|
40 | // THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR IMPLIED |
---|
41 | // WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF |
---|
42 | // MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO |
---|
43 | // EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, |
---|
44 | // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
---|
45 | // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; |
---|
46 | // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, |
---|
47 | // WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR |
---|
48 | // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF |
---|
49 | // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
---|
50 | // |
---|
51 | // Copyright (c) 2006 The Regents of the University of California. |
---|
52 | // All rights reserved. |
---|
53 | // UCRL-CODE-224807 |
---|
54 | // |
---|
55 | // $Id: G4LLNLFission.hh,v 1.2 2007/05/30 00:43:20 dennis Exp $ |
---|
56 | // |
---|
57 | |
---|
58 | // |
---|
59 | // This class is a copy of Fission.hh, made for use with Geant4. |
---|
60 | // |
---|
61 | |
---|
62 | extern void genspfissevt_(G4int *isotope, G4double *time); |
---|
63 | /* |
---|
64 | * This function is called to trigger a spontaneous fission. |
---|
65 | * Multiple neutrons and photons are generated and stored |
---|
66 | * in a stack along with their energies, directions and |
---|
67 | * emission times. |
---|
68 | * The arguments of this function are |
---|
69 | * isotope: 94239 for Pu-239 for instance |
---|
70 | * time: the time of the spontaneous fission |
---|
71 | */ |
---|
72 | |
---|
73 | extern void genfissevt_(G4int *isotope, G4double *time, G4double *nubar, G4double *eng); |
---|
74 | /* |
---|
75 | * This function is called to trigger a neutron-induced fission. |
---|
76 | * Multiple neutrons and photons are generated and stored |
---|
77 | * in a stack along with their energies, directions and |
---|
78 | * emission times. In addition to the arguments above, this |
---|
79 | * function needs |
---|
80 | * nubar: user-specified average number of neutrons emitted |
---|
81 | * per fission (e.g. as tabulated in the cross-section |
---|
82 | * libraries used by the particle transport code) |
---|
83 | * eng: energy of the neutron inducing fission |
---|
84 | */ |
---|
85 | |
---|
86 | extern G4int getnnu_(); |
---|
87 | /* |
---|
88 | * This function returns the number of neutrons emitted by the |
---|
89 | * fission, -1 if there is no neutron data for that isotope in |
---|
90 | * the fission library. |
---|
91 | */ |
---|
92 | |
---|
93 | extern G4int getpnu_(); |
---|
94 | /* |
---|
95 | * This function returns the number of photons emitted by the |
---|
96 | * fission, -1 if there is no photon data for that isotope in |
---|
97 | * the fission library. |
---|
98 | */ |
---|
99 | |
---|
100 | extern G4double getneng_(G4int *index); |
---|
101 | /* |
---|
102 | * Given the index of the emitted neutron, this function returns |
---|
103 | * its energy, -1 if index isout of range. |
---|
104 | */ |
---|
105 | |
---|
106 | extern G4double getnvel_(G4int *index); |
---|
107 | /* |
---|
108 | * Given the index of the emitted neutron, this function returns |
---|
109 | * the amplitude of its velocity, -1 if index is out of range. |
---|
110 | */ |
---|
111 | |
---|
112 | extern G4double getndircosu_(G4int *index); |
---|
113 | extern G4double getndircosv_(G4int *index); |
---|
114 | extern G4double getndircosw_(G4int *index); |
---|
115 | /* |
---|
116 | * Given the index of the emitted neutron, this function returns |
---|
117 | * the direction cosines of its velocity vector on the x, y and z |
---|
118 | * axes. |
---|
119 | */ |
---|
120 | |
---|
121 | extern G4double getpeng_(G4int *index); |
---|
122 | /* |
---|
123 | * Given the index of the emitted photon, this function returns |
---|
124 | * its energy, -1 if index is out of range. |
---|
125 | */ |
---|
126 | |
---|
127 | extern G4double getpvel_(G4int *index); |
---|
128 | /* |
---|
129 | * Given the index of the emitted photon, this function returns |
---|
130 | * the amplitude of its velocity, -1 if index is out of range. |
---|
131 | */ |
---|
132 | |
---|
133 | extern G4double getpdircosu_(G4int *index); |
---|
134 | extern G4double getpdircosv_(G4int *index); |
---|
135 | extern G4double getpdircosw_(G4int *index); |
---|
136 | /* |
---|
137 | * Given the index of the emitted photon, this function returns |
---|
138 | * the direction cosines of its velocity. |
---|
139 | */ |
---|
140 | |
---|
141 | extern G4double getnage_(G4int *index); |
---|
142 | /* |
---|
143 | * Given the index of the emitted neutron, this function returns |
---|
144 | * its age, -1 if index is out of range. |
---|
145 | * This age will be different from the time specified |
---|
146 | * in generateFissionEvent and generateSpontaneousFissionEvent |
---|
147 | * for non-prompt neutrons, i.e. delayed neutrons. |
---|
148 | */ |
---|
149 | |
---|
150 | extern G4double getpage_(G4int *index); |
---|
151 | /* |
---|
152 | * Given the index of the emitted photon, this function returns |
---|
153 | * its age, -1 of index is out of range. |
---|
154 | * This age will be different from the time specified |
---|
155 | * in generateFissionEvent and generateSpontaneousFissionEvent |
---|
156 | * for photons that are emitted by beta-decay of the fission |
---|
157 | * fragments. |
---|
158 | */ |
---|
159 | |
---|
160 | extern void setdelay_(G4int *delay); |
---|
161 | /* |
---|
162 | * This function is called to enable delayed neutrons and photons |
---|
163 | * Input |
---|
164 | * delay: |
---|
165 | * 0 (default) for strictly prompt neutrons and |
---|
166 | * photons |
---|
167 | * 1 (n/a) for prompt neutrons, prompt and delayed |
---|
168 | * photons |
---|
169 | * 2 (n/a) for prompt and delayed neutrons, prompt |
---|
170 | * photons |
---|
171 | * 3 (n/a) for prompt and delayed neutrons, prompt |
---|
172 | * and delayed photons |
---|
173 | */ |
---|
174 | |
---|
175 | extern void setcorrel_(G4int *correlation); |
---|
176 | /* |
---|
177 | * This function is called to set the type of neutron photon correlation |
---|
178 | * Input |
---|
179 | * correlation: |
---|
180 | * 0 (default) for no correlation between neutrons and |
---|
181 | * photons |
---|
182 | * 1 (n/a) for number correlation between neutrons and |
---|
183 | * photons |
---|
184 | * 2 (n/a) for number and energy correlation between |
---|
185 | * neutrons and photons |
---|
186 | */ |
---|
187 | |
---|
188 | extern void setnudist_(G4int *nudist); |
---|
189 | /* |
---|
190 | * This function is called to set the data to be sampled for the neutron |
---|
191 | * number distributions in induced fissions |
---|
192 | * Input |
---|
193 | * nudist: |
---|
194 | * 0 to use the fit to the Zucker and Holden tabulated |
---|
195 | * P(nu) distributions as a function of energy for |
---|
196 | * U235, U238 and Pu239. Terrell for other isotopes. |
---|
197 | * 1 to use fits to the Zucker and Holden tabulated |
---|
198 | * P(nu) distribution as a function of energy for |
---|
199 | * U238 and Pu239, and a fit to the Zucker and Holden |
---|
200 | * data as well as the Gwin, Spencer and Ingle data |
---|
201 | * (at thermal energies) as a function of energy for |
---|
202 | * U235. Terrell for other isotopes. |
---|
203 | * 2 (default) to use the fit to the Zucker and Holden |
---|
204 | * tabulated P(nu) distributions as a function of nubar. |
---|
205 | * The U238 fit is used for the U232, U234, U236 and |
---|
206 | * U238 isotopes, the U235 fit for U233 and U235, the |
---|
207 | * Pu239 fit for Pu239 and Pu241. Terrell for other |
---|
208 | * isotopes. |
---|
209 | */ |
---|
210 | |
---|
211 | |
---|
212 | extern void setcf252_(G4int *ndist, G4int *neng); |
---|
213 | /* |
---|
214 | * This function is called to set the data to be sampled for the |
---|
215 | * (a) Cf252 spontaneous fission number distribution, and |
---|
216 | * (b) Cf252 spontaneous fission neutron energy spectrum |
---|
217 | * Input |
---|
218 | * ndist: |
---|
219 | * 0 (default) to sample the number of neutrons from the |
---|
220 | * tabulated data measured by Spencer |
---|
221 | * 1 to sample the number of neutrons from Boldeman's data |
---|
222 | * neng: |
---|
223 | * 0 to sample the spontaneous fission neutron energy from |
---|
224 | * Mannhart corrected Maxwellian spectrum |
---|
225 | * 1 to sample the spontaneous fission neutron energy from |
---|
226 | * Madland-Nix theoretical spectrum |
---|
227 | * 2 to sample the spontaneous fission neutron energy from |
---|
228 | * the Froehner Watt spectrum |
---|
229 | */ |
---|
230 | |
---|
231 | extern void setrngf_(G4float (*funcptr) (void)); |
---|
232 | /* |
---|
233 | * This function sets the random number generator to the user-defined |
---|
234 | * one specified in the argument. If either setrngf_ or setrngd_ are |
---|
235 | * not specified, the default system call srand48 will be called. |
---|
236 | * Input |
---|
237 | * funcptr: |
---|
238 | * a random number generator function that returns a |
---|
239 | * variable of type G4float |
---|
240 | */ |
---|
241 | |
---|
242 | extern void setrngd_(G4double (*funcptr) (void)); |
---|
243 | /* |
---|
244 | * This function sets the random number generator to the user-defined |
---|
245 | * one specified in the argument. If either setrngf_ or setrngd_ are |
---|
246 | * not specified, the default system call srand48 will be called. |
---|
247 | * Input |
---|
248 | * funcptr: |
---|
249 | * a random number generator function that returns a |
---|
250 | * variable of type G4double |
---|
251 | */ |
---|
252 | // } |
---|