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
00031
00032
00033
00034
00035 #include "G4FissionParameters.hh"
00036 #include "G4SystemOfUnits.hh"
00037 #include "G4HadronicException.hh"
00038
00039
00040 const G4double G4FissionParameters::A1 = 134.0;
00041 const G4double G4FissionParameters::A2 = 141.0;
00042
00043 G4FissionParameters::G4FissionParameters(G4int A, G4int Z, G4double ExEnergy,
00044 G4double FissionBarrier)
00045 {
00046 G4double U = ExEnergy;
00047
00048 As = A/2.0;
00049
00050 if (A <= 235) Sigma2 = 5.6;
00051 else Sigma2 = 5.6 + 0.096*(A-235);
00052
00053 Sigma1 = 0.5*Sigma2;
00054
00055 SigmaS = std::exp(0.00553*U/MeV + 2.1386);
00056
00057
00058
00059
00060
00061 SigmaS*=0.8;
00062
00063
00064 G4double FasymAsym = 2.0*std::exp(-((A2-As)*(A2-As))/(2.0*Sigma2*Sigma2)) +
00065 std::exp(-((A1-As)*(A1-As))/(2.0*Sigma1*Sigma1));
00066
00067 G4double FsymA1A2 = std::exp(-((As-(A1+A2)/2.0)*(As-(A1+A2)/2.0))
00068 /(2.0*SigmaS*SigmaS));
00069
00070
00071 G4double wa = 0.0;
00072 w = 0.0;
00073 if (Z >= 90) {
00074 if (U <= 16.25) wa = std::exp(0.5385*U/MeV-9.9564);
00075 else wa = std::exp(0.09197*U/MeV-2.7003);
00076 } else if (Z == 89) {
00077 wa = std::exp(0.09197*U-1.0808);
00078 } else if (Z >= 82) {
00079 G4double X = FissionBarrier - 7.5*MeV;
00080 if (X < 0.0) X = 0.0;
00081 wa = std::exp(0.09197*(U-X)/MeV-1.0808);
00082 } else {
00083 w = 1001.0;
00084 }
00085
00086 if (w == 0.0) {
00087 G4double w1 = std::max(1.03*wa - FasymAsym, 0.0001);
00088 G4double w2 = std::max(1.0 - FsymA1A2*wa, 0.0001);
00089
00090 w = w1/w2;
00091
00092 if (82 <= Z && Z < 89 && A < 227) w *= std::exp(0.3*(227-A));
00093 }
00094
00095 }
00096
00097 G4FissionParameters::~G4FissionParameters()
00098 {}
00099