Geant4_10
G4EvaporationLevelDensityParameter.cc
Go to the documentation of this file.
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 // $Id: G4EvaporationLevelDensityParameter.cc 68724 2013-04-05 09:26:32Z gcosmo $
27 //
28 // Hadronic Process: Nuclear De-excitations
29 // by V. Lara (Oct 1998)
30 //
31 // Modified:
32 // 21.03.2013 V.Ivanchenko redesigned and comment out unused part
33 
35 
37 {}
38 
40 {}
41 
42 /*
43 
44 #include "G4ShellCorrections.hh"
45 #include "G4SystemOfUnits.hh"
46 #include "G4HadronicException.hh"
47 
48 // Those values are from table 3 in
49 // A.S. Iljinov et al. Nucl Phys A543 (1992) 517-557
50 // Table 3. alpha, beta and gamma for Cameron Shell corrections
51 // whithout collective effects. f-factor = 2.31.
52 
53 //JMQ 17-04-08 these are not used at present in G4Evaporation
54 const G4double G4EvaporationLevelDensityParameter::ConstEvapLevelDensityParameter = 0.125/MeV;
55 //const G4double G4EvaporationLevelDensityParameter::ConstEvapLevelDensityParameter= 0.0769231/MeV;
56 const G4double G4EvaporationLevelDensityParameter::alpha = 0.072/MeV;
57 const G4double G4EvaporationLevelDensityParameter::beta = 0.257/MeV;
58 const G4double G4EvaporationLevelDensityParameter::gamma = 0.059/MeV;
59 const G4double G4EvaporationLevelDensityParameter::Bs = 1.0;
60 
61 G4EvaporationLevelDensityParameter::G4EvaporationLevelDensityParameter()
62 {
63  // SPtr = G4ShellCorrections::GetInstance()->GetCameronTruranHilfShellCorrections();
64 }
65 
66 G4EvaporationLevelDensityParameter::~G4EvaporationLevelDensityParameter()
67 {}
68 
69 G4double
70 G4EvaporationLevelDensityParameter::LevelDensityParameter(G4int A, G4int, G4double) const
71 //JMQ (Apr .08) this is the method used in G4Evaporation
72 {
73 
74  //JMQ 25/04/08 a=A/10 according to original Gudima's prescription
75  G4double a=static_cast<G4double>(A)/10.;
76  return a;
77 //
78 
79 // G4int N = A - Z;
80 
81  // Asymptotic Level Density Parameter
82 // G4double AsymptoticLDP = (alpha*static_cast<G4double>(A) + beta*std::pow(static_cast<G4double>(A),2./3.)*Bs)/MeV;
83 
84  // Shape of the LDP U dependence
85 // G4double exponent = -gamma*U;
86 // G4double f = 1.;
87 // if (exponent > -300.) f -= std::exp(exponent);
88 
89  // Level Density Parameter
90 // G4double a = AsymptoticLDP*(1. + ShellCorrection(Z,N)*f/U);
91 // return a;
92 }
93 */