Geant4 Cross Reference |
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 // Hadronic Process: Nuclear De-excitations 28 // by V. Lara 29 // 30 // Modified: 31 // 21.03.2013 V.Ivanchenko redesigned and cleaned up 32 33 #include "G4CameronTruranHilfShellCorrections.hh" 34 #include <CLHEP/Units/SystemOfUnits.h> 35 36 // Data comes from: 37 // J.W. Truran, A.G.W. Cameron, and E. Hilf, 38 // Proc. Int. Conf. on the Properties of Nuclei Far From the Beta-Stability, 39 // Leysin, Switzerland, August 31 - September 4, 1970, Vol.1, p. 275 40 // S(Z) 41 G4double G4CameronTruranHilfShellCorrections::ShellZTable[] = 42 { // 93 from Z = 10 to Z = 102 43 2.349, 1.936, 1.596, 1.061, 0.341,-0.040, 0.565, 1.065, 1.536, 44 1.972, 1.855, 2.043, 1.931, 1.652, 1.347, 0.973, 0.579, 0.159, 45 -0.487,-0.192, 0.443, 0.932, 1.387, 1.810, 1.969, 2.067, 2.064, 46 1.825, 1.539, 1.251, 0.957, 1.128, 1.007, 0.603, 0.013,-0.635, 47 -1.258,-1.905,-2.562,-3.266,-4.099,-3.615,-3.171,-2.814,-2.337, 48 -1.778,-1.220,-0.694,-0.181, 0.323, 0.624, 0.841, 0.904, 0.906, 49 0.930, 0.919, 0.934, 0.941, 0.978, 0.982, 1.083, 1.201, 1.281, 50 1.189, 0.963, 0.781, 0.738, 0.696, 0.119,-0.619,-1.265,-1.898, 51 -2.431,-1.326,-0.268, 0.737, 1.451, 2.138, 2.307, 2.221, 2.041, 52 1.827, 1.239, 0.747, 0.214,-0.263,-0.778,-1.272,-1.800,-2.302, 53 -2.846,-3.499,-3.042 54 }; 55 // S(N) 56 G4double G4CameronTruranHilfShellCorrections::ShellNTable[] = 57 { // 146 from N = 10 to N = 155 58 2.439, 1.829, 1.419, 0.746,-0.082,-0.832,-0.960,-1.006,-1.045, 59 -1.114,-0.900,-0.081, 0.334, 0.064,-0.639,-1.363,-2.138,-2.987, 60 -4.042,-4.001,-3.582,-3.120,-2.677,-2.259,-1.778,-1.315,-0.944, 61 -0.599,-0.285,-0.020, 0.121, 0.140, 0.149,-0.001,-0.230,-0.604, 62 -1.010,-1.570,-2.466,-3.489,-4.552,-4.214,-3.375,-2.526,-1.725, 63 -0.923,-0.164, 0.601, 1.316, 1.947, 2.482, 2.971, 3.398, 3.737, 64 3.979, 4.183, 4.374, 4.517, 4.605, 4.539, 4.375, 4.043, 3.672, 65 3.250, 2.776, 2.254, 1.715, 1.151, 0.463,-0.237,-1.031,-1.850, 66 -2.722,-1.663,-0.724, 0.035, 0.786, 1.587, 2.145, 2.669, 2.680, 67 2.488, 2.243, 1.969, 1.778, 1.663, 1.487, 1.325, 1.148, 0.962, 68 0.843, 0.727, 0.574, 0.436, 0.320, 0.264, 0.397, 0.507, 0.405, 69 0.346, 0.369, 0.397, 0.403, 0.379, 0.184,-0.226,-0.737,-1.305, 70 -1.950,-2.565,-3.126,-3.721,-4.393,-5.082,-5.921,-6.712,-6.853, 71 -5.592,-4.413,-3.333,-2.413,-1.582,-0.966,-0.421,-0.123, 0.228, 72 0.543, 0.874, 1.059, 1.181, 1.186, 1.029, 1.029, 1.153, 1.227, 73 1.330, 1.449, 1.596, 1.712, 1.851, 1.949, 2.044, 2.155, 2.307, 74 2.621, 3.096 75 }; 76 77 G4CameronTruranHilfShellCorrections::G4CameronTruranHilfShellCorrections() 78 { 79 for(size_t i=0; i<ZTableSize; ++i) { ShellZTable[i] *= CLHEP::MeV; } 80 for(size_t i=0; i<NTableSize; ++i) { ShellNTable[i] *= CLHEP::MeV; } 81 } 82 83