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