| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 /// \file Run.cc
27 /// \brief Implementation of the Run class
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
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........oooOO0OOooo........oooOO0OOooo......
43
44 Run::Run(DetectorConstruction* det) : fDetector(det) {}
45
46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47
48 void Run::SetPrimary(G4ParticleDefinition* particle, G4double energy)
49 {
50 fParticle = particle;
51 fEkin = energy;
52 }
53
54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55
56 void Run::CountProcesses(const G4VProcess* process)
57 {
58 if (process == nullptr) return;
59 G4String procName = process->GetProcessName();
60 std::map<G4String, G4int>::iterator it = fProcCounter.find(procName);
61 if (it == fProcCounter.end()) {
62 fProcCounter[procName] = 1;
63 }
64 else {
65 fProcCounter[procName]++;
66 }
67 }
68
69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
70
71 void Run::ParticleCount(G4String name, G4double Ekin)
72 {
73 std::map<G4String, ParticleData>::iterator it = fParticleDataMap.find(name);
74 if (it == fParticleDataMap.end()) {
75 fParticleDataMap[name] = ParticleData(1, Ekin, Ekin, Ekin);
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........oooOO0OOooo........oooOO0OOooo......
90
91 void Run::SumTrackLength(G4int nstep1, G4int nstep2, G4double trackl1, G4double trackl2,
92 G4double time1, G4double time2)
93 {
94 fNbStep1 += nstep1;
95 fNbStep2 += nstep2;
96 fTrackLen1 += trackl1;
97 fTrackLen2 += trackl2;
98 fTime1 += time1;
99 fTime2 += time2;
100 }
101
102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
103
104 void Run::Merge(const G4Run* run)
105 {
106 const Run* localRun = static_cast<const Run*>(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 itp;
124 for (itp = localRun->fProcCounter.begin(); itp != localRun->fProcCounter.end(); ++itp) {
125 G4String procName = itp->first;
126 G4int localCount = itp->second;
127 if (fProcCounter.find(procName) == fProcCounter.end()) {
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_iterator itn;
137 for (itn = localRun->fParticleDataMap.begin(); itn != localRun->fParticleDataMap.end(); ++itn) {
138 G4String name = itn->first;
139 const ParticleData& localData = itn->second;
140 if (fParticleDataMap.find(name) == fParticleDataMap.end()) {
141 fParticleDataMap[name] =
142 ParticleData(localData.fCount, localData.fEmean, localData.fEmin, localData.fEmax);
143 }
144 else {
145 ParticleData& data = fParticleDataMap[name];
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........oooOO0OOooo........oooOO0OOooo......
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->GetMaterial();
168 G4double density = material->GetDensity();
169
170 G4String Particle = fParticle->GetParticleName();
171 G4cout << "\n The run is " << numberOfEvent << " " << Particle << " of "
172 << G4BestUnit(fEkin, "Energy") << " through "
173 << G4BestUnit(0.5 * (fDetector->GetSize()), "Length") << " of " << material->GetName()
174 << " (density: " << G4BestUnit(density, "Volumic Mass") << ")" << G4endl;
175
176 if (numberOfEvent == 0) {
177 G4cout.precision(dfprec);
178 return;
179 }
180
181 // frequency of processes
182 //
183 G4cout << "\n Process calls frequency :" << G4endl;
184 G4int survive = 0;
185 std::map<G4String, G4int>::iterator it;
186 for (it = fProcCounter.begin(); it != fProcCounter.end(); it++) {
187 G4String procName = it->first;
188 G4int count = it->second;
189 G4cout << "\t" << procName << "= " << count;
190 if (procName == "Transportation") survive = count;
191 }
192 G4cout << G4endl;
193
194 if (survive > 0) {
195 G4cout << "\n Nb of incident particles surviving after "
196 << G4BestUnit(0.5 * (fDetector->GetSize()), "Length") << " of "
197 << fDetector->GetMaterial()->GetName() << " : " << survive << G4endl;
198 }
199
200 // total track length of incident neutron
201 //
202 G4cout << "\n Parcours of incident neutron:";
203
204 G4double meanCollision1 = (G4double)fNbStep1 / numberOfEvent;
205 G4double meanCollision2 = (G4double)fNbStep2 / numberOfEvent;
206 G4double meanCollisTota = meanCollision1 + meanCollision2;
207
208 G4cout << "\n nb of collisions E>1*eV= " << meanCollision1
209 << " E<1*eV= " << meanCollision2 << " total= " << meanCollisTota;
210
211 G4double meanTrackLen1 = fTrackLen1 / numberOfEvent;
212 G4double meanTrackLen2 = fTrackLen2 / numberOfEvent;
213 G4double meanTrackLtot = meanTrackLen1 + meanTrackLen2;
214
215 G4cout << "\n track length E>1*eV= " << G4BestUnit(meanTrackLen1, "Length")
216 << " E<1*eV= " << G4BestUnit(meanTrackLen2, "Length")
217 << " total= " << G4BestUnit(meanTrackLtot, "Length");
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= " << G4BestUnit(meanTime1, "Time")
224 << " E<1*eV= " << G4BestUnit(meanTime2, "Time")
225 << " total= " << G4BestUnit(meanTimeTo, "Time") << G4endl;
226
227 // particles count
228 //
229 G4cout << "\n List of generated particles:" << G4endl;
230
231 std::map<G4String, ParticleData>::iterator itn;
232 for (itn = fParticleDataMap.begin(); itn != fParticleDataMap.end(); 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 << ": " << std::setw(7) << count
241 << " Emean = " << std::setw(wid) << G4BestUnit(eMean, "Energy") << "\t( "
242 << G4BestUnit(eMin, "Energy") << " --> " << G4BestUnit(eMax, "Energy") << ")" << G4endl;
243 }
244
245 // normalize histograms
246 ////G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
247 ////G4double factor = 1./numberOfEvent;
248 ////analysisManager->ScaleH1(3,factor);
249
250 // remove all contents in fProcCounter, fCount
251 fProcCounter.clear();
252 fParticleDataMap.clear();
253
254 // restore default format
255 G4cout.precision(dfprec);
256 }
257
258 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
259