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 // Hadrontherapy advanced example for Geant4 27 // See more at: https://twiki.cern.ch/twiki/bi 28 29 #include "G4SystemOfUnits.hh" 30 #include "G4Event.hh" 31 #include "G4EventManager.hh" 32 #include "G4HCofThisEvent.hh" 33 #include "G4VHitsCollection.hh" 34 #include "G4SDManager.hh" 35 #include "G4VVisManager.hh" 36 37 #include "HadrontherapyEventAction.hh" 38 #include "HadrontherapyDetectorHit.hh" 39 #include "HadrontherapyDetectorSD.hh" 40 #include "HadrontherapyDetectorConstruction.hh 41 #include "HadrontherapyMatrix.hh" 42 #include "HadrontherapyEventActionMessenger.hh 43 44 ////////////////////////////////////////////// 45 HadrontherapyEventAction::HadrontherapyEventAc 46 drawFlag("all" ),printModulo(10), pointerEvent 47 { 48 hitsCollectionID = -1; 49 pointerEventMessenger = new HadrontherapyE 50 } 51 52 ////////////////////////////////////////////// 53 HadrontherapyEventAction::~HadrontherapyEventA 54 { 55 delete pointerEventMessenger; 56 } 57 58 ////////////////////////////////////////////// 59 void HadrontherapyEventAction::BeginOfEventAct 60 { 61 G4SDManager* pSDManager = G4SDManager::Get 62 if(hitsCollectionID == -1) 63 hitsCollectionID = pSDManager -> GetCo 64 } 65 66 ////////////////////////////////////////////// 67 void HadrontherapyEventAction::EndOfEventActio 68 { 69 if(hitsCollectionID < 0) 70 return; 71 G4HCofThisEvent* HCE = evt -> GetHCofThisE 72 73 // Clear voxels hit list 74 HadrontherapyMatrix* matrix = Hadrontherap 75 if (matrix) matrix -> ClearHitTrack(); 76 77 if(HCE) 78 { 79 HadrontherapyDetectorHitsCollection* C 80 if(CHC) 81 { 82 if(matrix) 83 { 84 // Fill the matrix with the in 85 // in the detector at the end 86 87 size_t HitCount = CHC -> entri 88 for (size_t h=0; h<HitCount; h 89 { 90 G4int i = ((*CHC)[h]) -> G 91 G4int j = ((*CHC)[h]) -> G 92 G4int k = ((*CHC)[h]) -> G 93 G4double energyDeposit = ( 94 matrix -> Fill(i, j, k, en 95 } 96 } 97 } 98 } 99 } 100