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