47 G4cout<<
"Regular X-ray TR radiator EM process is called"<<
G4endl ;
72 G4double result,
sum = 0.,
tmp, cof1, cof2, cofMin, cofPHC, theta2, theta2k;
78 sigma = 0.5*(aMa + bMb);
97 if (cofMin > kMin) kMin++;
116 for( k = kMin; k <=
kMax; k++ )
119 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
121 if( k == kMin && kMin ==
G4int(cofMin) )
123 sum += 0.5*std::sin(tmp)*std::sin(tmp)*
std::abs(k-cofMin)/result;
127 sum += std::sin(tmp)*std::sin(tmp)*
std::abs(k-cofMin)/result;
129 theta2k = std::sqrt(theta2*
std::abs(k-cofMin));
135 G4cout<<k<<
" "<<theta2k<<
" "<<std::sin(tmp)*std::sin(tmp)*
std::abs(k-cofMin)/result
139 result = 2*( cof1 + cof2 )*( cof1 + cof2 )*sum/
energy;
195 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
200 result = 2.0*std::real(R);