| 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 // $Id: exampleB4a.cc 95508 2016-02-12 13:52:06Z gcosmo $
26 // 27 //
27 /// \file exampleB4a.cc 28 /// \file exampleB4a.cc
28 /// \brief Main program of the B4a example 29 /// \brief Main program of the B4a example
29 30
30 #include "ActionInitialization.hh" << 31 #include "B4DetectorConstruction.hh"
31 #include "DetectorConstruction.hh" << 32 #include "B4aActionInitialization.hh"
32 #include "FTFP_BERT.hh" << 33
>> 34 #ifdef G4MULTITHREADED
>> 35 #include "G4MTRunManager.hh"
>> 36 #else
>> 37 #include "G4RunManager.hh"
>> 38 #endif
33 39
34 #include "G4RunManagerFactory.hh" <<
35 #include "G4SteppingVerbose.hh" <<
36 #include "G4UIExecutive.hh" <<
37 #include "G4UIcommand.hh" <<
38 #include "G4UImanager.hh" 40 #include "G4UImanager.hh"
>> 41 #include "G4UIcommand.hh"
>> 42 #include "FTFP_BERT.hh"
>> 43
>> 44 #include "Randomize.hh"
>> 45
39 #include "G4VisExecutive.hh" 46 #include "G4VisExecutive.hh"
40 // #include "Randomize.hh" << 47 #include "G4UIExecutive.hh"
41 48
42 //....oooOO0OOooo........oooOO0OOooo........oo 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43 50
44 namespace << 51 namespace {
45 { << 52 void PrintUsage() {
46 void PrintUsage() << 53 G4cerr << " Usage: " << G4endl;
47 { << 54 G4cerr << " exampleB4a [-m macro ] [-u UIsession] [-t nThreads]" << G4endl;
48 G4cerr << " Usage: " << G4endl; << 55 G4cerr << " note: -t option is available only for multi-threaded mode."
49 G4cerr << " exampleB4a [-m macro ] [-u UIses << 56 << G4endl;
50 G4cerr << " note: -t option is available o << 57 }
51 } 58 }
52 } // namespace <<
53 59
54 //....oooOO0OOooo........oooOO0OOooo........oo 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55 61
56 int main(int argc, char** argv) << 62 int main(int argc,char** argv)
57 { 63 {
58 // Evaluate arguments 64 // Evaluate arguments
59 // 65 //
60 if (argc > 7) { << 66 if ( argc > 7 ) {
61 PrintUsage(); 67 PrintUsage();
62 return 1; 68 return 1;
63 } 69 }
64 << 70
65 G4String macro; 71 G4String macro;
66 G4String session; 72 G4String session;
67 G4bool verboseBestUnits = true; <<
68 #ifdef G4MULTITHREADED 73 #ifdef G4MULTITHREADED
69 G4int nThreads = 0; 74 G4int nThreads = 0;
70 #endif 75 #endif
71 for (G4int i = 1; i < argc; i = i + 2) { << 76 for ( G4int i=1; i<argc; i=i+2 ) {
72 if (G4String(argv[i]) == "-m") << 77 if ( G4String(argv[i]) == "-m" ) macro = argv[i+1];
73 macro = argv[i + 1]; << 78 else if ( G4String(argv[i]) == "-u" ) session = argv[i+1];
74 else if (G4String(argv[i]) == "-u") <<
75 session = argv[i + 1]; <<
76 #ifdef G4MULTITHREADED 79 #ifdef G4MULTITHREADED
77 else if (G4String(argv[i]) == "-t") { << 80 else if ( G4String(argv[i]) == "-t" ) {
78 nThreads = G4UIcommand::ConvertToInt(arg << 81 nThreads = G4UIcommand::ConvertToInt(argv[i+1]);
79 } 82 }
80 #endif 83 #endif
81 else if (G4String(argv[i]) == "-vDefault") <<
82 verboseBestUnits = false; <<
83 --i; // this option is not followed wit <<
84 } <<
85 else { 84 else {
86 PrintUsage(); 85 PrintUsage();
87 return 1; 86 return 1;
88 } 87 }
89 } << 88 }
90 << 89
91 // Detect interactive mode (if no macro prov 90 // Detect interactive mode (if no macro provided) and define UI session
92 // 91 //
93 G4UIExecutive* ui = nullptr; << 92 G4UIExecutive* ui = 0;
94 if (!macro.size()) { << 93 if ( ! macro.size() ) {
95 ui = new G4UIExecutive(argc, argv, session 94 ui = new G4UIExecutive(argc, argv, session);
96 } 95 }
97 96
98 // Optionally: choose a different Random eng << 97 // Choose the Random engine
99 // G4Random::setTheEngine(new CLHEP::MTwistE << 98 //
100 << 99 G4Random::setTheEngine(new CLHEP::RanecuEngine);
101 // Use G4SteppingVerboseWithUnits << 100
102 if (verboseBestUnits) { <<
103 G4int precision = 4; <<
104 G4SteppingVerbose::UseBestUnit(precision); <<
105 } <<
106 <<
107 // Construct the default run manager 101 // Construct the default run manager
108 // 102 //
109 auto runManager = G4RunManagerFactory::Creat <<
110 #ifdef G4MULTITHREADED 103 #ifdef G4MULTITHREADED
111 if (nThreads > 0) { << 104 G4MTRunManager * runManager = new G4MTRunManager;
>> 105 if ( nThreads > 0 ) {
112 runManager->SetNumberOfThreads(nThreads); 106 runManager->SetNumberOfThreads(nThreads);
113 } << 107 }
>> 108 #else
>> 109 G4RunManager * runManager = new G4RunManager;
114 #endif 110 #endif
115 111
116 // Set mandatory initialization classes 112 // Set mandatory initialization classes
117 // 113 //
118 auto detConstruction = new B4::DetectorConst << 114 B4DetectorConstruction* detConstruction = new B4DetectorConstruction();
119 runManager->SetUserInitialization(detConstru 115 runManager->SetUserInitialization(detConstruction);
120 116
121 auto physicsList = new FTFP_BERT; << 117 G4VModularPhysicsList* physicsList = new FTFP_BERT;
122 runManager->SetUserInitialization(physicsLis 118 runManager->SetUserInitialization(physicsList);
123 << 119
124 auto actionInitialization = new B4a::ActionI << 120 B4aActionInitialization* actionInitialization
>> 121 = new B4aActionInitialization(detConstruction);
125 runManager->SetUserInitialization(actionInit 122 runManager->SetUserInitialization(actionInitialization);
126 << 123
127 // Initialize visualization 124 // Initialize visualization
128 // 125 //
129 auto visManager = new G4VisExecutive; << 126 G4VisManager* visManager = new G4VisExecutive;
130 // G4VisExecutive can take a verbosity argum 127 // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
131 // auto visManager = new G4VisExecutive("Qui << 128 // G4VisManager* visManager = new G4VisExecutive("Quiet");
132 visManager->Initialize(); 129 visManager->Initialize();
133 130
134 // Get the pointer to the User Interface man 131 // Get the pointer to the User Interface manager
135 auto UImanager = G4UImanager::GetUIpointer() << 132 G4UImanager* UImanager = G4UImanager::GetUIpointer();
136 133
137 // Process macro or start UI session 134 // Process macro or start UI session
138 // 135 //
139 if (macro.size()) { << 136 if ( macro.size() ) {
140 // batch mode 137 // batch mode
141 G4String command = "/control/execute "; 138 G4String command = "/control/execute ";
142 UImanager->ApplyCommand(command + macro); << 139 UImanager->ApplyCommand(command+macro);
143 } 140 }
144 else { << 141 else {
145 // interactive mode : define UI session 142 // interactive mode : define UI session
146 UImanager->ApplyCommand("/control/execute 143 UImanager->ApplyCommand("/control/execute init_vis.mac");
147 if (ui->IsGUI()) { 144 if (ui->IsGUI()) {
148 UImanager->ApplyCommand("/control/execut 145 UImanager->ApplyCommand("/control/execute gui.mac");
149 } 146 }
150 ui->SessionStart(); 147 ui->SessionStart();
151 delete ui; 148 delete ui;
152 } 149 }
153 150
154 // Job termination 151 // Job termination
155 // Free the store: user actions, physics_lis 152 // Free the store: user actions, physics_list and detector_description are
156 // owned and deleted by the run manager, so << 153 // owned and deleted by the run manager, so they should not be deleted
157 // in the main() program ! 154 // in the main() program !
158 155
159 delete visManager; 156 delete visManager;
160 delete runManager; 157 delete runManager;
161 } 158 }
162 159
163 //....oooOO0OOooo........oooOO0OOooo........oo 160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
164 161