B1ConRunAction Class Reference

#include <B1ConRunAction.hh>

Inheritance diagram for B1ConRunAction:

Public Member Functions
	B1ConRunAction ()
virtual	~B1ConRunAction ()
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

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 47 of file B1ConRunAction.hh.

Constructor & Destructor Documentation

B1ConRunAction::B1ConRunAction

(

)

Definition at line 44 of file B1ConRunAction.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);        

}

B1ConRunAction::~B1ConRunAction ( )

virtual

Definition at line 63 of file B1ConRunAction.cc.

{}

Member Function Documentation

void B1ConRunAction::BeginOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 75 of file B1ConRunAction.cc.

References G4cout, G4endl, G4Run::GetRunID(), G4RunManager::GetRunManager(), G4UserRunAction::IsMaster(), and G4RunManager::SetRandomNumberStore().

{ 
  G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;

  //inform the runManager to save random number seed
  G4RunManager::GetRunManager()->SetRandomNumberStore(false);
  
  if (IsMaster()) {
     fdose_tally = new G4ConvergenceTester("DOSE_TALLY");
     //fdose_tally = new G4ConvergenceTester();
  }

}

void B1ConRunAction::EndOfRunAction ( const G4Run *  aRun )

virtual

Reimplemented from G4UserRunAction.

Definition at line 91 of file B1ConRunAction.cc.

References G4ConvergenceTester::AddScore(), G4BestUnit, G4cout, G4endl, B1Run::GetEdep(), B1Run::GetEdep2(), B1ConRun::GetEdepPerEvent(), G4LogicalVolume::GetMass(), B1ConRun::GetNumberOfEvent(), G4Run::GetNumberOfEvent(), G4ParticleGun::GetParticleDefinition(), G4ParticleGun::GetParticleEnergy(), B1PrimaryGeneratorAction::GetParticleGun(), G4ParticleDefinition::GetParticleName(), G4RunManager::GetRunManager(), B1DetectorConstruction::GetScoringVolume(), G4RunManager::GetUserDetectorConstruction(), G4RunManager::GetUserPrimaryGeneratorAction(), python.hepunit::gray, G4UserRunAction::IsMaster(), G4ConvergenceTester::ShowHistory(), and G4ConvergenceTester::ShowResult().

{
  G4int nofEvents = aRun->GetNumberOfEvent();
  if (nofEvents == 0) return;
  
  const B1ConRun* b1ConRun = static_cast<const B1ConRun*>(aRun);

  // Compute dose
  //
  G4double edep  = b1ConRun->GetEdep();
  G4double edep2 = b1ConRun->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())
  {

     for ( G4int i = 0 ; i != b1ConRun->GetNumberOfEvent(); i++ ) {
        G4double aDose = b1ConRun->GetEdepPerEvent(i)/mass/gray; 
        fdose_tally->AddScore( aDose ); 
     }

     fdose_tally->ShowResult();
     fdose_tally->ShowHistory();
     delete fdose_tally;

    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 * B1ConRunAction::GenerateRun ( void  )

virtual

Reimplemented from G4UserRunAction.

Definition at line 68 of file B1ConRunAction.cc.

{ 
  return new B1ConRun; 
}

---

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

• B1ConRunAction.hh
• B1ConRunAction.cc