66 G4cout <<
"Screened Rutherford Elastic model is constructed "
96 G4cout <<
"Calling G4DNAScreenedRutherfordElasticModel::Initialise()"
103 G4Exception (
"*** WARNING: the G4DNAScreenedRutherfordElasticModel is not "
104 "intented to be used with another particle than the electron",
111 G4Exception(
"*** WARNING: the G4DNAScreenedRutherfordElasticModel class is "
112 "not validated below 9 eV",
118 G4Exception(
"*** WARNING: the G4DNAScreenedRutherfordElasticModel class is "
119 "not validated above 1 MeV",
127 G4cout <<
"Screened Rutherford elastic model is initialized " <<
G4endl
204 G4cout <<
"Calling CrossSectionPerVolume() of "
205 "G4DNAScreenedRutherfordElasticModel"
220 sigma =
pi * crossSection / (n * (n + 1.));
226 G4cout <<
"__________________________________" <<
G4endl;
227 G4cout <<
"=== G4DNAScreenedRutherfordElasticModel - XS INFO START"
229 G4cout <<
"=== Kinetic energy(eV)=" << ekin/
eV
230 <<
" particle : " << particleDefinition->GetParticleName()
232 G4cout <<
"=== Cross section per water molecule (cm^2)=" << sigma/
cm/
cm
234 G4cout <<
"=== Cross section per water molecule (cm^-1)="
235 << sigma*waterDensity/(1./
cm) << G4endl;
236 G4cout <<
"=== G4DNAScreenedRutherfordElasticModel - XS INFO END"
241 return sigma*waterDensity;
287 * std::pow(z, 2. / 3.);
294 if (denominator > 0.) value = numerator / denominator;
311 G4cout <<
"Calling SampleSecondaries() of "
312 "G4DNAScreenedRutherfordElasticModel"
324 {
G4cout <<
"---> Using Brenner & Zaider model" <<
G4endl;}
333 {
G4cout <<
"---> Using Screened Rutherford model" <<
G4endl;}
345 G4double xDir = std::sqrt(1. - cosTheta*cosTheta);
347 xDir *= std::cos(phi);
348 yDir *= std::sin(phi);
350 G4ThreeVector zPrimeVers((xDir*xVers + yDir*yVers + cosTheta*zVers));
401 / (1. / (4. * gamma * gamma) + beta
402 / ((2. + 2. * delta) * (2. + 2. *
delta)));
413 leftDenominator = (1. + 2.*gamma - cosTheta);
414 rightDenominator = (1. + 2.*delta + cosTheta);
415 if ( (leftDenominator * rightDenominator) != 0. )
417 fCosTheta = oneOverMax * (1./(leftDenominator*leftDenominator)
418 + beta/(rightDenominator*rightDenominator));
477 G4double S = 2.0 / (a1 * a2) + 2.0 * beta / (b1 * b2); S = 1.0 /
S;
480 G4double A = S * (b1 - beta * a2) + cp * a2 * b1;
481 G4double B = S * (b1 * b2 + beta * a1 * a2) - cp * a2 * b1 * c1;
482 G4double C = S * (b * b1 + beta * a * a2) - cp * a2 * b1 * c2;
485 return (-1.0 * B + std::sqrt(B * B - 4.0 * A * C)) / (2.0 * A);
533 std::vector<G4double>& vec)
540 size_t size = vec.size();
585 fCosTheta = (1 + 2.*n - cosTheta);
586 if (fCosTheta !=0.) fCosTheta = oneOverMax / (fCosTheta*fCosTheta);
618 G4double numerator = cp * (1.0 + 2.0 *
n) - n;
620 return numerator / denominator;