ICRP145_HumanPhantoms/src/TETRunAction.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 //  TETRunAction.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 "TETRunAction.hh"
0034
0035 TETRunAction::TETRunAction(TETModelImport* _tetData, G4String _output)
0036 :fTetData(_tetData), fRun(nullptr), fNumOfEvent(0), fRunID(0), fOutputFile(_output)
0037 {}
0038
0039 G4Run* TETRunAction::GenerateRun()
0040 {
0041  // generate run
0042  fRun = new TETRun();
0043  return fRun;
0044 }
0045
0046 void TETRunAction::BeginOfRunAction(const G4Run* aRun)
0047 {
0048  // print the progress at the interval of 10%
0049  fNumOfEvent=aRun->GetNumberOfEventToBeProcessed();
0050  G4RunManager::GetRunManager()->SetPrintProgress(int(fNumOfEvent*0.1));
0051 }
0052
0053 void TETRunAction::EndOfRunAction(const G4Run* aRun)
0054 {
0055  // print the result only in the Master
0056  if(!isMaster) return;
0057
0058  // get the run ID
0059  fRunID = aRun->GetRunID();
0060
0061  // Print the run result by G4cout and std::ofstream
0062  //
0063  // print by G4cout
0064  PrintResult(G4cout);
0065
0066  // print by std::ofstream
0067  std::ofstream ofs(fOutputFile.c_str());
0068  PrintResult(ofs);
0069  ofs.close();
0070 }
0071
0072 void TETRunAction::PrintResult(std::ostream &out)
0073 {
0074  // Print run result
0075  //
0076  using namespace std;
0077  EDEPMAP edepMap = *fRun->GetEdepMap();
0078
0079     out << G4endl
0080      << "=====================================================================" << G4endl
0081      << " Run #" << fRunID << " / Number of event processed : "<< fNumOfEvent    << G4endl
0082      << "=====================================================================" << G4endl
0083      << "organ ID| "
0084      << setw(19) << "Organ Mass (g)"
0085          << setw(19) << "Dose (Gy/source)"
0086      << setw(19) << "Relative Error" << G4endl;
0087
0088     out.precision(3);
0089     auto massMap = fTetData->GetMassMap();
0090     for(auto itr : massMap){
0091         G4double meanDose    = edepMap[itr.first].first  / itr.second / fNumOfEvent;
0092         G4double squareDose = edepMap[itr.first].second / (itr.second*itr.second);
0093         G4double variance    = ((squareDose/fNumOfEvent) - (meanDose*meanDose))/fNumOfEvent;
0094         G4double relativeE(1);
0095         if(meanDose > 0) relativeE   = sqrt(variance)/meanDose;
0096
0097         out << setw(8)  << itr.first << "| "
0098             << setw(19) << fixed      << itr.second/g;
0099         out << setw(19) << scientific << meanDose/(joule/kg);
0100         out << setw(19) << fixed      << relativeE << G4endl;
0101     }
0102     out << "=====================================================================" << G4endl << G4endl;
0103 }