64 G4cout <<
"Screened Rutherford Elastic model is constructed " <<
G4endl<<
"Energy range: "
89 G4cout <<
"Calling G4DNAScreenedRutherfordElasticModel::Initialise()"
96 G4cout <<
"G4DNAScreenedRutherfordElasticModel: low energy limit increased from " <<
103 G4cout <<
"G4DNAScreenedRutherfordElasticModel: high energy limit decreased from " <<
150 G4cout <<
"Screened Rutherford elastic model is initialized " << G4endl
176 G4cout <<
"Calling CrossSectionPerVolume() of G4DNAScreenedRutherfordElasticModel"
185 if(waterDensity!= 0.0)
197 sigma =
pi * crossSection / (n * (n + 1.));
202 G4cout <<
"__________________________________" <<
G4endl;
203 G4cout <<
"=== G4DNAScreenedRutherfordElasticModel - XS INFO START" <<
G4endl;
205 G4cout <<
"=== Cross section per water molecule (cm^2)=" << sigma/
cm/
cm <<
G4endl;
206 G4cout <<
"=== Cross section per water molecule (cm^-1)=" << sigma*waterDensity/(1./
cm) << G4endl;
208 G4cout <<
"=== G4DNAScreenedRutherfordElasticModel - XS INFO END" <<
G4endl;
213 return sigma*waterDensity;
230 G4double length = (e_squared * (k + electron_mass_c2))
231 / (4 *
pi * epsilon0 * k * (k + 2 * electron_mass_c2));
232 G4double cross = z * (z + 1) * length * length;
258 G4double numerator = (alpha_1 + beta_1 * std::log(k /
eV)) * constK
259 * std::pow(z, 2. / 3.);
261 k /= electron_mass_c2;
266 if (denominator > 0.) value = numerator / denominator;
283 G4cout <<
"Calling SampleSecondaries() of G4DNAScreenedRutherfordElasticModel"
320 G4double xDir = std::sqrt(1. - cosTheta*cosTheta);
322 xDir *= std::cos(phi);
323 yDir *= std::sin(phi);
325 G4ThreeVector zPrimeVers((xDir*xVers + yDir*yVers + cosTheta*zVers));
377 / (1. / (4. * gamma * gamma) + beta
378 / ((2. + 2. * delta) * (2. + 2. * delta)));
389 leftDenominator = (1. + 2.*gamma - cosTheta);
390 rightDenominator = (1. + 2.*delta + cosTheta);
391 if ( (leftDenominator * rightDenominator) != 0. )
393 fCosTheta = oneOverMax * (1./(leftDenominator*leftDenominator) + beta/(rightDenominator*rightDenominator));
418 for (
G4int i=0; i<iMax; i++)
420 cosTheta = -1 + i*2./(iMax-1);
421 leftDenominator = (1. + 2.*gamma - cosTheta);
422 rightDenominator = (1. + 2.*delta + cosTheta);
423 if ( (leftDenominator * rightDenominator) != 0. )
425 cumul = cumul + (1./(leftDenominator*leftDenominator) + beta/(rightDenominator*rightDenominator));
430 for (
G4int i=0; i<iMax; i++)
432 cosTheta = -1 + i*2./(iMax-1);
433 leftDenominator = (1. + 2.*gamma - cosTheta);
434 rightDenominator = (1. + 2.*delta + cosTheta);
435 if (cumul !=0 && (leftDenominator * rightDenominator) != 0.)
436 value = value + (1./(leftDenominator*leftDenominator) + beta/(rightDenominator*rightDenominator)) / cumul;
437 if (random < value)
break;
457 size_t size = vec.size();
501 fCosTheta = (1 + 2.*n - cosTheta);
502 if (fCosTheta !=0.) fCosTheta = oneOverMax / (fCosTheta*fCosTheta);
525 for (
G4int i=0; i<iMax; i++)
527 cosTheta = -1 + i*2./(iMax-1);
528 fCosTheta = (1 + 2.*n - cosTheta);
529 if (fCosTheta !=0.) cumul = cumul + 1./(fCosTheta*fCosTheta);
533 for (
G4int i=0; i<iMax; i++)
535 cosTheta = -1 + i*2./(iMax-1);
536 fCosTheta = (1 + 2.*n - cosTheta);
537 if (cumul !=0.) value = value + (1./(fCosTheta*fCosTheta)) / cumul;
538 if (random < value)
break;
G4double LowEnergyLimit() const
G4double CalculatePolynomial(G4double k, std::vector< G4double > &vec)
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4double HighEnergyLimit() const
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
G4double RutherfordCrossSection(G4double energy, G4double z)
std::vector< G4double > gamma035_10Coeff
std::vector< G4double > betaCoeff
std::vector< G4double > gamma100_200Coeff
std::vector< G4double > gamma10_100Coeff
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
const G4String & GetParticleName() const
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetHighEnergyLimit(G4double)
const std::vector< G4double > * fpWaterDensity
G4GLOB_DLL std::ostream G4cout
const std::vector< double > * GetNumMolPerVolTableFor(const G4Material *) const
const G4ThreeVector & GetMomentumDirection() const
G4double BrennerZaiderRandomizeCosTheta(G4double k)
std::vector< G4double > deltaCoeff
static G4DNAMolecularMaterial * Instance()
G4double ScreeningFactor(G4double energy, G4double z)
void SetProposedKineticEnergy(G4double proposedKinEnergy)
G4DNAScreenedRutherfordElasticModel(const G4ParticleDefinition *p=0, const G4String &nam="DNAScreenedRutherfordElasticModel")
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
void ProposeTrackStatus(G4TrackStatus status)
G4double intermediateEnergyLimit
virtual G4double CrossSectionPerVolume(const G4Material *material, const G4ParticleDefinition *p, G4double ekin, G4double emin, G4double emax)
void SetLowEnergyLimit(G4double)
G4ParticleChangeForGamma * fParticleChangeForGamma
virtual ~G4DNAScreenedRutherfordElasticModel()
G4ParticleChangeForGamma * GetParticleChangeForGamma()
G4double ScreenedRutherfordRandomizeCosTheta(G4double k, G4double z)