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 exoticphysics/monopole/monopole.cc 26 /// \file exoticphysics/monopole/monopole.cc 27 /// \brief Main program of the exoticphysics/m 27 /// \brief Main program of the exoticphysics/monopole example 28 // 28 // 29 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 32 33 #include "ActionInitialization.hh" << 33 #include "G4Types.hh" 34 #include "DetectorConstruction.hh" << 35 34 >> 35 #ifdef G4MULTITHREADED >> 36 #include "G4MTRunManager.hh" >> 37 #else >> 38 #include "G4RunManager.hh" >> 39 #endif >> 40 >> 41 #include "G4UImanager.hh" >> 42 #include "Randomize.hh" >> 43 >> 44 #include "G4VisExecutive.hh" >> 45 #include "G4UIExecutive.hh" >> 46 >> 47 #include "DetectorConstruction.hh" 36 #include "G4MonopolePhysics.hh" 48 #include "G4MonopolePhysics.hh" 37 #include "G4PhysListFactory.hh" 49 #include "G4PhysListFactory.hh" 38 #include "G4RunManagerFactory.hh" << 39 #include "G4Types.hh" << 40 #include "G4UIExecutive.hh" << 41 #include "G4UImanager.hh" << 42 #include "G4VModularPhysicsList.hh" 50 #include "G4VModularPhysicsList.hh" 43 #include "G4VisExecutive.hh" << 51 #include "ActionInitialization.hh" 44 #include "Randomize.hh" << 52 45 #include "globals.hh" 53 #include "globals.hh" 46 54 47 //....oooOO0OOooo........oooOO0OOooo........oo 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 48 56 49 namespace << 57 namespace { 50 { << 58 void PrintUsage() { 51 void PrintUsage() << 59 G4cerr 52 { << 60 << " Usage: " << G4endl 53 G4cerr << " Usage: " << G4endl << " monopole << 61 << " exampleB4a [-m macro ] [-s setupMonopole] [-t nThreads]" << G4endl 54 << G4endl << " Note: " << G4endl << 62 << " Note: " << G4endl 55 << " -s should be followed by a co << 63 << " -s should be followed by a composed string, eg. \'1 0 100 GeV\'" << G4endl 56 << " -t option is for multi-thread << 64 << " -t option is available only for multi-threaded mode." << G4endl >> 65 << G4endl; >> 66 } 57 } 67 } 58 } // namespace << 59 68 60 //....oooOO0OOooo........oooOO0OOooo........oo 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 61 70 62 int main(int argc, char** argv) << 71 int main(int argc,char** argv) 63 { 72 { 64 // Evaluate arguments 73 // Evaluate arguments 65 // 74 // 66 if (argc > 7) { << 75 if ( argc > 7 ) { 67 PrintUsage(); 76 PrintUsage(); 68 return 1; 77 return 1; 69 } 78 } 70 79 71 G4String macro; 80 G4String macro; 72 G4String setupMonopole; 81 G4String setupMonopole; 73 G4int nThreads = 1; << 82 #ifdef G4MULTITHREADED 74 for (G4int i = 1; i < argc; i = i + 2) { << 83 G4int nThreads = 0; 75 if (G4String(argv[i]) == "-m") << 84 #endif 76 macro = argv[i + 1]; << 85 for ( G4int i=1; i<argc; i=i+2 ) { 77 else if (G4String(argv[i]) == "-s") << 86 if ( G4String(argv[i]) == "-m" ) macro = argv[i+1]; 78 setupMonopole = argv[i + 1]; << 87 else if ( G4String(argv[i]) == "-s" ) setupMonopole = argv[i+1]; 79 else if (G4String(argv[i]) == "-t") { << 88 #ifdef G4MULTITHREADED 80 nThreads = G4UIcommand::ConvertToInt(arg << 89 else if ( G4String(argv[i]) == "-t" ) { >> 90 nThreads = G4UIcommand::ConvertToInt(argv[i+1]); 81 } 91 } >> 92 #endif 82 else { 93 else { 83 PrintUsage(); 94 PrintUsage(); 84 return 1; 95 return 1; 85 } 96 } 86 } << 97 } 87 << 98 88 // Construct the default run manager << 99 // Instantiate G4UIExecutive if interactive mode 89 auto* runManager = G4RunManagerFactory::Crea << 90 if (nThreads > 0) { << 91 runManager->SetNumberOfThreads(nThreads); << 92 } << 93 G4cout << "===== Example is started with " < << 94 << " threads =====" << G4endl; << 95 << 96 // Instantiate G4UIExecutive if interactive << 97 G4UIExecutive* ui = nullptr; 100 G4UIExecutive* ui = nullptr; 98 if (macro.empty()) { << 101 if ( ! macro.size() ) { 99 ui = new G4UIExecutive(argc, argv); 102 ui = new G4UIExecutive(argc, argv); 100 } 103 } 101 104 102 // get the pointer to the User Interface man << 105 //choose the Random engine 103 G4UImanager* UImanager = G4UImanager::GetUIp << 106 CLHEP::HepRandom::setTheEngine(new CLHEP::MixMaxRng); >> 107 >> 108 // Construct the default run manager >> 109 #ifdef G4MULTITHREADED >> 110 G4MTRunManager* runManager = new G4MTRunManager; >> 111 if ( nThreads > 0 ) { >> 112 runManager->SetNumberOfThreads(nThreads); >> 113 } >> 114 G4cout << "===== Monopole is started with " >> 115 << runManager->GetNumberOfThreads() << " threads =====" << G4endl; >> 116 #else >> 117 G4RunManager* runManager = new G4RunManager(); >> 118 #endif 104 119 105 // create physicsList << 120 //create physicsList 106 // Physics List is defined via environment << 121 // Physics List is defined via environment variable PHYSLIST 107 G4PhysListFactory factory; 122 G4PhysListFactory factory; 108 G4VModularPhysicsList* phys = factory.GetRef 123 G4VModularPhysicsList* phys = factory.GetReferencePhysList("FTFP_BERT"); 109 124 110 // monopole physics is added 125 // monopole physics is added 111 G4MonopolePhysics* theMonopole = new G4Monop << 126 G4MonopolePhysics * theMonopole = new G4MonopolePhysics(); 112 127 >> 128 //get the pointer to the User Interface manager >> 129 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 113 // Setup monopole 130 // Setup monopole 114 if (setupMonopole.size()) { << 131 if ( setupMonopole.size() ) { 115 UImanager->ApplyCommand("/control/verbose 132 UImanager->ApplyCommand("/control/verbose 1"); 116 UImanager->ApplyCommand("/monopole/setup " 133 UImanager->ApplyCommand("/monopole/setup " + setupMonopole); 117 } 134 } 118 135 119 // regsiter monopole physics 136 // regsiter monopole physics 120 phys->RegisterPhysics(theMonopole); 137 phys->RegisterPhysics(theMonopole); 121 << 122 runManager->SetUserInitialization(phys); 138 runManager->SetUserInitialization(phys); 123 139 124 // visualization manager 140 // visualization manager 125 G4VisManager* visManager = nullptr; << 141 G4VisManager* visManager = new G4VisExecutive(); >> 142 visManager->Initialize(); 126 143 127 // set detector construction 144 // set detector construction 128 DetectorConstruction* det = new DetectorCons 145 DetectorConstruction* det = new DetectorConstruction(); 129 runManager->SetUserInitialization(det); 146 runManager->SetUserInitialization(det); 130 147 131 // set user action classes 148 // set user action classes 132 runManager->SetUserInitialization(new Action 149 runManager->SetUserInitialization(new ActionInitialization(det)); 133 150 134 // Process macro or start UI session 151 // Process macro or start UI session 135 // 152 // 136 if (macro.size()) { << 153 if ( macro.size() ) { 137 // batch mode 154 // batch mode 138 G4String command = "/control/execute "; 155 G4String command = "/control/execute "; 139 UImanager->ApplyCommand(command + macro); << 156 UImanager->ApplyCommand(command+macro); 140 } 157 } 141 else { << 158 else { 142 // interactive mode : define UI session 159 // interactive mode : define UI session 143 visManager = new G4VisExecutive(); << 144 visManager->Initialize(); << 145 UImanager->ApplyCommand("/control/execute 160 UImanager->ApplyCommand("/control/execute init_vis.mac"); 146 ui->SessionStart(); 161 ui->SessionStart(); 147 delete ui; 162 delete ui; 148 } 163 } 149 164 150 delete visManager; 165 delete visManager; 151 166 152 // job termination 167 // job termination 153 delete runManager; 168 delete runManager; 154 169 155 return 0; 170 return 0; 156 } 171 } 157 172 158 //....oooOO0OOooo........oooOO0OOooo........oo 173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 159 174