RunAction Class Reference

#include <RunAction.hh>

Inheritance diagram for RunAction:

[legend]

Collaboration diagram for RunAction:

[legend]

List of all members.

Public Member Functions
	RunAction (DetectorConstruction *)
	~RunAction ()
void	BeginOfRunAction (const G4Run *)
void	EndOfRunAction (const G4Run *)
void	SumEvents (G4double e)
void	SumDepthDose (G4double, G4double)
void	WriteDepthDose (G4String)
Private Attributes
DetectorConstruction *	m_detector
G4double	m_sumEdep
G4double	m_sumEdep2
G4int	m_nbEntry
std::vector< G4double >	m_depthDose
G4double	m_astronautHeight
G4double	m_sliceWidth
G4int	m_nbSlices

---

Detailed Description

Definition at line 19 of file RunAction.hh.

---

Constructor & Destructor Documentation

RunAction::RunAction

(

DetectorConstruction *

det

)

Definition at line 17 of file RunAction.cc.

:m_detector(det)
{ }

RunAction::~RunAction

(

)

Definition at line 23 of file RunAction.cc.

{ }

---

Member Function Documentation

void RunAction::BeginOfRunAction

(

const G4Run *

run

)

Definition at line 28 of file RunAction.cc.

{  
  G4cout << "### Run " << run->GetRunID() << " start." << G4endl;
  
  //initialize total energy deposit
  //
  m_sumEdep = m_sumEdep2 = 0.;
  m_nbEntry = 0;
  
  //construct vector of m_depthDose distribution
  //
  m_nbSlices = 100;
  m_astronautHeight = m_detector->GetAstronautHeight();
  m_sliceWidth = m_astronautHeight/m_nbSlices;
  m_depthDose.resize(m_nbSlices, 0.);
}

void RunAction::EndOfRunAction

(

const G4Run *

)

Definition at line 65 of file RunAction.cc.

{ 
  G4double meanEd = 0., meanEd2 = 0., rmsEd = 0., error_rel_sum = 0.;
  if (m_nbEntry) { 
    meanEd = m_sumEdep/m_nbEntry; meanEd2 = m_sumEdep2/m_nbEntry;
    G4double variance = meanEd2 - meanEd*meanEd;
    if (variance > 0.) rmsEd = std::sqrt(variance);
    error_rel_sum = rmsEd/(meanEd*std::sqrt(double(m_nbEntry)));
  }
   
  // choose printing format
  std::ios::fmtflags mode = G4cout.flags();
  G4cout.setf(std::ios::fixed,std::ios::floatfield);
  G4int prec = G4cout.precision(3);   
  
  //print total energy deposit
  //
  G4cout
    << "\n ---> Nb of non-empty events = " << m_nbEntry
    << "\n ---> Mean Edep per event    = " << G4BestUnit(meanEd,"Energy")     
    << "\n ---> Total Energy deposited = " << G4BestUnit(m_sumEdep,"Energy")
    << " +- " << G4BestUnit(m_sumEdep*error_rel_sum,"Energy")
    << "  (-> " << 100*error_rel_sum << " %)" 
    << G4endl;
    
  //get mass of astronaut
  //
  G4double mass = m_detector->GetAstronaut()->GetLogicalVolume()->GetMass();
  
  //compute dose
  //
  G4double dose = m_sumEdep/mass;
  G4cout.setf(mode,std::ios::floatfield);
  G4cout.precision(5);     
  G4cout << " ---> Dose = " << G4BestUnit(dose,"Dose") 
         << " +- " << G4BestUnit(dose*error_rel_sum,"Dose") << "\n" 
         << G4endl;
  
  // restaure default formats
  G4cout.setf(mode,std::ios::floatfield);
  G4cout.precision(prec);
        
  // convert m_depthDose array from energy to dose
  //
  for (G4int k=0; k<m_nbSlices; k++)  m_depthDose[k] /= mass;
  
  // write m_depthDose array on a file
  //
  WriteDepthDose("depthDose");
  
  // restaure default formats
  G4cout.setf(mode,std::ios::floatfield);
  G4cout.precision(prec);    
}

void RunAction::SumEvents

(

G4double

e

)

[inline]

Definition at line 29 of file RunAction.hh.

Referenced by EventAction::EndOfEventAction().

{ m_sumEdep += e; m_sumEdep2 += e*e; m_nbEntry++; };

void RunAction::SumDepthDose	(	G4double	zLocal,
		G4double	edep
	)

Definition at line 47 of file RunAction.cc.

Referenced by SteppingAction::UserSteppingAction().

{
  // check coherence
  if (zLocal > m_astronautHeight) {
    G4cout << "\n --> warning from RunAction::SumDepthDose() : "
           << " zLocal beyond astronaut height : " << zLocal/mm << G4endl;
    return;
  }
  
  // compute slice number               
  G4int k = int(zLocal/m_sliceWidth);
  
  // sum edep
  m_depthDose[k] += edep; 
}

void RunAction::WriteDepthDose

(

G4String

name

)

Definition at line 124 of file RunAction.cc.

Referenced by EndOfRunAction().

{ 
  G4String fileName = name + ".ascii";
  std::ofstream File(fileName, std::ios::out);
  
  std::ios::fmtflags mode = File.flags();  
  File.setf( std::ios::scientific, std::ios::floatfield );
  G4int prec = File.precision(4);
      
  File << "  Longitudinal depth dose distribution \n " 
       << "\n  zLocal (mm)\t  m_depthDose(gray) \n" << G4endl;
           
  G4double zLocal;    
  for (G4int k=0; k<m_nbSlices; k++) {
     zLocal = (k+0.5)*m_sliceWidth;
     File << "  "   << zLocal/mm 
          << "\t  " << m_depthDose[k]/gray << G4endl;
  }
  
  // restaure default formats
  File.setf(mode,std::ios::floatfield);
  File.precision(prec);         
}

---

Member Data Documentation

DetectorConstruction* RunAction::m_detector [private]

Definition at line 35 of file RunAction.hh.

Referenced by BeginOfRunAction(), and EndOfRunAction().

G4double RunAction::m_sumEdep [private]

Definition at line 37 of file RunAction.hh.

Referenced by BeginOfRunAction(), EndOfRunAction(), and SumEvents().

G4double RunAction::m_sumEdep2 [private]

Definition at line 37 of file RunAction.hh.

Referenced by BeginOfRunAction(), EndOfRunAction(), and SumEvents().

G4int RunAction::m_nbEntry [private]

Definition at line 38 of file RunAction.hh.

Referenced by BeginOfRunAction(), EndOfRunAction(), and SumEvents().

std::vector<G4double> RunAction::m_depthDose [private]

Definition at line 40 of file RunAction.hh.

Referenced by BeginOfRunAction(), EndOfRunAction(), SumDepthDose(), and WriteDepthDose().

G4double RunAction::m_astronautHeight [private]

Definition at line 41 of file RunAction.hh.

Referenced by BeginOfRunAction(), and SumDepthDose().

G4double RunAction::m_sliceWidth [private]

Definition at line 42 of file RunAction.hh.

Referenced by BeginOfRunAction(), SumDepthDose(), and WriteDepthDose().

G4int RunAction::m_nbSlices [private]

Definition at line 43 of file RunAction.hh.

Referenced by BeginOfRunAction(), EndOfRunAction(), and WriteDepthDose().

---

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

• include/RunAction.hh
• src/RunAction.cc