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 58 #include <cmath> 59 #include "G4Pow.hh" 60 #include "G4fissionEvent.hh" 61 62 G4int G4fissionEvent::G4SmpNuDistDataPu239(G4d 63 64 /* 65 Description 66 Sample Number of Neutrons from fission in 67 Zucker and Holden's tabulated data for Pu- 68 */ 69 70 /* 71 Input 72 erg - incident neutron energy 73 Output 74 G4SmpNuDistDataPu239 - sampled multiplici 75 76 */ 77 78 G4double cpnu; 79 G4double pnu[8] = {0.0, 0.0, 0.0, 0.0, 0.0, 80 G4double eng; 81 G4double r; 82 83 /* 84 Check if energy is within the range of exper 85 */ 86 if (erg > 10) eng=10.; 87 else eng=erg; 88 89 r=fisslibrng(); 90 91 /* 92 Pu-239 nu distribution 93 */ 94 G4Pow* Pow=G4Pow::GetInstance(); 95 if (eng <= 5.0) pnu[0] = 0.0108826e0 - 0.002 96 - 6.5e-4*Pow->powN(en 97 - 7.93e-5*Pow->powN(e 98 if (eng > 5 && eng <= 10) pnu[0] = 0.078606e 99 + 1.42034e- 100 + 1.34512e- 101 if (r <= pnu[0]) return 0; 102 103 104 if (eng <= 5.0) pnu[1] = 0.0994916e0 - 0.019 105 - 0.00236583e0*Pow->p 106 - 4.14016667e-4*Pow-> 107 if (eng > 5 && eng <= 10) pnu[1] = 0.10052e0 108 + 3.78355e- 109 + 1.95458e- 110 cpnu=pnu[0]+pnu[1]; 111 if (r <= cpnu) return 1; 112 113 114 if (eng <= 5.0) pnu[2] = 0.2748898e0 - 0.015 115 - 0.00749681e0*Pow->p 116 - 3.13041667e-4*Pow-> 117 if (eng > 5 && eng <= 10) pnu[2] = 0.282487e 118 - 1.16895e- 119 - 6.41257e- 120 cpnu=cpnu+pnu[2]; 121 if (r <= cpnu) return 2; 122 123 if (eng <= 5.0) pnu[3] = 0.3269196e0 + 0.004 124 - 0.00189322e0*Pow->p 125 + 1.18466667e-4*Pow-> 126 if (eng > 5 && eng <= 10) pnu[3] = 0.329058e 127 - 3.06402e- 128 - 1.50875e- 129 cpnu=cpnu+pnu[3]; 130 if (r <= cpnu) return 3; 131 132 if (eng <= 5.0) pnu[4] = 0.2046061e0 + 0.026 133 + 0.0041514e0*Pow->po 134 + 5.0325e-4*Pow->powN 135 if (eng > 5 && eng <= 10) pnu[4] = 0.18992e0 136 - 7.06316e- 137 - 4.71791e- 138 cpnu=cpnu+pnu[4]; 139 if (r <= cpnu) return 4; 140 141 if (eng <= 5.0) pnu[5] = 0.0726834e0 + 0.001 142 + 0.007572e0*Pow->pow 143 + 2.3545e-4*Pow->powN 144 if (eng > 5 && eng <= 10) pnu[5] = 0.0779212 145 + 6.68583e- 146 + 4.88625e- 147 cpnu=cpnu+pnu[5]; 148 if (r <= cpnu) return 5; 149 150 if (eng <= 5.0) pnu[6] = 0.0097282e0 + 0.004 151 + 0.00115294e0*Pow->p 152 + 6.00083333e-5*Pow-> 153 if (eng > 5 && eng <= 10) pnu[6] = 7.85432e- 154 - 2.03705e- 155 - 4.24164e- 156 cpnu=cpnu+pnu[6]; 157 if (r <= cpnu) return 6; 158 159 if (eng <= 5.0) pnu[7] = 6.301e-4 + 1.106666 160 + 4.28016667e-4*Pow-> 161 - 4.31666667e-6*Pow-> 162 if (eng > 5 && eng <= 10) pnu[7] = 1.5323e-3 163 + 8.01017e- 164 + 4.38333e- 165 cpnu=cpnu+pnu[7]; 166 if (r <= cpnu) return 7; 167 else return 8; 168 } 169