| 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 // This example is provided by the Geant4-DNA 26 // This example is provided by the Geant4-DNA collaboration
27 // chem6 example is derived from chem4 and che 27 // chem6 example is derived from chem4 and chem5 examples
28 // 28 //
29 // Any report or published results obtained us 29 // Any report or published results obtained using the Geant4-DNA software
30 // shall cite the following Geant4-DNA collabo 30 // shall cite the following Geant4-DNA collaboration publication:
31 // J. Appl. Phys. 125 (2019) 104301 31 // J. Appl. Phys. 125 (2019) 104301
32 // Med. Phys. 45 (2018) e722-e739 32 // Med. Phys. 45 (2018) e722-e739
33 // J. Comput. Phys. 274 (2014) 841-882 33 // J. Comput. Phys. 274 (2014) 841-882
34 // Med. Phys. 37 (2010) 4692-4708 34 // Med. Phys. 37 (2010) 4692-4708
35 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 35 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157-178
36 // The Geant4-DNA web site is available at htt 36 // The Geant4-DNA web site is available at http://geant4-dna.org
37 // 37 //
38 // Authors: W. G. Shin and S. Incerti (CENBG, 38 // Authors: W. G. Shin and S. Incerti (CENBG, France)
39 // 39 //
40 // $Id$ 40 // $Id$
41 // 41 //
42 /// \file RunAction.cc 42 /// \file RunAction.cc
43 /// \brief Implementation of the RunAction cla 43 /// \brief Implementation of the RunAction class
44 44
45 #include "RunAction.hh" 45 #include "RunAction.hh"
46 <<
47 #include "Run.hh" 46 #include "Run.hh"
48 47
49 #include "G4Run.hh" 48 #include "G4Run.hh"
50 #include "G4RunManager.hh" 49 #include "G4RunManager.hh"
51 #include "G4SystemOfUnits.hh" <<
52 #include "G4UnitsTable.hh" 50 #include "G4UnitsTable.hh"
>> 51 #include "G4SystemOfUnits.hh"
53 52
54 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
55 54
56 extern std::ofstream out; 55 extern std::ofstream out;
57 56
58 RunAction::RunAction() : G4UserRunAction() {} << 57 RunAction::RunAction()
>> 58 : G4UserRunAction()
>> 59 {
>> 60 }
59 61
60 //....oooOO0OOooo........oooOO0OOooo........oo 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
61 63
62 RunAction::~RunAction() {} << 64 RunAction::~RunAction()
>> 65 {
>> 66 }
63 67
64 //....oooOO0OOooo........oooOO0OOooo........oo 68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
65 69
66 G4Run* RunAction::GenerateRun() 70 G4Run* RunAction::GenerateRun()
67 { 71 {
68 Run* run = new Run(); 72 Run* run = new Run();
69 return run; 73 return run;
70 } 74 }
71 75
72 //....oooOO0OOooo........oooOO0OOooo........oo 76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
73 77
74 void RunAction::BeginOfRunAction(const G4Run* 78 void RunAction::BeginOfRunAction(const G4Run* run)
75 { 79 {
76 G4cout << "### Run " << run->GetRunID() << " 80 G4cout << "### Run " << run->GetRunID() << " starts." << G4endl;
77 81
78 // informs the runManager to save random num 82 // informs the runManager to save random number seed
79 G4RunManager::GetRunManager()->SetRandomNumb 83 G4RunManager::GetRunManager()->SetRandomNumberStore(false);
80 } 84 }
81 85
82 //....oooOO0OOooo........oooOO0OOooo........oo 86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
83 87
84 void RunAction::EndOfRunAction(const G4Run* ru 88 void RunAction::EndOfRunAction(const G4Run* run)
85 { 89 {
86 G4int nofEvents = run->GetNumberOfEvent(); 90 G4int nofEvents = run->GetNumberOfEvent();
87 if (nofEvents == 0) return; 91 if (nofEvents == 0) return;
88 92
89 // results 93 // results
90 // 94 //
91 const Run* chem6Run = static_cast<const Run* 95 const Run* chem6Run = static_cast<const Run*>(run);
92 G4double sumDose = chem6Run->GetSumDose(); << 96 G4double sumDose = chem6Run->GetSumDose();
93 97
94 // print 98 // print
95 // 99 //
96 if (IsMaster()) { << 100 if (IsMaster())
97 G4cout << G4endl << "--------------------E << 101 {
98 << " The run has " << nofEvents << << 102
>> 103 G4cout
>> 104 << G4endl
>> 105 << "--------------------End of Global Run-----------------------"
>> 106 << G4endl
>> 107 << " The run has " << nofEvents << " events "
>> 108 << G4endl;
99 109
100 ScoreSpecies* masterScorer = dynamic_cast< << 110 ScoreSpecies* masterScorer=
>> 111 dynamic_cast<ScoreSpecies*>(chem6Run->GetPrimitiveScorer());
101 112
102 G4cout << "Number of events recorded by th 113 G4cout << "Number of events recorded by the species scorer = "
103 << masterScorer->GetNumberOfRecorde << 114 << masterScorer->GetNumberOfRecordedEvents()
>> 115 << G4endl;
104 116
105 // LET 117 // LET
106 Run* aRun = (Run*)run; 118 Run* aRun = (Run*)run;
107 G4THitsMap<G4double>* totLET = aRun->GetLE 119 G4THitsMap<G4double>* totLET = aRun->GetLET();
108 G4int nOfEvent = totLET->entries(); 120 G4int nOfEvent = totLET->entries();
109 G4double LET_mean = 0; 121 G4double LET_mean = 0;
110 G4double LET_square = 0; 122 G4double LET_square = 0;
111 for (G4int i = 0; i < nOfEvent; i++) { << 123 for(G4int i=0;i<nOfEvent;i++){
112 G4double* LET = (*totLET)[i]; 124 G4double* LET = (*totLET)[i];
113 if (!LET) continue; << 125 if(!LET) continue;
114 LET_mean += *LET; 126 LET_mean += *LET;
115 LET_square += (*LET) * (*LET); << 127 LET_square += (*LET)*(*LET);
116 } 128 }
117 LET_mean /= nOfEvent; 129 LET_mean /= nOfEvent;
118 LET_square = std::sqrt(LET_square / nOfEve << 130 LET_square = std::sqrt(LET_square/nOfEvent - std::pow(LET_mean,2));
119 131
120 if (nOfEvent > 1) { << 132 if(nOfEvent > 1)
121 out << std::setw(12) << "LET" << std::se << 133 {
122 << std::setw(12) << LET_square / (nO << 134 out<<std::setw(12)<<"LET"<<std::setw(12)<<LET_mean
123 } << 135 <<std::setw(12)<<"LET_SD"<<std::setw(12)<<LET_square/(nOfEvent-1)<<'\n';
124 else { << 136 }else
125 out << std::setw(12) << "LET" << std::se << 137 {
126 << std::setw(12) << LET_square << '\ << 138 out<<std::setw(12)<<"LET"<<std::setw(12)<<LET_mean
>> 139 <<std::setw(12)<<"LET_SD"<<std::setw(12)<<LET_square<<'\n';
127 } 140 }
128 141
129 masterScorer->OutputAndClear(); 142 masterScorer->OutputAndClear();
130 143
131 out << '\n'; << 144 out<<'\n';
>> 145
132 } 146 }
133 else { << 147 else
134 G4cout << G4endl << "--------------------E << 148 {
135 << " The run has " << nofEvents << << 149 G4cout
>> 150 << G4endl
>> 151 << "--------------------End of Local Run------------------------"
>> 152 << G4endl
>> 153 << " The run has " << nofEvents << " events"
>> 154 << G4endl;
136 } 155 }
137 156
138 G4cout << " Total energy deposited in the wo << 157 G4cout
139 << " -------------------------------- << 158 << " Total energy deposited in the world volume : " << sumDose/eV << " eV"
>> 159 << G4endl
>> 160 << " ------------------------------------------------------------"
>> 161 << G4endl
>> 162 << G4endl;
140 } 163 }
141 164
142 //....oooOO0OOooo........oooOO0OOooo........oo 165 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
143 166