| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 /// \file scavenger/src/Run.cc 26 /// \file scavenger/src/Run.cc
27 /// \brief Implementation of the scavenger::Ru 27 /// \brief Implementation of the scavenger::Run class
28 28
29 #include "Run.hh" 29 #include "Run.hh"
30 <<
31 #include "ScoreSpecies.hh" <<
32 <<
33 #include "G4Event.hh" <<
34 #include "G4HCofThisEvent.hh" <<
35 #include "G4RunManager.hh" 30 #include "G4RunManager.hh"
>> 31 #include "G4Event.hh"
>> 32 #include "ScoreSpecies.hh"
36 #include "G4SDManager.hh" 33 #include "G4SDManager.hh"
>> 34 #include "G4HCofThisEvent.hh"
37 #include "G4THitsMap.hh" 35 #include "G4THitsMap.hh"
38 36
39 namespace scavenger 37 namespace scavenger
40 { 38 {
41 39
42 //....oooOO0OOooo........oooOO0OOooo........oo 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
43 41
44 Run::Run() : G4Run() << 42 Run::Run()
45 { << 43 : G4Run(){
46 G4MultiFunctionalDetector* mfdet = dynamic_c << 44 G4MultiFunctionalDetector *mfdet = dynamic_cast<G4MultiFunctionalDetector *>
47 G4SDManager::GetSDMpointer()->FindSensitiv << 45 (G4SDManager::GetSDMpointer()->FindSensitiveDetector("mfDetector"));
48 G4int CollectionIDspecies = G4SDManager::Get << 46 G4int CollectionIDspecies =
>> 47 G4SDManager::GetSDMpointer()->GetCollectionID("mfDetector/Species");
49 fScorerRun = mfdet->GetPrimitive(CollectionI 48 fScorerRun = mfdet->GetPrimitive(CollectionIDspecies);
50 } 49 }
51 50
52 //....oooOO0OOooo........oooOO0OOooo........oo 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
53 52
54 void Run::RecordEvent(const G4Event* event) << 53 void Run::RecordEvent(const G4Event *event) {
55 { <<
56 if (event->IsAborted()) { 54 if (event->IsAborted()) {
57 return; 55 return;
58 } 56 }
59 G4int CollectionID = G4SDManager::GetSDMpoin << 57 G4int CollectionID =
60 // Hits collections << 58 G4SDManager::GetSDMpointer()->GetCollectionID("mfDetector/Species");
61 G4HCofThisEvent* pHCE = event->GetHCofThisEv << 59 //Hits collections
>> 60 G4HCofThisEvent *pHCE = event->GetHCofThisEvent();
62 if (!pHCE) { 61 if (!pHCE) {
63 return; 62 return;
64 } 63 }
65 auto evtMap = dynamic_cast<G4THitsMap<G4doub << 64 auto evtMap = dynamic_cast<G4THitsMap<G4double> *>(pHCE->GetHC(CollectionID));
66 auto map = evtMap->GetMap(); 65 auto map = evtMap->GetMap();
67 for (const auto& it : *map) { << 66 for (const auto &it : *map) {
68 G4double edep = *(it.second); 67 G4double edep = *(it.second);
69 fSumEne += edep; 68 fSumEne += edep;
70 } 69 }
71 G4Run::RecordEvent(event); 70 G4Run::RecordEvent(event);
72 } 71 }
73 72
74 //....oooOO0OOooo........oooOO0OOooo........oo 73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
75 74
76 void Run::Merge(const G4Run* aRun) << 75 void Run::Merge(const G4Run *aRun) {
77 { <<
78 if (aRun == this) { 76 if (aRun == this) {
79 return; 77 return;
80 } 78 }
81 const auto localRun = dynamic_cast<const Run << 79 const auto localRun = dynamic_cast<const Run *>(aRun);
82 fSumEne += localRun->fSumEne; 80 fSumEne += localRun->fSumEne;
83 auto masterScorer = dynamic_cast<ScoreSpecie << 81 auto masterScorer = dynamic_cast<ScoreSpecies *>(this->fScorerRun);
84 auto localScorer = dynamic_cast<ScoreSpecies << 82 auto localScorer = dynamic_cast<ScoreSpecies *>(localRun->fScorerRun);
85 masterScorer->AbsorbResultsFromWorkerScorer( 83 masterScorer->AbsorbResultsFromWorkerScorer(localScorer);
86 G4Run::Merge(aRun); 84 G4Run::Merge(aRun);
87 } 85 }
88 86
89 //....oooOO0OOooo........oooOO0OOooo........oo 87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
90 88
91 } // namespace scavenger << 89 }