Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/fission/src/G4SmpNuDistDataPu239_241_MC.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/models/fission/src/G4SmpNuDistDataPu239_241_MC.cc (Version 11.3.0) and /processes/hadronic/models/fission/src/G4SmpNuDistDataPu239_241_MC.cc (Version 11.0.p4)


  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 //                                                 56 //
 57                                                    57 
 58 #include <cmath>                                   58 #include <cmath>
 59 #include "G4fissionEvent.hh"                       59 #include "G4fissionEvent.hh"
 60                                                    60 
 61 G4int G4fissionEvent::G4SmpNuDistDataPu239_241     61 G4int G4fissionEvent::G4SmpNuDistDataPu239_241_MC(G4double nubar) {
 62                                                    62 
 63 /*                                                 63 /*
 64   Description                                      64   Description
 65     Sample Number of Neutrons from fission in      65     Sample Number of Neutrons from fission in Pu-239 and Pu-241 using 
 66     Zucker and Holden's tabulated data for Pu-     66     Zucker and Holden's tabulated data for Pu-239
 67     The 11 P(nu) distributions are given as a      67     The 11 P(nu) distributions are given as a function of nubar, 
 68     the average number of neutrons from induce     68     the average number of neutrons from induced fission for the 
 69     11 different energies (0 to 10 MeV), based     69     11 different energies (0 to 10 MeV), based on the Pu-239 data 
 70     from Zucker and Holden.                        70     from Zucker and Holden.
 71 */                                                 71 */
 72                                                    72 
 73 /*                                                 73 /*
 74   Input                                            74   Input
 75     nubar    - average number of neutrons per      75     nubar    - average number of neutrons per fission
 76   Output                                           76   Output
 77     G4SmpNuDistDataPu239_241_MC  - sampled mul     77     G4SmpNuDistDataPu239_241_MC  - sampled multiplicity
 78                                                    78     
 79 */                                                 79 */
 80                                                    80 
 81   static G4double Pu239nu [11] [9] = {             81   static G4double Pu239nu [11] [9] = {
 82      {.0108826, .0994916, .2748898, .3269196,      82      {.0108826, .0994916, .2748898, .3269196, .2046061, .0726834, .0097282, .0006301, .0001685},
 83      {.0084842, .0790030, .2536175, .3289870,      83      {.0084842, .0790030, .2536175, .3289870, .2328111, .0800161, .0155581, .0011760, .0003469},
 84      {.0062555, .0611921, .2265608, .3260637,      84      {.0062555, .0611921, .2265608, .3260637, .2588354, .0956070, .0224705, .0025946, .0005205},
 85      {.0045860, .0477879, .1983002, .3184667,      85      {.0045860, .0477879, .1983002, .3184667, .2792811, .1158950, .0301128, .0048471, .0007233},
 86      {.0032908, .0374390, .1704196, .3071862,      86      {.0032908, .0374390, .1704196, .3071862, .2948565, .1392594, .0386738, .0078701, .0010046},
 87      {.0022750, .0291416, .1437645, .2928006,      87      {.0022750, .0291416, .1437645, .2928006, .3063902, .1641647, .0484343, .0116151, .0014149},
 88      {.0014893, .0222369, .1190439, .2756297,      88      {.0014893, .0222369, .1190439, .2756297, .3144908, .1892897, .0597353, .0160828, .0029917},
 89      {.0009061, .0163528, .0968110, .2558524,      89      {.0009061, .0163528, .0968110, .2558524, .3194566, .2134888, .0729739, .0213339, .0020017},
 90      {.0004647, .0113283, .0775201, .2335926,      90      {.0004647, .0113283, .0775201, .2335926, .3213289, .2356614, .0886183, .0274895, .0039531},
 91      {.0002800, .0071460, .0615577, .2089810,      91      {.0002800, .0071460, .0615577, .2089810, .3200121, .2545846, .1072344, .0347255, .0054786},
 92      {.0002064, .0038856, .0492548, .1822078,      92      {.0002064, .0038856, .0492548, .1822078, .3154159, .2687282, .1295143, .0432654, .0075217}
 93     };                                             93     };
 94   static G4double Pu239nubar [11] = {              94   static G4double Pu239nubar [11] = {
 95       2.8760000,                                   95       2.8760000,
 96       3.0088800,                                   96       3.0088800,
 97       3.1628300,                                   97       3.1628300,
 98       3.3167800,                                   98       3.3167800,
 99       3.4707300,                                   99       3.4707300,
100       3.6246800,                                  100       3.6246800,
101       3.7786300,                                  101       3.7786300,
102       3.9325800,                                  102       3.9325800,
103       4.0865300,                                  103       4.0865300,
104       4.2404900,                                  104       4.2404900,
105       4.3944400                                   105       4.3944400
106     };                                            106     };
107   G4double fraction, r, cum;                      107   G4double fraction, r, cum;
108   G4int engind, nu;                               108   G4int engind, nu;
109                                                   109 
110 /*                                                110 /* 
111   Check if nubar is within the range of experi    111   Check if nubar is within the range of experimental values
112 */                                                112 */
113   if(nubar >= Pu239nubar[0] && nubar <= Pu239n    113   if(nubar >= Pu239nubar[0] && nubar <= Pu239nubar[10]) {
114 /*                                                114 /*
115      Use Zucker and Holden Data                   115      Use Zucker and Holden Data
116 */                                                116 */
117      engind = 1;                                  117      engind = 1;
118      while (nubar > Pu239nubar[engind]){ engin    118      while (nubar > Pu239nubar[engind]){ engind++;}
119      // Loop checking, 11.03.2015, T. Koi         119      // Loop checking, 11.03.2015, T. Koi
120      fraction = (nubar-Pu239nubar[engind-1])/(    120      fraction = (nubar-Pu239nubar[engind-1])/(Pu239nubar[engind]-Pu239nubar[engind-1]);
121      if(fisslibrng() > fraction) engind--;        121      if(fisslibrng() > fraction) engind--;
122                                                   122 
123      r = fisslibrng();                            123      r = fisslibrng();
124      nu = 0;                                      124      nu = 0;
125      cum = Pu239nu[engind][0];                    125      cum = Pu239nu[engind][0];
126      while (r > cum && nu < 8){                   126      while (r > cum && nu < 8){ 
127      // Loop checking, 11.03.2015, T. Koi         127      // Loop checking, 11.03.2015, T. Koi
128        nu++;                                      128        nu++;
129        cum += Pu239nu[engind][nu];                129        cum += Pu239nu[engind][nu];
130      }                                            130      }
131      return nu;                                   131      return nu;
132   } else {                                        132   } else {
133 /*                                                133 /*
134      Use Terrell's formula                        134      Use Terrell's formula
135 */                                                135 */
136      return (G4int) G4SmpTerrell(nubar);          136      return (G4int) G4SmpTerrell(nubar);
137   }                                               137   }
138 }                                                 138 }
139                                                   139