Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm10/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 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 /// \file electromagnetic/TestEm10/src/RunAction.cc
 27 /// \brief Implementation of the RunAction class
 28 //
 29 //
 30 //
 31 //
 32 
 33 #include "RunAction.hh"
 34 
 35 #include "RunMessenger.hh"
 36 
 37 #include "G4AnalysisManager.hh"
 38 #include "G4PhysicalConstants.hh"
 39 #include "G4Run.hh"
 40 #include "G4SystemOfUnits.hh"
 41 #include "G4ios.hh"
 42 #include "Randomize.hh"
 43 
 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 45 
 46 RunAction::RunAction() : G4UserRunAction(), fRunMessenger(0), fRndmFreq(1)
 47 {
 48   fRunMessenger = new RunMessenger(this);
 49 
 50   BookHisto();
 51 }
 52 
 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 54 
 55 RunAction::~RunAction()
 56 {
 57   delete fRunMessenger;
 58 }
 59 
 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 61 
 62 void RunAction::BookHisto()
 63 {
 64   // Create or get analysis manager
 65   // The choice of analysis technology is done via selection of a namespace
 66   // in HistoManager.hh
 67   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 68   analysisManager->SetDefaultFileType("root");
 69   analysisManager->SetFileName("testem10");
 70   analysisManager->SetVerboseLevel(1);
 71   analysisManager->SetFirstHistoId(1);  // start histogram numbering from 1
 72   analysisManager->SetActivation(true);  // enable inactivation of histograms
 73 
 74   // Define histograms start values
 75   const G4int kMaxHisto = 5;
 76   const G4String id[] = {"1", "2", "3", "4", "5"};
 77   const G4String title[] = {"Edep",  // 1
 78                             "XTR Gamma spectrum",  // 2
 79                             "Secondary Gamma spectrum",  // 3
 80                             "Secondary e- spectrum",  // 4
 81                             "Edep.old"};  // 5
 82   // Default values (to be reset via /analysis/h1/set command)
 83   G4int nbins = 100;
 84   G4double vmin = 0.;
 85   G4double vmax = 100.;
 86 
 87   // Create all histograms as inactivated
 88   // as we have not yet set nbins, vmin, vmax
 89   for (G4int k = 0; k < kMaxHisto; k++) {
 90     G4int ih = analysisManager->CreateH1(id[k], title[k], nbins, vmin, vmax);
 91     analysisManager->SetH1Activation(ih, false);
 92   }
 93 }
 94 
 95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 96 
 97 void RunAction::BeginOfRunAction(const G4Run* aRun)
 98 {
 99   G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
100 
101   // save Rndm status
102   if (fRndmFreq > 0) {
103     CLHEP::HepRandom::showEngineStatus();
104     CLHEP::HepRandom::saveEngineStatus("beginOfRun.rndm");
105   }
106 
107   //  histograms
108   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
109   if (analysisManager->IsActive()) {
110     analysisManager->OpenFile();
111   }
112 }
113 
114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
115 
116 void RunAction::EndOfRunAction(const G4Run* run)
117 {
118   // print run statisctics
119   G4int nofEvents = run->GetNumberOfEvent();
120   G4cout << " ================== run summary =====================" << G4endl;
121   G4cout << " End of Run TotNbofEvents = " << nofEvents << G4endl;
122 
123   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
124   if (analysisManager->GetH1(1)) {
125     G4cout << " Mean energy deposit in absorber = " << analysisManager->GetH1(1)->mean() / MeV
126            << " +-" << analysisManager->GetH1(1)->rms() / MeV << "  MeV " << G4endl;
127   }
128   if (analysisManager->GetH1(2)) {
129     G4cout << " Total number of XTR gammas = " << analysisManager->GetH1(2)->entries() << G4endl;
130   }
131   if (analysisManager->GetH1(3)) {
132     G4cout << " Total number of all gammas = " << analysisManager->GetH1(3)->entries() << G4endl;
133   }
134 
135   // save Rndm status
136   if (fRndmFreq == 1) {
137     CLHEP::HepRandom::showEngineStatus();
138     CLHEP::HepRandom::saveEngineStatus("endOfRun.rndm");
139   }
140 
141   // save histograms
142   if (analysisManager->IsActive()) {
143     analysisManager->Write();
144     analysisManager->CloseFile();
145   }
146 }
147 
148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
149