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: G4SmpSpNugDistData.cc 69860 2013-05-16 14:39:02Z gcosmo $ 00057 // 00058 00059 #include "G4fissionEvent.hh" 00060 00061 #define nSPfissg 30 00062 00063 G4int G4fissionEvent::G4SmpSpNugDistData(G4int isotope) { 00064 00065 /* 00066 Description 00067 Sample Number of Photons from spontaneous fission in 00068 (a) Cf-252 using the double Poisson model from Brunson; 00069 (b) Th-232, 00070 U-232, U-233, U-234, U-235, U-236, U-238*, 00071 Np-237, 00072 Pu-238*, Pu-239, Pu-240*, Pu-241, Pu-242*, 00073 Am-241, 00074 Cm-242*, Cm-244*, 00075 Bk-249 00076 using the negative binomial distribution based on the 00077 spontaneous fission neutron nubar from Ensslin's 00078 tabulated data or Holden and Zucker's tabulated data 00079 (for isotopes denoted with asterix *). 00080 */ 00081 00082 /* 00083 Input 00084 iso - isotope 00085 Output 00086 G4SmpSpNugDistData - sampled multiplicity 00087 -1 if there is no multiplicity data for that isotope 00088 */ 00089 00090 G4int i; 00091 G4double sum, nubar; 00092 G4double r; 00093 00094 static G4double Cf252spdist [nSPfissg] = { 00095 5.162699e-4,3.742057e-3,1.360482e-2,3.312786e-2,6.090540e-2, 00096 9.043537e-2,1.133984e-1,1.240985e-1,1.216759e-1,1.092255e-1, 00097 9.137106e-2,7.219960e-2,5.438060e-2,3.923091e-2,2.714690e-2, 00098 1.800781e-2,1.143520e-2,6.942099e-3,4.025720e-3,2.229510e-3, 00099 1.179602e-3,5.966936e-4,2.888766e-4,1.340137e-4,5.965291e-5, 00100 2.551191e-5,1.049692e-5,4.160575e-6,1.590596e-6,0.000000e+0 00101 }; 00102 00103 /* 00104 sample the spontaneous fission photon number distribution 00105 */ 00106 nubar=0.; 00107 if (isotope == 98252) { 00108 // Cf-252 using the G4double Poisson model from Brunson; 00109 r=fisslibrng(); 00110 00111 sum = 0.; 00112 for (i = 0; i < nSPfissg-1; i++) { 00113 sum = sum + Cf252spdist[i]; 00114 if (r <= sum || Cf252spdist[i+1] == 0.) return i; 00115 } 00116 } else if (isotope == 92238) { 00117 /* 00118 using the spontaneous fission nubar from 00119 Holden and Zucker's tabulated data 00120 */ 00121 nubar = 1.9900002; 00122 } else if (isotope == 94240) { 00123 nubar = 2.1540006; 00124 } else if (isotope == 94242) { 00125 nubar = 2.1489998; 00126 } else if (isotope == 96242) { 00127 nubar = 2.54; 00128 } else if (isotope == 96244) { 00129 nubar = 2.7200005; 00130 } else if (isotope == 94238) { 00131 nubar = 2.2100301; 00132 } 00133 00134 if (nubar != 0.) { 00135 return G4SmpNugDist(isotope, nubar); 00136 } else { 00137 /* 00138 using the spontaneous fission nubar from 00139 N. Ensslin, et.al., "Application Guide to Neutron 00140 Multiplicity Counting," LA-13422-M (November 1998) 00141 */ 00142 nubar = G4SmpSpNubarData(isotope); 00143 if (nubar != -1.) { 00144 return G4SmpNugDist(isotope, nubar); 00145 } else { 00146 // There is no nubar information for that isotope, return -1, 00147 // meaning no data available for that isotope 00148 return -1; 00149 } 00150 } 00151 }