Geant4 Cross Reference |
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 // 27 // ------------------------------------------------------------ 28 // GEANT 4 class implementation file 29 // CERN Geneva Switzerland 30 // 31 // 32 // ------------ GammaRayTelRunAction ------ 33 // by R.Giannitrapani, F.Longo & G.Santin (13 nov 2000) 34 // 18.11.2001 G.Santin 35 // - Modified the analysis management according to the new design 36 // 37 // ************************************************************ 38 39 #include "GammaRayTelRunAction.hh" 40 #include "GammaRayTelAnalysis.hh" 41 42 #include <cstdlib> 43 44 #include "G4ios.hh" 45 #include "G4Run.hh" 46 #include "G4Threading.hh" 47 #include "G4UImanager.hh" 48 #include "G4VVisManager.hh" 49 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 51 52 GammaRayTelRunAction::GammaRayTelRunAction() { 53 } 54 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 56 57 GammaRayTelRunAction::~GammaRayTelRunAction() { 58 } 59 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 61 62 void GammaRayTelRunAction::BeginOfRunAction(const G4Run *run) { 63 fRunID = run->GetRunID(); 64 65 // Master mode or sequential 66 if (IsMaster()) { 67 G4cout << "### Run " << run->GetRunID() << " starts (master)." << G4endl; 68 } else { 69 G4cout << "### Run " << run->GetRunID() << " starts (worker)." << G4endl; 70 } 71 72 // Prepare the visualization 73 if (G4VVisManager::GetConcreteInstance() != nullptr) { 74 auto *userInterface = G4UImanager::GetUIpointer(); 75 userInterface->ApplyCommand("/vis/scene/notifyHandlers"); 76 } 77 78 // If analysis is used reset the histograms 79 auto *analysis = GammaRayTelAnalysis::getInstance(); 80 // analysis->BeginOfRun(run->GetRunID()); 81 analysis->BeginOfRun(); 82 } 83 84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 85 86 void GammaRayTelRunAction::EndOfRunAction(const G4Run *run) { 87 G4cout << "End of Run " << run->GetRunID() << G4endl; 88 89 // Close the file with the hits information 90 #ifdef G4STORE_DATA 91 if (outFile != nullptr) { 92 G4cout << "File " << fileName << G4endl; 93 outFile->close(); 94 delete outFile; 95 outFile = nullptr; 96 } 97 #endif 98 99 // If analysis is used, print out the histograms 100 auto *analysis = GammaRayTelAnalysis::getInstance(); 101 analysis->EndOfRun(); 102 } 103 104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 105 106 auto GammaRayTelRunAction::GetOutputFile() -> std::ofstream* { 107 if (outFile == nullptr) { 108 OpenFile(); 109 } 110 return outFile; 111 } 112 113 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 114 115 void GammaRayTelRunAction::OpenFile() { 116 // Open the file for the tracks of this run 117 #ifdef G4STORE_DATA 118 // check that we are in a worker: returns -1 in a master and -2 in sequential 119 // one file per thread is produced 120 // Tracks_runR.N.dat, where R = run number, N = thread ID 121 std::stringstream name; 122 123 if (G4Threading::G4GetThreadId() >= 0) { 124 name << "Tracks_run" << fRunID << "." << G4Threading::G4GetThreadId() << ".dat"; 125 } else { 126 name << "Tracks_run" << fRunID << "." << G4Threading::G4GetThreadId() << ".dat"; 127 } 128 129 if (outFile == nullptr) { 130 outFile = new std::ofstream; 131 outFile->open(name.str()); 132 fileName = G4String(name.str()); 133 } 134 G4cout << "Open file: " << fileName << G4endl; 135 #endif 136 } 137