B1RunAction Class Reference

#include <B1RunAction.hh>

Inheritance diagram for B1RunAction:

Public Member Functions
	B1RunAction ()
virtual	~B1RunAction ()
virtual G4Run *	GenerateRun ()
virtual void	BeginOfRunAction (const G4Run *)
virtual void	EndOfRunAction (const G4Run *)
	B1RunAction ()
virtual	~B1RunAction ()
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

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::B1RunAction

(

)

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

References python.hepunit::gray.

: G4UserRunAction()
{ 
  // 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);        
}

B1RunAction::~B1RunAction ( )

virtual

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

{}

B1RunAction::B1RunAction

(

)

virtual B1RunAction::~B1RunAction ( )

virtual

Member Function Documentation

virtual void B1RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

void B1RunAction::BeginOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

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

References G4RunManager::GetRunManager(), and G4RunManager::SetRandomNumberStore().

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

virtual void B1RunAction::EndOfRunAction ( const G4Run *  )

virtual

Reimplemented from G4UserRunAction.

void B1RunAction::EndOfRunAction ( const G4Run *  run )

virtual

Reimplemented from G4UserRunAction.

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

References G4BestUnit, G4cout, G4endl, B1Run::GetEdep(), B1Run::GetEdep2(), G4LogicalVolume::GetMass(), G4Run::GetNumberOfEvent(), G4ParticleGun::GetParticleDefinition(), G4ParticleGun::GetParticleEnergy(), B1PrimaryGeneratorAction::GetParticleGun(), G4ParticleDefinition::GetParticleName(), G4RunManager::GetRunManager(), B1DetectorConstruction::GetScoringVolume(), G4RunManager::GetUserDetectorConstruction(), G4RunManager::GetUserPrimaryGeneratorAction(), and G4UserRunAction::IsMaster().

{
  G4int nofEvents = run->GetNumberOfEvent();
  if (nofEvents == 0) return;
  
  const B1Run* b1Run = static_cast<const B1Run*>(run);

  // Compute dose
  //
  G4double edep  = b1Run->GetEdep();
  G4double edep2 = b1Run->GetEdep2();
  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
     << "\n--------------------End of Global Run-----------------------";
  }
  else {
    G4cout
     << "\n--------------------End of Local Run------------------------";
  }
  
  G4cout
     << "\n The run consists of " << nofEvents << " "<< runCondition
     << "\n Dose in scoring volume : " 
     << G4BestUnit(dose,"Dose") << " +- " << G4BestUnit(rmsDose,"Dose")
     << "\n------------------------------------------------------------\n"
     << G4endl;
}

G4Run * B1RunAction::GenerateRun ( void  )

virtual

Reimplemented from G4UserRunAction.

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

{
  return new B1Run; 
}

---

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

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