ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4SmpNuDistDataPu239.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4SmpNuDistDataPu239.cc
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 //
56 //
57 
58 #include <cmath>
59 #include "G4Pow.hh"
60 #include "G4fissionEvent.hh"
61 
62 G4int G4fissionEvent::G4SmpNuDistDataPu239(G4double erg) {
63 
64 /*
65  Description
66  Sample Number of Neutrons from fission in Pu-239 using
67  Zucker and Holden's tabulated data for Pu-239
68 */
69 
70 /*
71  Input
72  erg - incident neutron energy
73  Output
74  G4SmpNuDistDataPu239 - sampled multiplicity
75 
76 */
77 
78  G4double cpnu;
79  G4double pnu[8] = {0.0, 0.0, 0.0, 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 experimental values
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.00207694e0*eng
96  - 6.5e-4*Pow->powN(eng,2) + 4.023e-4*Pow->powN(eng,3)
97  - 7.93e-5*Pow->powN(eng,4) + 5.53666667e-6*Pow->powN(eng,5);
98  if (eng > 5 && eng <= 10) pnu[0] = 0.078606e0 - 5.17531e-2*eng
99  + 1.42034e-2*Pow->powN(eng,2) - 1.96292e-3*Pow->powN(eng,3)
100  + 1.34512e-4*Pow->powN(eng,4) - 3.63416e-6*Pow->powN(eng,5);
101  if (r <= pnu[0]) return 0;
102 
103 
104  if (eng <= 5.0) pnu[1] = 0.0994916e0 - 0.01979542e0*eng
105  - 0.00236583e0*Pow->powN(eng,2) + 0.0020581e0*Pow->powN(eng,3)
106  - 4.14016667e-4*Pow->powN(eng,4) + 2.85666667e-5*Pow->powN(eng,5);
107  if (eng > 5 && eng <= 10) pnu[1] = 0.10052e0 - 2.61361e-2*eng
108  + 3.78355e-3*Pow->powN(eng,2) - 3.70667e-4*Pow->powN(eng,3)
109  + 1.95458e-5*Pow->powN(eng,4) - 3.87499e-7*Pow->powN(eng,5);
110  cpnu=pnu[0]+pnu[1];
111  if (r <= cpnu) return 1;
112 
113 
114  if (eng <= 5.0) pnu[2] = 0.2748898e0 - 0.01565248e0*eng
115  - 0.00749681e0*Pow->powN(eng,2) + 0.00217121e0*Pow->powN(eng,3)
116  - 3.13041667e-4*Pow->powN(eng,4) + 1.88183333e-5*Pow->powN(eng,5);
117  if (eng > 5 && eng <= 10) pnu[2] = 0.282487e0 - 0.0261342e0*eng
118  - 1.16895e-3*Pow->powN(eng,2) + 1.9888e-4*Pow->powN(eng,3)
119  - 6.41257e-6*Pow->powN(eng,4) + 1.02502e-7*Pow->powN(eng,5);
120  cpnu=cpnu+pnu[2];
121  if (r <= cpnu) return 2;
122 
123  if (eng <= 5.0) pnu[3] = 0.3269196e0 + 0.00428312e0*eng
124  - 0.00189322e0*Pow->powN(eng,2) - 4.31925001e-4*Pow->powN(eng,3)
125  + 1.18466667e-4*Pow->powN(eng,4) - 9.04166668e-6*Pow->powN(eng,5);
126  if (eng > 5 && eng <= 10) pnu[3] = 0.329058e0 + 4.023e-3*eng
127  - 3.06402e-3*Pow->powN(eng,2) + 2.2628e-4*Pow->powN(eng,3)
128  - 1.50875e-5*Pow->powN(eng,4) + 4.39168e-7*Pow->powN(eng,5);
129  cpnu=cpnu+pnu[3];
130  if (r <= cpnu) return 3;
131 
132  if (eng <= 5.0) pnu[4] = 0.2046061e0 + 0.02633899e0*eng
133  + 0.0041514e0*Pow->powN(eng,2) - 0.00275542e0*Pow->powN(eng,3)
134  + 5.0325e-4*Pow->powN(eng,4) - 3.32158333e-5*Pow->powN(eng,5);
135  if (eng > 5 && eng <= 10) pnu[4] = 0.18992e0 + 4.55188e-2*eng
136  - 7.06316e-3*Pow->powN(eng,2) + 7.29916e-4*Pow->powN(eng,3)
137  - 4.71791e-5*Pow->powN(eng,4) + 1.185e-6*Pow->powN(eng,5);
138  cpnu=cpnu+pnu[4];
139  if (r <= cpnu) return 4;
140 
141  if (eng <= 5.0) pnu[5] = 0.0726834e0 + 0.00116043e0*eng
142  + 0.007572e0*Pow->powN(eng,2) - 0.00161972e0*Pow->powN(eng,3)
143  + 2.3545e-4*Pow->powN(eng,4) - 1.546e-5*Pow->powN(eng,5);
144  if (eng > 5 && eng <= 10) pnu[5] = 0.0779212e0 - 1.35849e-3*eng
145  + 6.68583e-3*Pow->powN(eng,2) - 7.98649e-4*Pow->powN(eng,3)
146  + 4.88625e-5*Pow->powN(eng,4) - 1.54167e-6*Pow->powN(eng,5);
147  cpnu=cpnu+pnu[5];
148  if (r <= cpnu) return 5;
149 
150  if (eng <= 5.0) pnu[6] = 0.0097282e0 + 0.00494589e0*eng
151  + 0.00115294e0*Pow->powN(eng,2) - 3.25191667e-4*Pow->powN(eng,3)
152  + 6.00083333e-5*Pow->powN(eng,4) - 3.745e-6*Pow->powN(eng,5);
153  if (eng > 5 && eng <= 10) pnu[6] = 7.85432e-3 + 7.33182e-3*eng
154  - 2.03705e-4*Pow->powN(eng,2) + 8.73787e-5*Pow->powN(eng,3)
155  - 4.24164e-6*Pow->powN(eng,4) + 2.37499e-7*Pow->powN(eng,5);
156  cpnu=cpnu+pnu[6];
157  if (r <= cpnu) return 6;
158 
159  if (eng <= 5.0) pnu[7] = 6.301e-4 + 1.10666667e-4*eng
160  + 4.28016667e-4*Pow->powN(eng,2) + 1.12041667e-5*Pow->powN(eng,3)
161  - 4.31666667e-6*Pow->powN(eng,4) + 3.29166667e-7*Pow->powN(eng,5);
162  if (eng > 5 && eng <= 10) pnu[7] = 1.5323e-3 - 7.91857e-4*eng
163  + 8.01017e-4*Pow->powN(eng,2) - 6.82833e-5*Pow->powN(eng,3)
164  + 4.38333e-6*Pow->powN(eng,4) - 6.0e-8*Pow->powN(eng,5);
165  cpnu=cpnu+pnu[7];
166  if (r <= cpnu) return 7;
167  else return 8;
168 }