| 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 // Hadrontherapy advanced example for Geant4 26 // Hadrontherapy advanced example for Geant4
27 // See more at: https://twiki.cern.ch/twiki/bi 27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
28 28
29 #include "HadrontherapyRunAction.hh" 29 #include "HadrontherapyRunAction.hh"
30 #include "HadrontherapyEventAction.hh" 30 #include "HadrontherapyEventAction.hh"
31 #include "G4AnalysisManager.hh" << 31 #include "HadrontherapyAnalysis.hh"
32 #include "G4Run.hh" 32 #include "G4Run.hh"
33 #include "G4RunManager.hh" 33 #include "G4RunManager.hh"
34 #include "G4UImanager.hh" 34 #include "G4UImanager.hh"
35 #include "G4ios.hh" 35 #include "G4ios.hh"
36 #include "HadrontherapyDetectorConstruction.hh 36 #include "HadrontherapyDetectorConstruction.hh"
37 #include "G4SDManager.hh" 37 #include "G4SDManager.hh"
38 #include "G4Timer.hh" 38 #include "G4Timer.hh"
39 #include "HadrontherapyRunAction.hh" 39 #include "HadrontherapyRunAction.hh"
40 #include "HadrontherapyMatrix.hh" 40 #include "HadrontherapyMatrix.hh"
41 #include "HadrontherapyRBE.hh" 41 #include "HadrontherapyRBE.hh"
42 #include "G4UnitsTable.hh" 42 #include "G4UnitsTable.hh"
43 #include "G4SystemOfUnits.hh" 43 #include "G4SystemOfUnits.hh"
44 44
45 #include <G4AccumulableManager.hh> 45 #include <G4AccumulableManager.hh>
46 46
47 ////////////////////////////////////////////// 47 /////////////////////////////////////////////////////////////////////////////
48 HadrontherapyRunAction::HadrontherapyRunAction 48 HadrontherapyRunAction::HadrontherapyRunAction()
49 { 49 {
50 G4AccumulableManager* accumulableManager = 50 G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
51 accumulableManager->Register(&fRBEAccumula << 51 accumulableManager->RegisterAccumulable(&fRBEAccumulable);
>> 52
>> 53 // Create analysis manager
>> 54 // The choice of analysis technology is done via selectin of a namespace
>> 55 // in Analysis.hh
>> 56 auto analysisManager =G4AnalysisManager::Instance();
>> 57 G4cout << "Using " << analysisManager -> GetType() << G4endl;
>> 58
>> 59 analysisManager->SetVerboseLevel(1);
>> 60 analysisManager->SetFirstHistoId(1);
>> 61
>> 62 // Comment out the following line to generate an N-tuple
>> 63 analysisManager-> SetFirstNtupleId(2);
>> 64
>> 65 // Creating the histograms of primary kinetic
>> 66 // energy (Ekin) and of the energy deposited (Edep)
>> 67 // in the first voxel/slice of the water phantom
>> 68 analysisManager -> CreateH1("Ekin","Ekin the voxel", 400,20*MeV, 60*MeV);
>> 69 analysisManager -> CreateH1("Edep","Edep the voxel", 200, -10, 10*MeV);
>> 70
>> 71 // Example of how to create an Ntuple (comment-out, if needed)
>> 72 //analysisManager->CreateNtuple("NYUPLA", "Edep and TrackL");
>> 73 //analysisManager->CreateNtupleDColumn("Ekin");
52 74
53 } 75 }
54 76
55 ////////////////////////////////////////////// 77 /////////////////////////////////////////////////////////////////////////////
56 HadrontherapyRunAction::~HadrontherapyRunActio 78 HadrontherapyRunAction::~HadrontherapyRunAction()
57 { 79 {
>> 80 delete G4AnalysisManager::Instance();
58 } 81 }
59 82
60 ////////////////////////////////////////////// 83 /////////////////////////////////////////////////////////////////////////////
61 void HadrontherapyRunAction::BeginOfRunAction( 84 void HadrontherapyRunAction::BeginOfRunAction(const G4Run* aRun)
62 { << 85 {
>> 86
>> 87 // Get analysis manager
>> 88 auto analysisManager = G4AnalysisManager::Instance();
>> 89
>> 90 // Open an output file
>> 91 //
>> 92 G4String fileName = "Hadrontherapy";
>> 93 analysisManager->OpenFile(fileName);
>> 94
63 G4AccumulableManager* accumulableManager = 95 G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
64 accumulableManager->Reset(); 96 accumulableManager->Reset();
65 97
66 G4RunManager::GetRunManager()-> SetRandomN 98 G4RunManager::GetRunManager()-> SetRandomNumberStore(true);
67 G4cout << "Run " << aRun -> GetRunID() << 99 G4cout << "Run " << aRun -> GetRunID() << " starts ..." << G4endl;
68 100
69 HadrontherapyRBE *rbe = HadrontherapyRBE:: 101 HadrontherapyRBE *rbe = HadrontherapyRBE::GetInstance();
70 if (rbe->IsCalculationEnabled() && IsMaste 102 if (rbe->IsCalculationEnabled() && IsMaster() && rbe->GetVerboseLevel() > 0)
71 { 103 {
72 rbe->PrintParameters(); 104 rbe->PrintParameters();
73 } 105 }
74 106
75 electromagnetic = 0; 107 electromagnetic = 0;
76 hadronic = 0; 108 hadronic = 0;
77 } 109 }
78 110
79 ////////////////////////////////////////////// 111 /////////////////////////////////////////////////////////////////////////////
80 void HadrontherapyRunAction::EndOfRunAction(co 112 void HadrontherapyRunAction::EndOfRunAction(const G4Run*)
81 { 113 {
82 auto analysisManager = G4AnalysisManager:: 114 auto analysisManager = G4AnalysisManager::Instance();
83 << 115
>> 116 //G4cout << " Summary of Run " << aRun -> GetRunID() <<" :"<< G4endl;
>> 117 //G4cout << "Number of electromagnetic processes of primary particles in the phantom:"
>> 118 // << electromagnetic << G4endl;
>> 119 //G4cout << "Number of hadronic processes of primary particles in the phantom:"
>> 120 // << hadronic << G4endl;
84 G4AccumulableManager* accumulableManager = 121 G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
85 accumulableManager->Merge(); 122 accumulableManager->Merge();
86 123
87 // Tell the RBE class what we have accumul 124 // Tell the RBE class what we have accumulated...
88 HadrontherapyRBE *rbe = HadrontherapyRBE:: 125 HadrontherapyRBE *rbe = HadrontherapyRBE::GetInstance();
89 if (rbe->IsCalculationEnabled() && IsMaste 126 if (rbe->IsCalculationEnabled() && IsMaster())
90 { 127 {
91 if (rbe->IsAccumulationEnabled()) 128 if (rbe->IsAccumulationEnabled())
92 { 129 {
93 rbe->AddAlphaNumerator(fRBEAccumul 130 rbe->AddAlphaNumerator(fRBEAccumulable.GetAlphaNumerator());
94 rbe->AddBetaNumerator(fRBEAccumula 131 rbe->AddBetaNumerator(fRBEAccumulable.GetBetaNumerator());
95 rbe->AddDenominator(fRBEAccumulabl 132 rbe->AddDenominator(fRBEAccumulable.GetDenominator());
96 rbe->AddEnergyDeposit(fRBEAccumula 133 rbe->AddEnergyDeposit(fRBEAccumulable.GetEnergyDeposit());
97 } 134 }
98 else 135 else
99 { 136 {
100 rbe->SetAlphaNumerator(fRBEAccumul 137 rbe->SetAlphaNumerator(fRBEAccumulable.GetAlphaNumerator());
101 rbe->SetBetaNumerator(fRBEAccumula 138 rbe->SetBetaNumerator(fRBEAccumulable.GetBetaNumerator());
102 rbe->SetDenominator(fRBEAccumulabl 139 rbe->SetDenominator(fRBEAccumulable.GetDenominator());
103 rbe->SetEnergyDeposit(fRBEAccumula 140 rbe->SetEnergyDeposit(fRBEAccumulable.GetEnergyDeposit());
104 } 141 }
105 142
106 rbe->StoreAlphaAndBeta(); 143 rbe->StoreAlphaAndBeta();
107 rbe->StoreRBE(); 144 rbe->StoreRBE();
108 } 145 }
109 146
110 if (analysisManager->IsOpenFile()) { << 147 analysisManager->Write();
111 analysisManager->Write(); << 148 analysisManager->CloseFile();
112 analysisManager->CloseFile(); << 149
113 } <<
114 } 150 }
115 ////////////////////////////////////////////// 151 /////////////////////////////////////////////////////////////////////////////
116 void HadrontherapyRunAction::AddEMProcess() 152 void HadrontherapyRunAction::AddEMProcess()
117 { 153 {
118 electromagnetic += 1; 154 electromagnetic += 1;
119 } 155 }
120 156
121 ////////////////////////////////////////////// 157 /////////////////////////////////////////////////////////////////////////////
122 void HadrontherapyRunAction::AddHadronicProces 158 void HadrontherapyRunAction::AddHadronicProcess()
123 { 159 {
124 hadronic += 1; 160 hadronic += 1;
125 } 161 }
126 162