57{
58
59
60 if ( argc > 7 ) {
62 return 1;
63 }
64
67 G4bool verboseBestUnits = true;
68#ifdef G4MULTITHREADED
69 G4int nThreads = 0;
70#endif
71 for ( G4int i=1; i<argc; i=i+2 ) {
72 if (
G4String(argv[i]) ==
"-m" ) macro = argv[i+1];
73 else if (
G4String(argv[i]) ==
"-u" ) session = argv[i+1];
74#ifdef G4MULTITHREADED
75 else if (
G4String(argv[i]) ==
"-t" ) {
76 nThreads = G4UIcommand::ConvertToInt(argv[i+1]);
77 }
78#endif
79 else if (
G4String(argv[i]) ==
"-vDefault" ) {
80 verboseBestUnits = false;
81 --i;
82 }
83 else {
85 return 1;
86 }
87 }
88
89
90
92 if ( ! macro.size() ) {
94 }
95
96
97
98
99
100 if ( verboseBestUnits ) {
101 G4int precision = 4;
102 G4SteppingVerbose::UseBestUnit(precision);
103 }
104
105
106
107 auto runManager =
108 G4RunManagerFactory::CreateRunManager(G4RunManagerType::Default);
109#ifdef G4MULTITHREADED
110 if ( nThreads > 0 ) {
111 runManager->SetNumberOfThreads(nThreads);
112 }
113#endif
114
115
116
118 runManager->SetUserInitialization(detConstruction);
119
121 runManager->SetUserInitialization(physicsList);
122
124 runManager->SetUserInitialization(actionInitialization);
125
126
128
129
130 visManager->Initialize();
131
132
133 auto UImanager = G4UImanager::GetUIpointer();
134
135
136
137 if ( macro.size() ) {
138
139 G4String command =
"/control/execute ";
140 UImanager->ApplyCommand(command+macro);
141 }
142 else {
143
144 UImanager->ApplyCommand("/control/execute init_vis.mac");
145 if (ui->IsGUI()) {
146 UImanager->ApplyCommand("/control/execute gui.mac");
147 }
148 ui->SessionStart();
149 delete ui;
150 }
151
152
153
154
155
156
157 delete visManager;
158 delete runManager;
159}
Action initialization class.
Detector construction class to define materials and geometry.