46 s200keV(0.), s1keV(0.),
49 cmin(0.), clow(0.), chigh(0.)
55 G4cout <<
" G4RToEConvForGamma::G4RToEConvForGamma() ";
86 for (
size_t ibin=0; ibin<size_t(
TotBin); ibin++) {
88 for (
size_t iel=0; iel<size_t(NumEl); iel++) {
89 G4int IndEl = (*elementVector)[iel]->GetIndex();
90 SIGMA += atomicNumDensityVector[iel]*
91 (*((*aCrossSectionTable)[IndEl]))[ibin];
94 absorptionLengthVector->
PutValue(ibin, 5./SIGMA);
95 if (absorptionLengthMax < 5./SIGMA ) absorptionLengthMax = 5./SIGMA;
121 s200keV = (0.2651-0.1501*Zlog+0.02283*Zlogsquare)*Zsquare;
122 tmin = (0.552+218.5/Z+557.17/Zsquare)*
MeV;
123 smin = (0.01239+0.005585*Zlog-0.000923*Zlogsquare)*std::exp(1.5*Zlog);
125 tlow = 0.2*std::exp(-7.355/std::sqrt(Z))*
MeV;
127 s1keV = 300.*Zsquare;
130 chigh=(7.55e-5-0.0542e-5*
Z)*Zsquare*Z/std::log(t100MeV/
tmin);
135 if ( KineticEnergy<
tlow ) {
136 if(KineticEnergy<t1keV) xs =
slow*std::exp(
clow*std::log(
tlow/t1keV));
137 else xs =
slow*std::exp(
clow*std::log(
tlow/KineticEnergy));
139 }
else if ( KineticEnergy<t200keV ) {
141 * std::exp(0.042*
Z*std::log(t200keV/KineticEnergy)*std::log(t200keV/KineticEnergy));
143 }
else if( KineticEnergy<
tmin ){
145 * std::exp(
cmin*std::log(
tmin/KineticEnergy)*std::log(
tmin/KineticEnergy));