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.2.p1)


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