Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/dnadamage2/src/RunAction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/medical/dna/dnadamage2/src/RunAction.cc (Version 11.3.0) and /examples/extended/medical/dna/dnadamage2/src/RunAction.cc (Version 11.2.2)


  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 // This example is provided by the Geant4-DNA      26 // This example is provided by the Geant4-DNA collaboration
 27 // dnadamage3 example is derived from the chem     27 // dnadamage3 example is derived from the chem6 example
 28 // chem6 example authors: W. G. Shin and S. In     28 // chem6 example authors: W. G. Shin and S. Incerti (CENBG, France)
 29 //                                                 29 //
 30 // Any report or published results obtained us     30 // Any report or published results obtained using the Geant4-DNA software
 31 // shall cite the following Geant4-DNA collabo     31 // shall cite the following Geant4-DNA collaboration publication:
 32 // J. Appl. Phys. 125 (2019) 104301                32 // J. Appl. Phys. 125 (2019) 104301
 33 // Med. Phys. 45 (2018) e722-e739                  33 // Med. Phys. 45 (2018) e722-e739
 34 // J. Comput. Phys. 274 (2014) 841-882             34 // J. Comput. Phys. 274 (2014) 841-882
 35 // Med. Phys. 37 (2010) 4692-4708                  35 // Med. Phys. 37 (2010) 4692-4708
 36 // Int. J. Model. Simul. Sci. Comput. 1 (2010)     36 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157-178
 37 // The Geant4-DNA web site is available at htt     37 // The Geant4-DNA web site is available at http://geant4-dna.org
 38 //                                                 38 //
 39 // Authors: J. Naoki D. Kondo (UCSF, US)           39 // Authors: J. Naoki D. Kondo (UCSF, US)
 40 //          J. Ramos-Mendez and B. Faddegon (U     40 //          J. Ramos-Mendez and B. Faddegon (UCSF, US)
 41 //                                                 41 //
 42 //                                                 42 //
 43 /// \file RunAction.cc                             43 /// \file RunAction.cc
 44 /// \brief Implementation of the RunAction cla     44 /// \brief Implementation of the RunAction class
 45                                                    45 
 46 #include "RunAction.hh"                            46 #include "RunAction.hh"
 47                                                << 
 48 #include "Run.hh"                                  47 #include "Run.hh"
 49                                                    48 
 50 #include "G4Run.hh"                                49 #include "G4Run.hh"
 51 #include "G4RunManager.hh"                         50 #include "G4RunManager.hh"
 52 #include "G4SystemOfUnits.hh"                  << 
 53 #include "G4UnitsTable.hh"                         51 #include "G4UnitsTable.hh"
                                                   >>  52 #include "G4SystemOfUnits.hh"
 54                                                    53 
 55 //....oooOO0OOooo........oooOO0OOooo........oo     54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
 56                                                    55 
 57 RunAction::RunAction() : G4UserRunAction() {}  <<  56 RunAction::RunAction()
                                                   >>  57  : G4UserRunAction()
                                                   >>  58 {
                                                   >>  59 }
 58                                                    60 
 59 //....oooOO0OOooo........oooOO0OOooo........oo     61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
 60                                                    62 
 61 G4Run* RunAction::GenerateRun()                    63 G4Run* RunAction::GenerateRun()
 62 {                                                  64 {
 63   Run* run = new Run();                            65   Run* run = new Run();
 64   return run;                                      66   return run;
 65 }                                                  67 }
 66                                                    68 
 67 //....oooOO0OOooo........oooOO0OOooo........oo     69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
 68                                                    70 
 69 void RunAction::BeginOfRunAction(const G4Run*)     71 void RunAction::BeginOfRunAction(const G4Run*)
 70 {                                                  72 {
 71   // informs the runManager to save random num     73   // informs the runManager to save random number seed
 72   G4RunManager::GetRunManager()->SetRandomNumb     74   G4RunManager::GetRunManager()->SetRandomNumberStore(false);
 73 }                                                  75 }
 74                                                    76 
 75 //....oooOO0OOooo........oooOO0OOooo........oo     77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
 76                                                    78 
 77 void RunAction::EndOfRunAction(const G4Run* ru     79 void RunAction::EndOfRunAction(const G4Run* run)
 78 {                                                  80 {
 79   G4int nofEvents = run->GetNumberOfEvent();       81   G4int nofEvents = run->GetNumberOfEvent();
 80   if (nofEvents == 0) return;                      82   if (nofEvents == 0) return;
 81                                                    83 
 82   // results                                       84   // results
 83   //                                               85   //
 84   const Run* dnadamage3Run = static_cast<const     86   const Run* dnadamage3Run = static_cast<const Run*>(run);
 85   G4double sumDose = dnadamage3Run->GetSumDose <<  87   G4double sumDose   = dnadamage3Run->GetSumDose();
 86                                                    88 
 87   // print                                         89   // print
 88   //                                               90   //
 89   if (IsMaster()) {                            <<  91   if (IsMaster())
 90     G4cout << G4endl << "--------------------E <<  92   {
 91            << "  The run has " << nofEvents << <<  93 
                                                   >>  94     G4cout
                                                   >>  95       << G4endl
                                                   >>  96       << "--------------------End of Global Run-----------------------"
                                                   >>  97       << G4endl
                                                   >>  98       << "  The run has " << nofEvents << " events "
                                                   >>  99       << G4endl;
 92                                                   100 
 93     ScoreSpecies* masterScorer = dynamic_cast< << 101     ScoreSpecies* masterScorer=
                                                   >> 102       dynamic_cast<ScoreSpecies*>(dnadamage3Run->GetPrimitiveScorer());
 94                                                   103 
 95     ScoreStrandBreaks* masterSBScorer =        << 104     ScoreStrandBreaks* masterSBScorer=
 96       dynamic_cast<ScoreStrandBreaks*>(dnadama    105       dynamic_cast<ScoreStrandBreaks*>(dnadamage3Run->GetSBScorer());
 97                                                   106 
 98     G4cout << "Number of events recorded by th    107     G4cout << "Number of events recorded by the species scorer = "
 99            << masterScorer->GetNumberOfRecorde << 108       << masterScorer->GetNumberOfRecordedEvents()
                                                   >> 109       << G4endl;
100                                                   110 
101     // LET                                        111     // LET
102     Run* aRun = (Run*)run;                        112     Run* aRun = (Run*)run;
103     G4THitsMap<G4double>* totLET = aRun->GetLE    113     G4THitsMap<G4double>* totLET = aRun->GetLET();
104     G4int nOfEvent = totLET->entries();           114     G4int nOfEvent = totLET->entries();
105     G4double LET_mean = 0;                        115     G4double LET_mean = 0;
106     G4double LET_square = 0;                      116     G4double LET_square = 0;
107     for (G4int i = 0; i < nOfEvent; i++) {     << 117     for(G4int i=0;i<nOfEvent;i++){
108       G4double* LET = (*totLET)[i];               118       G4double* LET = (*totLET)[i];
109       if (!LET) continue;                      << 119       if(!LET) continue;
110       LET_mean += *LET;                           120       LET_mean += *LET;
111       LET_square += (*LET) * (*LET);           << 121       LET_square += (*LET)*(*LET);
112     }                                             122     }
113     LET_mean /= nOfEvent;                         123     LET_mean /= nOfEvent;
114     LET_square = std::sqrt(LET_square / nOfEve << 124     LET_square = std::sqrt(LET_square/nOfEvent - std::pow(LET_mean,2));
115                                                   125 
116     masterScorer->OutputAndClear();               126     masterScorer->OutputAndClear();
117     masterSBScorer->OutputAndClear(LET_mean, L << 127     masterSBScorer->OutputAndClear(LET_mean,LET_square);
                                                   >> 128 
118   }                                               129   }
119   else {                                       << 130   else
120     G4cout << G4endl << "--------------------E << 131   {
121            << "  The run has " << nofEvents << << 132     G4cout
                                                   >> 133       << G4endl
                                                   >> 134       << "--------------------End of Local Run------------------------"
                                                   >> 135       << G4endl
                                                   >> 136       << "  The run has " << nofEvents << " events"
                                                   >> 137       << G4endl;
122   }                                               138   }
123                                                   139 
124   G4cout << " Total energy deposited in the wo << 140   G4cout
125          << " -------------------------------- << 141     << " Total energy deposited in the world volume : "
                                                   >> 142     << sumDose/eV << " eV"
                                                   >> 143     << G4endl
                                                   >> 144     << " ------------------------------------------------------------"
                                                   >> 145     << G4endl
                                                   >> 146     << G4endl;
126 }                                                 147 }
127                                                   148 
128 //....oooOO0OOooo........oooOO0OOooo........oo    149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
129                                                   150