#include <G4DNAEventScheduler.hh>

Inheritance diagram for G4DNAEventScheduler:

Public Types
using	MapCounter = std::map< MolType, G4int >

using	MapList = std::map< MolType, size_t >

using	MolType = const G4MolecularConfiguration *

Public Member Functions
void	AddTimeToRecord (const G4double &time)

void	ClearAndReChargeCounter ()

	G4DNAEventScheduler (const G4DNABoundingBox &boundingBox, G4int pixel)

	G4DNAEventScheduler (const G4DNAEventScheduler &)=delete

std::map< G4double, MapCounter >	GetCounterMap () const

G4double	GetEndTime () const

G4DNAMesh *	GetMesh () const

G4int	GetPixels () const

G4double	GetStartTime () const

G4double	GetTimeStep () const

G4int	GetVerbose () const

void	Initialize ()

void	InitializeInMesh ()

G4DNAEventScheduler &	operator= (const G4DNAEventScheduler &right)=delete

void	PrintRecordTime ()

void	RecordTime ()

void	Reset ()

void	ResetInMesh ()

void	ReVoxelizing (G4int)

void	Run ()

void	RunInMesh ()

void	SetChangeMesh (G4bool change)

void	SetEndTime (const G4double &)

void	SetMaxNbSteps (G4int)

void	SetStartTime (G4double time)

void	SetUserMeshAction (std::unique_ptr< G4UserMeshAction >)

void	SetVerbose (G4int verbose)

void	Stepping ()

void	Stop ()

void	Voxelizing ()

	~G4DNAEventScheduler () override

Static Public Member Functions
static G4bool	CheckingReactionRadius (G4double resolution)

Private Attributes
std::map< G4double, MapCounter >	fCounterMap

G4double	fEndTime

G4double	fGlobalTime

G4bool	fInitialized

G4double	fInitialPixels

G4bool	fIsChangeMesh

G4double	fJumpingNumber

std::set< G4double >::iterator	fLastRecoredTime

G4int	fMaxStep

std::unique_ptr< G4DNAEventSet >	fpEventSet

std::unique_ptr< G4DNAGillespieDirectMethod >	fpGillespieReaction

G4int	fPixel

std::unique_ptr< G4DNAMesh >	fpMesh

std::unique_ptr< G4DNAUpdateSystemModel >	fpUpdateSystem

std::unique_ptr< G4UserMeshAction >	fpUserMeshAction

G4double	fReactionNumber

G4bool	fRunning

G4bool	fSetChangeMesh

G4double	fStartTime

G4int	fStepNumber

G4int	fStepNumberInMesh

G4double	fTimeStep

std::set< G4double >	fTimeToRecord

G4int	fVerbose

Detailed Description

Definition at line 52 of file G4DNAEventScheduler.hh.

Member Typedef Documentation

◆ MapCounter

using G4DNAEventScheduler::MapCounter = std::map<MolType, G4int>

Definition at line 57 of file G4DNAEventScheduler.hh.

◆ MapList

using G4DNAEventScheduler::MapList = std::map<MolType, size_t>

Definition at line 56 of file G4DNAEventScheduler.hh.

◆ MolType

using G4DNAEventScheduler::MolType = const G4MolecularConfiguration*

Definition at line 55 of file G4DNAEventScheduler.hh.

Constructor & Destructor Documentation

◆ G4DNAEventScheduler() [1/2]

G4DNAEventScheduler::G4DNAEventScheduler	(	const G4DNABoundingBox &	boundingBox,
		G4int	pixel
	)

Definition at line 40 of file G4DNAEventScheduler.cc.

  : IEventScheduler()
  , fVerbose(0)
  , fInitialized(false)
  , fStartTime(1 * ps)
  , fEndTime(10000 * s)
  , fStepNumber(0)
  , fMaxStep(INT_MAX)
  , fRunning(true)
  , fTimeStep(DBL_MAX)
  , fGlobalTime(fStartTime)
  , fJumpingNumber(0)
  , fReactionNumber(0)
  , fPixel(pixel)
  , fIsChangeMesh(false)
  , fSetChangeMesh(true)
  , fStepNumberInMesh(0)
  , fInitialPixels(fPixel)
  , fpMesh(new G4DNAMesh(boundingBox, fPixel))
  , fpGillespieReaction(new G4DNAGillespieDirectMethod())
  , fpEventSet(new G4DNAEventSet())
  , fpUpdateSystem(new G4DNAUpdateSystemModel())
  , fpUserMeshAction(nullptr)
{
  if(!CheckingReactionRadius(fpMesh->GetResolution()))
  {
    G4String WarMessage = "resolution is not good : " +
                          std::to_string(fpMesh->GetResolution() / nm);
    G4Exception("G4DNAEventScheduler::InitializeInMesh()", "WrongResolution",
                JustWarning, WarMessage);
  }
}

References CheckingReactionRadius(), fpMesh, G4Exception(), JustWarning, and nm.

◆ ~G4DNAEventScheduler()

G4DNAEventScheduler::~G4DNAEventScheduler ( )

overridedefault

◆ G4DNAEventScheduler() [2/2]

G4DNAEventScheduler::G4DNAEventScheduler ( const G4DNAEventScheduler & )

delete

Member Function Documentation

◆ AddTimeToRecord()

void G4DNAEventScheduler::AddTimeToRecord ( const G4double & time )

Definition at line 125 of file G4DNAEventScheduler.cc.

{
  if(fTimeToRecord.find(time) == fTimeToRecord.end())
  {
    fTimeToRecord.insert(time);
  }
}

References fTimeToRecord.

◆ CheckingReactionRadius()

G4bool G4DNAEventScheduler::CheckingReactionRadius ( G4double resolution )

static

Definition at line 660 of file G4DNAEventScheduler.cc.

{
  auto pMolecularReactionTable = G4DNAMolecularReactionTable::Instance();
  auto reactionDataList = pMolecularReactionTable->GetVectorOfReactionData();
  if(reactionDataList.empty())
  {
    G4cout << "reactionDataList.empty()" << G4endl;
    return true;
  }
  else
  {
    for(auto it : reactionDataList)
    {
      if(it->GetEffectiveReactionRadius() >= resolution / CLHEP::pi)
      {
        G4cout << it->GetReactant1()->GetName() << " + "
               << it->GetReactant2()->GetName() << G4endl;
        G4cout << "G4DNAEventScheduler::ReactionRadius : "
               << G4BestUnit(it->GetEffectiveReactionRadius(), "Length")
               << G4endl;
        G4cout << "resolution : " << G4BestUnit(resolution, "Length") << G4endl;
        return false;
      }
    }
    return true;
  }
}

References G4BestUnit, G4cout, G4endl, G4DNAMolecularReactionTable::Instance(), and CLHEP::pi.

Referenced by G4DNAEventScheduler().

◆ ClearAndReChargeCounter()

void G4DNAEventScheduler::ClearAndReChargeCounter ( )

Definition at line 74 of file G4DNAEventScheduler.cc.

{
  fCounterMap.clear();
  if(fTimeToRecord.empty())
  {
    G4cout << "fTimeToRecord is empty " << G4endl;
  }
  fLastRecoredTime = fTimeToRecord.begin();
 
  if(G4VMoleculeCounter::Instance()->InUse())  // copy from MoleculeCounter
  {
    G4MoleculeCounter::RecordedMolecules species;
    species = G4MoleculeCounter::Instance()->GetRecordedMolecules();
    if(species.get() == nullptr)
    {
      return;
    }
    else if(species->empty())
    {
      G4MoleculeCounter::Instance()->ResetCounter();
      return;
    }
    for(auto time_mol : fTimeToRecord)
    {
      if(time_mol > fStartTime)
      {
        continue;
      }
 
      for(auto molecule : *species)
      {
        G4int n_mol = G4MoleculeCounter::Instance()->GetNMoleculesAtTime(
          molecule, time_mol);
 
        if(n_mol < 0)
        {
          G4cerr << "G4DNAEventScheduler::ClearAndReChargeCounter() ::N "
                    "molecules not valid < 0 "
                 << G4endl;
          G4Exception("", "N<0", FatalException, "");
        }
        fCounterMap[time_mol][molecule] = n_mol;
      }
 
      fLastRecoredTime++;
    }
    G4MoleculeCounter::Instance()->ResetCounter();  // reset
    G4MoleculeCounter::Instance()->Use(false);      // no more used
  }
}

References FatalException, fCounterMap, fLastRecoredTime, fStartTime, fTimeToRecord, G4cerr, G4cout, G4endl, G4Exception(), G4MoleculeCounter::GetNMoleculesAtTime(), G4MoleculeCounter::GetRecordedMolecules(), G4MoleculeCounter::Instance(), G4VMoleculeCounter::Instance(), G4MoleculeCounter::ResetCounter(), and G4VMoleculeCounter::Use().

Referenced by Initialize().

◆ GetCounterMap()

std::map< G4double, G4DNAEventScheduler::MapCounter > G4DNAEventScheduler::GetCounterMap ( ) const

Definition at line 645 of file G4DNAEventScheduler.cc.

{
  return fCounterMap;
}

References fCounterMap.

◆ GetEndTime()

G4double G4DNAEventScheduler::GetEndTime ( ) const

Definition at line 320 of file G4DNAEventScheduler.cc.

320{ return fEndTime; }

References fEndTime.

◆ GetMesh()

G4DNAMesh * G4DNAEventScheduler::GetMesh ( ) const

Definition at line 656 of file G4DNAEventScheduler.cc.

656{ return fpMesh.get(); }

References fpMesh.

◆ GetPixels()

G4int G4DNAEventScheduler::GetPixels ( ) const

Definition at line 658 of file G4DNAEventScheduler.cc.

658{ return fPixel; }

References fPixel.

◆ GetStartTime()

G4double G4DNAEventScheduler::GetStartTime ( ) const

Definition at line 318 of file G4DNAEventScheduler.cc.

318{ return fStartTime; }

References fStartTime.

◆ GetTimeStep()

G4double G4DNAEventScheduler::GetTimeStep ( ) const

Definition at line 322 of file G4DNAEventScheduler.cc.

322{ return fTimeStep; }

References fTimeStep.

◆ GetVerbose()

G4int G4DNAEventScheduler::GetVerbose ( ) const

inline

Definition at line 324 of file G4DNAEventScheduler.cc.

324{ return fVerbose; }

References fVerbose.

◆ Initialize()

void G4DNAEventScheduler::Initialize ( )

Definition at line 242 of file G4DNAEventScheduler.cc.

{
  if(!fInitialized)
  {
    fPixel = fInitialPixels;
    fpMesh = std::make_unique<G4DNAMesh>(fpMesh->GetBoundingBox(), fPixel);
 
    // Scavenger();
 
    auto pScavengerMaterial = dynamic_cast<G4DNAScavengerMaterial*>(
      G4Scheduler::Instance()->GetScavengerMaterial());
    if(pScavengerMaterial == nullptr)
    {
      G4cout << "pScavengerMaterial == nullptr" << G4endl;
    }
    else
    {
      if(fVerbose > 1)
      {
        pScavengerMaterial->PrintInfo();
      }
    }
 
    Voxelizing();
    fpGillespieReaction->SetVoxelMesh(*fpMesh);
    fpGillespieReaction->SetEventSet(fpEventSet.get());
    fpGillespieReaction->SetTimeStep(
      0);  // reset fTimeStep = 0 in fpGillespieReaction
    fpGillespieReaction->Initialize();
    fpUpdateSystem->SetMesh(fpMesh.get());
    ClearAndReChargeCounter();
    fInitialized = true;
  }
 
  if(fVerbose > 0)
  {
    fpUpdateSystem->SetVerbose(1);
  }
 
  if(fVerbose > 2)
  {
    fpMesh->PrintMesh();
  }
}

References ClearAndReChargeCounter(), fInitialized, fInitialPixels, fpEventSet, fpGillespieReaction, fPixel, fpMesh, fpUpdateSystem, fVerbose, G4cout, G4endl, G4Scheduler::GetScavengerMaterial(), G4Scheduler::Instance(), and Voxelizing().

◆ InitializeInMesh()

void G4DNAEventScheduler::InitializeInMesh ( )

Definition at line 286 of file G4DNAEventScheduler.cc.

{
  if(fPixel <= 1)
  {
    fRunning = false;
    return;
  }
  // TEST /3
  ReVoxelizing(fPixel / 2);  //
  // ReVoxelizing(fPixel/3);//
 
  fpGillespieReaction->SetVoxelMesh(*fpMesh);
  fpUpdateSystem->SetMesh(fpMesh.get());
  fpGillespieReaction->Initialize();
}

References fpGillespieReaction, fPixel, fpMesh, fpUpdateSystem, fRunning, and ReVoxelizing().

Referenced by RunInMesh().

◆ operator=()

G4DNAEventScheduler & G4DNAEventScheduler::operator= ( const G4DNAEventScheduler & right )

delete

◆ PrintRecordTime()

void G4DNAEventScheduler::PrintRecordTime ( )

Definition at line 621 of file G4DNAEventScheduler.cc.

{
  G4cout << "fCounterMap.size : " << fCounterMap.size() << G4endl;
 
  for(const auto& i : fCounterMap)
  {
    auto map   = i.second;
    auto begin = map.begin();  //
    auto end   = map.end();    //
    for(; begin != end; begin++)
    {
      auto molecule = begin->first;
      auto number   = begin->second;
      if(number == 0)
      {
        continue;
      }
      G4cout << "molecule : " << molecule->GetName() << " number : " << number
             << G4endl;
    }
  }
}

References fCounterMap, G4cout, G4endl, and anonymous_namespace{G4QuasiElRatios.cc}::map.

Referenced by RecordTime().

◆ RecordTime()

void G4DNAEventScheduler::RecordTime ( )

Definition at line 577 of file G4DNAEventScheduler.cc.

{
  auto recordTime = *fLastRecoredTime;
  if(fGlobalTime >= recordTime && fCounterMap[recordTime].empty())
  {
    auto begin = fpMesh->begin();
    auto end   = fpMesh->end();
    for(; begin != end; begin++)
    {
      auto node = begin->second;
      if(node == nullptr) {
        continue;
      }
      for(const auto& it : node->GetMapList())
      {
        fCounterMap[recordTime][it.first] += it.second;
      }
    }
    fLastRecoredTime++;
 
#ifdef DEBUG
    PrintRecordTime();
    G4MoleculeTable* pMoleculeTable = G4MoleculeTable::Instance();
    auto iter = pMoleculeTable->GetConfigurationIterator();
    iter.reset();
    while(iter())
    {
      auto conf = iter.value();
 
      G4cout << "GlobalTime : " << G4BestUnit(fGlobalTime, "Time")
             << "  recordTime : " << G4BestUnit(recordTime, "Time") << "  "
             << conf->GetName()
             << "  number : " << fCounterMap[recordTime][conf]
             << "  MoleculeCounter : "
             << G4MoleculeCounter::Instance()->GetCurrentNumberOf(conf)
             << G4endl;
 
      assert(G4MoleculeCounter::Instance()->GetCurrentNumberOf(conf) ==
             fCounterMap[recordTime][conf]);
    }
#endif
  }
}

References anonymous_namespace{G4MTcoutDestination.cc}::empty, fCounterMap, fGlobalTime, fLastRecoredTime, fpMesh, G4BestUnit, G4cout, G4endl, G4MoleculeTable::GetConfigurationIterator(), G4MoleculeCounter::Instance(), G4MoleculeTable::Instance(), PrintRecordTime(), and G4MoleculeIterator< MOLECULE >::reset().

Referenced by Stepping().

◆ Reset()

void G4DNAEventScheduler::Reset ( )

Definition at line 216 of file G4DNAEventScheduler.cc.

{
  // find another solultion
  fGlobalTime = fEndTime;
 
  //
  // RecordTime();//Last register for counter
 
  if(fVerbose > 0)
  {
    G4cout << "End Processing and reset Gird, ScavengerTable, EventSet for new "
              "simulation!!!!"
           << G4endl;
  }
  fInitialized    = false;
  fTimeStep       = 0;
  fStepNumber     = 0;
  fGlobalTime     = fStartTime;
  fRunning        = true;
  fReactionNumber = 0;
  fJumpingNumber  = 0;
 
  fpEventSet->RemoveEventSet();
  fpMesh->Reset();
}

References fEndTime, fGlobalTime, fInitialized, fJumpingNumber, fpEventSet, fpMesh, fReactionNumber, fRunning, fStartTime, fStepNumber, fTimeStep, fVerbose, G4cout, and G4endl.

Referenced by Run().

◆ ResetInMesh()

void G4DNAEventScheduler::ResetInMesh ( )

Definition at line 302 of file G4DNAEventScheduler.cc.

{
  if(fVerbose > 0)
  {
    G4cout
      << "*** End Processing In Mesh and reset Mesh, EventSet for new Mesh!!!!"
      << G4endl;
  }
  fpEventSet->RemoveEventSet();
  fInitialized      = false;
  fIsChangeMesh     = false;
  fReactionNumber   = 0;
  fJumpingNumber    = 0;
  fStepNumberInMesh = 0;
}

References fInitialized, fIsChangeMesh, fJumpingNumber, fpEventSet, fReactionNumber, fStepNumberInMesh, fVerbose, G4cout, and G4endl.

Referenced by RunInMesh().

◆ ReVoxelizing()

void G4DNAEventScheduler::ReVoxelizing ( G4int pixel )

Definition at line 176 of file G4DNAEventScheduler.cc.

{
  fPixel       = pixel;
  auto newMesh = new G4DNAMesh(fpMesh->GetBoundingBox(), fPixel);
 
  auto begin = fpMesh->begin();
  auto end   = fpMesh->end();
  std::map<G4DNAMesh::Key, MapList> TrackKeyMap;
  for(; begin != end; begin++)
  {
    auto index  = fpMesh->GetIndex(begin->first);
    auto newKey = newMesh->GetKey(fpMesh->GetIndex(index, fPixel));
    auto node   = begin->second;
    // if (node == nullptr) continue;
    if(TrackKeyMap.find(newKey) == TrackKeyMap.end())
    {
      TrackKeyMap[newKey] = node->GetMapList();
    }
    else
    {
      for(const auto& it : node->GetMapList())
      {
        TrackKeyMap[newKey][it.first] += it.second;
      }
      if(fVerbose > 1)
      {
        G4cout << "key : " << begin->first << " index : " << index
               << " new index : " << fpMesh->GetIndex(index, fPixel)
               << " new key : " << newKey
               << " number: " << node->GetMapList().begin()->second << G4endl;
      }
    }
  }
  fpMesh.reset(newMesh);
 
  for(auto& it : TrackKeyMap)
  {
    fpMesh->SetVoxelMapList(it.first, std::move(it.second));
  }
}

References fPixel, fpMesh, fVerbose, G4cout, and G4endl.

Referenced by InitializeInMesh().

◆ Run()

void G4DNAEventScheduler::Run ( )

Definition at line 334 of file G4DNAEventScheduler.cc.

{
  G4Timer localtimer;
  if(fVerbose > 0)
  {
    localtimer.Start();
    G4cout << "***G4DNAEventScheduler::Run*** for Pixel : " << fPixel << G4endl;
  }
  while(fEndTime > fGlobalTime && fRunning)
  {
    RunInMesh();
  }
  if(fVerbose > 0)
  {
    if(!fRunning)
    {
      G4cout << " StepNumber(" << fStepNumber << ") = MaxStep(" << fMaxStep
             << ")" << G4endl;
    }
    else if(fEndTime <= fGlobalTime)
    {
      G4cout << " GlobalTime(" << fGlobalTime << ") > EndTime(" << fEndTime
             << ")"
             << " StepNumber : " << fStepNumber << G4endl;
    }
    localtimer.Stop();
    G4cout << "***G4DNAEventScheduler::Ending::"
           << G4BestUnit(fGlobalTime, "Time")
           << " Events left : " << fpEventSet->size() << G4endl;
    if(fVerbose > 1) {
      fpMesh->PrintMesh();
    }
    G4cout << " Computing Time : " << localtimer << G4endl;
  }
  Reset();
}

References fEndTime, fGlobalTime, fMaxStep, fpEventSet, fPixel, fpMesh, fRunning, fStepNumber, fVerbose, G4BestUnit, G4cout, G4endl, Reset(), RunInMesh(), G4Timer::Start(), and G4Timer::Stop().

◆ RunInMesh()

void G4DNAEventScheduler::RunInMesh ( )

Definition at line 371 of file G4DNAEventScheduler.cc.

{
  if(!fInitialized)
  {
    InitializeInMesh();
  }
  G4Timer localtimerInMesh;
  // if (fVerbose > 0)
  {
    localtimerInMesh.Start();
    G4double C = 20;
    G4double D = G4MoleculeTable::Instance()
                 ->GetConfiguration("H2O2")
                 ->GetDiffusionCoefficient();
    G4double transferTime = std::pow(fpMesh->GetResolution(), 2) * C / (6 * D);
    G4cout << "***G4DNAEventScheduler::RunInMesh*** for Pixel : " << fPixel
           << "  transferTime : " << G4BestUnit(transferTime, "Time") << G4endl;
    G4cout << "  resolution : " << G4BestUnit(fpMesh->GetResolution(), "Length")
           << G4endl;
  }
 
  if(fVerbose > 2)
  {
    fpMesh->PrintMesh();
  }
 
  if(fpUserMeshAction != nullptr)
  {
    fpUserMeshAction->BeginOfMesh(fpMesh.get(), fGlobalTime);
  }
 
  // if diffusive jumping is avaiable, EventSet is never empty
  while(!fpEventSet->Empty() && !fIsChangeMesh && fEndTime > fGlobalTime)
  {
    Stepping();
    fGlobalTime = fTimeStep + fStartTime;
 
    if(fpUserMeshAction != nullptr)
    {
      fpUserMeshAction->InMesh(fpMesh.get(), fGlobalTime);
    }
 
    if(fVerbose > 2)
    {
      G4cout << "fGlobalTime : " << G4BestUnit(fGlobalTime, "Time")
             << " fTimeStep : " << G4BestUnit(fTimeStep, "Time") << G4endl;
    }
 
    G4double C = 20;
    G4double D = G4MoleculeTable::Instance()
                 ->GetConfiguration("H2O2")
                 ->GetDiffusionCoefficient();
    if(D == 0)
    {
      G4ExceptionDescription exceptionDescription;
      exceptionDescription << "D == 0";
      G4Exception("G4DNAEventScheduler::RunInMesh", "G4DNAEventScheduler001",
                  FatalErrorInArgument, exceptionDescription);
    }
    G4double transferTime = std::pow(fpMesh->GetResolution(), 2) * C / (6 * D);
 
    // if(fStepNumberInMesh > 40000 && fPixel != 1)
    if(transferTime < fTimeStep &&
       fPixel != 1)  // dont change Mesj if fPixel == 1
    {
      if(fVerbose > 1)
      {
        G4cout << " Pixels : " << fPixel << "  resolution : "
               << G4BestUnit(fpMesh->GetResolution(), "Length")
               << "  fStepNumberInMesh : " << fStepNumberInMesh
               << " at fGlobalTime : " << G4BestUnit(fGlobalTime, "Time")
               << " at fTimeStep : " << G4BestUnit(fTimeStep, "Time")
               << "  fReactionNumber : " << fReactionNumber
               << " transferTime : " << G4BestUnit(transferTime, "Time")
               << G4endl;
      }
      if(fSetChangeMesh)
      {
        fIsChangeMesh = true;
      }
    }
  }
 
  if(fVerbose > 1)
  {
    localtimerInMesh.Stop();
    G4cout << "***G4DNAEventScheduler::Ending::"
           << G4BestUnit(fGlobalTime, "Time")
           << " Event left : " << fpEventSet->size() << G4endl;
    G4cout << " Computing Time : " << localtimerInMesh << " Due to : ";
    if(fpEventSet->Empty())
    {
      G4cout << "EventSet is Empty" << G4endl;
    }
    else if(fIsChangeMesh)
    {
      G4cout << "Changing Mesh from : " << fPixel
             << " pixels to : " << fPixel / 2 << " pixels" << G4endl;
      G4cout << "Info : ReactionNumber : " << fReactionNumber
             << "   JumpingNumber : " << fJumpingNumber << G4endl;
    }
    else if(fEndTime > fGlobalTime)
    {
      G4cout << " GlobalTime(" << fGlobalTime << ") > EndTime(" << fEndTime
             << ")"
             << " StepNumber : " << fStepNumber << G4endl;
    }
    if(fVerbose > 2)
    {
      fpMesh->PrintMesh();
    }
    G4cout << G4endl;
  }
 
  if(fpUserMeshAction != nullptr)
  {
    fpUserMeshAction->EndOfMesh(fpMesh.get(), fGlobalTime);
  }
  ResetInMesh();
}

References C(), D(), FatalErrorInArgument, fEndTime, fGlobalTime, fInitialized, fIsChangeMesh, fJumpingNumber, fpEventSet, fPixel, fpMesh, fpUserMeshAction, fReactionNumber, fSetChangeMesh, fStartTime, fStepNumber, fStepNumberInMesh, fTimeStep, fVerbose, G4BestUnit, G4cout, G4endl, G4Exception(), G4MoleculeTable::GetConfiguration(), G4MolecularConfiguration::GetDiffusionCoefficient(), InitializeInMesh(), G4MoleculeTable::Instance(), ResetInMesh(), G4Timer::Start(), Stepping(), and G4Timer::Stop().

Referenced by Run().

◆ SetChangeMesh()

void G4DNAEventScheduler::SetChangeMesh ( G4bool change )

inline

Definition at line 75 of file G4DNAEventScheduler.hh.

75{ fSetChangeMesh = change; }

References fSetChangeMesh.

◆ SetEndTime()

void G4DNAEventScheduler::SetEndTime ( const G4double & endTime )

Definition at line 572 of file G4DNAEventScheduler.cc.

{
  fEndTime = endTime;
}

References fEndTime.

◆ SetMaxNbSteps()

void G4DNAEventScheduler::SetMaxNbSteps ( G4int max )

Definition at line 326 of file G4DNAEventScheduler.cc.

326{ fMaxStep = max; }

G4INCL::Math::max

T max(const T t1, const T t2)

brief Return the largest of the two arguments

Definition: G4INCLGlobals.hh:112

References fMaxStep, and G4INCL::Math::max().

◆ SetStartTime()

void G4DNAEventScheduler::SetStartTime ( G4double time )

Definition at line 328 of file G4DNAEventScheduler.cc.

{
  fStartTime = time;
}

References fStartTime.

◆ SetUserMeshAction()

void G4DNAEventScheduler::SetUserMeshAction ( std::unique_ptr< G4UserMeshAction > pUserMeshAction )

Definition at line 650 of file G4DNAEventScheduler.cc.

{
  fpUserMeshAction = std::move(pUserMeshAction);
}

References fpUserMeshAction.

◆ SetVerbose()

void G4DNAEventScheduler::SetVerbose ( G4int verbose )

inline

Definition at line 72 of file G4DNAEventScheduler.hh.

72{ fVerbose = verbose; }

References fVerbose.

◆ Stepping()

void G4DNAEventScheduler::Stepping ( )

Definition at line 492 of file G4DNAEventScheduler.cc.

{
  fStepNumber < fMaxStep ? fStepNumber++ : fRunning = false;
 
  if(fpEventSet->size() > fpMesh->size())
  {
    G4ExceptionDescription exceptionDescription;
    exceptionDescription << "fpEventSet->size() > fpMesh->size()";
    G4Exception("G4DNAEventScheduler::Stepping", "G4DNAEventScheduler002",
                FatalErrorInArgument, exceptionDescription);
  };
 
  auto selected   = fpEventSet->begin();
  const auto& key = (*selected)->GetKey();
  auto index      = fpMesh->GetIndex(key);
 
  if(fVerbose > 1)
  {
    G4cout << "G4DNAEventScheduler::Stepping()*********************************"
              "*******"
           << G4endl;
    (*selected)->PrintEvent();
  }
 
  // get selected time step
  fTimeStep = (*selected)->GetTime();
 
  // selected data
  auto pJumping  = (*selected)->GetJumpingData();
  auto pReaction = (*selected)->GetReactionData();
 
  fpUpdateSystem->SetGlobalTime(fTimeStep +
                                fStartTime);  // this is just for printing
 
  if(pJumping == nullptr && pReaction == nullptr)
  {
    G4ExceptionDescription exceptionDescription;
    exceptionDescription << "pJumping == nullptr && pReaction == nullptr";
    G4Exception("G4DNAEventScheduler::Stepping", "G4DNAEventScheduler003",
                FatalErrorInArgument, exceptionDescription);
  }
 
  fpGillespieReaction->SetTimeStep(fTimeStep);
 
  if(pJumping == nullptr)
  {
    fpUpdateSystem->UpdateSystem(index, *pReaction);
 
    fpEventSet->RemoveEvent(selected);
    // create new event
    fpGillespieReaction->CreateEvent(key);
    fReactionNumber++;
 
    // recordTime in reaction
    RecordTime();
  }
  else if(pReaction == nullptr)
  {
    // dont change this
    fpUpdateSystem->UpdateSystem(index, *pJumping);
    auto jumpingKey = fpMesh->GetKey(pJumping->second);
    fpEventSet->RemoveEvent(selected);
 
    // create new event
    // should create Jumping before key
    fpGillespieReaction->CreateEvent(jumpingKey);
    fpGillespieReaction->CreateEvent(key);
 
    fJumpingNumber++;
  }
  if(fVerbose > 1)
  {
    G4cout << "G4DNAEventScheduler::Stepping::end "
              "Print***********************************"
           << G4endl;
    G4cout << G4endl;
  }
  fStepNumberInMesh++;
}

References FatalErrorInArgument, fJumpingNumber, fMaxStep, fpEventSet, fpGillespieReaction, fpMesh, fpUpdateSystem, fReactionNumber, fRunning, fStartTime, fStepNumber, fStepNumberInMesh, fTimeStep, fVerbose, G4cout, G4endl, G4Exception(), and RecordTime().

Referenced by RunInMesh().

◆ Stop()

void G4DNAEventScheduler::Stop ( )

Definition at line 333 of file G4DNAEventScheduler.cc.

333{ fRunning = false; }

References fRunning.

◆ Voxelizing()

void G4DNAEventScheduler::Voxelizing ( )

Definition at line 135 of file G4DNAEventScheduler.cc.

{
  auto pMainList = G4ITTrackHolder::Instance()->GetMainList();
  std::map<G4DNAMesh::Key, MapList> TrackKeyMap;
  for(auto track : *pMainList)
  {
    auto molType = GetMolecule(track)->GetMolecularConfiguration();
 
    auto pScavengerMaterial = dynamic_cast<G4DNAScavengerMaterial*>(
      G4Scheduler::Instance()->GetScavengerMaterial());
    if(pScavengerMaterial != nullptr &&
       pScavengerMaterial->find(molType))  // avoid voxelize the scavenger
    {
      continue;
    }
 
    auto key = fpMesh->GetKey(track->GetPosition());
    if(TrackKeyMap.find(key) != TrackKeyMap.end())
    {
      std::map<MolType, size_t>& TrackTypeMap = TrackKeyMap[key];
      if(TrackTypeMap.find(molType) != TrackTypeMap.end())
      {
        TrackTypeMap[molType]++;
      }
      else
      {
        TrackTypeMap[molType] = 1;
      }
    }
    else
    {
      TrackKeyMap[key][molType] = 1;
    }
  }
 
  for(auto& it : TrackKeyMap)
  {
    fpMesh->SetVoxelMapList(it.first, std::move(it.second));
  }
}

References fpMesh, G4ITTrackHolder::GetMainList(), G4Molecule::GetMolecularConfiguration(), GetMolecule(), G4Scheduler::GetScavengerMaterial(), G4ITTrackHolder::Instance(), and G4Scheduler::Instance().

Referenced by Initialize().

Field Documentation

◆ fCounterMap

std::map<G4double , MapCounter> G4DNAEventScheduler::fCounterMap

private

Definition at line 120 of file G4DNAEventScheduler.hh.

Referenced by ClearAndReChargeCounter(), GetCounterMap(), PrintRecordTime(), and RecordTime().

◆ fEndTime

G4double G4DNAEventScheduler::fEndTime

private

Definition at line 99 of file G4DNAEventScheduler.hh.

Referenced by GetEndTime(), Reset(), Run(), RunInMesh(), and SetEndTime().

◆ fGlobalTime

G4double G4DNAEventScheduler::fGlobalTime

private

Definition at line 104 of file G4DNAEventScheduler.hh.

Referenced by RecordTime(), Reset(), Run(), and RunInMesh().

◆ fInitialized

G4bool G4DNAEventScheduler::fInitialized

private

Definition at line 97 of file G4DNAEventScheduler.hh.

Referenced by Initialize(), Reset(), ResetInMesh(), and RunInMesh().

◆ fInitialPixels

G4double G4DNAEventScheduler::fInitialPixels

private

Definition at line 111 of file G4DNAEventScheduler.hh.

Referenced by Initialize().

◆ fIsChangeMesh

G4bool G4DNAEventScheduler::fIsChangeMesh

private

Definition at line 108 of file G4DNAEventScheduler.hh.

Referenced by ResetInMesh(), and RunInMesh().

◆ fJumpingNumber

G4double G4DNAEventScheduler::fJumpingNumber

private

Definition at line 105 of file G4DNAEventScheduler.hh.

Referenced by Reset(), ResetInMesh(), RunInMesh(), and Stepping().

◆ fLastRecoredTime

std::set<G4double>::iterator G4DNAEventScheduler::fLastRecoredTime

private

Definition at line 122 of file G4DNAEventScheduler.hh.

Referenced by ClearAndReChargeCounter(), and RecordTime().

◆ fMaxStep

G4int G4DNAEventScheduler::fMaxStep

private

Definition at line 101 of file G4DNAEventScheduler.hh.

Referenced by Run(), SetMaxNbSteps(), and Stepping().

◆ fpEventSet

std::unique_ptr<G4DNAEventSet> G4DNAEventScheduler::fpEventSet

private

Definition at line 115 of file G4DNAEventScheduler.hh.

Referenced by Initialize(), Reset(), ResetInMesh(), Run(), RunInMesh(), and Stepping().

◆ fpGillespieReaction

std::unique_ptr<G4DNAGillespieDirectMethod> G4DNAEventScheduler::fpGillespieReaction

private

Definition at line 114 of file G4DNAEventScheduler.hh.

Referenced by Initialize(), InitializeInMesh(), and Stepping().

◆ fPixel

G4int G4DNAEventScheduler::fPixel

private

Definition at line 107 of file G4DNAEventScheduler.hh.

Referenced by GetPixels(), Initialize(), InitializeInMesh(), ReVoxelizing(), Run(), and RunInMesh().

◆ fpMesh

std::unique_ptr<G4DNAMesh> G4DNAEventScheduler::fpMesh

private

Definition at line 113 of file G4DNAEventScheduler.hh.

Referenced by G4DNAEventScheduler(), GetMesh(), Initialize(), InitializeInMesh(), RecordTime(), Reset(), ReVoxelizing(), Run(), RunInMesh(), Stepping(), and Voxelizing().

◆ fpUpdateSystem

std::unique_ptr<G4DNAUpdateSystemModel> G4DNAEventScheduler::fpUpdateSystem

private

Definition at line 116 of file G4DNAEventScheduler.hh.

Referenced by Initialize(), InitializeInMesh(), and Stepping().

◆ fpUserMeshAction

std::unique_ptr<G4UserMeshAction> G4DNAEventScheduler::fpUserMeshAction

private

Definition at line 117 of file G4DNAEventScheduler.hh.

Referenced by RunInMesh(), and SetUserMeshAction().

◆ fReactionNumber

G4double G4DNAEventScheduler::fReactionNumber

private

Definition at line 106 of file G4DNAEventScheduler.hh.

Referenced by Reset(), ResetInMesh(), RunInMesh(), and Stepping().

◆ fRunning

G4bool G4DNAEventScheduler::fRunning

private

Definition at line 102 of file G4DNAEventScheduler.hh.

Referenced by InitializeInMesh(), Reset(), Run(), Stepping(), and Stop().

◆ fSetChangeMesh

G4bool G4DNAEventScheduler::fSetChangeMesh

private

Definition at line 109 of file G4DNAEventScheduler.hh.

Referenced by RunInMesh(), and SetChangeMesh().

◆ fStartTime

G4double G4DNAEventScheduler::fStartTime

private

Definition at line 98 of file G4DNAEventScheduler.hh.

Referenced by ClearAndReChargeCounter(), GetStartTime(), Reset(), RunInMesh(), SetStartTime(), and Stepping().

◆ fStepNumber

G4int G4DNAEventScheduler::fStepNumber

private

Definition at line 100 of file G4DNAEventScheduler.hh.

Referenced by Reset(), Run(), RunInMesh(), and Stepping().

◆ fStepNumberInMesh

G4int G4DNAEventScheduler::fStepNumberInMesh

private

Definition at line 110 of file G4DNAEventScheduler.hh.

Referenced by ResetInMesh(), RunInMesh(), and Stepping().

◆ fTimeStep

G4double G4DNAEventScheduler::fTimeStep

private

Definition at line 103 of file G4DNAEventScheduler.hh.

Referenced by GetTimeStep(), Reset(), RunInMesh(), and Stepping().

◆ fTimeToRecord

std::set<G4double> G4DNAEventScheduler::fTimeToRecord

private

Definition at line 121 of file G4DNAEventScheduler.hh.

Referenced by AddTimeToRecord(), and ClearAndReChargeCounter().

◆ fVerbose

G4int G4DNAEventScheduler::fVerbose

private

Definition at line 96 of file G4DNAEventScheduler.hh.

Referenced by GetVerbose(), Initialize(), Reset(), ResetInMesh(), ReVoxelizing(), Run(), RunInMesh(), SetVerbose(), and Stepping().

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

source/processes/electromagnetic/dna/models/include/G4DNAEventScheduler.hh
source/processes/electromagnetic/dna/models/src/G4DNAEventScheduler.cc

Public Types

Public Member Functions

Static Public Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ MapCounter

◆ MapList

◆ MolType

Constructor & Destructor Documentation

◆ G4DNAEventScheduler() [1/2]

◆ ~G4DNAEventScheduler()

◆ G4DNAEventScheduler() [2/2]

Member Function Documentation

◆ AddTimeToRecord()

◆ CheckingReactionRadius()

◆ ClearAndReChargeCounter()

◆ GetCounterMap()

◆ GetEndTime()

◆ GetMesh()

◆ GetPixels()

◆ GetStartTime()

◆ GetTimeStep()

◆ GetVerbose()

◆ Initialize()

◆ InitializeInMesh()

◆ operator=()

◆ PrintRecordTime()

◆ RecordTime()

◆ Reset()

◆ ResetInMesh()

◆ ReVoxelizing()

◆ Run()

◆ RunInMesh()

◆ SetChangeMesh()

◆ SetEndTime()

◆ SetMaxNbSteps()

◆ SetStartTime()

◆ SetUserMeshAction()

◆ SetVerbose()

◆ Stepping()

◆ Stop()

◆ Voxelizing()

Field Documentation

◆ fCounterMap

◆ fEndTime

◆ fGlobalTime

◆ fInitialized

◆ fInitialPixels

◆ fIsChangeMesh

◆ fJumpingNumber

◆ fLastRecoredTime

◆ fMaxStep

◆ fpEventSet

◆ fpGillespieReaction

◆ fPixel

◆ fpMesh

◆ fpUpdateSystem

◆ fpUserMeshAction

◆ fReactionNumber

◆ fRunning

◆ fSetChangeMesh

◆ fStartTime

◆ fStepNumber

◆ fStepNumberInMesh

◆ fTimeStep

◆ fTimeToRecord

◆ fVerbose