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