Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/human_phantom/src/G4HumanPhantomEventAction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/advanced/human_phantom/src/G4HumanPhantomEventAction.cc (Version 11.3.0) and /examples/advanced/human_phantom/src/G4HumanPhantomEventAction.cc (Version 9.6.p1)


  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 // Previous authors: G. Guerrieri, S. Guatelli <<  26 // Authors: S. Guatelli and M. G. Pia, INFN Genova, Italy
 27 // Authors (since 2007): S. Guatelli, Universi <<  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.
 28 //                                                 33 // 
 29 #include "G4HumanPhantomEventAction.hh"            34 #include "G4HumanPhantomEventAction.hh"
 30 #include "G4HumanPhantomHit.hh"                    35 #include "G4HumanPhantomHit.hh"
 31 #include "G4Event.hh"                              36 #include "G4Event.hh"
 32 #include "G4EventManager.hh"                       37 #include "G4EventManager.hh"
 33 #include "G4ios.hh"                                38 #include "G4ios.hh"
 34 #include "G4SDManager.hh"                          39 #include "G4SDManager.hh"
 35 #include "G4UnitsTable.hh"                         40 #include "G4UnitsTable.hh"
 36 #include "G4HumanPhantomRunAction.hh"              41 #include "G4HumanPhantomRunAction.hh"
 37 #include "G4RunManager.hh"                         42 #include "G4RunManager.hh"
 38                                                    43 
 39 G4HumanPhantomEventAction::G4HumanPhantomEvent     44 G4HumanPhantomEventAction::G4HumanPhantomEventAction():
 40 fHitCollectionID(-1)                           <<  45 hitCollectionID(-1)
 41 {                                                  46 { 
 42                                                    47  
 43 }                                                  48 }
                                                   >>  49  
                                                   >>  50 G4HumanPhantomEventAction::~G4HumanPhantomEventAction()
                                                   >>  51 {
                                                   >>  52 }
 44                                                    53 
 45 void G4HumanPhantomEventAction::BeginOfEventAc     54 void G4HumanPhantomEventAction::BeginOfEventAction(const G4Event*)
 46 {                                                  55 {
 47  fEnergyTotal["logicalHead"]=0.;               <<  56  energyTotal["logicalHead"]=0.;
 48  fEnergyTotal["logicalTrunk"]=0.;              <<  57  energyTotal["logicalTrunk"]=0.;
 49  fEnergyTotal["logicalLeftLeg"]=0.;            <<  58  energyTotal["logicalLeftLeg"]=0.;
 50  fEnergyTotal["logicalRightLeg"]=0.;           <<  59  energyTotal["logicalRightLeg"]=0.;
 51  fEnergyTotal["logicalSkull"]=0.;              <<  60  energyTotal["logicalBrain"]=0.;
 52  fEnergyTotal["logicalLeftArmBone"]=0.;        <<  61  energyTotal["logicalLeftArmBone"]=0.;
 53  fEnergyTotal["logicalRightArmBone"]=0.;       <<  62  energyTotal["logicalRightArmBone"]=0.;
 54  fEnergyTotal["logicalUpperSpine"]=0.;         <<  63  energyTotal["logicalSkull"]=0.;
 55  fEnergyTotal["logicalMiddleLowerSpine"]=0.;   <<  64  energyTotal["logicalUpperSpine"]=0.;
 56  fEnergyTotal["logicalPelvis"]=0.;             <<  65  energyTotal["logicalMiddleLowerSpine"]=0.;
 57  fEnergyTotal["logicalRibCage"]=0.;            <<  66  energyTotal["logicalPelvis"]=0.;
 58  fEnergyTotal["logicalLeftClavicle"]=0.;       <<  67  energyTotal["logicalStomach"]=0.;
 59  fEnergyTotal["logicalRightClavicle"]=0.;      <<  68  energyTotal["logicalUpperLargeIntestine"]=0.;
 60  fEnergyTotal["logicalLeftLegBone"]=0.;        <<  69  energyTotal["logicalLowerLargeIntestine"]=0.;
 61  fEnergyTotal["logicalRightLegBone"]=0.;       <<  70  energyTotal["logicalRibCage"]=0.;
 62  fEnergyTotal["logicalLeftScapula"]=0.;        <<  71  energyTotal["logicalSpleen"]=0.;
 63  fEnergyTotal["logicalRightScapula"]=0.;       <<  72  energyTotal["logicalPancreas"]=0.;
 64  fEnergyTotal["logicalHeart"]=0.;              <<  73  energyTotal["logicalLeftKidney"]=0.;
 65  fEnergyTotal["logicalThyroid"]=0.;            <<  74  energyTotal["logicalRightKidney"]=0.;
 66  fEnergyTotal["logicalThymus"]=0.;             <<  75  energyTotal["logicalUrinaryBladder"]=0.;
 67  fEnergyTotal["logicalMaleGenitalia"]=0.;      <<  76  energyTotal["logicalUterus"]=0.;
 68  fEnergyTotal["logicalBrain"]=0.;              <<  77  energyTotal["logicalLeftLung"]=0.;
 69  fEnergyTotal["logicalStomach"]=0.;            <<  78  energyTotal["logicalRightLung"]=0.;
 70  fEnergyTotal["logicalUpperLargeIntestine"]=0. <<  79  energyTotal["logicalLeftOvary"]=0.;
 71  fEnergyTotal["logicalLowerLargeIntestine"]=0. <<  80  energyTotal["logicalRightOvary"]=0.;
 72  fEnergyTotal["logicalSmallIntestine"]=0;      <<  81  energyTotal["logicalLeftLegBone"]=0.;
 73  fEnergyTotal["logicalSpleen"]=0.;             <<  82  energyTotal["logicalRightLegBone"]=0.;
 74  fEnergyTotal["logicalPancreas"]=0.;           <<  83  energyTotal["logicalLeftBreast"]=0.;
 75  fEnergyTotal["logicalLeftKidney"]=0.;         <<  84  energyTotal["logicalRightBreast"]=0.; 
 76  fEnergyTotal["logicalRightKidney"]=0.;        <<  85  energyTotal["logicalLeftScapula"]=0.; 
 77  fEnergyTotal["logicalUrinaryBladder"]=0.;     <<  86  energyTotal["logicalRightScapula"]=0.;
 78  fEnergyTotal["logicalUterus"]=0.;             <<  87  energyTotal["logicalLeftAdrenal"]=0.; 
 79  fEnergyTotal["logicalLeftLung"]=0.;           <<  88  energyTotal["logicalRightAdrenal"]=0.;
 80  fEnergyTotal["logicalRightLung"]=0.;          << 
 81  fEnergyTotal["logicalLeftOvary"]=0.;          << 
 82  fEnergyTotal["logicalRightOvary"]=0.;         << 
 83  fEnergyTotal["logicalLeftTeste"]=0;           << 
 84  fEnergyTotal["logicalRightTeste"]=0;          << 
 85  fEnergyTotal["logicalLeftBreast"]=0.;         << 
 86  fEnergyTotal["logicalRightBreast"]=0.;        << 
 87  fEnergyTotal["logicalLeftAdrenal"]=0.;        << 
 88  fEnergyTotal["logicalRightAdrenal"]=0.;       << 
 89                                                << 
 90  G4SDManager * SDman = G4SDManager::GetSDMpoin     89  G4SDManager * SDman = G4SDManager::GetSDMpointer();  
 91                                                    90 
 92   if (fHitCollectionID==-1) {                  <<  91   if (hitCollectionID==-1) {
 93     fHitCollectionID = SDman->GetCollectionID( <<  92     hitCollectionID = SDman->GetCollectionID("HumanPhantomCollection");
 94   }                                                93   }
 95 }                                                  94 }
 96                                                    95  
 97 void G4HumanPhantomEventAction::EndOfEventActi     96 void G4HumanPhantomEventAction::EndOfEventAction(const G4Event* evt)
 98 {                                                  97 {  
 99                                                    98 
100  G4HCofThisEvent* HCE = evt->GetHCofThisEvent(     99  G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
101                                                   100  
102  G4HumanPhantomHitsCollection* HC = nullptr;   << 101  G4HumanPhantomHitsCollection* HC = 0;
103                                                   102 
104  if (HCE)                                         103  if (HCE)
105      HC = (G4HumanPhantomHitsCollection*)(HCE- << 104      HC = (G4HumanPhantomHitsCollection*)(HCE->GetHC(hitCollectionID));
106                                                   105 
107  if (HC)                                          106  if (HC)
108   {                                               107   {
109     G4int hitNumber = HC->entries();              108     G4int hitNumber = HC->entries();
110           G4double edep =0;                       109           G4double edep =0;
111           G4String bodyPart;                      110           G4String bodyPart;
112           for (G4int i=0;i<hitNumber;i++)         111           for (G4int i=0;i<hitNumber;i++) 
113       {                                           112       {
114         edep = (*HC)[i]->GetEdep();               113         edep = (*HC)[i]->GetEdep();
115         bodyPart = (*HC)[i]->GetBodyPartID();     114         bodyPart = (*HC)[i]->GetBodyPartID();
116               Fill(bodyPart, edep);               115               Fill(bodyPart, edep);
117         //        if(edep !=0.) G4cout << body    116         //        if(edep !=0.) G4cout << bodyPart <<": "<< edep/MeV << G4endl; 
118       }                                           117       }
119   }                                               118   }
120                                                   119 
121  totalEventEnergyDeposit();                       120  totalEventEnergyDeposit();
122 }                                                 121 }
123                                                   122 
124 void G4HumanPhantomEventAction:: Fill(G4String    123 void G4HumanPhantomEventAction:: Fill(G4String bName, 
125               G4double energyDeposit)             124               G4double energyDeposit)
126                                                   125 
127 {                                                 126 {
128  fEnergyTotal[bName] += energyDeposit;         << 127  energyTotal[bName] += energyDeposit;
129 }                                                 128 }
130                                                   129 
131 void G4HumanPhantomEventAction::totalEventEner    130 void G4HumanPhantomEventAction::totalEventEnergyDeposit() 
132 {                                                 131 {
133                                                   132 
134  G4RunManager* runManager = G4RunManager::GetR    133  G4RunManager* runManager = G4RunManager::GetRunManager();
135  auto* pointerRun = (G4HumanPhantomRunAction*) << 134  G4HumanPhantomRunAction* pointerRun = (G4HumanPhantomRunAction*)(runManager->GetUserRunAction());
136                                                   135 
137  std::map<std::string,G4double>::iterator i =  << 136  std::map<std::string,G4double>::iterator i = energyTotal.begin();
138   std::map<std::string,G4double>::iterator end << 137   std::map<std::string,G4double>::iterator end = energyTotal.end();
139                                                   138 
140   while(i!=end)                                   139   while(i!=end)
141     {                                             140     {
142                                                   141 
143       G4String bodypart = i->first;               142       G4String bodypart = i->first;
144       G4double energyDep = i->second;             143       G4double energyDep = i->second;
145                                                   144       
146       if(energyDep != 0.)                         145       if(energyDep != 0.)
147   {                                               146   {
148     pointerRun->Fill(bodypart, energyDep);        147     pointerRun->Fill(bodypart, energyDep);
149   }                                               148   }
150       i++;                                        149       i++;
151     }                                             150     }
152                                                   151   
153 }                                                 152 }
154                                                   153