73 lowEnergyLimit(250*
eV),
74 highEnergyLimit(100*
GeV),
75 intrinsicLowEnergyLimit(10*
eV),
76 intrinsicHighEnergyLimit(100*
GeV)
81 G4Exception(
"G4LowEnergyRayleigh::G4LowEnergyRayleigh()",
83 "Energy limit outside intrinsic process validity range!");
89 G4String formFactorFile =
"rayl/re-ff-";
116 G4String crossSectionFile =
"rayl/re-cs-";
162 fcostheta = ( 1. + cosTheta*cosTheta)/2.;
165 G4double sinThetaHalf = std::sqrt((1. - cosTheta) / 2.);
166 x = sinThetaHalf / (wlPhoton/
cm);
172 sinTheta = std::sqrt(1. - cosTheta*cosTheta);
173 gReject = dataFormFactor * dataFormFactor;
175 }
while( gReject < randomFormFactor);
179 G4double dirX = sinTheta*std::cos(phi);
180 G4double dirY = sinTheta*std::sin(phi);
186 photonDirection1.
rotateUz(photonDirection0);
207 size_t materialIndex = couple->
GetIndex();