G4DNAWaterDissociationDisplacer Class Reference

#include <G4DNAWaterDissociationDisplacer.hh>

Inheritance diagram for G4DNAWaterDissociationDisplacer:

Public Member Functions
	G4DNAWaterDissociationDisplacer ()
virtual	~G4DNAWaterDissociationDisplacer ()
std::vector< G4ThreeVector >	GetProductsDisplacement (const G4MolecularDissociationChannel *) const override
G4ThreeVector	GetMotherMoleculeDisplacement (const G4MolecularDissociationChannel *) const override
G4ThreeVector	radialDistributionOfElectron () const
G4ThreeVector	radialDistributionOfProducts (G4double r_rms) const
Public Member Functions inherited from G4VMolecularDissociationDisplacer
virtual	~G4VMolecularDissociationDisplacer ()=default
virtual std::vector< G4ThreeVector >	GetProductsDisplacement (const G4MolecularDissociationChannel *) const =0
virtual G4ThreeVector	GetMotherMoleculeDisplacement (const G4MolecularDissociationChannel *) const =0
void	SetVerbose (G4int verbose)

Static Public Member Functions
static G4double	ElectronProbaDistribution (G4double r)

Additional Inherited Members
Protected Member Functions inherited from G4VMolecularDissociationDisplacer
	G4VMolecularDissociationDisplacer ()=default
Protected Attributes inherited from G4VMolecularDissociationDisplacer
G4int	fVerbose

Detailed Description

Definition at line 57 of file G4DNAWaterDissociationDisplacer.hh.

Constructor & Destructor Documentation

G4DNAWaterDissociationDisplacer()

G4DNAWaterDissociationDisplacer::G4DNAWaterDissociationDisplacer

(

)

~G4DNAWaterDissociationDisplacer()

G4DNAWaterDissociationDisplacer::~G4DNAWaterDissociationDisplacer ( )

virtual

Definition at line 156 of file G4DNAWaterDissociationDisplacer.cc.

{
    ;
}

Member Function Documentation

ElectronProbaDistribution()

static G4double G4DNAWaterDissociationDisplacer::ElectronProbaDistribution ( G4double  r )

static

GetMotherMoleculeDisplacement()

G4ThreeVector G4DNAWaterDissociationDisplacer::GetMotherMoleculeDisplacement ( const G4MolecularDissociationChannel *  theDecayChannel ) const

overridevirtual

Implements G4VMolecularDissociationDisplacer.

Definition at line 164 of file G4DNAWaterDissociationDisplacer.cc.

{
    G4int decayType = theDecayChannel->GetDisplacementType();
    G4double RMSMotherMoleculeDisplacement(0.);
 
    switch (decayType)
    {
        case Ionisation_DissociationDecay:
            RMSMotherMoleculeDisplacement = 2.0 * nanometer;
            break;
        case A1B1_DissociationDecay:
            RMSMotherMoleculeDisplacement = 0. * nanometer;
            break;
        case B1A1_DissociationDecay:
            RMSMotherMoleculeDisplacement = 0. * nanometer;
            break;
        case AutoIonisation:
            RMSMotherMoleculeDisplacement = 2.0 * nanometer;
            break;
        case DissociativeAttachment:
            RMSMotherMoleculeDisplacement = 0. * nanometer;
            break;
    }
 
    if (RMSMotherMoleculeDisplacement == 0)
    {
        return G4ThreeVector(0, 0, 0);
    }
    auto RandDirection =
            radialDistributionOfProducts(RMSMotherMoleculeDisplacement);
 
    return RandDirection;
}

GetProductsDisplacement()

vector< G4ThreeVector > G4DNAWaterDissociationDisplacer::GetProductsDisplacement ( const G4MolecularDissociationChannel *  pDecayChannel ) const

overridevirtual

Implements G4VMolecularDissociationDisplacer.

Definition at line 203 of file G4DNAWaterDissociationDisplacer.cc.

{
    G4int nbProducts = pDecayChannel->GetNbProducts();
    vector<G4ThreeVector> theProductDisplacementVector(nbProducts);
 
    typedef map<const G4MoleculeDefinition*, G4double> RMSmap;
    RMSmap theRMSmap;
 
    G4int decayType = pDecayChannel->GetDisplacementType();
 
    switch (decayType)
    {
        case Ionisation_DissociationDecay:
        {
            if (fVerbose)
            {
                G4cout << "Ionisation_DissociationDecay" << G4endl;
                G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
            }
            G4double RdmValue = G4UniformRand();
 
            if (RdmValue < 0.5)
            {
                // H3O
                theRMSmap[G4H3O::Definition()] = 0. * nanometer;
                // OH
                theRMSmap[G4OH::Definition()] = 0.8 * nanometer;
            }
            else
            {
                // H3O
                theRMSmap[G4H3O::Definition()] = 0.8 * nanometer;
                // OH
                theRMSmap[G4OH::Definition()] = 0. * nanometer;
            }
 
            for (int i = 0; i < nbProducts; i++)
            {
                auto pProduct = pDecayChannel->GetProduct(i);
                G4double theRMSDisplacement = theRMSmap[pProduct->GetDefinition()];
 
                if (theRMSDisplacement == 0.)
                {
                    theProductDisplacementVector[i] = G4ThreeVector();
                }
                else
                {
                    auto RandDirection = radialDistributionOfProducts(theRMSDisplacement);
                    theProductDisplacementVector[i] = RandDirection;
                }
            }
            break;
        }
        case A1B1_DissociationDecay:
        {
            if (fVerbose)
            {
                G4cout << "A1B1_DissociationDecay" << G4endl;
                G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
            }
            G4double theRMSDisplacement = 2.4 * nanometer;
            auto RandDirection =
                    radialDistributionOfProducts(theRMSDisplacement);
 
            for (G4int i = 0; i < nbProducts; i++)
            {
                auto pProduct = pDecayChannel->GetProduct(i);
 
                if (pProduct->GetDefinition() == G4OH::Definition())
                {
                    theProductDisplacementVector[i] = -1. / 18. * RandDirection;
                }
                else if (pProduct->GetDefinition() == G4Hydrogen::Definition())
                {
                    theProductDisplacementVector[i] = +17. / 18. * RandDirection;
                }
            }
            break;
        }
        case B1A1_DissociationDecay:
        {
            if (fVerbose)
            {
                G4cout << "B1A1_DissociationDecay" << G4endl;
                G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
            }
 
            G4double theRMSDisplacement = 0.8 * nanometer;
            auto RandDirection =
                    radialDistributionOfProducts(theRMSDisplacement);
 
            G4int NbOfOH = 0;
            for (G4int i = 0; i < nbProducts; ++i)
            {
                auto pProduct = pDecayChannel->GetProduct(i);
                if (pProduct->GetDefinition() == G4H2::Definition())
                {
                    // H2
                    theProductDisplacementVector[i] = -2. / 18. * RandDirection;
                }
                else if (pProduct->GetDefinition() == G4OH::Definition())
                {
                    // OH
                    G4ThreeVector OxygenDisplacement = +16. / 18. * RandDirection;
                    G4double OHRMSDisplacement = 1.1 * nanometer;
 
                    auto OHDisplacement =
                            radialDistributionOfProducts(OHRMSDisplacement);
 
                    if (NbOfOH == 0)
                    {
                        OHDisplacement = 0.5 * OHDisplacement;
                    }
                    else
                    {
                        OHDisplacement = -0.5 * OHDisplacement;
                    }
 
                    theProductDisplacementVector[i] =
                            OHDisplacement + OxygenDisplacement;
 
                    ++NbOfOH;
                }
            }
            break;
        }
        case AutoIonisation:
        {
            if (fVerbose)
            {
                G4cout << "AutoIonisation" << G4endl;
                G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
            }
 
            G4double RdmValue = G4UniformRand();
 
            if (RdmValue < 0.5)
            {
                // H3O
                theRMSmap[G4H3O::Definition()] = 0. * nanometer;
                // OH
                theRMSmap[G4OH::Definition()] = 0.8 * nanometer;
            }
            else
            {
                // H3O
                theRMSmap[G4H3O::Definition()] = 0.8 * nanometer;
                // OH
                theRMSmap[G4OH::Definition()] = 0. * nanometer;
            }
 
            for (G4int i = 0; i < nbProducts; i++)
            {
                auto pProduct = pDecayChannel->GetProduct(i);
                auto theRMSDisplacement = theRMSmap[pProduct->GetDefinition()];
 
                if (theRMSDisplacement == 0)
                {
                    theProductDisplacementVector[i] = G4ThreeVector();
                }
                else
                {
                    auto RandDirection =
                            radialDistributionOfProducts(theRMSDisplacement);
                    theProductDisplacementVector[i] = RandDirection;
                }
                if (pProduct->GetDefinition() == G4Electron_aq::Definition())
                {
                    theProductDisplacementVector[i] = radialDistributionOfElectron();
                }
            }
            break;
        }
        case DissociativeAttachment:
        {
            if (fVerbose)
            {
                G4cout << "DissociativeAttachment" << G4endl;
                G4cout << "Channel's name: " << pDecayChannel->GetName() << G4endl;
            }
            G4double theRMSDisplacement = 0.8 * nanometer;
            auto RandDirection =
                    radialDistributionOfProducts(theRMSDisplacement);
 
            G4int NbOfOH = 0;
            for (G4int i = 0; i < nbProducts; ++i)
            {
                auto pProduct = pDecayChannel->GetProduct(i);
                if (pProduct->GetDefinition() == G4H2::Definition())
                {
                    // In the paper of Kreipl (2009)
                    // theProductDisplacementVector[i] = -2. / 18. * RandDirection;
 
                    // Based on momentum conservation
                    theProductDisplacementVector[i] = -16. / 18. * RandDirection;
                }
                else if (pProduct->GetDefinition() == G4OH::Definition())
                {
                    // In the paper of Kreipl (2009)
                    // G4ThreeVector OxygenDisplacement = +16. / 18. * RandDirection;
 
                    // Based on momentum conservation
                    G4ThreeVector OxygenDisplacement = +2. / 18. * RandDirection;
                    G4double OHRMSDisplacement = 1.1 * nanometer;
 
                    auto OHDisplacement =
                            radialDistributionOfProducts(OHRMSDisplacement);
 
                    if (NbOfOH == 0)
                    {
                        OHDisplacement = 0.5 * OHDisplacement;
                    }
                    else
                    {
                        OHDisplacement = -0.5 * OHDisplacement;
                    }
                    theProductDisplacementVector[i] = OHDisplacement +
                                                      OxygenDisplacement;
                    ++NbOfOH;
                }
            }
            break;
        }
    }
    return theProductDisplacementVector;
}

radialDistributionOfElectron()

G4ThreeVector G4DNAWaterDissociationDisplacer::radialDistributionOfElectron

(

)

const

Definition at line 448 of file G4DNAWaterDissociationDisplacer.cc.

{/*
    G4double rand_value = G4UniformRand();
    size_t nBins = fElectronThermalization.size();
    size_t bin = size_t(floor(rand_value * nBins));
    size_t bin_p1 = min(bin + 1, nBins - 1);
 
    return (fElectronThermalization[bin] * (1. - rand_value)
            + fElectronThermalization[bin_p1] * rand_value) *
           G4RandomDirection();*/
 
    G4ThreeVector pdf = G4ThreeVector(0,0,0);
 
    if(dnaSubType == fRitchie1994eSolvation) DNA::Penetration::Ritchie1994::GetPenetration(ke,pdf);
    else if(dnaSubType == fTerrisol1990eSolvation) DNA::Penetration::Terrisol1990::GetPenetration(ke,pdf);
    else if(dnaSubType == fMeesungnoensolid2002eSolvation) DNA::Penetration::Meesungnoen2002_amorphous::GetPenetration(ke,pdf);
    else if(dnaSubType == fKreipl2009eSolvation) DNA::Penetration::Kreipl2009::GetPenetration(ke,pdf);
    else DNA::Penetration::Meesungnoen2002::GetPenetration(ke,pdf);
 
    return pdf;
}

Referenced by GetProductsDisplacement().

radialDistributionOfProducts()

G4ThreeVector G4DNAWaterDissociationDisplacer::radialDistributionOfProducts

(

G4double

r_rms

)

const

Definition at line 434 of file G4DNAWaterDissociationDisplacer.cc.

{
    static const double inverse_sqrt_3 = 1. / sqrt(3.);
    double sigma = Rrms * inverse_sqrt_3;
    double x = G4RandGauss::shoot(0., sigma);
    double y = G4RandGauss::shoot(0., sigma);
    double z = G4RandGauss::shoot(0., sigma);
    return G4ThreeVector(x, y, z);
}

Referenced by GetMotherMoleculeDisplacement(), and GetProductsDisplacement().

---

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

• G4DNAWaterDissociationDisplacer.hh
• G4DNAWaterDissociationDisplacer.cc