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 79 of file G4DNAOneStepThermalizationModel.hh.

Member Function Documentation

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

static

Definition at line 132 of file G4DNAOneStepThermalizationModel.cc.

                                                     {
    G4double k_eV = energy/eV;
    if(k_eV < 0.2) return 1e-3*CLHEP::nanometer;
    // rare events:
    //  prevent H2O and secondary electron to be at the spot
    
    if(k_eV == 9.) return gStdDev_T1990[10];
    // TODO if k_eV > 9

    size_t lowBin, upBin;
    
    if(k_eV >= 1.){
      lowBin=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];
    
    // G4cout << lowE << " " << upE << G4endl;
    
    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;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

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

static

Definition at line 178 of file G4DNAOneStepThermalizationModel.cc.

                                                                {
    double sigma3D=Get3DStdDeviation(energy);
    // G4cout << "sigma3D = " << sigma3D/CLHEP::nanometer << G4endl;
    
    static constexpr double factor = 2.20496999539;
    // 1./(3. - 8./CLHEP::pi);
    
    double sigma1D = sqrt(pow(sigma3D, 2.)*factor);
    
    // G4cout << "sigma1D = " << sigma1D/CLHEP::nanometer << G4endl;

    double x = G4RandGauss::shoot(0.,sigma1D);
    double y = G4RandGauss::shoot(0.,sigma1D);
    double z = G4RandGauss::shoot(0.,sigma1D);
    displacement=G4ThreeVector(x,y,z);
    // G4cout << "displacement[nm]: "
    //        << displacement.mag()/CLHEP::nanometer << G4endl;
  }

Here is the call graph for this function:

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

static

Definition at line 168 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;
  }

Here is the call graph for this function:

Member Data Documentation

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 86 of file G4DNAOneStepThermalizationModel.hh.

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 87 of file G4DNAOneStepThermalizationModel.hh.

---

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

• geant4.10.03.p01/source/processes/electromagnetic/dna/models/include/G4DNAOneStepThermalizationModel.hh
• geant4.10.03.p01/source/processes/electromagnetic/dna/models/src/G4DNAOneStepThermalizationModel.cc