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 #include "FlashSensitiveDetector.hh" 27 28 #include "G4HCofThisEvent.hh" 29 #include "G4SDManager.hh" 30 31 FlashSensitiveDetector::FlashSensitiveDetector 32 { 33 34 } 35 36 FlashSensitiveDetector::~FlashSensitiveDetecto 37 { 38 39 } 40 41 42 G4bool FlashSensitiveDetector::ProcessHits(G4S 43 44 { 45 G4double edep = step -> GetTotalEnergyDepo 46 47 if (edep == 0) return false; 48 49 const G4VTouchable *touchable = step -> Ge 50 G4VPhysicalVolume *physVol = touchable -> 51 G4ThreeVector posDetector = physVol -> Get 52 53 G4Track *track = step -> GetTrack(); 54 55 auto material = track -> GetMaterial(); 56 auto density = material -> GetDensity(); 57 G4double cavityVolume = (track->GetVolume( 58 G4double massOfCavity = cavityVolume * den 59 G4double dose = edep / massOfCavity / gray 60 edep = edep / MeV; 61 62 63 64 G4AnalysisManager *fAnalysisManager = G4An 65 fAnalysisManager -> FillNtupleDColumn(0, 0 66 fAnalysisManager -> FillNtupleDColumn(0, 1 67 fAnalysisManager -> FillNtupleDColumn(0, 2 68 fAnalysisManager -> FillNtupleDColumn(0, 3 69 fAnalysisManager -> FillNtupleDColumn(0, 4 70 fAnalysisManager -> FillNtupleIColumn(0, 5 71 fAnalysisManager -> FillNtupleIColumn(0, 6 72 fAnalysisManager -> FillNtupleSColumn(0, 7 73 fAnalysisManager -> AddNtupleRow(0); 74 75 return true; 76 } 77