G4SmpNuDistDataU232_234_236_238_MC.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //
00027 // This software was developed by Lawrence Livermore National Laboratory.
00028 //
00029 // Redistribution and use in source and binary forms, with or without
00030 // modification, are permitted provided that the following conditions are met:
00031 //
00032 // 1. Redistributions of source code must retain the above copyright notice,
00033 //   this list of conditions and the following disclaimer.
00034 // 2. Redistributions in binary form must reproduce the above copyright notice,
00035 //   this list of conditions and the following disclaimer in the documentation
00036 //   and/or other materials provided with the distribution.
00037 // 3. The name of the author may not be used to endorse or promote products
00038 //   derived from this software without specific prior written permission.
00039 //
00040 // THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR IMPLIED
00041 // WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
00042 // MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO
00043 // EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
00044 // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00045 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
00046 // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
00047 // WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
00048 // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
00049 // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00050 //
00051 // Copyright (c) 2006 The Regents of the University of California.
00052 // All rights reserved.
00053 // UCRL-CODE-224807
00054 //
00055 //
00056 // $Id: G4SmpNuDistDataU232_234_236_238_MC.cc 69860 2013-05-16 14:39:02Z gcosmo $
00057 //
00058 
00059 #include <cmath>
00060 #include "G4fissionEvent.hh"
00061 
00062 G4int G4fissionEvent::G4SmpNuDistDataU232_234_236_238_MC(G4double nubar) {
00063 
00064 /*
00065   Description
00066     Sample Number of Neutrons from fission in U-232, U-234, U-236, 
00067     and U-238 using Zucker and Holden's tabulated data for U-238
00068     The 11 P(nu) distributions are given as a function of nubar, 
00069     the average number of neutrons from induced fission for the 
00070     11 different energies (0 to 10 MeV), based on the U-238 data 
00071     from Zucker and Holden.
00072 */
00073 
00074 /*
00075   Input
00076     nubar    - average number of neutrons per fission
00077   Output
00078     G4SmpNuDistDataU232_234_236_238_MC  - sampled multiplicity
00079     
00080 */
00081 
00082   static G4double U238nu [11] [9] = {
00083      {.0396484, .2529541, .2939544, .2644470, .1111758, .0312261, .0059347, .0005436, .0001158},
00084      {.0299076, .2043215, .2995886, .2914889, .1301480, .0363119, .0073638, .0006947, .0001751},
00085      {.0226651, .1624020, .2957263, .3119098, .1528786, .0434233, .0097473, .0009318, .0003159},
00086      {.0170253, .1272992, .2840540, .3260192, .1779579, .0526575, .0130997, .0013467, .0005405},
00087      {.0124932, .0984797, .2661875, .3344938, .2040116, .0640468, .0173837, .0020308, .0008730},
00088      {.0088167, .0751744, .2436570, .3379711, .2297901, .0775971, .0225619, .0030689, .0013626},
00089      {.0058736, .0565985, .2179252, .3368863, .2541575, .0933127, .0286200, .0045431, .0031316},
00090      {.0035997, .0420460, .1904095, .3314575, .2760413, .1112075, .0355683, .0065387, .0031316},
00091      {.0019495, .0309087, .1625055, .3217392, .2943792, .1313074, .0434347, .0091474, .0046284},
00092      {.0008767, .0226587, .1356058, .3076919, .3080816, .1536446, .0522549, .0124682, .0067176},
00093      {.0003271, .0168184, .1111114, .2892434, .3160166, .1782484, .0620617, .0166066, .0095665}
00094     };
00095   static G4double U238nubar [11] = {
00096       2.2753781,
00097       2.4305631,
00098       2.5857481,
00099       2.7409331,
00100       2.8961181,
00101       3.0513031,
00102       3.2064881,
00103       3.3616731,
00104       3.5168581,
00105       3.6720432,
00106       3.8272281
00107     };
00108   G4double fraction, r, cum;
00109   G4int engind, nu;
00110 
00111 /* 
00112   Check if nubar is within the range of experimental values
00113 */
00114   if(nubar >= U238nubar[0] && nubar <= U238nubar[10]) {
00115 /*
00116      Use Zucker and Holden Data
00117 */
00118      engind = 1;
00119      while (nubar > U238nubar[engind]){ engind++;}
00120      fraction = (nubar-U238nubar[engind-1])/(U238nubar[engind]-U238nubar[engind-1]);
00121      if(fisslibrng() > fraction) engind--;
00122 
00123      r = fisslibrng();
00124      nu = 0;
00125      cum = U238nu[engind][0];
00126      while (r > cum && nu < 8){ 
00127        nu++;
00128        cum += U238nu[engind][nu];
00129      }
00130      return nu;
00131   } else {
00132 /*
00133      Use Terrell's formula
00134 */
00135      return (G4int) G4SmpTerrell(nubar);
00136   }
00137 }

Generated on Mon May 27 17:49:51 2013 for Geant4 by  doxygen 1.4.7