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 #include "G4InuclSpecialFunctions.hh"
00035 #include "G4CascadeInterpolator.hh"
00036
00037 void
00038 G4InuclSpecialFunctions::paraMaker(G4double Z,
00039 std::pair<std::vector<G4double>, std::vector<G4double> >& parms) {
00040 G4int verboseLevel(0);
00041
00042 if (verboseLevel > 3) {
00043 G4cout << " >>> G4InuclSpecialFunctions::paraMaker" << G4endl;
00044 }
00045
00046
00047
00048
00049 static const G4double Z1[5] = {10.0, 20.0, 30.0, 50.0, 70.0};
00050 static const G4double AP[5] = {0.42, 0.58, 0.68, 0.77, 0.80};
00051 static const G4double CP[5] = {0.50, 0.28, 0.20, 0.15, 0.10};
00052 static const G4double AA[5] = {0.68, 0.82, 0.91, 0.97, 0.98};
00053 static const G4double CA[5] = {0.10, 0.10, 0.10, 0.08, 0.06};
00054
00055
00056 std::vector<G4double>& AK = parms.first;
00057 AK.resize(6,0.);
00058
00059 std::vector<G4double>& CPA = parms.second;
00060 CPA.resize(6,0.);
00061
00062 AK[0] = 0.0;
00063 CPA[0] = 0.0;
00064
00065 static G4CascadeInterpolator<5> interp(Z1, false);
00066 AK[1] = interp.interpolate(Z, AP);
00067 AK[5] = interp.interpolate(Z, AA);
00068 CPA[1] = interp.interpolate(Z, CP);
00069 CPA[5] = interp.interpolate(Z, CA);
00070
00071 AK[2] = AK[1] + 0.06;
00072 AK[3] = AK[1] + 0.12;
00073 AK[4] = AK[5] - 0.06;
00074
00075 CPA[2] = CPA[1] * 0.5;
00076 CPA[3] = CPA[1] / 3.0;
00077 CPA[4] = 4.0 * CPA[5] / 3.0;
00078
00079 return;
00080 }
00081
00082 void
00083 G4InuclSpecialFunctions::paraMakerTruncated(G4double Z,
00084 std::pair<G4double,G4double>& parms) {
00085 G4int verboseLevel(0);
00086
00087 if (verboseLevel > 3) {
00088 G4cout << " >>> G4InuclSpecialFunctions::paraMakerTruncated" << G4endl;
00089 }
00090
00091
00092 static const G4double Z1[5] = {10.0, 20.0, 30.0, 50.0, 70.0};
00093 static const G4double AP[5] = {0.42, 0.58, 0.68, 0.77, 0.8};
00094 static const G4double CP[5] = {0.5, 0.28, 0.2, 0.15, 0.1};
00095
00096
00097 G4double& AK2=parms.first;
00098 G4double& CP2=parms.second;
00099
00100 static G4CascadeInterpolator<5> interp(Z1, false);
00101 AK2 = interp.interpolate(Z, AP);
00102 CP2 = interp.interpolate(Z, CP);
00103
00104 return;
00105 }