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