| 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 OpNovice/OpNovice.cc 26 /// \file OpNovice/OpNovice.cc
27 /// \brief Main program of the OpNovice exampl 27 /// \brief Main program of the OpNovice example
28 // 28 //
29 // 29 //
30 // 30 //
31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 // 33 //
34 // Description: Test of Continuous Process G4C 34 // Description: Test of Continuous Process G4Cerenkov
35 // and RestDiscrete Process G4Sci 35 // and RestDiscrete Process G4Scintillation
36 // -- Generation Cerenkov Photons 36 // -- Generation Cerenkov Photons --
37 // -- Generation Scintillation Ph 37 // -- Generation Scintillation Photons --
38 // -- Transport of optical Photon 38 // -- Transport of optical Photons --
39 // Version: 5.0 39 // Version: 5.0
40 // Created: 1996-04-30 40 // Created: 1996-04-30
41 // Author: Juliet Armstrong 41 // Author: Juliet Armstrong
42 // << 42 // mail: gum@triumf.ca
>> 43 //
43 //....oooOO0OOooo........oooOO0OOooo........oo 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 45
45 #include "OpNoviceDetectorConstruction.hh" << 46 #ifdef G4MULTITHREADED
46 #ifdef GEANT4_USE_GDML << 47 #include "G4MTRunManager.hh"
47 # include "OpNoviceGDMLDetectorConstruction.h << 48 #else
>> 49 #include "G4RunManager.hh"
48 #endif 50 #endif
49 #include "FTFP_BERT.hh" <<
50 #include "OpNoviceActionInitialization.hh" <<
51 51
52 #include "G4EmStandardPhysics_option4.hh" <<
53 #include "G4OpticalPhysics.hh" <<
54 #include "G4RunManagerFactory.hh" <<
55 #include "G4Types.hh" <<
56 #include "G4UIExecutive.hh" <<
57 #include "G4UImanager.hh" 52 #include "G4UImanager.hh"
>> 53
>> 54 #include "OpNovicePhysicsList.hh"
>> 55 #include "OpNoviceDetectorConstruction.hh"
>> 56
>> 57 #include "OpNoviceActionInitialization.hh"
>> 58
>> 59 #ifdef G4VIS_USE
58 #include "G4VisExecutive.hh" 60 #include "G4VisExecutive.hh"
>> 61 #endif
59 62
60 //....oooOO0OOooo........oooOO0OOooo........oo << 63 #ifdef G4UI_USE
61 namespace << 64 #include "G4UIExecutive.hh"
62 { <<
63 void PrintUsage() <<
64 { <<
65 G4cerr << " Usage: " << G4endl; <<
66 #ifdef GEANT4_USE_GDML <<
67 G4cerr << " OpNovice [-g gdmlfile] [-m macro <<
68 "nThreads] [-r seed] " <<
69 << G4endl; <<
70 #else <<
71 G4cerr << " OpNovice [-m macro ] [-u UIsess <<
72 #endif 65 #endif
73 G4cerr << " note: -t option is available o << 66
>> 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 68 namespace {
>> 69 void PrintUsage() {
>> 70 G4cerr << " Usage: " << G4endl;
>> 71 G4cerr << " OpNovice [-m macro ] [-u UIsession] [-t nThreads] [-r seed] "
>> 72 << G4endl;
>> 73 G4cerr << " note: -t option is available only for multi-threaded mode."
>> 74 << G4endl;
>> 75 }
74 } 76 }
75 } // namespace <<
76 77
77 //....oooOO0OOooo........oooOO0OOooo........oo 78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78 79
79 int main(int argc, char** argv) << 80 int main(int argc,char** argv)
80 { 81 {
81 // Evaluate arguments 82 // Evaluate arguments
82 // 83 //
83 if (argc > 9) { << 84 if ( argc > 9 ) {
84 PrintUsage(); 85 PrintUsage();
85 return 1; 86 return 1;
86 } 87 }
87 G4String gdmlfile; << 88
88 G4String macro; 89 G4String macro;
89 G4String session; 90 G4String session;
90 #ifdef G4MULTITHREADED 91 #ifdef G4MULTITHREADED
91 G4int nThreads = 0; 92 G4int nThreads = 0;
92 #endif 93 #endif
93 94
94 G4long myseed = 345354; 95 G4long myseed = 345354;
95 for (G4int i = 1; i < argc; i = i + 2) { << 96 for ( G4int i=1; i<argc; i=i+2 ) {
96 if (G4String(argv[i]) == "-g") << 97 if ( G4String(argv[i]) == "-m" ) macro = argv[i+1];
97 gdmlfile = argv[i + 1]; << 98 else if ( G4String(argv[i]) == "-u" ) session = argv[i+1];
98 else if (G4String(argv[i]) == "-m") << 99 else if ( G4String(argv[i]) == "-r" ) myseed = atoi(argv[i+1]);
99 macro = argv[i + 1]; <<
100 else if (G4String(argv[i]) == "-u") <<
101 session = argv[i + 1]; <<
102 else if (G4String(argv[i]) == "-r") <<
103 myseed = atoi(argv[i + 1]); <<
104 #ifdef G4MULTITHREADED 100 #ifdef G4MULTITHREADED
105 else if (G4String(argv[i]) == "-t") { << 101 else if ( G4String(argv[i]) == "-t" ) {
106 nThreads = G4UIcommand::ConvertToInt(arg << 102 nThreads = G4UIcommand::ConvertToInt(argv[i+1]);
107 } 103 }
108 #endif 104 #endif
109 else { 105 else {
110 PrintUsage(); 106 PrintUsage();
111 return 1; 107 return 1;
112 } 108 }
113 } 109 }
114 110
115 // Instantiate G4UIExecutive if interactive << 111 // Choose the Random engine
116 G4UIExecutive* ui = nullptr; << 112 //
117 if (macro.size() == 0) { << 113 G4Random::setTheEngine(new CLHEP::RanecuEngine);
118 ui = new G4UIExecutive(argc, argv); <<
119 } <<
120 114
121 // Construct the default run manager 115 // Construct the default run manager
122 auto runManager = G4RunManagerFactory::Creat << 116 //
123 #ifdef G4MULTITHREADED 117 #ifdef G4MULTITHREADED
124 if (nThreads > 0) runManager->SetNumberOfThr << 118 G4MTRunManager * runManager = new G4MTRunManager;
>> 119 if ( nThreads > 0 ) runManager->SetNumberOfThreads(nThreads);
>> 120 #else
>> 121 G4RunManager * runManager = new G4RunManager;
125 #endif 122 #endif
126 123
127 // Seed the random number generator manually 124 // Seed the random number generator manually
128 G4Random::setTheSeed(myseed); 125 G4Random::setTheSeed(myseed);
129 126
130 // Set mandatory initialization classes 127 // Set mandatory initialization classes
131 // 128 //
132 // Detector construction 129 // Detector construction
133 if (gdmlfile != "") { << 130 runManager-> SetUserInitialization(new OpNoviceDetectorConstruction());
134 #ifdef GEANT4_USE_GDML <<
135 runManager->SetUserInitialization(new OpNo <<
136 #else <<
137 G4cout << "Error! Input gdml file specifie <<
138 << "built with gdml support." << G4 <<
139 return 1; <<
140 #endif <<
141 } <<
142 else { <<
143 runManager->SetUserInitialization(new OpNo <<
144 } <<
145 // Physics list 131 // Physics list
146 G4VModularPhysicsList* physicsList = new FTF << 132 runManager-> SetUserInitialization(new OpNovicePhysicsList());
147 physicsList->ReplacePhysics(new G4EmStandard << 133 // User action initialization
148 auto opticalPhysics = new G4OpticalPhysics() <<
149 physicsList->RegisterPhysics(opticalPhysics) <<
150 runManager->SetUserInitialization(physicsLis <<
151 <<
152 runManager->SetUserInitialization(new OpNovi 134 runManager->SetUserInitialization(new OpNoviceActionInitialization());
153 135
154 G4VisManager* visManager = new G4VisExecutiv << 136 // Initialize G4 kernel
155 visManager->Initialize(); << 137 //
>> 138 runManager->Initialize();
156 139
157 G4UImanager* UImanager = G4UImanager::GetUIp << 140 #ifdef G4VIS_USE
>> 141 // Initialize visualization
>> 142 //
>> 143 G4VisManager* visManager = new G4VisExecutive;
>> 144 // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
>> 145 // G4VisManager* visManager = new G4VisExecutive("Quiet");
>> 146 visManager->Initialize();
>> 147 #endif
158 148
159 if (macro.size()) { << 149 // Get the pointer to the User Interface manager
160 G4String command = "/control/execute "; << 150 //
161 UImanager->ApplyCommand(command + macro); << 151 G4UImanager* UImanager = G4UImanager::GetUIpointer();
>> 152
>> 153 if ( macro.size() ) {
>> 154 // Batch mode
>> 155 G4String command = "/control/execute ";
>> 156 UImanager->ApplyCommand(command+macro);
162 } 157 }
163 else // Define UI session for interactive m << 158 else // Define UI session for interactive mode
164 { 159 {
165 UImanager->ApplyCommand("/control/execute << 160 #ifdef G4UI_USE
166 if (ui->IsGUI()) UImanager->ApplyCommand(" << 161 G4UIExecutive * ui = new G4UIExecutive(argc,argv,session);
167 ui->SessionStart(); << 162 #ifdef G4VIS_USE
168 delete ui; << 163 UImanager->ApplyCommand("/control/execute vis.mac");
>> 164 #else
>> 165 UImanager->ApplyCommand("/control/execute OpNovice.in");
>> 166 #endif
>> 167 if (ui->IsGUI())
>> 168 UImanager->ApplyCommand("/control/execute gui.mac");
>> 169 ui->SessionStart();
>> 170 delete ui;
>> 171 #endif
169 } 172 }
170 173
>> 174 // Job termination
>> 175 // Free the store: user actions, physics_list and detector_description are
>> 176 // owned and deleted by the run manager, so they should not
>> 177 // be deleted in the main() program !
>> 178
>> 179 #ifdef G4VIS_USE
171 delete visManager; 180 delete visManager;
>> 181 #endif
172 delete runManager; 182 delete runManager;
173 183
174 return 0; 184 return 0;
175 } 185 }
176 186
177 //....oooOO0OOooo........oooOO0OOooo........oo 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178 188