#include <G4DNATotallyDiffusionControlled.hh>

Inheritance diagram for G4DNATotallyDiffusionControlled:

Public Member Functions
	G4DNATotallyDiffusionControlled ()

	G4DNATotallyDiffusionControlled (const G4DNATotallyDiffusionControlled &other)=delete

G4bool	GeminateRecombinationProbability (const G4MolecularConfiguration pConfMolA, const G4MolecularConfiguration pConfMolB) override

G4double	GetTimeToEncounter (const G4Track &trackA, const G4Track &trackB) override

G4DNATotallyDiffusionControlled &	operator= (const G4DNATotallyDiffusionControlled &other)=delete

	~G4DNATotallyDiffusionControlled () override

Private Member Functions
G4double	GetDiffusionCoefficient (const G4MolecularConfiguration pMA, const G4MolecularConfiguration pMB)

Detailed Description

Definition at line 39 of file G4DNATotallyDiffusionControlled.hh.

Constructor & Destructor Documentation

◆ G4DNATotallyDiffusionControlled() [1/2]

G4DNATotallyDiffusionControlled::G4DNATotallyDiffusionControlled ( )

Definition at line 42 of file G4DNATotallyDiffusionControlled.cc.

43 : G4VReactionType()

44{}

G4VReactionType::G4VReactionType

G4VReactionType()=default

◆ ~G4DNATotallyDiffusionControlled()

G4DNATotallyDiffusionControlled::~G4DNATotallyDiffusionControlled ( )

overridedefault

◆ G4DNATotallyDiffusionControlled() [2/2]

G4DNATotallyDiffusionControlled::G4DNATotallyDiffusionControlled ( const G4DNATotallyDiffusionControlled & other )

delete

Member Function Documentation

◆ GeminateRecombinationProbability()

G4bool G4DNATotallyDiffusionControlled::GeminateRecombinationProbability	(	const G4MolecularConfiguration *	pConfMolA,
		const G4MolecularConfiguration *	pConfMolB
	)

overridevirtual

Implements G4VReactionType.

Definition at line 85 of file G4DNATotallyDiffusionControlled.cc.

{
    if(pMolA->GetDefinition() == G4Electron_aq::Definition() ||
       pMolA->GetDefinition() == G4Hydrogen::Definition())
    {
        G4bool spinA;
        G4bool spinB;
        spinA = G4UniformRand() < 0.5;
        if(spinA &&
          (pMolB->GetDefinition() == G4Electron_aq::Definition() ||
           pMolB->GetDefinition() ==  G4Hydrogen::Definition()))
        {
            spinB = G4UniformRand() < 0.5;
            if( !spinB )
            {
                return true;
            }
        }
        return false;
    }
    return true;
}

References G4Electron_aq::Definition(), G4Hydrogen::Definition(), G4UniformRand, and G4MolecularConfiguration::GetDefinition().

◆ GetDiffusionCoefficient()

G4double G4DNATotallyDiffusionControlled::GetDiffusionCoefficient	(	const G4MolecularConfiguration *	pMA,
		const G4MolecularConfiguration *	pMB
	)

private

Definition at line 111 of file G4DNATotallyDiffusionControlled.cc.

{
    G4double D;
     if(mA == mB)
     {    
         D = (mA->GetDiffusionCoefficient());
     }
     else
    {
        D = (mA->GetDiffusionCoefficient() +
             mB->GetDiffusionCoefficient());//
    }
    return D;
}

References D(), and G4MolecularConfiguration::GetDiffusionCoefficient().

Referenced by GetTimeToEncounter().

◆ GetTimeToEncounter()

G4double G4DNATotallyDiffusionControlled::GetTimeToEncounter	(	const G4Track &	trackA,
		const G4Track &	trackB
	)

overridevirtual

Implements G4VReactionType.

Definition at line 47 of file G4DNATotallyDiffusionControlled.cc.

{
    auto pMolConfA = GetMolecule(trackA)->GetMolecularConfiguration();
    auto pMolConfB = GetMolecule(trackB)->GetMolecularConfiguration();
 
    G4double D = GetDiffusionCoefficient(pMolConfA, pMolConfB);
    auto reactionData = G4DNAMolecularReactionTable::Instance()
    ->GetReactionData(pMolConfA, pMolConfB);
    G4double kobs = reactionData->GetObservedReactionRateConstant();
    G4double distance = (trackA.GetPosition() - trackB.GetPosition()).mag();
    G4double Reff  = kobs / ( 4 * CLHEP::pi * D * Avogadro );
 
    if( distance < Reff )
    {
         G4ExceptionDescription exceptionDescription;
         exceptionDescription << "distance = "<< distance
                              << " is uncorrected with "
                              <<" Reff = "<<Reff
                              <<" for : "<<pMolConfA->GetName()
                              <<" and "<<pMolConfB->GetName();
         G4Exception("G4DNATotallyDiffusionControlled"
                     "::GetTimeToEncounter()", "G4DNATotallyDiffusionControlled02",
         FatalException, exceptionDescription);
    }
    
    G4double Winf = Reff / distance;
    G4double U = G4UniformRand();
    G4double irt = -1.0 * ps;
 
    if ( U < Winf )
    {
        G4double d = ( distance - Reff ) /
                    G4ErrorFunction::erfcInv( U / Winf );
        irt = ( 1.0 / ( 4 * D ) ) * d * d;
    }
    return irt;
}

References source.hepunit::Avogadro, D(), G4ErrorFunction::erfcInv(), FatalException, G4Exception(), G4UniformRand, GetDiffusionCoefficient(), G4Molecule::GetMolecularConfiguration(), GetMolecule(), G4DNAMolecularReactionData::GetObservedReactionRateConstant(), G4Track::GetPosition(), G4DNAMolecularReactionTable::GetReactionData(), G4DNAMolecularReactionTable::Instance(), CLHEP::pi, and ps.

◆ operator=()

G4DNATotallyDiffusionControlled & G4DNATotallyDiffusionControlled::operator= ( const G4DNATotallyDiffusionControlled & other )