00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030 #include "G4QMDParameters.hh"
00031 #include "G4PhysicalConstants.hh"
00032
00033 G4QMDParameters* G4QMDParameters::parameters = NULL;
00034
00035 G4QMDParameters::G4QMDParameters()
00036 {
00037
00038 wl = 2.0;
00039 hbc = 0.1973;
00040
00041
00042 cpw = 1.0 / 2.0 / wl;
00043
00044 cph = 2.0 * wl / (hbc*hbc);
00045
00046 cpc = 4.0;
00047
00048 epsx = -20.0 ;
00049
00050 rho0 = 0.168;
00051
00052
00053 G4double rpot = 1.0/3.0;
00054
00055 G4double ebinm = -16.0;
00056 G4double ebin = ebinm * 0.001;
00057
00058 G4double pfer = hbc * std::pow ( 3./2. *pi*pi * rho0 , 1./3. );
00059
00060 G4double rmass = 0.938;
00061
00062 G4double efer = pfer*pfer / 2. / rmass;
00063
00064 G4double t3 = 8. / 3. / rpot / std::pow( rho0 , ( 1.+rpot ) ) * ( efer / 5. - ebin );
00065
00066 G4double t0 = -16./15. * efer / rho0 - ( 1.+rpot ) * t3 * std::pow( rho0 , rpot );
00067
00068
00069 G4double aaa = 3./4. * t0 * rho0;
00070 G4double bbb = 3./8. * t3 * ( 2.+rpot ) * std::pow( rho0 , ( 1.+rpot ) );
00071 G4double esymm = 25 * 0.001;
00072
00073 gamm = rpot + 1.0;
00074
00075
00076 c0 = aaa / ( rho0 * std::pow( 4 * pi * wl , 1.5 ) * 2.0 );
00077
00078 c3 = bbb / ( std::pow( rho0 , gamm ) * std::pow ( (4.0*pi*wl) , (1.5*gamm) ) * ( gamm+1.0) );
00079
00080 cs = esymm / ( rho0 * std::pow( (4.0*pi*wl) , 1.5 ) * 2.0 );
00081
00082 G4double ccoul = 0.001439767;
00083 cl = ccoul/2.0 * 1;
00084
00085
00086
00087
00088
00089 cdp = 1.0 / std::pow ( ( 4.0 * pi * wl ) , 1.5 );
00090 c0p = c0 * 2.0;
00091 c3p = c3 * ( gamm + 1.0 );
00092 csp = cs * 2.0;
00093 clp = cl * 2.0;
00094
00095 }
00096
00097
00098
00099 G4QMDParameters::~G4QMDParameters()
00100 {
00101 ;
00102 }
00103