| 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 // 26 //
27 /// \file B3/B3b/src/Run.cc 27 /// \file B3/B3b/src/Run.cc
28 /// \brief Implementation of the B3b::Run clas 28 /// \brief Implementation of the B3b::Run class
29 29
30 #include "Run.hh" 30 #include "Run.hh"
31 31
32 #include "G4Event.hh" <<
33 #include "G4HCofThisEvent.hh" <<
34 #include "G4RunManager.hh" 32 #include "G4RunManager.hh"
>> 33 #include "G4Event.hh"
>> 34
35 #include "G4SDManager.hh" 35 #include "G4SDManager.hh"
36 #include "G4SystemOfUnits.hh" << 36 #include "G4HCofThisEvent.hh"
37 #include "G4THitsMap.hh" 37 #include "G4THitsMap.hh"
>> 38 #include "G4SystemOfUnits.hh"
38 39
39 namespace B3b 40 namespace B3b
40 { 41 {
41 42
42 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 44
44 void Run::RecordEvent(const G4Event* event) 45 void Run::RecordEvent(const G4Event* event)
45 { 46 {
46 if (fCollID_cryst < 0) { << 47 if ( fCollID_cryst < 0 ) {
47 fCollID_cryst = G4SDManager::GetSDMpointer << 48 fCollID_cryst
48 // G4cout << " fCollID_cryst: " << fCollID << 49 = G4SDManager::GetSDMpointer()->GetCollectionID("crystal/edep");
>> 50 //G4cout << " fCollID_cryst: " << fCollID_cryst << G4endl;
49 } 51 }
50 52
51 if (fCollID_patient < 0) { << 53 if ( fCollID_patient < 0 ) {
52 fCollID_patient = G4SDManager::GetSDMpoint << 54 fCollID_patient
53 // G4cout << " fCollID_patient: " << fColl << 55 = G4SDManager::GetSDMpointer()->GetCollectionID("patient/dose");
>> 56 //G4cout << " fCollID_patient: " << fCollID_patient << G4endl;
54 } 57 }
55 58
56 G4int evtNb = event->GetEventID(); 59 G4int evtNb = event->GetEventID();
57 60
58 if (evtNb % fPrintModulo == 0) { << 61 if (evtNb%fPrintModulo == 0) {
59 G4cout << G4endl << "---> end of event: " 62 G4cout << G4endl << "---> end of event: " << evtNb << G4endl;
60 } 63 }
61 64
62 // Hits collections << 65 //Hits collections
63 // 66 //
64 G4HCofThisEvent* HCE = event->GetHCofThisEve 67 G4HCofThisEvent* HCE = event->GetHCofThisEvent();
65 if (!HCE) return; << 68 if(!HCE) return;
66 69
67 // Energy in crystals : identify 'good event << 70 //Energy in crystals : identify 'good events'
68 // 71 //
69 const G4double eThreshold = 500 * keV; << 72 const G4double eThreshold = 500*keV;
70 G4int nbOfFired = 0; 73 G4int nbOfFired = 0;
71 74
72 auto evtMap = static_cast<G4THitsMap<G4doubl 75 auto evtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID_cryst));
73 76
74 for (auto& mapElement : (*evtMap->GetMap())) << 77 std::map<G4int,G4double*>::iterator itr;
75 auto edep = *(mapElement.second); << 78 for (itr = evtMap->GetMap()->begin(); itr != evtMap->GetMap()->end(); itr++) {
>> 79 G4double edep = *(itr->second);
76 if (edep > eThreshold) nbOfFired++; 80 if (edep > eThreshold) nbOfFired++;
77 // auto copyNb = mapElement.first; << 81 ///G4int copyNb = (itr->first);
78 // G4cout << "\n cryst" << copyNb << ": " << 82 ///G4cout << G4endl << " cryst" << copyNb << ": " << edep/keV << " keV ";
79 } 83 }
80 if (nbOfFired == 2) fGoodEvents++; 84 if (nbOfFired == 2) fGoodEvents++;
81 85
82 // Dose deposit in patient << 86 //Dose deposit in patient
83 // 87 //
84 G4double dose = 0.; 88 G4double dose = 0.;
85 89
86 evtMap = static_cast<G4THitsMap<G4double>*>( 90 evtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID_patient));
87 91
88 for (auto& mapElement : (*evtMap->GetMap())) << 92 for (itr = evtMap->GetMap()->begin(); itr != evtMap->GetMap()->end(); itr++) {
89 dose += *(mapElement.second); << 93 ///G4int copyNb = (itr->first);
90 // auto copyNb = mapElement.first; << 94 dose = *(itr->second);
91 // G4cout << "\n patient" << copyNb << ": <<
92 } 95 }
93 fSumDose += dose; 96 fSumDose += dose;
94 fStatDose += dose; 97 fStatDose += dose;
95 98
96 G4Run::RecordEvent(event); 99 G4Run::RecordEvent(event);
97 } 100 }
98 101
99 //....oooOO0OOooo........oooOO0OOooo........oo 102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
100 103
101 void Run::Merge(const G4Run* run) << 104 void Run::Merge(const G4Run* aRun)
102 { 105 {
103 const Run* localRun = static_cast<const Run* << 106 const Run* localRun = static_cast<const Run*>(aRun);
104 fGoodEvents += localRun->fGoodEvents; 107 fGoodEvents += localRun->fGoodEvents;
105 fSumDose += localRun->fSumDose; << 108 fSumDose += localRun->fSumDose;
106 fStatDose += localRun->fStatDose; << 109 fStatDose += localRun->fStatDose;
107 G4Run::Merge(run); << 110 G4Run::Merge(aRun);
108 } 111 }
109 112
110 //....oooOO0OOooo........oooOO0OOooo........oo 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
111 114
112 } // namespace B3b << 115 }
>> 116
113 117