Geant4 Cross Reference |
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 // 27 /// file: 28 /// brief: 29 #include "EventAction.hh" 30 31 #include "AnalysisManager.hh" 32 #include "ChromosomeHit.hh" 33 #include "ChromosomeMapper.hh" 34 #include "DNAGeometry.hh" 35 #include "DNAHit.hh" 36 #include "DetectorConstruction.hh" 37 38 #include "G4RunManager.hh" 39 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 41 42 EventAction::EventAction(AnalysisManager* man) : fAnalysisManager(man) {} 43 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 void EventAction::Initialize() 47 { 48 // Get DNA Geometry 49 const auto* det = dynamic_cast<const DetectorConstruction*>( 50 G4RunManager::GetRunManager()->GetUserDetectorConstruction()); 51 fDNAGeometry = det->GetDNAGeometry(); 52 53 // Prepare Chromo Hits Map 54 std::vector<G4String> keys = fDNAGeometry->GetChromosomeMapper()->GetChromosomeKeys(); 55 for (const auto& key : keys) { 56 uint32_t key_i = G4::hashing::crc32::Hash(key); 57 fChromoHitMap[key_i] = nullptr; 58 } 59 fInitialized = true; 60 } 61 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 63 64 void EventAction::BeginOfEventAction(const G4Event*) 65 { 66 if (!fInitialized) { 67 Initialize(); 68 } 69 70 fEdepCell = 0; 71 fTraLenCell = 0; 72 fTraLenChro = 0; 73 74 // Prepare DNA Hits vector 75 fDNAHits.reserve(10000); 76 77 std::vector<G4String> keys = fDNAGeometry->GetChromosomeMapper()->GetChromosomeKeys(); 78 // ChromosomeHit* theHit; 79 for (const auto& it : keys) { 80 uint32_t key_i = G4::hashing::crc32::Hash(it); 81 fChromoHitMap[key_i] = new ChromosomeHit(it); 82 } 83 } 84 85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 86 87 void EventAction::EndOfEventAction(const G4Event*) 88 { 89 // post-analysis 90 fAnalysisManager->ProcessPrimary(fprimstoppos, fTraLenCell, fTraLenChro); 91 fAnalysisManager->ProcessDNAHitsVector(fDNAHits); 92 fAnalysisManager->ProcessChromosomeHitMap(fChromoHitMap); 93 fAnalysisManager->ProcessCellEdep(fEdepCell); 94 95 // destruction ///////////////////////////////////////////////////////////// 96 for (auto& it : fChromoHitMap) { 97 delete it.second; 98 } 99 100 for (auto& fDNAHit : fDNAHits) { 101 delete fDNAHit; 102 } 103 fDNAHits.clear(); 104 } 105 106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 107 108 void EventAction::AddChromosomeEdep(const G4String& key, const G4double& e, const G4bool& isDNA) 109 { 110 uint32_t key_i = G4::hashing::crc32::Hash(key); 111 auto it = fChromoHitMap.find(key_i); 112 if (it != fChromoHitMap.end()) { 113 if (isDNA) { 114 it->second->AddDNAEdep(e); 115 } 116 else { 117 it->second->AddChromosomeEdep(e); 118 } 119 } 120 else { 121 G4ExceptionDescription errmsg; 122 errmsg << "Energy deposit in unknown chromosome: " << key << G4endl; 123 G4Exception("EventAction::AddChromosomeEdep", "ERR_UNKNOWN_CHROMOSOME", FatalException, errmsg); 124 } 125 } 126 // 2600754154 127 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 128 129 void EventAction::AddCellEdep(const G4double& edep) 130 { 131 fEdepCell += edep; 132 } 133 134 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 135 136 void EventAction::AddTrackLengthCell(const G4double& tl) 137 { 138 fTraLenCell += tl; 139 } 140 141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 142 143 void EventAction::AddTrackLengthChro(const G4double& tl) 144 { 145 fTraLenChro += tl; 146 } 147 148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 149