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 medical/dna/range/src/Run.cc 27 /// \brief Implementation of the Run class 28 // 29 // $Id: Run.cc 71376 2013-06-14 07:44:50Z maire $ 30 // 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 33 34 #include "Run.hh" 35 36 #include "DetectorConstruction.hh" 37 #include "HistoManager.hh" 38 #include "PrimaryGeneratorAction.hh" 39 40 #include "G4Material.hh" 41 #include "G4MultiFunctionalDetector.hh" 42 #include "G4SDManager.hh" 43 #include "G4SystemOfUnits.hh" 44 #include "G4UnitsTable.hh" 45 #include "G4VPrimitiveScorer.hh" 46 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 48 49 Run::Run(const DetectorConstruction* /*detector*/) 50 : G4Run(), fParticle(0), fEkin(0.), fEdeposit(0.), fCollName(), fRunMap(0) 51 52 {} 53 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 55 56 Run::~Run() {} 57 void Run::RecordEvent(const G4Event* /*aEvent*/) 58 { 59 numberOfEvent++; // This is an original line. 60 } 61 G4THitsMap<G4double>* Run::GetHitsMap(const G4String& detName, const G4String& colName) 62 { 63 G4String fullName = detName + "/" + colName; 64 return GetHitsMap(fullName); 65 } 66 G4THitsMap<G4double>* Run::GetHitsMap(const G4String& fullName) 67 { 68 if (fCollName == fullName) { 69 return fRunMap; 70 } 71 return NULL; 72 } 73 void Run::DumpAllScorer() 74 { 75 // - Number of HitsMap in this RUN. 76 // - GetHitsMap and dump values. 77 78 G4THitsMap<G4double>* RunMap = GetHitsMap(); 79 if (RunMap) { 80 G4cout << " PrimitiveScorer RUN " << RunMap->GetSDname() << "," << RunMap->GetName() << G4endl; 81 G4cout << " Number of entries " << RunMap->entries() << G4endl; 82 std::map<G4int, G4double*>::iterator itr = RunMap->GetMap()->begin(); 83 for (; itr != RunMap->GetMap()->end(); itr++) { 84 G4cout << " copy no.: " << itr->first << " Run Value : " << *(itr->second) << G4endl; 85 } 86 } 87 } 88 89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 90 91 void Run::SetPrimary(G4ParticleDefinition* particle, G4double energy) 92 { 93 fParticle = particle; 94 fEkin = energy; 95 } 96 97 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 98 99 void Run::Merge(const G4Run* run) 100 { 101 const Run* localRun = static_cast<const Run*>(run); 102 103 fParticle = localRun->fParticle; 104 fEkin = localRun->fEkin; 105 106 G4Run::Merge(run); 107 } 108 109 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 110 111 void Run::EndOfRun() 112 { 113 std::ios::fmtflags mode = G4cout.flags(); 114 G4cout.setf(std::ios::fixed, std::ios::floatfield); 115 G4int prec = G4cout.precision(2); 116 117 G4cout.setf(mode, std::ios::floatfield); 118 G4cout.precision(prec); 119 } 120 121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 122