B1RunAction Class Reference

#include <B1RunAction.hh>

Inheritance diagram for B1RunAction:

Collaboration diagram for B1RunAction:

Public Member Functions
	B1RunAction ()
virtual	~B1RunAction ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
void	AddEdep (G4double edep)
	B1RunAction ()
virtual	~B1RunAction ()
virtual G4Run *	GenerateRun ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
Public Member Functions inherited from G4UserRunAction
	G4UserRunAction ()
virtual	~G4UserRunAction ()
virtual void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

Additional Inherited Members
Protected Attributes inherited from G4UserRunAction
G4bool	isMaster

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 46 of file B1RunAction.hh.

Constructor & Destructor Documentation

B1RunAction::B1RunAction

(

)

Definition at line 46 of file B1RunAction.cc.

: G4UserRunAction(),
  fEdep(0.),
  fEdep2(0.)
{ 
  // 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); 
}

Here is the call graph for this function:

B1RunAction::~B1RunAction ( )

virtual

Definition at line 71 of file B1RunAction.cc.

{}

B1RunAction::B1RunAction

(

)

virtual B1RunAction::~B1RunAction ( )

virtual

Member Function Documentation

void B1RunAction::AddEdep

(

G4double

edep

)

Definition at line 156 of file B1RunAction.cc.

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

Here is the caller graph for this function:

virtual void B1RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

void B1RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 76 of file B1RunAction.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();

}

Here is the call graph for this function:

virtual void B1RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

void B1RunAction::EndOfRunAction ( const G4Run *  run )

virtual

Reimplemented from G4UserRunAction.

Definition at line 89 of file B1RunAction.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 B1DetectorConstruction* detectorConstruction
   = static_cast<const B1DetectorConstruction*>
     (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
  G4double mass = detectorConstruction->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 B1PrimaryGeneratorAction* generatorAction
   = static_cast<const B1PrimaryGeneratorAction*>
     (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;
}

Here is the call graph for this function:

G4Run * B1RunAction::GenerateRun ( void  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 67 of file B1RunAction.cc.

{
  return new B1Run; 
}

---

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

• source/geant4.10.03.p03/examples/basic/B1/include/B1RunAction.hh
• source/geant4.10.03.p03/examples/extended/visualization/userVisAction/include/B1RunAction.hh
• source/geant4.10.03.p03/examples/basic/B1/src/B1RunAction.cc
• source/geant4.10.03.p03/examples/extended/visualization/userVisAction/src/B1RunAction.cc