73{
74 static G4String regName[3] = {
"Calor-A",
"Calor-B",
"Calor-C"};
75
77
78 if(IsMaster()) {
79
80 G4cout
81 << "############################################################" << G4endl;
82 G4cout
83 << " Run Summary - Number of events : " << theRun->GetNumberOfEvent()
84 << G4endl;
85 G4cout
86 << "############################################################" << G4endl;
87
88 G4double nEvt = (G4double)(theRun->GetNumberOfEvent());
89 for(size_t i=0;i<3;i++)
90 {
91 size_t ih1 = 2*i;
92 size_t ih2 = 2*i+1;
93
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")
101 << G4endl;
102 G4cout << " Energy deposition in an event :" << G4endl << " "
105 << 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 "
117 << G4endl << " "
118 << " e- in Absorber "
121 << G4endl << " "
122 << " e+ in Absorber "
125 << G4endl;
126 G4cout << " Total track length of e+/e- in an event :" << G4endl
127 << " "
130 << G4endl;
131 G4cout << " Total number of steps of e+/e- in an event :" << G4endl
132 << " "
133 <<
" Absorber " << (theRun->
GetNStep(ih1))/nEvt
134 <<
" Gap " << (theRun->
GetNStep(ih2))/nEvt
135 << G4endl;
136 G4cout
137 << "------------------------------------------------------------"
138 << G4endl;
139 G4cout << "Scores in parallel geometry" << G4endl;
140 G4cout
141 << "layer eDep/evt nGamma/evt nElec/evt nPosi/evt stpLen/evt nStep/evt"
142 << G4endl;
143 for(size_t k=0;k<20;k++)
144 {
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;
148 }
149 G4cout << G4endl;
150 }
151 G4cout
152 << "############################################################"
153 << G4endl;
154 }
155 } else {
156
157 G4cout << "CPU Time spent by each region" << G4endl;
161
162 }
163}
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