| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26
27 #include "G4NuclearShellModelDensity.hh"
28 #include "G4PhysicalConstants.hh"
29 #include "G4SystemOfUnits.hh"
30 #include "G4Exp.hh"
31 #include "G4Log.hh"
32 #include "G4Pow.hh"
33
34 G4NuclearShellModelDensity::G4NuclearShellMode
35 : theA(anA)//, theZ(aZ)
36 {
37 const G4double r0sq=0.8133*fermi*fermi
38 theRsquare= r0sq * G4Pow::GetInstance()->Z23
39 G4double x = 1./(pi*theRsquare);
40 Setrho0(x*std::sqrt(x));
41 }
42
43 G4NuclearShellModelDensity::~G4NuclearShellMod
44
45 G4double G4NuclearShellModelDensity::GetRelati
46 {
47 return G4Exp(-1*aPosition.mag2()/theRsquare)
48 }
49
50 G4double G4NuclearShellModelDensity::GetRadius
51 {
52
53 return (maxRelativeDensity>0 && maxRelati
54 std::sqrt(theRsquare * G4Log(1/ma
55 }
56
57 G4double G4NuclearShellModelDensity::GetDeri
58 {
59 return -2* aPosition.mag() / theRsquare *
60 }
61