Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/human_phantom/src/G4HumanPhantomRunAction.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 ]

  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 // Previous authors: G. Guerrieri, S. Guatelli and M. G. Pia, INFN Genova, Italy
 27 // Authors (since 2007): S. Guatelli, University of Wollongong, Australia
 28 //
 29 
 30 #include "G4HumanPhantomRunAction.hh"
 31 #include "G4HumanPhantomAnalysisManager.hh"
 32 #include "G4SystemOfUnits.hh"
 33 #include "G4UnitsTable.hh"
 34 #include "G4ios.hh"
 35 #include "G4Run.hh"
 36 
 37 G4HumanPhantomRunAction::G4HumanPhantomRunAction(G4HumanPhantomAnalysisManager* analysis) :fAnalysisMan(analysis)
 38  { }
 39 
 40 G4HumanPhantomRunAction::~G4HumanPhantomRunAction()
 41 {}
 42 
 43 void G4HumanPhantomRunAction::BeginOfRunAction(const G4Run* aRun)
 44 {
 45  G4int run_number = aRun->GetRunID();
 46  G4cout << "### Run " << run_number << " start." << G4endl;
 47  fEnergyTotal["logicalHead"]=0.;
 48  fEnergyTotal["logicalTrunk"]=0.;
 49  fEnergyTotal["logicalLeftLeg"]=0.;
 50  fEnergyTotal["logicalRightLeg"]=0.;
 51  fEnergyTotal["logicalSkull"]=0.;
 52  fEnergyTotal["logicalLeftArmBone"]=0.;
 53  fEnergyTotal["logicalRightArmBone"]=0.;
 54  fEnergyTotal["logicalUpperSpine"]=0.;
 55  fEnergyTotal["logicalMiddleLowerSpine"]=0.;
 56  fEnergyTotal["logicalPelvis"]=0.;
 57  fEnergyTotal["logicalRibCage"]=0.;
 58  fEnergyTotal["logicalLeftClavicle"]=0.;
 59  fEnergyTotal["logicalRightClavicle"]=0.;
 60  fEnergyTotal["logicalLeftLegBone"]=0.;
 61  fEnergyTotal["logicalRightLegBone"]=0.;
 62  fEnergyTotal["logicalLeftScapula"]=0.; 
 63  fEnergyTotal["logicalRightScapula"]=0.;
 64  fEnergyTotal["logicalHeart"]=0.;
 65  fEnergyTotal["logicalThyroid"]=0.;
 66  fEnergyTotal["logicalThymus"]=0.;
 67  fEnergyTotal["logicalMaleGenitalia"]=0.;
 68  fEnergyTotal["logicalBrain"]=0.;
 69  fEnergyTotal["logicalStomach"]=0.;
 70  fEnergyTotal["logicalUpperLargeIntestine"]=0.;
 71  fEnergyTotal["logicalLowerLargeIntestine"]=0.;
 72  fEnergyTotal["logicalSmallIntestine"]=0;
 73  fEnergyTotal["logicalSpleen"]=0.;
 74  fEnergyTotal["logicalPancreas"]=0.;
 75  fEnergyTotal["logicalLeftKidney"]=0.;
 76  fEnergyTotal["logicalRightKidney"]=0.;
 77  fEnergyTotal["logicalUrinaryBladder"]=0.;
 78  fEnergyTotal["logicalUterus"]=0.;
 79  fEnergyTotal["logicalLeftLung"]=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  // Create ROOT file, histograms and ntuple
 91  fAnalysisMan -> book();
 92 }
 93 
 94 void G4HumanPhantomRunAction::EndOfRunAction(const G4Run* aRun)
 95 {
 96   G4cout << "Number of events = " << aRun->GetNumberOfEvent() << G4endl;
 97   totalRunEnergyDeposit();
 98 
 99 // Close the output ROOT file with the results
100    fAnalysisMan -> save(); 
101 }
102 
103 void G4HumanPhantomRunAction::Fill(G4String bName, 
104               G4double energyDeposit)
105 {
106  fEnergyTotal[bName] += energyDeposit;
107 }
108 
109 void G4HumanPhantomRunAction::totalRunEnergyDeposit() 
110 {
111  std::map<std::string,G4double>::iterator i = fEnergyTotal.begin();
112   std::map<std::string,G4double>::iterator end = fEnergyTotal.end();
113 
114   G4double totalEnergyDepositInPhantom =0.;
115   G4int k=0;
116   while(i!=end)
117     {
118       G4String bodypart = i->first;
119       G4double energyDep = i->second;
120       //  if(energyDep != 0.)
121       //  {
122      
123       G4cout << "Energy Total in Run:" << bodypart << ", ID: "  << k 
124              << ", Energy Deposition (MeV): "  
125        << energyDep/MeV
126        << G4endl;
127 
128      // Fill Ntuple
129      fAnalysisMan -> FillNtupleWithEnergyDeposition(k, energyDep/MeV);
130 
131      i++;
132      k++;
133 
134      totalEnergyDepositInPhantom += energyDep;
135     }
136   
137  G4cout << "Total Energy deposit in the body is: " 
138   << totalEnergyDepositInPhantom/MeV << " MeV" <<G4endl; 
139 }
140