34#include "G4RegionStore.hh"
36#include "G4ProductionCuts.hh"
38#include "G4UnitsTable.hh"
40#include "G4VSteppingVerbose.hh"
74 static G4String regName[3] = {
"Calor-A",
"Calor-B",
"Calor-C"};
81 <<
"############################################################" << G4endl;
83 <<
" Run Summary - Number of events : " << theRun->GetNumberOfEvent()
86 <<
"############################################################" << G4endl;
88 G4double nEvt = (G4double)(theRun->GetNumberOfEvent());
89 for(
size_t i=0;i<3;i++)
94 G4Region* region = G4RegionStore::GetInstance()->GetRegion(regName[i]);
96 G4cout <<
"Region : " << region->GetName() << G4endl;
97 G4cout <<
" Production thresholds :" << G4endl <<
" "
98 <<
" gamma " <<
G4BestUnit(cuts->GetProductionCut(
"gamma"),
"Length")
99 <<
" e- " <<
G4BestUnit(cuts->GetProductionCut(
"e-"),
"Length")
100 <<
" e+ " <<
G4BestUnit(cuts->GetProductionCut(
"e+"),
"Length")
102 G4cout <<
" Energy deposition in an event :" << G4endl <<
" "
106 G4cout <<
" Number of secondaries in an event :" << G4endl <<
" "
107 <<
" gamma in Absorber " << (theRun->
GetNGamma(ih1))/nEvt
108 <<
" in Gap " << (theRun->
GetNGamma(ih2))/nEvt << G4endl <<
" "
109 <<
" e- in Absorber " << (theRun->
GetNElectron(ih1))/nEvt
110 <<
" in Gap " << (theRun->
GetNElectron(ih2))/nEvt << G4endl <<
" "
111 <<
" e+ in Absorber " << (theRun->
GetNPositron(ih1))/nEvt
112 <<
" in Gap " << (theRun->
GetNPositron(ih2))/nEvt << G4endl;
113 G4cout <<
" Minimum kinetic energy of generated secondaries :" << G4endl
114 <<
" " <<
" gamma in Absorber "
118 <<
" e- in Absorber "
122 <<
" e+ in Absorber "
126 G4cout <<
" Total track length of e+/e- in an event :" << G4endl
131 G4cout <<
" Total number of steps of e+/e- in an event :" << G4endl
133 <<
" Absorber " << (theRun->
GetNStep(ih1))/nEvt
134 <<
" Gap " << (theRun->
GetNStep(ih2))/nEvt
137 <<
"------------------------------------------------------------"
139 G4cout <<
"Scores in parallel geometry" << G4endl;
141 <<
"layer eDep/evt nGamma/evt nElec/evt nPosi/evt stpLen/evt nStep/evt"
143 for(
size_t k=0;k<20;k++)
145 G4cout << std::setw(8) << k;
146 for(
size_t j=0;j<6;j++) {
147 G4cout <<
" " << std::setw(10) << (theRun->
GetParaValue(i,j,k))/nEvt;
152 <<
"############################################################"
157 G4cout <<
"CPU Time spent by each region" << G4endl;
Definition of the RE06RunAction class.
Definition of the RE06Run class.
Definition of the RE06SteppingVerbose class.
virtual void BeginOfRunAction(const G4Run *)
virtual void EndOfRunAction(const G4Run *)
G4double GetTotalE(G4int i) const
G4double GetEMinPositron(G4int i) const
G4double GetEMinGamma(G4int i) const
G4double GetNElectron(G4int i) const
G4double GetNGamma(G4int i) const
G4double GetTotalL(G4int i) const
G4double GetParaValue(G4int i, G4int j, G4int k) const
G4double GetNPositron(G4int i) const
G4double GetNStep(G4int i) const
G4double GetEMinElectron(G4int i) const