Geant4  10.02
G4HumanPhantomRunAction.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // Authors: S. Guatelli , M. G. Pia, INFN Genova and F. Ambroglini INFN Perugia, Italy
27 //
28 // Based on code developed by the undergraduate student G. Guerrieri
29 // Note: this is a preliminary beta-version of the code; an improved
30 // version will be distributed in the next Geant4 public release, compliant
31 // with the design in a forthcoming publication, and subject to a
32 // design and code review.
33 //
34 
37 #include "G4SystemOfUnits.hh"
38 #include "G4UnitsTable.hh"
39 #include "G4ios.hh"
40 #include "G4Run.hh"
41 
42 #ifdef ANALYSIS_USE
44  :analysisMan(analysis)
45  {
46  }
47 #else
49  {}
50 #endif
51 
53 {
54 }
55 
57 {
58  G4int run_number = aRun->GetRunID();
59  G4cout << "### Run " << run_number << " start." << G4endl;
60  energyTotal["logicalHead"]=0.;
61  energyTotal["logicalTrunk"]=0.;
62  energyTotal["logicalLeftLeg"]=0.;
63  energyTotal["logicalRightLeg"]=0.;
64  energyTotal["logicalSkull"]=0.;
65  energyTotal["logicalLeftArmBone"]=0.;
66  energyTotal["logicalRightArmBone"]=0.;
67  energyTotal["logicalUpperSpine"]=0.;
68  energyTotal["logicalMiddleLowerSpine"]=0.;
69  energyTotal["logicalPelvis"]=0.;
70  energyTotal["logicalRibCage"]=0.;
71  energyTotal["logicalLeftClavicle"]=0.;
72  energyTotal["logicalRightClavicle"]=0.;
73  energyTotal["logicalLeftLegBone"]=0.;
74  energyTotal["logicalRightLegBone"]=0.;
75  energyTotal["logicalLeftScapula"]=0.;
76  energyTotal["logicalRightScapula"]=0.;
77  energyTotal["logicalHeart"]=0.;
78  energyTotal["logicalThyroid"]=0.;
79  energyTotal["logicalThymus"]=0.;
80  energyTotal["logicalMaleGenitalia"]=0.;
81  energyTotal["logicalBrain"]=0.;
82  energyTotal["logicalStomach"]=0.;
83  energyTotal["logicalUpperLargeIntestine"]=0.;
84  energyTotal["logicalLowerLargeIntestine"]=0.;
85  energyTotal["logicalSmallIntestine"]=0;
86  energyTotal["logicalSpleen"]=0.;
87  energyTotal["logicalPancreas"]=0.;
88  energyTotal["logicalLeftKidney"]=0.;
89  energyTotal["logicalRightKidney"]=0.;
90  energyTotal["logicalUrinaryBladder"]=0.;
91  energyTotal["logicalUterus"]=0.;
92  energyTotal["logicalLeftLung"]=0.;
93  energyTotal["logicalRightLung"]=0.;
94  energyTotal["logicalLeftOvary"]=0.;
95  energyTotal["logicalRightOvary"]=0.;
96  energyTotal["logicalLeftTeste"]=0;
97  energyTotal["logicalRightTeste"]=0;
98  energyTotal["logicalLeftBreast"]=0.;
99  energyTotal["logicalRightBreast"]=0.;
100  energyTotal["logicalLeftAdrenal"]=0.;
101  energyTotal["logicalRightAdrenal"]=0.;
102 
103 #ifdef ANALYSIS_USE
104  // Create ROOT file, histograms and ntuple
105  analysisMan -> book();
106 #endif
107 }
108 
110 {
111  G4cout << "Number of events = " << aRun->GetNumberOfEvent() << G4endl;
113 
114 #ifdef ANALYSIS_USE
115 // Close the output ROOT file with the results
116  analysisMan -> save();
117 #endif
118 }
119 
120 
122  G4double energyDeposit)
123 
124 {
125  energyTotal[bName] += energyDeposit;
126 }
127 
129 {
130  std::map<std::string,G4double>::iterator i = energyTotal.begin();
131  std::map<std::string,G4double>::iterator end = energyTotal.end();
132 
133  G4double totalEnergyDepositInPhantom =0.;
134  G4int k=0;
135  while(i!=end)
136  {
137  G4String bodypart = i->first;
138  G4double energyDep = i->second;
139  // if(energyDep != 0.)
140  // {
141 
142  G4cout << "Energy Total in Run:" << bodypart << ", ID: " << k
143  << ", Energy Deposition (MeV): "
144  << energyDep/MeV
145  << G4endl;
146 
147 #ifdef ANALYSIS_USE
148  // Fill Ntuple
149  analysisMan -> FillNtupleWithEnergyDeposition(k, energyDep/MeV);
150 #endif
151  i++;
152  k++;
153 
154  totalEnergyDepositInPhantom += energyDep;
155  }
156 
157 
158  G4cout << "Total Energy deposit in the body is: "
159  << totalEnergyDepositInPhantom/MeV << " MeV" <<G4endl;
160 }
static const double MeV
Definition: G4SIunits.hh:211
G4int first(char) const
std::map< std::string, G4double > energyTotal
int G4int
Definition: G4Types.hh:78
void Fill(G4String bodypartName, G4double energyDeposit)
void BeginOfRunAction(const G4Run *)
G4GLOB_DLL std::ostream G4cout
G4int GetNumberOfEvent() const
Definition: G4Run.hh:79
G4int GetRunID() const
Definition: G4Run.hh:76
Definition: G4Run.hh:46
void EndOfRunAction(const G4Run *)
#define G4endl
Definition: G4ios.hh:61
G4HumanPhantomRunAction(G4HumanPhantomAnalysisManager *analysis)
double G4double
Definition: G4Types.hh:76