48 G4cout<<
"Regular transparent X-ray TR radiator EM process is called"<<
G4endl;
72 G4double result,
sum = 0.,
tmp, cof1, cof2, cofMin, cofPHC, theta2, theta2k ;
95 if (cofMin > kMin) kMin++;
114 for( k = kMin; k <=
kMax; k++ )
117 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
119 if( k == kMin && kMin ==
G4int(cofMin) )
121 sum += 0.5*std::sin(tmp)*std::sin(tmp)*
std::abs(k-cofMin)/result;
125 sum += std::sin(tmp)*std::sin(tmp)*
std::abs(k-cofMin)/result;
127 theta2k = std::sqrt(theta2*
std::abs(k-cofMin));
133 G4cout<<k<<
" "<<theta2k<<
" "<<std::sin(tmp)*std::sin(tmp)*
std::abs(k-cofMin)/result
137 result = 4.*( cof1 + cof2 )*( cof1 + cof2 )*sum/
energy;
198 G4double result, Qa, Qb,
Q, aZa, bZb, aMa, bMb,
D, sigma;
205 Qa = std::exp(-0.5*aMa);
206 Qb = std::exp(-0.5*bMb);
209 G4complex Ha( Qa*std::cos(aZa), -Qa*std::sin(aZa) );
210 G4complex Hb( Qb*std::cos(bZb), -Qb*std::sin(bZb) );
213 D = 1.0 /( (1 -
Q)*(1 - Q) +
214 4*Q*std::sin(0.5*(aZa + bZb))*std::sin(0.5*(aZa + bZb)) );
217 G4complex F2 = (1.0 - Ha)*(1.0 - Ha)*Hb*(1.0 - Hs)*(1.0 - Hs)
221 result = 2.0*std::real(R);