51 G4cout<<
"Regular transparent X-ray TR radiator EM process is called"<<
G4endl;
75 G4double result, sum = 0.,
tmp, cof1, cof2, cofMin, cofPHC,aMa, bMb, sigma;
105 kMin =
G4int(cofMin);
106 if (cofMin > kMin) kMin++;
121 for( k = kMin; k <= kMax; k++ )
124 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
126 if( k == kMin && kMin ==
G4int(cofMin) )
128 sum += 0.5*sin(tmp)*sin(tmp)*std::abs(k-cofMin)/result;
132 sum += sin(tmp)*sin(tmp)*std::abs(k-cofMin)/result;
136 result = 4.*( cof1 + cof2 )*( cof1 + cof2 )*sum/
energy;
137 result *= ( 1. - exp(-
fPlateNumber*sigma) )/( 1. - exp(-sigma) );
181 G4double result, Qa, Qb, Q, aZa, bZb, aMa, bMb, D, sigma;
192 G4complex Ha( Qa*cos(aZa), -Qa*sin(aZa) );
193 G4complex Hb( Qb*cos(bZb), -Qb*sin(bZb) );
196 D = 1.0 /( (1 - Q)*(1 - Q) +
197 4*Q*sin(0.5*(aZa + bZb))*sin(0.5*(aZa + bZb)) );
198 G4complex F1 = (1.0 - Ha)*(1.0 - Hb)*(1.0 - Hs)
200 G4complex F2 = (1.0 - Ha)*(1.0 - Ha)*Hb*(1.0 - Hs)*(1.0 - Hs)
204 result = 2.0*real(R);