| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 /// \file Run.cc
27 /// \brief Implementation of the Run class
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo
31 //....oooOO0OOooo........oooOO0OOooo........oo
32
33 #include "Run.hh"
34
35 #include "DetectorConstruction.hh"
36 #include "HistoManager.hh"
37 #include "PrimaryGeneratorAction.hh"
38
39 #include "G4SystemOfUnits.hh"
40 #include "G4UnitsTable.hh"
41
42 //....oooOO0OOooo........oooOO0OOooo........oo
43
44 Run::Run(DetectorConstruction* det) : fDetecto
45
46 //....oooOO0OOooo........oooOO0OOooo........oo
47
48 void Run::SetPrimary(G4ParticleDefinition* par
49 {
50 fParticle = particle;
51 fEkin = energy;
52 }
53
54 //....oooOO0OOooo........oooOO0OOooo........oo
55
56 void Run::CountProcesses(const G4VProcess* pro
57 {
58 if (process == nullptr) return;
59 G4String procName = process->GetProcessName(
60 std::map<G4String, G4int>::iterator it = fPr
61 if (it == fProcCounter.end()) {
62 fProcCounter[procName] = 1;
63 }
64 else {
65 fProcCounter[procName]++;
66 }
67 }
68
69 //....oooOO0OOooo........oooOO0OOooo........oo
70
71 void Run::ParticleCount(G4String name, G4doubl
72 {
73 std::map<G4String, ParticleData>::iterator i
74 if (it == fParticleDataMap.end()) {
75 fParticleDataMap[name] = ParticleData(1, E
76 }
77 else {
78 ParticleData& data = it->second;
79 data.fCount++;
80 data.fEmean += Ekin;
81 // update min max
82 G4double emin = data.fEmin;
83 if (Ekin < emin) data.fEmin = Ekin;
84 G4double emax = data.fEmax;
85 if (Ekin > emax) data.fEmax = Ekin;
86 }
87 }
88
89 //....oooOO0OOooo........oooOO0OOooo........oo
90
91 void Run::SumTrackLength(G4int nstep1, G4int n
92 G4double time1, G4dou
93 {
94 fNbStep1 += nstep1;
95 fNbStep2 += nstep2;
96 fTrackLen1 += trackl1;
97 fTrackLen2 += trackl2;
98 fTime1 += time1;
99 fTime2 += time2;
100 }
101
102 //....oooOO0OOooo........oooOO0OOooo........oo
103
104 void Run::Merge(const G4Run* run)
105 {
106 const Run* localRun = static_cast<const Run*
107
108 // primary particle info
109 //
110 fParticle = localRun->fParticle;
111 fEkin = localRun->fEkin;
112
113 // accumulate sums
114 //
115 fNbStep1 += localRun->fNbStep1;
116 fNbStep2 += localRun->fNbStep2;
117 fTrackLen1 += localRun->fTrackLen1;
118 fTrackLen2 += localRun->fTrackLen2;
119 fTime1 += localRun->fTime1;
120 fTime2 += localRun->fTime2;
121
122 // map: processes count
123 std::map<G4String, G4int>::const_iterator it
124 for (itp = localRun->fProcCounter.begin(); i
125 G4String procName = itp->first;
126 G4int localCount = itp->second;
127 if (fProcCounter.find(procName) == fProcCo
128 fProcCounter[procName] = localCount;
129 }
130 else {
131 fProcCounter[procName] += localCount;
132 }
133 }
134
135 // map: created particles count
136 std::map<G4String, ParticleData>::const_iter
137 for (itn = localRun->fParticleDataMap.begin(
138 G4String name = itn->first;
139 const ParticleData& localData = itn->secon
140 if (fParticleDataMap.find(name) == fPartic
141 fParticleDataMap[name] =
142 ParticleData(localData.fCount, localDa
143 }
144 else {
145 ParticleData& data = fParticleDataMap[na
146 data.fCount += localData.fCount;
147 data.fEmean += localData.fEmean;
148 G4double emin = localData.fEmin;
149 if (emin < data.fEmin) data.fEmin = emin
150 G4double emax = localData.fEmax;
151 if (emax > data.fEmax) data.fEmax = emax
152 }
153 }
154
155 G4Run::Merge(run);
156 }
157
158 //....oooOO0OOooo........oooOO0OOooo........oo
159
160 void Run::EndOfRun()
161 {
162 G4int prec = 5, wid = prec + 2;
163 G4int dfprec = G4cout.precision(prec);
164
165 // run condition
166 //
167 G4Material* material = fDetector->GetMateria
168 G4double density = material->GetDensity();
169
170 G4String Particle = fParticle->GetParticleNa
171 G4cout << "\n The run is " << numberOfEvent
172 << G4BestUnit(fEkin, "Energy") << " t
173 << G4BestUnit(0.5 * (fDetector->GetSi
174 << " (density: " << G4BestUnit(densit
175
176 if (numberOfEvent == 0) {
177 G4cout.precision(dfprec);
178 return;
179 }
180
181 // frequency of processes
182 //
183 G4cout << "\n Process calls frequency :" <<
184 G4int survive = 0;
185 std::map<G4String, G4int>::iterator it;
186 for (it = fProcCounter.begin(); it != fProcC
187 G4String procName = it->first;
188 G4int count = it->second;
189 G4cout << "\t" << procName << "= " << coun
190 if (procName == "Transportation") survive
191 }
192 G4cout << G4endl;
193
194 if (survive > 0) {
195 G4cout << "\n Nb of incident particles sur
196 << G4BestUnit(0.5 * (fDetector->Get
197 << fDetector->GetMaterial()->GetNam
198 }
199
200 // total track length of incident neutron
201 //
202 G4cout << "\n Parcours of incident neutron:"
203
204 G4double meanCollision1 = (G4double)fNbStep1
205 G4double meanCollision2 = (G4double)fNbStep2
206 G4double meanCollisTota = meanCollision1 + m
207
208 G4cout << "\n nb of collisions E>1*eV=
209 << " E<1*eV= " << meanCollision2
210
211 G4double meanTrackLen1 = fTrackLen1 / number
212 G4double meanTrackLen2 = fTrackLen2 / number
213 G4double meanTrackLtot = meanTrackLen1 + mea
214
215 G4cout << "\n track length E>1*eV=
216 << " E<1*eV= " << G4BestUnit(meanTra
217 << " total= " << G4BestUnit(meanTra
218
219 G4double meanTime1 = fTime1 / numberOfEvent;
220 G4double meanTime2 = fTime2 / numberOfEvent;
221 G4double meanTimeTo = meanTime1 + meanTime2;
222
223 G4cout << "\n time of flight E>1*eV=
224 << " E<1*eV= " << G4BestUnit(meanTim
225 << " total= " << G4BestUnit(meanTim
226
227 // particles count
228 //
229 G4cout << "\n List of generated particles:"
230
231 std::map<G4String, ParticleData>::iterator i
232 for (itn = fParticleDataMap.begin(); itn !=
233 G4String name = itn->first;
234 ParticleData data = itn->second;
235 G4int count = data.fCount;
236 G4double eMean = data.fEmean / count;
237 G4double eMin = data.fEmin;
238 G4double eMax = data.fEmax;
239
240 G4cout << " " << std::setw(13) << name <<
241 << " Emean = " << std::setw(wid) <
242 << G4BestUnit(eMin, "Energy") << "
243 }
244
245 // normalize histograms
246 ////G4AnalysisManager* analysisManager = G4A
247 ////G4double factor = 1./numberOfEvent;
248 ////analysisManager->ScaleH1(3,factor);
249
250 // remove all contents in fProcCounter, fCou
251 fProcCounter.clear();
252 fParticleDataMap.clear();
253
254 // restore default format
255 G4cout.precision(dfprec);
256 }
257
258 //....oooOO0OOooo........oooOO0OOooo........oo
259