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 (Oct 1998) 29 // 30 // Modified: 31 // 21.03.2013 V.Ivanchenko redesigned and comment out unused part 32 33 #include "G4EvaporationLevelDensityParameter.hh" 34 #include "G4NuclearLevelData.hh" 35 36 G4EvaporationLevelDensityParameter::G4EvaporationLevelDensityParameter() 37 { 38 fNucData = G4NuclearLevelData::GetInstance(); 39 } 40 41 G4EvaporationLevelDensityParameter::~G4EvaporationLevelDensityParameter() 42 {} 43 44 G4double 45 G4EvaporationLevelDensityParameter::LevelDensityParameter(G4int A, G4int Z, G4double U) const 46 { 47 return fNucData->GetLevelDensity(Z, A, U); 48 } 49 50 /* 51 #include "G4ShellCorrections.hh" 52 #include "G4SystemOfUnits.hh" 53 54 // Those values are from table 3 in 55 // A.S. Iljinov et al. Nucl Phys A543 (1992) 517-557 56 // Table 3. alpha, beta and gamma for Cameron Shell corrections 57 // whithout collective effects. f-factor = 2.31. 58 59 //JMQ 17-04-08 these are not used at present in G4Evaporation 60 const G4double 61 G4EvaporationLevelDensityParameter::ConstEvapLevelDensityParameter = 0.125/MeV; 62 const G4double 63 G4EvaporationLevelDensityParameter::ConstEvapLevelDensityParameter= 0.0769231/MeV; 64 const G4double G4EvaporationLevelDensityParameter::alpha = 0.072/MeV; 65 const G4double G4EvaporationLevelDensityParameter::beta = 0.257/MeV; 66 const G4double G4EvaporationLevelDensityParameter::gamma = 0.059/MeV; 67 const G4double G4EvaporationLevelDensityParameter::Bs = 1.0; 68 69 // Asymptotic Level Density Parameter 70 //G4double AsymptoticLDP = (alpha*A + beta*g4pow-Z23(A)*Bs)/MeV; 71 72 // Shape of the LDP U dependence 73 G4double exponent = -gamma*U; 74 G4double f = 1.; 75 if (exponent > -300.) f -= G4Exp(exponent); 76 G4double a = AsymptoticLDP*(1. + ShellCorrection(Z,N)*f/U); 77 return a; 78 */ 79