B1RunAction Class Reference

Run action class. More...

#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 *)

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

Constructor & Destructor Documentation

B1RunAction()

B1RunAction::B1RunAction

(

)

Definition at line 43 of file B1RunAction.cc.

: 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::~B1RunAction ( )

virtual

Definition at line 61 of file B1RunAction.cc.

{}

Member Function Documentation

GenerateRun()

G4Run * B1RunAction::GenerateRun ( )

virtual

Definition at line 66 of file B1RunAction.cc.

{
  return new B1Run; 
}

BeginOfRunAction()

void B1RunAction::BeginOfRunAction ( const G4Run *  )

virtual

Definition at line 73 of file B1RunAction.cc.

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

EndOfRunAction()

void B1RunAction::EndOfRunAction ( const G4Run *  run )

virtual

Definition at line 81 of file B1RunAction.cc.

{
  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;
}

---

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

• extended/visualization/userVisAction/include/B1RunAction.hh
• extended/visualization/userVisAction/src/B1RunAction.cc