DNA::Penetration::Terrisol1990 Struct Reference

#include <G4DNAOneStepThermalizationModel.hh>

Static Public Member Functions
static void	GetPenetration (G4double energy, G4ThreeVector &displacement)
static double	GetRmean (double energy)
static double	Get3DStdDeviation (double energy)

Static Public Attributes
static const double	gEnergies_T1990 [11]
static const double	gStdDev_T1990 [11]

Detailed Description

Definition at line 101 of file G4DNAOneStepThermalizationModel.hh.

Member Function Documentation

Get3DStdDeviation()

double DNA::Penetration::Terrisol1990::Get3DStdDeviation ( double  energy )

static

Definition at line 185 of file G4DNAOneStepThermalizationModel.cc.

                                                   {
  G4double k_eV = energy/eV;
  if(k_eV < 0.2){
    // rare events:
    //  prevent H2O and secondary electron to be at the spot
    return 1e-3*CLHEP::nanometer;
  }
  else if(k_eV == 9.){
    return gStdDev_T1990[10];
  }
  else if(k_eV > 9.){
    G4ExceptionDescription description;
    description << "Terrisol1990 is not tabulated for energies greater than 9eV";
    G4Exception("Terrisol1990::Get3DStdDeviation",
                "INVALID_ARGUMENT",
                FatalErrorInArgument,
                description);
  }
 
  size_t lowBin, upBin;
 
  if(k_eV >= 1.){
    lowBin=std::floor(k_eV)+1;
    upBin=std::min(lowBin+1, size_t(10));
  }
  else{
    auto it=std::lower_bound(&gEnergies_T1990[0],
                             &gEnergies_T1990[2],
                             k_eV);
    lowBin = it-&gEnergies_T1990[0];
    upBin = lowBin+1;
  }
 
  double lowE = gEnergies_T1990[lowBin];
  double upE = gEnergies_T1990[upBin];
 
  double lowS = gStdDev_T1990[lowBin];
  double upS = gStdDev_T1990[upBin];
 
  double tanA = (lowS-upS)/(lowE-upE);
  double sigma3D = lowS + (k_eV-lowE)*tanA;
  return sigma3D;
}

Referenced by GetPenetration(), and GetRmean().

GetPenetration()

void DNA::Penetration::Terrisol1990::GetPenetration ( G4double  energy, G4ThreeVector &  displacement  )

static

Definition at line 238 of file G4DNAOneStepThermalizationModel.cc.

                                                              {
  double sigma3D = Get3DStdDeviation(energy);
 
  static constexpr double factor = 2.20496999539; // = 1./(3. - 8./CLHEP::pi);
 
  double sigma1D = std::sqrt(std::pow(sigma3D, 2.)*factor);
 
  displacement = G4ThreeVector(G4RandGauss::shoot(0, sigma1D),
                               G4RandGauss::shoot(0, sigma1D),
                               G4RandGauss::shoot(0, sigma1D));
}

Referenced by G4DNAWaterDissociationDisplacer::radialDistributionOfElectron().

GetRmean()

double DNA::Penetration::Terrisol1990::GetRmean ( double  energy )

static

Definition at line 229 of file G4DNAOneStepThermalizationModel.cc.

                                          {
  double sigma3D=Get3DStdDeviation(energy);
 
  static constexpr double s2r=1.595769121605731; // = pow(2,3./2.)/sqrt(CLHEP::pi)
 
  double r_mean=sigma3D*s2r;
  return r_mean;
}

Member Data Documentation

gEnergies_T1990

const double DNA::Penetration::Terrisol1990::gEnergies_T1990

static

Initial value:

=
{ 0.2, 0.5, 1, 2, 3, 4, 5, 6, 7,
  
  8, 9}

Definition at line 108 of file G4DNAOneStepThermalizationModel.hh.

Referenced by Get3DStdDeviation().

gStdDev_T1990

const double DNA::Penetration::Terrisol1990::gStdDev_T1990

static

Initial value:

=
{ 17.68*CLHEP::angstrom,
  22.3*CLHEP::angstrom,
  28.49*CLHEP::angstrom,
  45.35*CLHEP::angstrom,
  70.03*CLHEP::angstrom,
  98.05*CLHEP::angstrom,
  120.56*CLHEP::angstrom,
  132.73*CLHEP::angstrom,
  142.60*CLHEP::angstrom,
  
  
  
  
  
  137.9*CLHEP::angstrom,
  120.7*CLHEP::angstrom
}

Definition at line 109 of file G4DNAOneStepThermalizationModel.hh.

Referenced by Get3DStdDeviation().

---

The documentation for this struct was generated from the following files:

• G4DNAOneStepThermalizationModel.hh
• G4DNAOneStepThermalizationModel.cc