ICRP145_HumanPhantoms/src/TETRun.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // TETRun.cc
0027 //
0028 // Author: Haegin Han
0029 // Reference: ICRP Publication 145. Ann. ICRP 49(3), 2020.
0030 // Geant4 Contributors: J. Allison and S. Guatelli
0031 //
0032
0033 #include "TETRun.hh"
0034
0035 TETRun::TETRun()
0036 :G4Run(), fCollID(-1)
0037 {}
0038
0039 TETRun::~TETRun()
0040 {
0041  fEdepMap.clear();
0042 }
0043
0044 void TETRun::RecordEvent(const G4Event* event)
0045 {
0046  if(fCollID<0)
0047   fCollID = G4SDManager::GetSDMpointer()->GetCollectionID("PhantomSD/eDep");
0048
0049  // Hits collections
0050  //
0051  G4HCofThisEvent* HCE = event->GetHCofThisEvent();
0052  if(!HCE) return;
0053
0054  auto* evtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID));
0055  // sum up the energy deposition and the square of it
0056  for (auto itr : *evtMap->GetMap()) {
0057         fEdepMap[itr.first].first  += *itr.second;                   //sum
0058         fEdepMap[itr.first].second += (*itr.second) * (*itr.second); //sum square
0059     }
0060 }
0061
0062 void TETRun::Merge(const G4Run* run)
0063 {
0064  // merge the data from each thread
0065  EDEPMAP localMap = static_cast<const TETRun*>(run)->fEdepMap;
0066
0067  for(auto itr : localMap){
0068      fEdepMap[itr.first].first  += itr.second.first;
0069          fEdepMap[itr.first].second += itr.second.second;
0070     }
0071
0072  G4Run::Merge(run);
0073 }
0074
0075
0076
0077
0078