Geant4  10.03
G4HumanPhantomEventAction.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 //
35 #include "G4HumanPhantomHit.hh"
36 #include "G4Event.hh"
37 #include "G4EventManager.hh"
38 #include "G4ios.hh"
39 #include "G4SDManager.hh"
40 #include "G4UnitsTable.hh"
42 #include "G4RunManager.hh"
43 
45 hitCollectionID(-1)
46 {
47 
48 }
49 
51 {
52 }
53 
55 {
56  energyTotal["logicalHead"]=0.;
57  energyTotal["logicalTrunk"]=0.;
58  energyTotal["logicalLeftLeg"]=0.;
59  energyTotal["logicalRightLeg"]=0.;
60  energyTotal["logicalSkull"]=0.;
61  energyTotal["logicalLeftArmBone"]=0.;
62  energyTotal["logicalRightArmBone"]=0.;
63  energyTotal["logicalUpperSpine"]=0.;
64  energyTotal["logicalMiddleLowerSpine"]=0.;
65  energyTotal["logicalPelvis"]=0.;
66  energyTotal["logicalRibCage"]=0.;
67  energyTotal["logicalLeftClavicle"]=0.;
68  energyTotal["logicalRightClavicle"]=0.;
69  energyTotal["logicalLeftLegBone"]=0.;
70  energyTotal["logicalRightLegBone"]=0.;
71  energyTotal["logicalLeftScapula"]=0.;
72  energyTotal["logicalRightScapula"]=0.;
73  energyTotal["logicalHeart"]=0.;
74  energyTotal["logicalThyroid"]=0.;
75  energyTotal["logicalThymus"]=0.;
76  energyTotal["logicalMaleGenitalia"]=0.;
77  energyTotal["logicalBrain"]=0.;
78  energyTotal["logicalStomach"]=0.;
79  energyTotal["logicalUpperLargeIntestine"]=0.;
80  energyTotal["logicalLowerLargeIntestine"]=0.;
81  energyTotal["logicalSmallIntestine"]=0;
82  energyTotal["logicalSpleen"]=0.;
83  energyTotal["logicalPancreas"]=0.;
84  energyTotal["logicalLeftKidney"]=0.;
85  energyTotal["logicalRightKidney"]=0.;
86  energyTotal["logicalUrinaryBladder"]=0.;
87  energyTotal["logicalUterus"]=0.;
88  energyTotal["logicalLeftLung"]=0.;
89  energyTotal["logicalRightLung"]=0.;
90  energyTotal["logicalLeftOvary"]=0.;
91  energyTotal["logicalRightOvary"]=0.;
92  energyTotal["logicalLeftTeste"]=0;
93  energyTotal["logicalRightTeste"]=0;
94  energyTotal["logicalLeftBreast"]=0.;
95  energyTotal["logicalRightBreast"]=0.;
96  energyTotal["logicalLeftAdrenal"]=0.;
97  energyTotal["logicalRightAdrenal"]=0.;
98 
100 
101  if (hitCollectionID==-1) {
102  hitCollectionID = SDman->GetCollectionID("HumanPhantomCollection");
103  }
104 }
105 
107 {
108 
109  G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
110 
112 
113  if (HCE)
115 
116  if (HC)
117  {
118  G4int hitNumber = HC->entries();
119  G4double edep =0;
120  G4String bodyPart;
121  for (G4int i=0;i<hitNumber;i++)
122  {
123  edep = (*HC)[i]->GetEdep();
124  bodyPart = (*HC)[i]->GetBodyPartID();
125  Fill(bodyPart, edep);
126  // if(edep !=0.) G4cout << bodyPart <<": "<< edep/MeV << G4endl;
127  }
128  }
129 
131 }
132 
134  G4double energyDeposit)
135 
136 {
137  energyTotal[bName] += energyDeposit;
138 }
139 
141 {
142 
144  G4HumanPhantomRunAction* pointerRun = (G4HumanPhantomRunAction*)(runManager->GetUserRunAction());
145 
146  std::map<std::string,G4double>::iterator i = energyTotal.begin();
147  std::map<std::string,G4double>::iterator end = energyTotal.end();
148 
149  while(i!=end)
150  {
151 
152  G4String bodypart = i->first;
153  G4double energyDep = i->second;
154 
155  if(energyDep != 0.)
156  {
157  pointerRun->Fill(bodypart, energyDep);
158  }
159  i++;
160  }
161 
162 }
G4VHitsCollection * GetHC(G4int i)
G4int GetCollectionID(G4String colName)
Definition: G4SDManager.cc:135
G4int first(char) const
G4int entries() const
int G4int
Definition: G4Types.hh:78
void Fill(G4String bodypartName, G4double energyDeposit)
void BeginOfEventAction(const G4Event *)
std::map< std::string, G4double > energyTotal
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:79
static G4SDManager * GetSDMpointer()
Definition: G4SDManager.cc:40
void Fill(G4String bodypartName, G4double energyDeposit)
G4HCofThisEvent * GetHCofThisEvent() const
Definition: G4Event.hh:185
const G4UserRunAction * GetUserRunAction() const
double G4double
Definition: G4Types.hh:76
void EndOfEventAction(const G4Event *)