Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // TETRun.cc 26 // TETRun.cc 27 // 27 // 28 // Author: Haegin Han 28 // Author: Haegin Han 29 // Reference: ICRP Publication 145. Ann. ICRP 29 // Reference: ICRP Publication 145. Ann. ICRP 49(3), 2020. 30 // Geant4 Contributors: J. Allison and S. Guat 30 // Geant4 Contributors: J. Allison and S. Guatelli 31 // 31 // 32 32 33 #include "TETRun.hh" 33 #include "TETRun.hh" 34 34 35 TETRun::TETRun() 35 TETRun::TETRun() 36 :G4Run(), fCollID(-1) << 36 :G4Run() 37 {} 37 {} 38 38 39 TETRun::~TETRun() 39 TETRun::~TETRun() 40 { 40 { 41 fEdepMap.clear(); 41 fEdepMap.clear(); 42 } 42 } 43 43 44 void TETRun::RecordEvent(const G4Event* event) 44 void TETRun::RecordEvent(const G4Event* event) 45 { 45 { 46 if(fCollID<0) << 46 auto fCollID 47 fCollID = G4SDManager::GetSDMpointer()->GetC << 47 = G4SDManager::GetSDMpointer()->GetCollectionID("PhantomSD/eDep"); 48 48 49 // Hits collections << 49 // Hits collections 50 // << 50 // 51 G4HCofThisEvent* HCE = event->GetHCofThisEven << 51 G4HCofThisEvent* HCE = event->GetHCofThisEvent(); 52 if(!HCE) return; << 52 if(!HCE) return; 53 << 53 54 auto* evtMap = static_cast<G4THitsMap<G4doubl << 54 //G4THitsMap<G4double>* evtMap = 55 // sum up the energy deposition and the squar << 55 // static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID)); 56 for (auto itr : *evtMap->GetMap()) { << 56 auto* evtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID)); >> 57 // sum up the energy deposition and the square of it >> 58 for (auto itr : *evtMap->GetMap()) { 57 fEdepMap[itr.first].first += *itr.second; 59 fEdepMap[itr.first].first += *itr.second; //sum 58 fEdepMap[itr.first].second += (*itr.second 60 fEdepMap[itr.first].second += (*itr.second) * (*itr.second); //sum square 59 } 61 } 60 } 62 } 61 63 62 void TETRun::Merge(const G4Run* run) 64 void TETRun::Merge(const G4Run* run) 63 { 65 { 64 // merge the data from each thread 66 // merge the data from each thread 65 EDEPMAP localMap = static_cast<const TETRun*> 67 EDEPMAP localMap = static_cast<const TETRun*>(run)->fEdepMap; 66 68 67 for(auto itr : localMap){ 69 for(auto itr : localMap){ 68 fEdepMap[itr.first].first += itr.second.fi 70 fEdepMap[itr.first].first += itr.second.first; 69 fEdepMap[itr.first].second += itr.sec 71 fEdepMap[itr.first].second += itr.second.second; 70 } 72 } 71 73 72 G4Run::Merge(run); << 74 G4Run::Merge(run); 73 } 75 } 74 76 75 77 76 78 77 79 78 80 79 81