Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 26 // 27 // This software was developed by Lawrence Liv 27 // This software was developed by Lawrence Livermore National Laboratory. 28 // 28 // 29 // Redistribution and use in source and binary 29 // Redistribution and use in source and binary forms, with or without 30 // modification, are permitted provided that t 30 // modification, are permitted provided that the following conditions are met: 31 // 31 // 32 // 1. Redistributions of source code must reta 32 // 1. Redistributions of source code must retain the above copyright notice, 33 // this list of conditions and the following 33 // this list of conditions and the following disclaimer. 34 // 2. Redistributions in binary form must repr 34 // 2. Redistributions in binary form must reproduce the above copyright notice, 35 // this list of conditions and the following 35 // this list of conditions and the following disclaimer in the documentation 36 // and/or other materials provided with the 36 // and/or other materials provided with the distribution. 37 // 3. The name of the author may not be used t 37 // 3. The name of the author may not be used to endorse or promote products 38 // derived from this software without specif 38 // derived from this software without specific prior written permission. 39 // 39 // 40 // THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``A 40 // THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR IMPLIED 41 // WARRANTIES, INCLUDING, BUT NOT LIMITED TO, 41 // WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF 42 // MERCHANTABILITY AND FITNESS FOR A PARTICULA 42 // MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO 43 // EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DI 43 // EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, 44 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGE 44 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 45 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES 45 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; 46 // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AN 46 // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, 47 // WHETHER IN CONTRACT, STRICT LIABILITY, OR T 47 // WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR 48 // OTHERWISE) ARISING IN ANY WAY OUT OF THE US 48 // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF 49 // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 49 // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 50 // 50 // 51 // Copyright (c) 2006 The Regents of the Unive 51 // Copyright (c) 2006 The Regents of the University of California. 52 // All rights reserved. 52 // All rights reserved. 53 // UCRL-CODE-224807 53 // UCRL-CODE-224807 54 // 54 // 55 // 55 // >> 56 // $Id: G4SmpNuDistDataPu239_241_MC.cc 68799 2013-04-05 13:29:46Z gcosmo $ 56 // 57 // 57 58 58 #include <cmath> 59 #include <cmath> 59 #include "G4fissionEvent.hh" 60 #include "G4fissionEvent.hh" 60 61 61 G4int G4fissionEvent::G4SmpNuDistDataPu239_241 62 G4int G4fissionEvent::G4SmpNuDistDataPu239_241_MC(G4double nubar) { 62 63 63 /* 64 /* 64 Description 65 Description 65 Sample Number of Neutrons from fission in 66 Sample Number of Neutrons from fission in Pu-239 and Pu-241 using 66 Zucker and Holden's tabulated data for Pu- 67 Zucker and Holden's tabulated data for Pu-239 67 The 11 P(nu) distributions are given as a 68 The 11 P(nu) distributions are given as a function of nubar, 68 the average number of neutrons from induce 69 the average number of neutrons from induced fission for the 69 11 different energies (0 to 10 MeV), based 70 11 different energies (0 to 10 MeV), based on the Pu-239 data 70 from Zucker and Holden. 71 from Zucker and Holden. 71 */ 72 */ 72 73 73 /* 74 /* 74 Input 75 Input 75 nubar - average number of neutrons per 76 nubar - average number of neutrons per fission 76 Output 77 Output 77 G4SmpNuDistDataPu239_241_MC - sampled mul 78 G4SmpNuDistDataPu239_241_MC - sampled multiplicity 78 79 79 */ 80 */ 80 81 81 static G4double Pu239nu [11] [9] = { 82 static G4double Pu239nu [11] [9] = { 82 {.0108826, .0994916, .2748898, .3269196, 83 {.0108826, .0994916, .2748898, .3269196, .2046061, .0726834, .0097282, .0006301, .0001685}, 83 {.0084842, .0790030, .2536175, .3289870, 84 {.0084842, .0790030, .2536175, .3289870, .2328111, .0800161, .0155581, .0011760, .0003469}, 84 {.0062555, .0611921, .2265608, .3260637, 85 {.0062555, .0611921, .2265608, .3260637, .2588354, .0956070, .0224705, .0025946, .0005205}, 85 {.0045860, .0477879, .1983002, .3184667, 86 {.0045860, .0477879, .1983002, .3184667, .2792811, .1158950, .0301128, .0048471, .0007233}, 86 {.0032908, .0374390, .1704196, .3071862, 87 {.0032908, .0374390, .1704196, .3071862, .2948565, .1392594, .0386738, .0078701, .0010046}, 87 {.0022750, .0291416, .1437645, .2928006, 88 {.0022750, .0291416, .1437645, .2928006, .3063902, .1641647, .0484343, .0116151, .0014149}, 88 {.0014893, .0222369, .1190439, .2756297, 89 {.0014893, .0222369, .1190439, .2756297, .3144908, .1892897, .0597353, .0160828, .0029917}, 89 {.0009061, .0163528, .0968110, .2558524, 90 {.0009061, .0163528, .0968110, .2558524, .3194566, .2134888, .0729739, .0213339, .0020017}, 90 {.0004647, .0113283, .0775201, .2335926, 91 {.0004647, .0113283, .0775201, .2335926, .3213289, .2356614, .0886183, .0274895, .0039531}, 91 {.0002800, .0071460, .0615577, .2089810, 92 {.0002800, .0071460, .0615577, .2089810, .3200121, .2545846, .1072344, .0347255, .0054786}, 92 {.0002064, .0038856, .0492548, .1822078, 93 {.0002064, .0038856, .0492548, .1822078, .3154159, .2687282, .1295143, .0432654, .0075217} 93 }; 94 }; 94 static G4double Pu239nubar [11] = { 95 static G4double Pu239nubar [11] = { 95 2.8760000, 96 2.8760000, 96 3.0088800, 97 3.0088800, 97 3.1628300, 98 3.1628300, 98 3.3167800, 99 3.3167800, 99 3.4707300, 100 3.4707300, 100 3.6246800, 101 3.6246800, 101 3.7786300, 102 3.7786300, 102 3.9325800, 103 3.9325800, 103 4.0865300, 104 4.0865300, 104 4.2404900, 105 4.2404900, 105 4.3944400 106 4.3944400 106 }; 107 }; 107 G4double fraction, r, cum; 108 G4double fraction, r, cum; 108 G4int engind, nu; 109 G4int engind, nu; 109 110 110 /* 111 /* 111 Check if nubar is within the range of experi 112 Check if nubar is within the range of experimental values 112 */ 113 */ 113 if(nubar >= Pu239nubar[0] && nubar <= Pu239n 114 if(nubar >= Pu239nubar[0] && nubar <= Pu239nubar[10]) { 114 /* 115 /* 115 Use Zucker and Holden Data 116 Use Zucker and Holden Data 116 */ 117 */ 117 engind = 1; 118 engind = 1; 118 while (nubar > Pu239nubar[engind]){ engin 119 while (nubar > Pu239nubar[engind]){ engind++;} 119 // Loop checking, 11.03.2015, T. Koi 120 // Loop checking, 11.03.2015, T. Koi 120 fraction = (nubar-Pu239nubar[engind-1])/( 121 fraction = (nubar-Pu239nubar[engind-1])/(Pu239nubar[engind]-Pu239nubar[engind-1]); 121 if(fisslibrng() > fraction) engind--; 122 if(fisslibrng() > fraction) engind--; 122 123 123 r = fisslibrng(); 124 r = fisslibrng(); 124 nu = 0; 125 nu = 0; 125 cum = Pu239nu[engind][0]; 126 cum = Pu239nu[engind][0]; 126 while (r > cum && nu < 8){ 127 while (r > cum && nu < 8){ 127 // Loop checking, 11.03.2015, T. Koi 128 // Loop checking, 11.03.2015, T. Koi 128 nu++; 129 nu++; 129 cum += Pu239nu[engind][nu]; 130 cum += Pu239nu[engind][nu]; 130 } 131 } 131 return nu; 132 return nu; 132 } else { 133 } else { 133 /* 134 /* 134 Use Terrell's formula 135 Use Terrell's formula 135 */ 136 */ 136 return (G4int) G4SmpTerrell(nubar); 137 return (G4int) G4SmpTerrell(nubar); 137 } 138 } 138 } 139 } 139 140