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