Run action class. More...

#include <Doxymodules_basic.h>

Inheritance diagram for B3a::RunAction:

Public Member Functions
	RunAction ()

	~RunAction () override=default

void	BeginOfRunAction (const G4Run *) override

void	EndOfRunAction (const G4Run *) override

void	CountEvent ()

void	SumDose (G4double dose)

Private Attributes
G4Accumulable< G4int >	fGoodEvents = 0

G4Accumulable< G4double >	fSumDose = 0.

Detailed Description

Run action class.

Definition at line 90 of file Doxymodules_basic.h.

Constructor & Destructor Documentation

◆ RunAction()

B3a::RunAction::RunAction ( )

Definition at line 46 of file RunAction.cc.

{
  //add new units for dose
  //
  const G4double milligray = 1.e-3*gray;
  const G4double microgray = 1.e-6*gray;
  const G4double nanogray  = 1.e-9*gray;
  const G4double picogray  = 1.e-12*gray;
 
  new G4UnitDefinition("milligray", "milliGy" , "Dose", milligray);
  new G4UnitDefinition("microgray", "microGy" , "Dose", microgray);
  new G4UnitDefinition("nanogray" , "nanoGy"  , "Dose", nanogray);
  new G4UnitDefinition("picogray" , "picoGy"  , "Dose", picogray);
 
  // Register accumulable to the accumulable manager
  G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
  accumulableManager->RegisterAccumulable(fGoodEvents);
  accumulableManager->RegisterAccumulable(fSumDose);
}

◆ ~RunAction()

B3a::RunAction::~RunAction ( )

overridedefault

Member Function Documentation

◆ BeginOfRunAction()

void B3a::RunAction::BeginOfRunAction ( const G4Run * run )

override

Definition at line 68 of file RunAction.cc.

{
  G4cout << "### Run " << run->GetRunID() << " start." << G4endl;
 
  // reset accumulables to their initial values
  G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
  accumulableManager->Reset();
 
  //inform the runManager to save random number seed
  G4RunManager::GetRunManager()->SetRandomNumberStore(false);
}

◆ EndOfRunAction()

void B3a::RunAction::EndOfRunAction ( const G4Run * run )

override

Definition at line 82 of file RunAction.cc.

{
  G4int nofEvents = run->GetNumberOfEvent();
  if (nofEvents == 0) return;
 
  // Merge accumulables
  G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
  accumulableManager->Merge();
 
  // Run conditions
  //  note: There is no primary generator action object for "master"
  //        run manager for multi-threaded mode.
  const auto generatorAction = static_cast<const PrimaryGeneratorAction*>(
    G4RunManager::GetRunManager()->GetUserPrimaryGeneratorAction());
  G4String partName;
  if (generatorAction)
  {
    G4ParticleDefinition* particle
      = generatorAction->GetParticleGun()->GetParticleDefinition();
    partName = particle->GetParticleName();
  }
 
  // Print results
  //
  if (IsMaster())
  {
    G4cout
     << G4endl
     << "--------------------End of Global Run-----------------------"
     << G4endl
     << "  The run was " << nofEvents << " events ";
  }
  else
  {
    G4cout
     << G4endl
     << "--------------------End of Local Run------------------------"
     << G4endl
     << "  The run was " << nofEvents << " "<< partName;
  }
  G4cout
     << "; Nb of 'good' e+ annihilations: " << fGoodEvents.GetValue()  << G4endl
     << " Total dose in patient : " << G4BestUnit(fSumDose.GetValue(),"Dose")
     << G4endl
     << "------------------------------------------------------------" << G4endl
     << G4endl;
}

◆ CountEvent()

void B3a::RunAction::CountEvent ( )

inline

Definition at line 51 of file RunAction.hh.

51{ fGoodEvents += 1; };

◆ SumDose()

void B3a::RunAction::SumDose ( G4double dose )

inline

Definition at line 52 of file RunAction.hh.

52{ fSumDose += dose; };

Member Data Documentation

◆ fGoodEvents

G4Accumulable<G4int> B3a::RunAction::fGoodEvents = 0

private

Definition at line 55 of file RunAction.hh.

◆ fSumDose

G4Accumulable<G4double> B3a::RunAction::fSumDose = 0.