G4DNAUpdateSystemModel Class Reference

#include <G4DNAUpdateSystemModel.hh>

Inheritance diagram for G4DNAUpdateSystemModel:

Public Types
using	Index = G4Voxel::Index
using	JumpingData = std::pair< MolType, Index >
using	MolType = const G4MolecularConfiguration *
using	ReactionData = const G4DNAMolecularReactionData

Public Member Functions
void	CreateMolecule (const Index &index, MolType)
	G4DNAUpdateSystemModel ()
void	JumpIn (const Index &index, MolType)
void	JumpTo (const Index &index, MolType type)
void	KillMolecule (const Index &index, MolType type)
void	SetGlobalTime (const G4double &globalTime)
void	SetMesh (G4DNAMesh *)
void	SetVerbose (G4int verbose)
void	UpdateSystem (const Index &index, const JumpingData &data)
void	UpdateSystem (const Index &index, const ReactionData &data)
	~G4DNAUpdateSystemModel () override=default

Private Attributes
G4double	fGlobalTime
G4DNAMesh *	fpMesh
G4int	fVerbose

Detailed Description

Definition at line 36 of file G4DNAUpdateSystemModel.hh.

Member Typedef Documentation

Index

using G4DNAUpdateSystemModel::Index = G4Voxel::Index

Definition at line 39 of file G4DNAUpdateSystemModel.hh.

JumpingData

using G4DNAUpdateSystemModel::JumpingData = std::pair<MolType, Index>

Definition at line 41 of file G4DNAUpdateSystemModel.hh.

MolType

using G4DNAUpdateSystemModel::MolType = const G4MolecularConfiguration*

Definition at line 40 of file G4DNAUpdateSystemModel.hh.

ReactionData

using G4DNAUpdateSystemModel::ReactionData = const G4DNAMolecularReactionData

Definition at line 42 of file G4DNAUpdateSystemModel.hh.

Constructor & Destructor Documentation

G4DNAUpdateSystemModel()

G4DNAUpdateSystemModel::G4DNAUpdateSystemModel

(

)

Definition at line 34 of file G4DNAUpdateSystemModel.cc.

  : G4VUpdateSystemModel()
  , fpMesh(nullptr)
  , fVerbose(0)
  , fGlobalTime(DBL_MAX)
{}

~G4DNAUpdateSystemModel()

G4DNAUpdateSystemModel::~G4DNAUpdateSystemModel ( )

overridedefault

Member Function Documentation

CreateMolecule()

void G4DNAUpdateSystemModel::CreateMolecule	(	const Index &	index,
		MolType	type
	)

Definition at line 105 of file G4DNAUpdateSystemModel.cc.

{
  // for scavenger
  auto pScavengerMaterial = dynamic_cast<G4DNAScavengerMaterial*>(
    G4Scheduler::Instance()->GetScavengerMaterial());
  if(pScavengerMaterial != nullptr && pScavengerMaterial->find(type))
  {
    pScavengerMaterial->AddNumberMoleculePerVolumeUnitForMaterialConf(
      type, fGlobalTime);
    return;
  }
  // for molecule
  auto& node = fpMesh->GetVoxelMapList(index);
  auto iter  = node.find(type);
  if(iter != node.end())
  {
    iter->second++;
  }
  else
  {
    node[type] = 1;
  }
 
  if(G4VMoleculeCounter::Instance()->InUse())
  {
    G4VMoleculeCounter::Instance()->AddAMoleculeAtTime(type, fGlobalTime);
  }
}

References G4VMoleculeCounter::AddAMoleculeAtTime(), G4DNAScavengerMaterial::AddNumberMoleculePerVolumeUnitForMaterialConf(), fGlobalTime, fpMesh, G4Scheduler::GetScavengerMaterial(), G4DNAMesh::GetVoxelMapList(), G4Scheduler::Instance(), G4VMoleculeCounter::Instance(), and G4VMoleculeCounter::InUse().

Referenced by UpdateSystem().

JumpIn()

void G4DNAUpdateSystemModel::JumpIn	(	const Index &	index,
		MolType	type
	)

Definition at line 134 of file G4DNAUpdateSystemModel.cc.

{
  // for molecule
  auto& node = fpMesh->GetVoxelMapList(index);
  auto iter  = node.find(type);
  if(iter != node.end())
  {
    iter->second++;
  }
  else
  {
    node[type] = 1;
  }
}

References fpMesh, and G4DNAMesh::GetVoxelMapList().

Referenced by UpdateSystem().

JumpTo()

void G4DNAUpdateSystemModel::JumpTo	(	const Index &	index,
		MolType	type
	)

Definition at line 81 of file G4DNAUpdateSystemModel.cc.

{
  auto& node = fpMesh->GetVoxelMapList(index);
 
  auto iter = node.find(type);
  if(iter != node.end())
  {
    if(iter->second <= 0)
    {
      G4cout << "G4DNAUpdateSystemModel::KillMolecule::molecule : "
             << type->GetName() << " index : " << index
             << " number : " << iter->second << G4endl;
      assert(false);
    }
    iter->second--;
  }
  else
  {
    G4cout << "index : " << index << " " << type->GetName() << G4endl;
    G4cout << "This molecule is not belong  particle-base" << G4endl;
    assert(false);
  }
}

References fpMesh, G4cout, G4endl, G4MolecularConfiguration::GetName(), and G4DNAMesh::GetVoxelMapList().

Referenced by UpdateSystem().

KillMolecule()

void G4DNAUpdateSystemModel::KillMolecule	(	const Index &	index,
		MolType	type
	)

Definition at line 42 of file G4DNAUpdateSystemModel.cc.

{
  // kill normal molecule
  auto& node = fpMesh->GetVoxelMapList(index);
  auto iter  = node.find(type);
  if(iter != node.end())
  {
    if(iter->second <= 0)
    {
      G4cout << "G4DNAUpdateSystemModel::KillMolecule::molecule : "
             << type->GetName() << " index : " << index
             << " number : " << iter->second << G4endl;
      assert(false);
    }
    iter->second--;
    if(G4VMoleculeCounter::Instance()->InUse())
    {
      G4VMoleculeCounter::Instance()->RemoveAMoleculeAtTime(type, fGlobalTime);
    }
  }
  else
  {
    auto pScavengerMaterial = dynamic_cast<G4DNAScavengerMaterial*>(
      G4Scheduler::Instance()->GetScavengerMaterial());
    if(pScavengerMaterial != nullptr)
    {
      pScavengerMaterial->ReduceNumberMoleculePerVolumeUnitForMaterialConf(
        type, fGlobalTime);
    }
    else
    {
      G4cout << "index : " << index << " " << type->GetName() << G4endl;
      G4cout << "This molecule is not belong scavengers or particle-base"
             << G4endl;
      assert(false);
    }
  }
}

References fGlobalTime, fpMesh, G4cout, G4endl, G4MolecularConfiguration::GetName(), G4Scheduler::GetScavengerMaterial(), G4DNAMesh::GetVoxelMapList(), G4Scheduler::Instance(), G4VMoleculeCounter::Instance(), G4DNAScavengerMaterial::ReduceNumberMoleculePerVolumeUnitForMaterialConf(), and G4VMoleculeCounter::RemoveAMoleculeAtTime().

Referenced by UpdateSystem().

SetGlobalTime()

void G4DNAUpdateSystemModel::SetGlobalTime ( const G4double &  globalTime )

inline

Definition at line 56 of file G4DNAUpdateSystemModel.hh.

{ fGlobalTime = globalTime; }

References fGlobalTime.

SetMesh()

void G4DNAUpdateSystemModel::SetMesh

(

G4DNAMesh *

pMesh

)

Definition at line 41 of file G4DNAUpdateSystemModel.cc.

{ fpMesh = pMesh; }

References fpMesh.

SetVerbose()

void G4DNAUpdateSystemModel::SetVerbose ( G4int  verbose )

inline

Definition at line 57 of file G4DNAUpdateSystemModel.hh.

{ fVerbose = verbose; }

References fVerbose.

UpdateSystem() [1/2]

void G4DNAUpdateSystemModel::UpdateSystem	(	const Index &	index,
		const JumpingData &	data
	)

Definition at line 260 of file G4DNAUpdateSystemModel.cc.

{
  auto reactant    = std::get<0>(data);
  auto JunpToIndex = std::get<1>(data);
#ifdef G4VERBOSE
  if(fVerbose > 1)
  {
    G4cout << "At time : " << std::setw(7) << G4BestUnit(fGlobalTime, "Time")
           << " Jumping : " << reactant->GetName() << " from " << index
           << " -> " << JunpToIndex << G4endl;
  }
#endif
  JumpTo(index, reactant);
  JumpIn(JunpToIndex, reactant);
}

References fGlobalTime, fVerbose, G4BestUnit, G4cout, G4endl, JumpIn(), and JumpTo().

UpdateSystem() [2/2]

void G4DNAUpdateSystemModel::UpdateSystem	(	const Index &	index,
		const ReactionData &	data
	)

Definition at line 149 of file G4DNAUpdateSystemModel.cc.

{
  auto reactant1 = data.GetReactant1();
  auto reactant2 = data.GetReactant2();
#ifdef G4VERBOSE
  if(fVerbose != 0)
  {
    G4cout << "At time : " << std::setw(7) << G4BestUnit(fGlobalTime, "Time")
           << " Reaction : " << reactant1->GetName() << " + "
           << reactant2->GetName() << " -> ";
  }
#endif
  const G4int nbProducts = data.GetNbProducts();
  if(nbProducts != 0)
  {
    for(size_t j = 0; j < (size_t) nbProducts; ++j)
    {
#ifdef G4VERBOSE
      if((fVerbose != 0) && j != 0)
      {
        G4cout << " + ";
      }
      if(fVerbose != 0)
      {
        G4cout << data.GetProduct(j)->GetName();
        // for test
        // G4cout<<"  fGlobalTime : "<<fGlobalTime;
        // end fortest
      }
#endif
      CreateMolecule(index, data.GetProduct(j));
      //#define DEBUG 1
 
#ifdef DEBUG
      if(G4MoleculeCounter::Instance()->InUse())
        if(fpMesh->GetNumberOfType(data.GetProduct(j)) !=
           G4MoleculeCounter::Instance()->GetCurrentNumberOf(
             data.GetProduct(j)))
        {
          G4cout << "*********G4DNAUpdateSystemModel::DEBUG::GetNumberOfType("
                 << data.GetProduct(j)->GetName()
                 << ") : " << fpMesh->GetNumberOfType(data.GetProduct(j))
                 << G4endl;
          G4cout << "G4MoleculeCounter::GetCurrentNumberOf ("
                 << data.GetProduct(j)->GetName() << ") : "
                 << G4MoleculeCounter::Instance()->GetCurrentNumberOf(
                      data.GetProduct(j))
                 << G4endl;
          G4MoleculeCounter::Instance()->Dump();
          throw;
        }
 
#endif
    }
  }
  else
  {
#ifdef G4VERBOSE
    if(fVerbose != 0)
    {
      G4cout << "No product";
      // for test
      // G4cout<<"  fGlobalTime : "<<fGlobalTime;
      // end fortest
    }
#endif
  }
#ifdef G4VERBOSE
  if(fVerbose != 0)
  {
    G4cout << G4endl;
  }
#endif
  KillMolecule(index, reactant1);
#ifdef DEBUG
  if(G4MoleculeCounter::Instance()->InUse())
    if(fpMesh->GetNumberOfType(reactant1) !=
       G4MoleculeCounter::Instance()->GetCurrentNumberOf(reactant1))
    {
      G4cout << "*********G4DNAUpdateSystemModel::DEBUG::GetNumberOfType("
             << reactant1->GetName()
             << ") : " << fpMesh->GetNumberOfType(reactant1) << G4endl;
      G4cout << "G4MoleculeCounter::GetCurrentNumberOf ("
             << reactant1->GetName() << ") : "
             << G4MoleculeCounter::Instance()->GetCurrentNumberOf(reactant1)
             << G4endl;
      G4MoleculeCounter::Instance()->Dump();
      throw;
    }
#endif
  KillMolecule(index, reactant2);
#ifdef DEBUG
 
  if(G4MoleculeCounter::Instance()->InUse())
    if(fpMesh->GetNumberOfType(reactant2) !=
       G4MoleculeCounter::Instance()->GetCurrentNumberOf(reactant2))
    {
      G4cout << "*********G4DNAUpdateSystemModel::DEBUG::GetNumberOfType("
             << reactant2->GetName()
             << ") : " << fpMesh->GetNumberOfType(reactant2) << G4endl;
      G4cout << "G4MoleculeCounter::GetCurrentNumberOf ("
             << reactant2->GetName() << ") : "
             << G4MoleculeCounter::Instance()->GetCurrentNumberOf(reactant2)
             << G4endl;
      G4MoleculeCounter::Instance()->Dump();
      throw;
    }
#endif
}

References CreateMolecule(), G4MoleculeCounter::Dump(), fGlobalTime, fpMesh, fVerbose, G4BestUnit, G4cout, G4endl, G4MolecularConfiguration::GetName(), G4DNAMolecularReactionData::GetNbProducts(), G4DNAMesh::GetNumberOfType(), G4DNAMolecularReactionData::GetProduct(), G4DNAMolecularReactionData::GetReactant1(), G4DNAMolecularReactionData::GetReactant2(), G4MoleculeCounter::Instance(), and KillMolecule().

Field Documentation

fGlobalTime

G4double G4DNAUpdateSystemModel::fGlobalTime

private

Definition at line 62 of file G4DNAUpdateSystemModel.hh.

Referenced by CreateMolecule(), KillMolecule(), SetGlobalTime(), and UpdateSystem().

fpMesh

G4DNAMesh* G4DNAUpdateSystemModel::fpMesh

private

Definition at line 60 of file G4DNAUpdateSystemModel.hh.

Referenced by CreateMolecule(), JumpIn(), JumpTo(), KillMolecule(), SetMesh(), and UpdateSystem().

fVerbose

G4int G4DNAUpdateSystemModel::fVerbose

private

Definition at line 61 of file G4DNAUpdateSystemModel.hh.

Referenced by SetVerbose(), and UpdateSystem().

---

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

• source/processes/electromagnetic/dna/models/include/G4DNAUpdateSystemModel.hh
• source/processes/electromagnetic/dna/models/src/G4DNAUpdateSystemModel.cc