| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //
28 //
29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31 #include "G4RunManagerFactory.hh"
32 #include "G4SteppingVerbose.hh"
33 #include "G4Types.hh"
34 #include "G4UIExecutive.hh"
35 #include "G4UImanager.hh"
36 #include "G4VisExecutive.hh"
37
38 #include "ActionInitialization.hh"
39 #include "DetectorConstruction.hh"
40 #include "PhysicsList.hh"
41
42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
43
44 int main(int argc, char** argv)
45 {
46 // detect interactive mode (if no arguments) and define UI session
47 G4UIExecutive* ui = nullptr;
48 if (argc == 1) ui = new G4UIExecutive(argc, argv);
49
50 // Use SteppingVerbose with Unit
51 G4int precision = 4;
52 G4SteppingVerbose::UseBestUnit(precision);
53
54 if (argc == 2) {
55 G4Exception("argument", "Fatal error in Argument", FatalErrorInArgument,
56 "Type of simulation must be precised: either ./stim_pixe_tomography -p "
57 "pixe3d.mac or ./stim_pixe_tomography -s pixe3d.mac");
58 }
59 G4bool isPIXE = false;
60 if (argc > 2) {
61 G4String s;
62 s = argv[1];
63 if (s == "-p")
64 isPIXE = true;
65 else if (s == "-s")
66 isPIXE = false;
67 else {
68 G4Exception("argument", "Fatal error in Argument", FatalErrorInArgument,
69 "Type of simulation must be precised: /stim_pixe_tomography -p "
70 "pixe3d.mac or ./stim_pixe_tomography -s pixe3d.mac");
71 }
72 }
73 G4int nThreads = 4;
74 if (argc == 4) nThreads = G4UIcommand::ConvertToInt(argv[3]);
75 auto* runManager = G4RunManagerFactory::CreateRunManager();
76 runManager->SetNumberOfThreads(nThreads);
77
78 // UserInitialization classes - mandatory
79 runManager->SetUserInitialization(new DetectorConstruction());
80 runManager->SetUserInitialization(new PhysicsList());
81 runManager->SetUserInitialization(new ActionInitialization(isPIXE));
82
83 // initialize visualization
84 auto visManager = new G4VisExecutive;
85 visManager->Initialize();
86
87 // get the pointer to the User Interface manager
88 auto UImanager = G4UImanager::GetUIpointer();
89
90 // User interactions
91 // Define (G)UI for interactive mode
92 if (ui) {
93 UImanager->ApplyCommand("/control/execute init_vis.mac");
94 ui->SessionStart();
95 delete ui;
96 }
97 else {
98 G4String command;
99 command = "/control/execute ";
100 G4String fileName;
101 fileName = argv[2];
102 UImanager->ApplyCommand(command + fileName);
103 }
104
105 // job termination
106 delete visManager;
107 delete runManager;
108 return 0;
109 }
110
111 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
112