Geant4 Cross Reference |
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 electromagnetic/TestEm4/src/RunActio 26 /// \file electromagnetic/TestEm4/src/RunAction.cc 27 /// \brief Implementation of the RunAction cla 27 /// \brief Implementation of the RunAction class 28 // 28 // 29 // << 29 // $Id: RunAction.cc 103468 2017-04-11 07:28:24Z gcosmo $ >> 30 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 33 33 #include "RunAction.hh" 34 #include "RunAction.hh" 34 35 35 #include "PrimaryGeneratorAction.hh" << 36 << 37 #include "G4AnalysisManager.hh" << 38 #include "G4ParticleGun.hh" << 39 #include "G4Run.hh" 36 #include "G4Run.hh" 40 #include "G4RunManager.hh" 37 #include "G4RunManager.hh" >> 38 41 #include "G4SystemOfUnits.hh" 39 #include "G4SystemOfUnits.hh" 42 #include "G4UnitsTable.hh" << 43 #include "Randomize.hh" 40 #include "Randomize.hh" 44 41 45 //....oooOO0OOooo........oooOO0OOooo........oo 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 46 43 47 RunAction::RunAction() 44 RunAction::RunAction() >> 45 : G4UserRunAction() 48 { 46 { 49 G4AnalysisManager* analysisManager = G4Analy 47 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); 50 analysisManager->SetDefaultFileType("root"); << 48 analysisManager->SetVerboseLevel(1); 51 analysisManager->SetVerboseLevel(1); << 52 analysisManager->SetFirstHistoId(1); 49 analysisManager->SetFirstHistoId(1); 53 << 50 54 // Creating histograms 51 // Creating histograms 55 // 52 // 56 analysisManager->CreateH1("1", "energy (MeV) << 53 analysisManager->CreateH1("1","energy (MeV) deposited in C6F6",100,0.,10.); >> 54 } >> 55 >> 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 57 >> 58 RunAction::~RunAction() >> 59 { >> 60 delete G4AnalysisManager::Instance(); 57 } 61 } 58 62 59 //....oooOO0OOooo........oooOO0OOooo........oo 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 60 64 61 void RunAction::BeginOfRunAction(const G4Run*) 65 void RunAction::BeginOfRunAction(const G4Run*) 62 { 66 { 63 // show Rndm status 67 // show Rndm status 64 if (isMaster) G4Random::showEngineStatus(); 68 if (isMaster) G4Random::showEngineStatus(); 65 69 66 // Get analysis manager << 70 >> 71 // Get analysis manager 67 G4AnalysisManager* analysisManager = G4Analy 72 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); 68 73 69 // Open an output file 74 // Open an output file 70 // 75 // 71 G4String fileName = "testem4"; 76 G4String fileName = "testem4"; 72 analysisManager->OpenFile(fileName); << 77 analysisManager->OpenFile(fileName); 73 } 78 } 74 79 75 //....oooOO0OOooo........oooOO0OOooo........oo 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 76 81 77 void RunAction::EndOfRunAction(const G4Run* ru << 82 void RunAction::EndOfRunAction(const G4Run*) 78 { << 83 { 79 G4int nofEvents = run->GetNumberOfEvent(); << 80 if (nofEvents == 0) return; << 81 << 82 G4AnalysisManager* analysisManager = G4Analy << 83 << 84 if (isMaster) { << 85 // Run conditions << 86 // note: There is no primary generator ac << 87 // run manager for multi-threaded m << 88 PrimaryGeneratorAction* generatorAction = << 89 (PrimaryGeneratorAction*)(G4RunManager:: << 90 G4String runCondition = " "; << 91 if (generatorAction) { << 92 G4ParticleGun* particleGun = generatorAc << 93 runCondition = particleGun->GetParticleD << 94 runCondition += " of "; << 95 G4double particleEnergy = particleGun->G << 96 runCondition += G4BestUnit(particleEnerg << 97 runCondition += " in a volume of C6F6 "; << 98 } << 99 << 100 G4cout << G4endl << "--------------------E << 101 G4cout << G4endl << " The run consists of << 102 if (analysisManager->GetH1(1)) { << 103 G4cout << " Energy deposit per event: m << 104 << G4BestUnit(analysisManager->Ge << 105 << " rms = " << G4BestUnit(ana << 106 << "----------------------------- << 107 } << 108 } << 109 << 110 // show Rndm status 84 // show Rndm status 111 if (isMaster) G4Random::showEngineStatus(); << 85 if (isMaster) G4Random::showEngineStatus(); 112 86 113 // save histograms << 87 //save histograms >> 88 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); 114 analysisManager->Write(); 89 analysisManager->Write(); 115 analysisManager->CloseFile(); 90 analysisManager->CloseFile(); 116 } 91 } 117 92 118 //....oooOO0OOooo........oooOO0OOooo........oo 93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 119 94