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 ()
void	SetMaster (G4bool val=true)
G4bool	IsMaster () const

Private Attributes
G4Parameter< G4double >	fEdep
G4Parameter< G4double >	fEdep2

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 basic/B1/include/B1RunAction.hh.

Constructor & Destructor Documentation

B1RunAction() [1/2]

B1RunAction::B1RunAction

(

)

Definition at line 46 of file basic/B1/src/B1RunAction.cc.

: G4UserRunAction(),
  fEdep("Edep", 0.),
  fEdep2("Edep2", 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 parameter to the parameter manager
  G4ParameterManager* parameterManager = G4ParameterManager::Instance();
  parameterManager->RegisterParameter(fEdep);
  parameterManager->RegisterParameter(fEdep2); 
}

Here is the call graph for this function:

~B1RunAction() [1/2]

B1RunAction::~B1RunAction ( )

virtual

Definition at line 71 of file basic/B1/src/B1RunAction.cc.

{}

B1RunAction() [2/2]

B1RunAction::B1RunAction

(

)

~B1RunAction() [2/2]

virtual B1RunAction::~B1RunAction ( )

virtual

Member Function Documentation

AddEdep()

void B1RunAction::AddEdep

(

G4double

edep

)

Definition at line 156 of file basic/B1/src/B1RunAction.cc.

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

Here is the caller graph for this function:

BeginOfRunAction() [1/2]

virtual void B1RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

BeginOfRunAction() [2/2]

void B1RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 76 of file basic/B1/src/B1RunAction.cc.

{ 
  // inform the runManager to save random number seed
  G4RunManager::GetRunManager()->SetRandomNumberStore(false);

  // reset parameters to their initial values
  G4ParameterManager* parameterManager = G4ParameterManager::Instance();
  parameterManager->Reset();

}

Here is the call graph for this function:

Here is the caller graph for this function:

EndOfRunAction() [1/2]

virtual void B1RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

EndOfRunAction() [2/2]

void B1RunAction::EndOfRunAction ( const G4Run *  run )

virtual

Reimplemented from G4UserRunAction.

Definition at line 89 of file basic/B1/src/B1RunAction.cc.

{
  G4int nofEvents = run->GetNumberOfEvent();
  if (nofEvents == 0) return;

  // Merge parameters 
  G4ParameterManager* parameterManager = G4ParameterManager::Instance();
  parameterManager->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:

Here is the caller graph for this function:

GenerateRun()

G4Run * B1RunAction::GenerateRun ( void  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 67 of file extended/visualization/userVisAction/src/B1RunAction.cc.

{
  return new B1Run; 
}

Here is the call graph for this function:

Member Data Documentation

fEdep

G4Parameter<G4double> B1RunAction::fEdep

private

Definition at line 59 of file basic/B1/include/B1RunAction.hh.

fEdep2

G4Parameter<G4double> B1RunAction::fEdep2

private

Definition at line 60 of file basic/B1/include/B1RunAction.hh.

---

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

• basic/B1/include/B1RunAction.hh
• basic/B1/src/B1RunAction.cc