| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 /// \file electromagnetic/TestEm5/TestEm5.cc <<
27 /// \brief Main program of the electromagnetic <<
28 // <<
29 // 23 //
>> 24 // $Id: TestEm5.cc,v 1.14 2005/12/06 11:46:46 gcosmo Exp $
>> 25 // GEANT4 tag $Name: geant4-08-00 $
30 // 26 //
31 //....oooOO0OOooo........oooOO0OOooo........oo 27 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oo 28 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 29
34 #include "ActionInitialization.hh" << 30 #include "G4RunManager.hh"
35 #include "DetectorConstruction.hh" <<
36 #include "PhysicsList.hh" <<
37 <<
38 #include "G4RunManagerFactory.hh" <<
39 #include "G4SteppingVerbose.hh" <<
40 #include "G4Types.hh" <<
41 #include "G4UIExecutive.hh" <<
42 #include "G4UImanager.hh" 31 #include "G4UImanager.hh"
43 #include "G4VisExecutive.hh" << 32 #include "G4UIterminal.hh"
>> 33 #include "G4UItcsh.hh"
44 #include "Randomize.hh" 34 #include "Randomize.hh"
45 35
>> 36 #ifdef G4VIS_USE
>> 37 #include "G4VisExecutive.hh"
>> 38 #endif
>> 39
>> 40 #include "DetectorConstruction.hh"
>> 41 #include "PhysicsList.hh"
>> 42 #include "PrimaryGeneratorAction.hh"
>> 43 #include "RunAction.hh"
>> 44 #include "EventAction.hh"
>> 45 #include "TrackingAction.hh"
>> 46 #include "SteppingAction.hh"
>> 47 #include "SteppingVerbose.hh"
>> 48 #include "StackingAction.hh"
>> 49 #include "HistoManager.hh"
>> 50
46 //....oooOO0OOooo........oooOO0OOooo........oo 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 52
48 int main(int argc, char** argv) << 53 int main(int argc,char** argv) {
49 { << 54
50 // detect interactive mode (if no arguments) << 55 //choose the Random engine
51 G4UIExecutive* ui = nullptr; << 56 CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);
52 if (argc == 1) ui = new G4UIExecutive(argc, << 57
53 << 58 //my Verbose output class
54 // Use SteppingVerbose with Unit << 59 G4VSteppingVerbose::SetInstance(new SteppingVerbose);
55 G4int precision = 4; << 60
56 G4SteppingVerbose::UseBestUnit(precision); << 61 // Construct the default run manager
57 << 62 G4RunManager * runManager = new G4RunManager;
58 // Creating run manager << 63
59 auto* runManager = G4RunManagerFactory::Crea <<
60 <<
61 if (argc == 3) { <<
62 G4int nThreads = G4UIcommand::ConvertToInt <<
63 runManager->SetNumberOfThreads(nThreads); <<
64 <<
65 G4cout << "===== TestEm5 is started with " <<
66 << " threads =====" << G4endl; <<
67 } <<
68 // set mandatory initialization classes 64 // set mandatory initialization classes
69 DetectorConstruction* detector = new Detecto << 65 DetectorConstruction* detector;
>> 66 detector = new DetectorConstruction;
70 runManager->SetUserInitialization(detector); 67 runManager->SetUserInitialization(detector);
71 runManager->SetUserInitialization(new Physic 68 runManager->SetUserInitialization(new PhysicsList());
72 69
73 // set user action classes << 70 #ifdef G4VIS_USE
74 runManager->SetUserInitialization(new Action << 71 // visualization manager
75 << 72 G4VisManager* visManager = new G4VisExecutive;
76 // initialize visualization << 73 visManager->Initialize();
77 G4VisManager* visManager = nullptr; << 74 #endif
78 75
79 // get the pointer to the User Interface man << 76 HistoManager* histo = new HistoManager();
80 G4UImanager* UImanager = G4UImanager::GetUIp <<
81 <<
82 if (ui) { <<
83 // interactive mode <<
84 visManager = new G4VisExecutive; <<
85 visManager->Initialize(); <<
86 ui->SessionStart(); <<
87 delete ui; <<
88 } <<
89 else { <<
90 // batch mode <<
91 G4String command = "/control/execute "; <<
92 G4String fileName = argv[1]; <<
93 UImanager->ApplyCommand(command + fileName <<
94 } <<
95 77
>> 78 // set user action classes
>> 79 //
>> 80 //primaryGenerator
>> 81 PrimaryGeneratorAction* primary = new PrimaryGeneratorAction(detector);
>> 82 runManager->SetUserAction(primary);
>> 83
>> 84 //runAction
>> 85 RunAction* runaction = new RunAction(detector,primary,histo);
>> 86 runManager->SetUserAction(runaction);
>> 87
>> 88 //eventAction
>> 89 EventAction* eventaction = new EventAction(runaction,histo);
>> 90 runManager->SetUserAction(eventaction);
>> 91
>> 92 //trackAction
>> 93 TrackingAction* trackingaction = new TrackingAction(detector, runaction,
>> 94 eventaction, histo);
>> 95 runManager->SetUserAction(trackingaction);
>> 96
>> 97 //stepAction
>> 98 SteppingAction* steppingaction = new SteppingAction(detector, runaction,
>> 99 eventaction, histo);
>> 100 runManager->SetUserAction(steppingaction);
>> 101
>> 102 //stackAction
>> 103 StackingAction* stackingaction = new StackingAction(runaction,
>> 104 eventaction, histo);
>> 105 runManager->SetUserAction(stackingaction);
>> 106
>> 107 // get the pointer to the User Interface manager
>> 108 G4UImanager* UI = G4UImanager::GetUIpointer();
>> 109
>> 110 if (argc==1) // Define UI terminal for interactive mode
>> 111 {
>> 112 G4UIsession * session = 0;
>> 113 #ifdef G4UI_USE_TCSH
>> 114 session = new G4UIterminal(new G4UItcsh);
>> 115 #else
>> 116 session = new G4UIterminal();
>> 117 #endif
>> 118 session->SessionStart();
>> 119 delete session;
>> 120 }
>> 121 else // Batch mode
>> 122 {
>> 123 G4String command = "/control/execute ";
>> 124 G4String fileName = argv[1];
>> 125 UI->ApplyCommand(command+fileName);
>> 126 }
>> 127
96 // job termination 128 // job termination
>> 129
>> 130 #ifdef G4VIS_USE
97 delete visManager; 131 delete visManager;
>> 132 #endif
>> 133
>> 134 delete histo;
98 delete runManager; 135 delete runManager;
>> 136
>> 137 return 0;
99 } 138 }
100 139
101 //....oooOO0OOooo........oooOO0OOooo........oo 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 141
102 142