| 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 | // }
|
|---|
