G4DNAWaterDissociationDisplacer Class Reference

#include <G4DNAWaterDissociationDisplacer.hh>

Inheritance diagram for G4DNAWaterDissociationDisplacer:

Collaboration 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 G4VMolecularDecayDisplacer
virtual	~G4VMolecularDecayDisplacer ()
void	SetVerbose (G4int verbose)

Static Public Member Functions
static G4double	ElectronProbaDistribution (G4double r)

Additional Inherited Members
Protected Member Functions inherited from G4VMolecularDecayDisplacer
	G4VMolecularDecayDisplacer ()
Protected Attributes inherited from G4VMolecularDecayDisplacer
G4int	fVerbose

Detailed Description

Definition at line 57 of file G4DNAWaterDissociationDisplacer.hh.

Constructor & Destructor Documentation

G4DNAWaterDissociationDisplacer::G4DNAWaterDissociationDisplacer

(

)

G4DNAWaterDissociationDisplacer::~G4DNAWaterDissociationDisplacer ( )

virtual

Definition at line 146 of file G4DNAWaterDissociationDisplacer.cc.

{;}

Member Function Documentation

static G4double G4DNAWaterDissociationDisplacer::ElectronProbaDistribution ( G4double  r )

static

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

overridevirtual

Implements G4VMolecularDecayDisplacer.

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

Here is the call graph for this function:

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

overridevirtual

Implements G4VMolecularDecayDisplacer.

Definition at line 190 of file G4DNAWaterDissociationDisplacer.cc.

{
  G4int nbProducts = theDecayChannel->GetNbProducts();
  vector<G4ThreeVector> theProductDisplacementVector(nbProducts);

  typedef map<const G4MoleculeDefinition*, G4double> RMSmap;
  RMSmap theRMSmap;

  G4int decayType = theDecayChannel->GetDisplacementType();
  
  switch(decayType){
    case Ionisation_DissociationDecay:
    {
      if (fVerbose){
        G4cout << "Ionisation_DissociationDecay" << G4endl;
        G4cout << "Channel's name: " << theDecayChannel->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++){
        G4MolecularConfiguration* product = theDecayChannel->GetProduct(i);
        G4double theRMSDisplacement = theRMSmap[product->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: " << theDecayChannel->GetName() << G4endl;
      }
      G4double theRMSDisplacement = 2.4 * nanometer;
      auto RandDirection =
      radialDistributionOfProducts(theRMSDisplacement);
      
      for(G4int i =0; i < nbProducts; i++){
        G4MolecularConfiguration* product = theDecayChannel->GetProduct(i);
        
        if(product->GetDefinition()== G4OH::Definition()){
          // OH
          theProductDisplacementVector[i] = -1./18.*RandDirection;
        }
        else if(product->GetDefinition() == G4Hydrogen::Definition()){
          // H
          theProductDisplacementVector[i] = +17./18.*RandDirection;
        }
      }
      break;
    }
    case B1A1_DissociationDecay:
    {
      if(fVerbose){
        G4cout<<"B1A1_DissociationDecay"<<G4endl;
        G4cout << "Channel's name: " << theDecayChannel->GetName() << G4endl;
      }
      
      G4double theRMSDisplacement = 0.8 * nanometer;
      auto RandDirection =
      radialDistributionOfProducts(theRMSDisplacement);
      
      G4int NbOfOH = 0;
      for(G4int i =0; i < nbProducts; ++i)
      {
        G4MolecularConfiguration* product = theDecayChannel->GetProduct(i);
        if(product->GetDefinition() == G4H2::Definition()){
          // H2
          theProductDisplacementVector[i] = -2./18.*RandDirection;
        }
        else if(product->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: " << theDecayChannel->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++){
        G4MolecularConfiguration* product = theDecayChannel->GetProduct(i);
        auto theRMSDisplacement = theRMSmap[product->GetDefinition()];
        
        if(theRMSDisplacement==0){
          theProductDisplacementVector[i] = G4ThreeVector();
        }
        else{
          auto RandDirection =
          radialDistributionOfProducts(theRMSDisplacement);
          theProductDisplacementVector[i] = RandDirection;
        }
        if(product->GetDefinition() == G4Electron_aq::Definition()){
          theProductDisplacementVector[i]=radialDistributionOfElectron();
        }
      }
      break;
    }
    case DissociativeAttachment:
    {
      if(fVerbose){
        G4cout<<"DissociativeAttachment"<<G4endl;
        G4cout << "Channel's name: " << theDecayChannel->GetName() << G4endl;
      }
      G4double theRMSDisplacement = 0.8 * nanometer;
      auto RandDirection =
      radialDistributionOfProducts(theRMSDisplacement);
      
      G4int NbOfOH = 0;
      for(G4int i =0; i < nbProducts; ++i){
        G4MolecularConfiguration* product = theDecayChannel->GetProduct(i);
        if(product->GetDefinition() == G4H2::Definition()){
          theProductDisplacementVector[i] = -2./18.*RandDirection;
        }
        else if(product->GetDefinition() == G4OH::Definition()){
          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;
    }
  }
  return theProductDisplacementVector;
}

Here is the call graph for this function:

G4ThreeVector G4DNAWaterDissociationDisplacer::radialDistributionOfElectron

(

)

const

Definition at line 400 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();
//   return fElectronThermalization[bin]*G4RandomDirection();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector G4DNAWaterDissociationDisplacer::radialDistributionOfProducts

(

G4double

r_rms

)

const

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

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

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