| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // This example is provided by the Geant4-DNA
27 // Any report or published results obtained us
28 // shall cite the following Geant4-DNA collabo
29 // Med. Phys. 37 (2010) 4692-4708
30 // The Geant4-DNA web site is available at htt
31 //
32 //
33 /// \file clustering.cc
34 /// \brief clustering example
35
36 #include "ActionInitialization.hh"
37 #include "CommandLineParser.hh"
38 #include "DetectorConstruction.hh"
39 #include "PhysicsList.hh"
40
41 #include "G4DNAChemistryManager.hh"
42 #include "G4RunManagerFactory.hh"
43 #include "G4UIExecutive.hh"
44 #include "G4UImanager.hh"
45 #include "G4VisExecutive.hh"
46
47 using namespace G4DNAPARSER;
48 CommandLineParser* parser(0);
49
50 //....oooOO0OOooo........oooOO0OOooo........oo
51
52 void Parse(int& argc, char** argv);
53
54 //....oooOO0OOooo........oooOO0OOooo........oo
55
56 int main(int argc, char** argv)
57 {
58 //////////
59 // Parse options given in commandLine
60 //
61 Parse(argc, argv);
62
63 //////////
64 // Construct the run manager according to wh
65 //
66 Command* commandLine(0);
67
68 auto* runManager = G4RunManagerFactory::Crea
69 if ((commandLine = parser->GetCommandIfActiv
70 int nThreads = 2;
71 if (commandLine->GetOption() == "NMAX") {
72 nThreads = G4Threading::G4GetNumberOfCor
73 }
74 else {
75 nThreads = G4UIcommand::ConvertToInt(com
76 }
77 runManager->SetNumberOfThreads(nThreads);
78
79 G4cout << "===== clustering is started wit
80 << " threads =====" << G4endl;
81 }
82
83 //////////
84 // Set mandatory user initialization classes
85 //
86 DetectorConstruction* detector = new Detecto
87 runManager->SetUserInitialization(new Physic
88 runManager->SetUserInitialization(detector);
89 runManager->SetUserInitialization(new Action
90
91 // Initialize G4 kernel
92 runManager->Initialize();
93
94 // Initialize visualization
95 G4VisManager* visManager = new G4VisExecutiv
96 visManager->Initialize();
97
98 // Get the pointer to the User Interface man
99 G4UImanager* UImanager = G4UImanager::GetUIp
100 G4UIExecutive* ui(0);
101
102 // interactive mode : define UI session
103 if ((commandLine = parser->GetCommandIfActiv
104 ui = new G4UIExecutive(argc, argv, command
105
106 if (ui->IsGUI()) UImanager->ApplyCommand("
107
108 if (parser->GetCommandIfActive("-novis") =
109 // visualization is used by default
110 {
111 if ((commandLine = parser->GetCommandIfA
112 // select a visualization driver if need
113 {
114 UImanager->ApplyCommand(G4String("/vis
115 }
116 else
117 // by default OGL is used
118 {
119 UImanager->ApplyCommand("/vis/open OGL
120 }
121 UImanager->ApplyCommand("/control/execut
122 }
123 }
124 else
125 // to be use visualization file (= store the
126 // an external file:
127 // ASCIITree ; DAWNFILE ; HepRepFile ; VRML
128 {
129 if ((commandLine = parser->GetCommandIfAct
130 UImanager->ApplyCommand(G4String("/vis/o
131 UImanager->ApplyCommand("/control/execut
132 }
133 }
134
135 if ((commandLine = parser->GetCommandIfActiv
136 G4String command = "/control/execute ";
137 UImanager->ApplyCommand(command + commandL
138 }
139 else {
140 UImanager->ApplyCommand("/control/execute
141 }
142
143 if ((commandLine = parser->GetCommandIfActiv
144 ui->SessionStart();
145 delete ui;
146 }
147
148 // Job termination
149 // Free the store: user actions, physics_lis
150 // owned and deleted by the run manager, so
151 // in the main() program !
152
153 delete visManager;
154 delete runManager;
155
156 CommandLineParser::DeleteInstance();
157
158 return 0;
159 }
160
161 //....oooOO0OOooo........oooOO0OOooo........oo
162
163 void Parse(int& argc, char** argv)
164 {
165 //////////
166 // Parse options given in commandLine
167 //
168 parser = CommandLineParser::GetParser();
169
170 parser->AddCommand("-gui", Command::OptionNo
171 "Select geant4 UI or just
172
173 parser->AddCommand("-mac", Command::WithOpti
174
175 // You cann your own command, as for instanc
176 // parser->AddCommand("-seed",
177 // Command::WithOption,
178 // "Give a seed value in
179 // it is then up to you to manage this optio
180
181 parser->AddCommand("-mt", Command::WithOptio
182 "Launch in MT mode (event
183 " NOT RECOMMENDED WITH CH
184 "2");
185
186 parser->AddCommand("-vis", Command::WithOpti
187 "OGL 600x600-0+0");
188
189 parser->AddCommand("-novis", Command::Withou
190
191 //////////
192 // If -h or --help is given in option : prin
193 //
194 if (parser->Parse(argc, argv) != 0) // help
195 {
196 // if you are using ROOT, create a TApplic
197 // to print the help from ROOT as well
198 CommandLineParser::DeleteInstance();
199 // Exit is kept
200 std::exit(0);
201 }
202
203 ///////////
204 // Kill application if wrong argument in com
205 //
206 if (parser->CheckIfNotHandledOptionsExists(a
207 // if you are using ROOT, you should initi
208 // before this condition
209 abort();
210 }
211 }
212