Run action class. More...

#include <Doxymodules_basic.h>

Inheritance diagram for B1::RunAction:

Public Member Functions
	RunAction ()

	~RunAction () override=default

void	BeginOfRunAction (const G4Run *) override

void	EndOfRunAction (const G4Run *) override

void	AddEdep (G4double edep)

Private Attributes
G4Accumulable< G4double >	fEdep = 0.

G4Accumulable< G4double >	fEdep2 = 0.

Detailed Description

Run action class.

In EndOfRunAction(), it calculates the dose in the selected volume from the energy deposit accumulated via stepping and event actions. The computed dose is then printed on the screen.

Definition at line 24 of file Doxymodules_basic.h.

Constructor & Destructor Documentation

◆ RunAction()

B1::RunAction::RunAction ( )

Definition at line 48 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(fEdep);
  accumulableManager->RegisterAccumulable(fEdep2);
}

◆ ~RunAction()

B1::RunAction::~RunAction ( )

overridedefault

Member Function Documentation

◆ BeginOfRunAction()

void B1::RunAction::BeginOfRunAction ( const G4Run * )

override

Definition at line 70 of file RunAction.cc.

{
  // inform the runManager to save random number seed
  G4RunManager::GetRunManager()->SetRandomNumberStore(false);
 
  // reset accumulables to their initial values
  G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
  accumulableManager->Reset();
 
}

◆ EndOfRunAction()

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

override

Definition at line 83 of file RunAction.cc.

{
  G4int nofEvents = run->GetNumberOfEvent();
  if (nofEvents == 0) return;
 
  // Merge accumulables
  G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
  accumulableManager->Merge();
 
  // Compute dose = total energy deposit in a run and its variance
  //
  G4double edep  = fEdep.GetValue();
  G4double edep2 = fEdep2.GetValue();
 
  G4double rms = edep2 - edep*edep/nofEvents;
  if (rms > 0.) rms = std::sqrt(rms); else rms = 0.;
 
  const auto detConstruction = static_cast<const DetectorConstruction*>(
    G4RunManager::GetRunManager()->GetUserDetectorConstruction());
  G4double mass = detConstruction->GetScoringVolume()->GetMass();
  G4double dose = edep/mass;
  G4double rmsDose = rms/mass;
 
  // 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 runCondition;
  if (generatorAction)
  {
    const G4ParticleGun* particleGun = generatorAction->GetParticleGun();
    runCondition += particleGun->GetParticleDefinition()->GetParticleName();
    runCondition += " of ";
    G4double particleEnergy = particleGun->GetParticleEnergy();
    runCondition += G4BestUnit(particleEnergy,"Energy");
  }
 
  // Print
  //
  if (IsMaster()) {
    G4cout
     << G4endl
     << "--------------------End of Global Run-----------------------";
  }
  else {
    G4cout
     << G4endl
     << "--------------------End of Local Run------------------------";
  }
 
  G4cout
     << G4endl
     << " The run consists of " << nofEvents << " "<< runCondition
     << G4endl
     << " Cumulated dose per run, in scoring volume : "
     << G4BestUnit(dose,"Dose") << " rms = " << G4BestUnit(rmsDose,"Dose")
     << G4endl
     << "------------------------------------------------------------"
     << G4endl
     << G4endl;
}

◆ AddEdep()

void B1::RunAction::AddEdep ( G4double edep )

Definition at line 148 of file RunAction.cc.

{
  fEdep  += edep;
  fEdep2 += edep*edep;
}

Member Data Documentation

◆ fEdep

G4Accumulable<G4double> B1::RunAction::fEdep = 0.

private

Definition at line 60 of file RunAction.hh.

◆ fEdep2

G4Accumulable<G4double> B1::RunAction::fEdep2 = 0.