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 /// \file field/field06/field06.cc 27 /// \file field/field06/field06.cc 28 /// \brief Main program of the field/field06 e 28 /// \brief Main program of the field/field06 example 29 // 29 // >> 30 // >> 31 // >> 32 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 35 33 #ifndef WIN32 36 #ifndef WIN32 34 # include <unistd.h> << 37 #include <unistd.h> 35 #endif 38 #endif 36 39 37 #include "F06ActionInitialization.hh" << 40 #include "G4Types.hh" 38 #include "F06DetectorConstruction.hh" << 39 #include "F06PhysicsList.hh" << 40 41 41 #include "G4RunManagerFactory.hh" << 42 #include "F06PhysicsList.hh" 42 #include "G4Transportation.hh" 43 #include "G4Transportation.hh" 43 #include "G4Types.hh" << 44 44 #include "G4UIExecutive.hh" << 45 #include "F06DetectorConstruction.hh" 45 #include "G4UIcommand.hh" << 46 #include "F06ActionInitialization.hh" >> 47 >> 48 #ifdef G4MULTITHREADED >> 49 #include "G4MTRunManager.hh" >> 50 #else >> 51 #include "G4RunManager.hh" >> 52 #endif >> 53 46 #include "G4UImanager.hh" 54 #include "G4UImanager.hh" 47 #include "G4VisExecutive.hh" << 55 #include "G4UIcommand.hh" >> 56 48 #include "Randomize.hh" 57 #include "Randomize.hh" 49 58 >> 59 #include "G4VisExecutive.hh" >> 60 #include "G4UIExecutive.hh" >> 61 50 //....oooOO0OOooo........oooOO0OOooo........oo 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 51 namespace << 63 namespace { 52 { << 64 void PrintUsage() { 53 void PrintUsage() << 65 G4cerr << " Usage: " << G4endl; 54 { << 66 G4cerr << " field06 [-m macro ] [-u UIsession] [-t nThreads] [-r randomSeed] " 55 G4cerr << " Usage: " << G4endl; << 67 << G4endl; 56 G4cerr << " field06 [-m macro ] [-u UIsessio << 68 G4cerr << " note: -t option is available only for multi-threaded mode." 57 G4cerr << " note: -t option is available o << 69 << G4endl; >> 70 } 58 } 71 } 59 } // namespace << 60 72 61 //....oooOO0OOooo........oooOO0OOooo........oo 73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 62 74 63 int main(int argc, char** argv) << 75 int main(int argc,char** argv) 64 { 76 { 65 // Evaluate arguments 77 // Evaluate arguments 66 // 78 // 67 if (argc > 9) { << 79 if ( argc > 9 ) { 68 PrintUsage(); 80 PrintUsage(); 69 return 1; 81 return 1; 70 } 82 } 71 83 72 G4String macro; 84 G4String macro; 73 G4String session; 85 G4String session; >> 86 #ifdef G4MULTITHREADED 74 G4int nThreads = 0; 87 G4int nThreads = 0; >> 88 #endif 75 89 76 G4long randomSeed = 1234; 90 G4long randomSeed = 1234; 77 for (G4int i = 1; i < argc; i = i + 2) { << 91 for ( G4int i=1; i<argc; i=i+2 ) { 78 if (G4String(argv[i]) == "-m") << 92 if ( G4String(argv[i]) == "-m" ) macro = argv[i+1]; 79 macro = argv[i + 1]; << 93 else if ( G4String(argv[i]) == "-u" ) session = argv[i+1]; 80 else if (G4String(argv[i]) == "-u") << 94 else if ( G4String(argv[i]) == "-r" ) randomSeed = atoi(argv[i+1]); 81 session = argv[i + 1]; << 95 #ifdef G4MULTITHREADED 82 else if (G4String(argv[i]) == "-r") << 96 else if ( G4String(argv[i]) == "-t" ) { 83 randomSeed = atoi(argv[i + 1]); << 97 nThreads = G4UIcommand::ConvertToInt(argv[i+1]); 84 else if (G4String(argv[i]) == "-t") { << 85 nThreads = G4UIcommand::ConvertToInt(arg << 86 } 98 } >> 99 #endif 87 else { 100 else { 88 PrintUsage(); 101 PrintUsage(); 89 return 1; 102 return 1; 90 } 103 } 91 } 104 } 92 105 93 // Instantiate G4UIExecutive for interactive 106 // Instantiate G4UIExecutive for interactive mode 94 G4UIExecutive* ui = nullptr; 107 G4UIExecutive* ui = nullptr; 95 if (macro.size() == 0) { << 108 if ( macro.size() == 0 ) { 96 ui = new G4UIExecutive(argc, argv, session 109 ui = new G4UIExecutive(argc, argv, session); 97 } 110 } 98 111 >> 112 // Choose the Random engine >> 113 // >> 114 G4Random::setTheEngine(new CLHEP::RanecuEngine); >> 115 99 // Construct the default run manager 116 // Construct the default run manager 100 // 117 // 101 auto* runManager = G4RunManagerFactory::Crea << 118 #ifdef G4MULTITHREADED 102 if (nThreads > 0) runManager->SetNumberOfThr << 119 G4MTRunManager * runManager = new G4MTRunManager; >> 120 if ( nThreads > 0 ) runManager->SetNumberOfThreads(nThreads); >> 121 #else >> 122 G4RunManager * runManager = new G4RunManager; >> 123 #endif 103 124 104 // Seed the random number generator manually 125 // Seed the random number generator manually 105 G4Random::setTheSeed(randomSeed); 126 G4Random::setTheSeed(randomSeed); 106 127 107 // Set mandatory initialization classes 128 // Set mandatory initialization classes 108 // 129 // 109 // Detector construction 130 // Detector construction 110 runManager->SetUserInitialization(new F06Det 131 runManager->SetUserInitialization(new F06DetectorConstruction()); 111 // Physics list 132 // Physics list 112 runManager->SetUserInitialization(new F06Phy 133 runManager->SetUserInitialization(new F06PhysicsList()); 113 134 114 // Ensure that Transportation considers grav 135 // Ensure that Transportation considers gravity fields. 115 G4Transportation::EnableGravity(true); 136 G4Transportation::EnableGravity(true); 116 137 117 // User action initialization 138 // User action initialization 118 runManager->SetUserInitialization(new F06Act 139 runManager->SetUserInitialization(new F06ActionInitialization()); 119 140 120 // Initialize visualization 141 // Initialize visualization 121 // 142 // 122 G4VisManager* visManager = new G4VisExecutiv 143 G4VisManager* visManager = new G4VisExecutive; 123 // G4VisExecutive can take a verbosity argum 144 // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance. 124 // G4VisManager* visManager = new G4VisExecu 145 // G4VisManager* visManager = new G4VisExecutive("Quiet"); 125 visManager->Initialize(); 146 visManager->Initialize(); 126 147 127 // Get the pointer to the User Interface man 148 // Get the pointer to the User Interface manager 128 // 149 // 129 G4UImanager* UImanager = G4UImanager::GetUIp 150 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 130 << 151 131 if (!ui) { << 152 if ( !ui ) { 132 // batch mode << 153 // batch mode 133 G4String command = "/control/execute "; << 154 G4String command = "/control/execute "; 134 UImanager->ApplyCommand(command + macro); << 155 UImanager->ApplyCommand(command+macro); 135 } << 156 } 136 else { << 157 else 137 UImanager->ApplyCommand("/control/execute << 158 { 138 if (ui->IsGUI()) UImanager->ApplyCommand(" << 159 UImanager->ApplyCommand("/control/execute init_vis.mac"); 139 ui->SessionStart(); << 160 if (ui->IsGUI()) 140 delete ui; << 161 UImanager->ApplyCommand("/control/execute gui.mac"); >> 162 ui->SessionStart(); >> 163 delete ui; 141 } 164 } 142 165 143 // Job termination 166 // Job termination 144 // Free the store: user actions, physics_lis 167 // Free the store: user actions, physics_list and detector_description are 145 // owned and deleted by the 168 // owned and deleted by the run manager, so they should not 146 // be deleted in the main() 169 // be deleted in the main() program ! 147 << 170 148 delete visManager; 171 delete visManager; 149 delete runManager; 172 delete runManager; 150 173 151 return 0; 174 return 0; 152 } 175 } 153 176 154 //....oooOO0OOooo........oooOO0OOooo........oo 177 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 155 178