G4DNAWaterDissociationDisplacer.cc File Reference

#include "G4DNAWaterDissociationDisplacer.hh"
#include "G4PhysicalConstants.hh"
#include "G4SystemOfUnits.hh"
#include "G4H2O.hh"
#include "G4H2.hh"
#include "G4Hydrogen.hh"
#include "G4OH.hh"
#include "G4H3O.hh"
#include "G4Electron_aq.hh"
#include "G4H2O2.hh"
#include "Randomize.hh"
#include "G4Molecule.hh"
#include "G4MolecularConfiguration.hh"
#include "G4RandomDirection.hh"
#include "G4Exp.hh"
#include "G4UnitsTable.hh"

Include dependency graph for G4DNAWaterDissociationDisplacer.cc:

Go to the source code of this file.

Functions
	G4CT_COUNT_IMPL (G4DNAWaterDissociationDisplacer, Ionisation_DissociationDecay) G4CT_COUNT_IMPL(G4DNAWaterDissociationDisplacer
A1B1_DissociationDecay	G4CT_COUNT_IMPL (G4DNAWaterDissociationDisplacer, B1A1_DissociationDecay) G4CT_COUNT_IMPL(G4DNAWaterDissociationDisplacer
A1B1_DissociationDecay AutoIonisation	G4CT_COUNT_IMPL (G4DNAWaterDissociationDisplacer, DissociativeAttachment) G4DNAWaterDissociationDisplacer

Function Documentation

G4CT_COUNT_IMPL	(	G4DNAWaterDissociationDisplacer	,
		Ionisation_DissociationDecay
	)

A1B1_DissociationDecay G4CT_COUNT_IMPL	(	G4DNAWaterDissociationDisplacer	,
		B1A1_DissociationDecay
	)

A1B1_DissociationDecay AutoIonisation G4CT_COUNT_IMPL	(	G4DNAWaterDissociationDisplacer	,
		DissociativeAttachment
	)

Definition at line 68 of file G4DNAWaterDissociationDisplacer.cc.

{
  return (2.*r+1.)*G4Exp(-2.*r);
}
#endif

//------------------------------------------------------------------------------
#ifdef _WATER_DISPLACER_USE_TERRISOL_
// This function is used to generate the following density distribution:
//   f(r) := 4*x^2/(sqrt(%pi)*(b)^3)*exp(-x^2/(b)^2)
// with b=27.22 for 7 eV
// reference
// Terrissol M, Beaudre A (1990) Simulation of space and time evolution
// of radiolytic species induced by electrons in water.
// Radiat Prot Dosimetry 31:171–175

G4double G4DNAWaterDissociationDisplacer::ElectronProbaDistribution(G4double r)
{
#define b 27.22 //*nanometer
  static constexpr double sqrt_pi=1.77245; // sqrt(CLHEP::pi);
  static constexpr double b_to3 = 20168.1; // pow(b,3.);
  static constexpr double b_to2 = 740.928; // pow(b,2.);
  static constexpr double inverse_b_to2 = 1./b_to2;

  static constexpr double main_factor=4./(sqrt_pi*b_to3);
  static constexpr double factorA=sqrt_pi*b_to3/4.;
  static constexpr double factorB=b_to2/2.;

  return (main_factor*
          (factorA*erf(r/b)
            - factorB*r*G4Exp(-pow(r,2.)*inverse_b_to2)));
}
#endif
//------------------------------------------------------------------------------

G4DNAWaterDissociationDisplacer::G4DNAWaterDissociationDisplacer() :
    G4VMolecularDecayDisplacer(),
#ifdef _WATER_DISPLACER_USE_KREIPL_
    fFastElectronDistrib(0., 5., 0.001)
#elif defined _WATER_DISPLACER_USE_TERRISOL_
    fFastElectronDistrib(0., 100., 0.001)
#endif
{
  fProba1DFunction =
      std::bind((G4double(*)(G4double))
                &G4DNAWaterDissociationDisplacer::ElectronProbaDistribution,
                std::placeholders::_1);

  size_t nBins = 500;
  fElectronThermalization.reserve(nBins);
  double eps = 1./((int)nBins);
  double proba = eps;

  fElectronThermalization.push_back(0.);

  for(size_t i = 1 ; i < nBins ; ++i){
    double r = fFastElectronDistrib.Revert(proba, fProba1DFunction);
    fElectronThermalization.push_back(r*nanometer);
    proba+=eps;
//  G4cout << G4BestUnit(r*nanometer, "Length") << G4endl;
  }
//   SetVerbose(1);
}

Here is the call graph for this function: