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