Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // Previous authors: G. Guerrieri, S. Guatelli 27 // Authors (since 2007): S. Guatelli, Universi 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::G4HumanPhantomRunActi 38 { } 39 40 G4HumanPhantomRunAction::~G4HumanPhantomRunAct 41 {} 42 43 void G4HumanPhantomRunAction::BeginOfRunAction 44 { 45 G4int run_number = aRun->GetRunID(); 46 G4cout << "### Run " << run_number << " start 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(c 95 { 96 G4cout << "Number of events = " << aRun->Get 97 totalRunEnergyDeposit(); 98 99 // Close the output ROOT file with the results 100 fAnalysisMan -> save(); 101 } 102 103 void G4HumanPhantomRunAction::Fill(G4String bN 104 G4double energyDeposit) 105 { 106 fEnergyTotal[bName] += energyDeposit; 107 } 108 109 void G4HumanPhantomRunAction::totalRunEnergyDe 110 { 111 std::map<std::string,G4double>::iterator i = 112 std::map<std::string,G4double>::iterator 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:" << body 124 << ", Energy Deposition (MeV): " 125 << energyDep/MeV 126 << G4endl; 127 128 // Fill Ntuple 129 fAnalysisMan -> FillNtupleWithEnergyDepos 130 131 i++; 132 k++; 133 134 totalEnergyDepositInPhantom += energyDep; 135 } 136 137 G4cout << "Total Energy deposit in the body i 138 << totalEnergyDepositInPhantom/MeV << " MeV" 139 } 140