Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/runAndEvent/RE02/src/RE02RunAction.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/runAndEvent/RE02/src/RE02RunAction.cc (Version 11.3.0) and /examples/extended/runAndEvent/RE02/src/RE02RunAction.cc (Version 9.6)


  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 /// \file runAndEvent/RE02/src/RE02RunAction.c     26 /// \file runAndEvent/RE02/src/RE02RunAction.cc
 27 /// \brief Implementation of the RE02RunAction     27 /// \brief Implementation of the RE02RunAction class
 28 //                                                 28 //
 29 //                                                 29 //
 30 //                                             <<  30 // $Id$
                                                   >>  31 // 
 31 #include "RE02RunAction.hh"                        32 #include "RE02RunAction.hh"
 32                                                << 
 33 #include "RE02Run.hh"                              33 #include "RE02Run.hh"
 34                                                    34 
 35 //-- In order to obtain detector information.      35 //-- In order to obtain detector information.
 36 #include "RE02DetectorConstruction.hh"         << 
 37                                                << 
 38 #include "G4RunManager.hh"                         36 #include "G4RunManager.hh"
 39 #include "G4SystemOfUnits.hh"                  <<  37 #include "RE02DetectorConstruction.hh"
 40 #include "G4THitsMap.hh"                           38 #include "G4THitsMap.hh"
                                                   >>  39 
 41 #include "G4UnitsTable.hh"                         40 #include "G4UnitsTable.hh"
                                                   >>  41 #include "G4SystemOfUnits.hh"    
 42                                                    42 
 43 #include <fstream>                                 43 #include <fstream>
 44 //============================================     44 //=======================================================================
 45 // RE02RunAction                                   45 // RE02RunAction
 46 //                                             <<  46 //  
 47 //                                                 47 //
 48 //                                                 48 //
 49 //============================================     49 //=======================================================================
 50                                                    50 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52 // Constructor                                     52 // Constructor
 53 RE02RunAction::RE02RunAction() : G4UserRunActi <<  53 RE02RunAction::RE02RunAction()
                                                   >>  54   : fNx(0), fNy(0), fNz(0)
 54 {                                                  55 {
 55   // - Prepare data member for RE02Run.            56   // - Prepare data member for RE02Run.
 56   //   vector represents a list of MultiFuncti     57   //   vector represents a list of MultiFunctionalDetector names.
 57   fSDName.push_back(G4String("PhantomSD"));        58   fSDName.push_back(G4String("PhantomSD"));
 58 }                                                  59 }
 59                                                    60 
 60 //....oooOO0OOooo........oooOO0OOooo........oo     61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 61 // Destructor.                                     62 // Destructor.
 62 RE02RunAction::~RE02RunAction()                    63 RE02RunAction::~RE02RunAction()
 63 {                                                  64 {
 64   fSDName.clear();                                 65   fSDName.clear();
 65 }                                                  66 }
 66                                                    67 
 67 //....oooOO0OOooo........oooOO0OOooo........oo     68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 68 //==                                           <<  69 //== 
 69 G4Run* RE02RunAction::GenerateRun()                70 G4Run* RE02RunAction::GenerateRun()
 70 {                                                  71 {
 71   // Generate new RUN object, which is special     72   // Generate new RUN object, which is specially
 72   // dedicated for MultiFunctionalDetector sch     73   // dedicated for MultiFunctionalDetector scheme.
 73   //  Detail description can be found in RE02R     74   //  Detail description can be found in RE02Run.hh/cc.
 74   return new RE02Run(fSDName);                     75   return new RE02Run(fSDName);
 75 }                                                  76 }
 76                                                    77 
 77 //....oooOO0OOooo........oooOO0OOooo........oo     78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 78 //==                                               79 //==
 79 void RE02RunAction::BeginOfRunAction(const G4R     80 void RE02RunAction::BeginOfRunAction(const G4Run* aRun)
 80 {                                                  81 {
 81   G4cout << "### Run " << aRun->GetRunID() <<      82   G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
 82 }                                                  83 }
 83                                                    84 
 84 //....oooOO0OOooo........oooOO0OOooo........oo     85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 85 //==                                           <<  86 //== 
 86 void RE02RunAction::EndOfRunAction(const G4Run     87 void RE02RunAction::EndOfRunAction(const G4Run* aRun)
 87 {                                                  88 {
 88   if (!IsMaster()) return;                     << 
 89                                                << 
 90   //- RE02Run object.                              89   //- RE02Run object.
 91   RE02Run* re02Run = (RE02Run*)aRun;               90   RE02Run* re02Run = (RE02Run*)aRun;
 92   //--- Dump all socred quantities involved in     91   //--- Dump all socred quantities involved in RE02Run.
 93   re02Run->DumpAllScorer();                        92   re02Run->DumpAllScorer();
 94   //---                                            93   //---
 95                                                    94 
 96   //                                               95   //
 97   //- water phantom (Detector) Information.        96   //- water phantom (Detector) Information.
 98   //-- Number of segments in the water phantom     97   //-- Number of segments in the water phantom.
 99   const RE02DetectorConstruction* detector =       98   const RE02DetectorConstruction* detector =
100     (const RE02DetectorConstruction*)(G4RunMan <<  99       (const RE02DetectorConstruction*)
101   detector->GetNumberOfSegmentsInPhantom(fNx,  << 100       (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
                                                   >> 101   detector->GetNumberOfSegmentsInPhantom(fNx,fNy,fNz); //Fill fNx,y,z.
102                                                   102 
103   //------------------------------------------    103   //---------------------------------------------
104   // Dump accumulated quantities for this RUN.    104   // Dump accumulated quantities for this RUN.
105   //  (Display only central region of x-y plan    105   //  (Display only central region of x-y plane)
106   //------------------------------------------    106   //---------------------------------------------
107   G4THitsMap<G4double>* totEdep = re02Run->Get << 107   G4THitsMap<G4double>* totEdep  = re02Run->GetHitsMap("PhantomSD/totalEDep");
108   G4THitsMap<G4double>* proEdep = re02Run->Get << 108   G4THitsMap<G4double>* proEdep  = re02Run->GetHitsMap("PhantomSD/protonEDep");
109   G4THitsMap<G4double>* proNstep = re02Run->Ge    109   G4THitsMap<G4double>* proNstep = re02Run->GetHitsMap("PhantomSD/protonNStep");
110   G4THitsMap<G4double>* passCFx = re02Run->Get << 110   G4THitsMap<G4double>* passCFx  =
111   G4THitsMap<G4double>* cFx = re02Run->GetHits << 111     re02Run->GetHitsMap("PhantomSD/chargedPassCellFlux");
112   G4THitsMap<G4double>* surfFx = re02Run->GetH << 112   G4THitsMap<G4double>*     cFx  =
113   G4THitsMap<G4double>* gCurr00 = re02Run->Get << 113     re02Run->GetHitsMap("PhantomSD/chargedCellFlux");
114   G4THitsMap<G4double>* gCurr01 = re02Run->Get << 114   G4THitsMap<G4double>*  surfFx  =
115   G4THitsMap<G4double>* gCurr02 = re02Run->Get << 115     re02Run->GetHitsMap("PhantomSD/chargedSurfFlux");
116   G4THitsMap<G4double>* gCurr03 = re02Run->Get << 116   G4THitsMap<G4double>* gCurr00  =
117                                                << 117     re02Run->GetHitsMap("PhantomSD/gammaSurfCurr000");
118   G4cout << "================================= << 118   G4THitsMap<G4double>* gCurr01  =
119   G4cout << " Number of event processed : " << << 119     re02Run->GetHitsMap("PhantomSD/gammaSurfCurr001");
120   G4cout << "================================= << 120   G4THitsMap<G4double>* gCurr02  =
121   G4cout << std::setw(8) << "#Z Cell#";        << 121     re02Run->GetHitsMap("PhantomSD/gammaSurfCurr002");
                                                   >> 122   G4THitsMap<G4double>* gCurr03  =
                                                   >> 123     re02Run->GetHitsMap("PhantomSD/gammaSurfCurr003");
                                                   >> 124 
                                                   >> 125   G4cout << "============================================================="
                                                   >> 126          << G4endl;
                                                   >> 127   G4cout << " Number of event processed : "<< aRun->GetNumberOfEvent() << G4endl;
                                                   >> 128   G4cout << "============================================================="
                                                   >> 129          << G4endl;
                                                   >> 130   G4cout << std::setw( 8) << "#Z Cell#";
122   G4cout << std::setw(16) << totEdep->GetName(    131   G4cout << std::setw(16) << totEdep->GetName();
123   G4cout << std::setw(16) << proEdep->GetName(    132   G4cout << std::setw(16) << proEdep->GetName();
124   G4cout << std::setw(12) << proNstep->GetName    133   G4cout << std::setw(12) << proNstep->GetName();
125   G4cout << std::setw(21) << passCFx->GetName(    134   G4cout << std::setw(21) << passCFx->GetName();
126   G4cout << std::setw(20) << cFx->GetName();      135   G4cout << std::setw(20) << cFx->GetName();
127   G4cout << std::setw(20) << surfFx->GetName()    136   G4cout << std::setw(20) << surfFx->GetName();
128   G4cout << std::setw(20) << gCurr00->GetName(    137   G4cout << std::setw(20) << gCurr00->GetName();
129   G4cout << std::setw(20) << gCurr01->GetName(    138   G4cout << std::setw(20) << gCurr01->GetName();
130   G4cout << std::setw(20) << gCurr02->GetName(    139   G4cout << std::setw(20) << gCurr02->GetName();
131   G4cout << std::setw(20) << gCurr03->GetName( << 140   G4cout << std::setw(20) << gCurr03->GetName()
132   G4int ix = fNx / 2;                          << 141          << G4endl;
133   G4int iy = fNy / 2;                          << 142   G4int ix = fNx/2;  
                                                   >> 143   G4int iy = fNy/2;
134   G4int iz;                                       144   G4int iz;
135   // G4double totE, proE, proN,pasCF,CF,surfF, << 145   //G4double totE, proE, proN,pasCF,CF,surfF,gCr0,gCr1,gCr2,gCr3;
136   for (iz = 0; iz < fNz; iz++) {               << 146   for ( iz = 0; iz < fNz; iz++){ 
137     G4double* totED = (*totEdep)[CopyNo(ix, iy << 147       G4double* totED = (*totEdep)[CopyNo(ix,iy,iz)];
138     G4double* proED = (*proEdep)[CopyNo(ix, iy << 148       G4double* proED = (*proEdep)[CopyNo(ix,iy,iz)];
139     G4double* proNS = (*proNstep)[CopyNo(ix, i << 149       G4double* proNS = (*proNstep)[CopyNo(ix,iy,iz)];
140     G4double* pasCF = (*passCFx)[CopyNo(ix, iy << 150       G4double* pasCF = (*passCFx)[CopyNo(ix,iy,iz)];
141     G4double* cF = (*cFx)[CopyNo(ix, iy, iz)]; << 151       G4double* cF    = (*cFx)[CopyNo(ix,iy,iz)];
142     G4double* sfx = (*surfFx)[CopyNo(ix, iy, i << 152       G4double* sfx   = (*surfFx)[CopyNo(ix,iy,iz)];
143     G4double* gcur0 = (*gCurr00)[CopyNo(ix, iy << 153       G4double* gcur0 = (*gCurr00)[CopyNo(ix,iy,iz)];
144     G4double* gcur1 = (*gCurr01)[CopyNo(ix, iy << 154       G4double* gcur1 = (*gCurr01)[CopyNo(ix,iy,iz)];
145     G4double* gcur2 = (*gCurr02)[CopyNo(ix, iy << 155       G4double* gcur2 = (*gCurr02)[CopyNo(ix,iy,iz)];
146     G4double* gcur3 = (*gCurr03)[CopyNo(ix, iy << 156       G4double* gcur3 = (*gCurr03)[CopyNo(ix,iy,iz)];
147     if (!totED) totED = new G4double(0.0);     << 157       if ( !totED ) totED = new G4double(0.0);
148     if (!proED) proED = new G4double(0.0);     << 158       if ( !proED ) proED = new G4double(0.0);
149     if (!proNS) proNS = new G4double(0.0);     << 159       if ( !proNS ) proNS = new G4double(0.0);
150     if (!pasCF) pasCF = new G4double(0.0);     << 160       if ( !pasCF ) pasCF = new G4double(0.0);
151     if (!cF) cF = new G4double(0.0);           << 161       if ( !cF    )    cF = new G4double(0.0);
152     if (!sfx) sfx = new G4double(0.0);         << 162       if ( !sfx   )   sfx = new G4double(0.0);
153     if (!gcur0) gcur0 = new G4double(0.0);     << 163       if ( !gcur0 ) gcur0 = new G4double(0.0);
154     if (!gcur1) gcur1 = new G4double(0.0);     << 164       if ( !gcur1 ) gcur1 = new G4double(0.0);
155     if (!gcur2) gcur2 = new G4double(0.0);     << 165       if ( !gcur2 ) gcur2 = new G4double(0.0);
156     if (!gcur3) gcur3 = new G4double(0.0);     << 166       if ( !gcur3 ) gcur3 = new G4double(0.0);
157     G4cout << std::setw(6) << iz << "  " << st << 167       G4cout << std::setw( 6) << iz << "  "
158            << std::setw(12) << G4BestUnit(*pro << 168              << std::setw(12) << G4BestUnit(*totED,"Energy")
159            << std::setw(13) << (*pasCF) * cm * << 169              << std::setw(12) << G4BestUnit(*proED,"Energy")
160            << " /cm2" << std::setw(15) << (*sf << 170              << std::setw(12) << (*proNS)        << "   "
161            << (*gcur0) * cm * cm << " /cm2" << << 171              << std::setw(13) << (*pasCF)*cm*cm  <<" /cm2"
162            << std::setw(15) << (*gcur2) * cm * << 172              << std::setw(15) << (*cF)*cm*cm     <<" /cm2"
163            << " /cm2" << G4endl;               << 173              << std::setw(15) << (*sfx)*cm*cm    <<" /cm2"
                                                   >> 174              << std::setw(15) << (*gcur0)*cm*cm  <<" /cm2"
                                                   >> 175              << std::setw(15) << (*gcur1)*cm*cm  <<" /cm2"
                                                   >> 176              << std::setw(15) << (*gcur2)*cm*cm  <<" /cm2"
                                                   >> 177              << std::setw(15) << (*gcur3)*cm*cm  <<" /cm2"
                                                   >> 178              << G4endl;
164   }                                               179   }
165   G4cout << "================================= << 180   G4cout << "============================================="<<G4endl;
                                                   >> 181 
                                                   >> 182   std::ofstream  file("totED.txt");
                                                   >> 183   for ( iz = 0; iz < fNz; iz++){   
                                                   >> 184     for ( iy = 0; iy < fNy; iy++){ 
                                                   >> 185       for ( ix = 0; ix < fNx; ix++){ 
                                                   >> 186         G4double* totED = (*totEdep)[CopyNo(ix,iy,iz)];
                                                   >> 187         if ( !totED ) totED = new G4double(0.0);
                                                   >> 188         file << ix << " "<<iy<<" "<<iz<<" "<< *totED/MeV << G4endl;
166                                                   189 
167   std::ofstream file("totED.txt");             << 
168   for (iz = 0; iz < fNz; iz++) {               << 
169     for (iy = 0; iy < fNy; iy++) {             << 
170       for (ix = 0; ix < fNx; ix++) {           << 
171         G4double* totED = (*totEdep)[CopyNo(ix << 
172         if (!totED) totED = new G4double(0.0); << 
173         file << ix << " " << iy << " " << iz < << 
174       }                                           190       }
175     }                                             191     }
176   }                                               192   }
177   file.close();                                   193   file.close();
178 }                                                 194 }
179 //                                                195 //
180 // --                                             196 // --
181                                                   197