| 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 // ------------------------------------------- 26 // --------------------------------------------------------------
27 // GEANT 4 - Brachytherapy exa 27 // GEANT 4 - Brachytherapy example
28 // ------------------------------------------- 28 // --------------------------------------------------------------
29 // 29 //
30 // Code developed currently by: 30 // Code developed currently by:
31 // S.Guatelli & D. Cutajar 31 // S.Guatelli & D. Cutajar
32 32
33 // 33 //
34 // ******************************* 34 // *******************************
35 // * * 35 // * *
36 // * Brachy.cc * 36 // * Brachy.cc *
37 // * * 37 // * *
38 // ******************************* 38 // *******************************
39 // 39 //
40 // 40 //
41 #include "G4Types.hh" 41 #include "G4Types.hh"
42 #include "G4RunManagerFactory.hh" 42 #include "G4RunManagerFactory.hh"
43 #include "G4UImanager.hh" 43 #include "G4UImanager.hh"
44 #include "G4UIExecutive.hh" 44 #include "G4UIExecutive.hh"
45 #include "BrachyActionInitialization.hh" 45 #include "BrachyActionInitialization.hh"
46 #include "G4VisExecutive.hh" 46 #include "G4VisExecutive.hh"
>> 47 #include "BrachyAnalysisManager.hh"
47 #include "BrachyDetectorConstruction.hh" 48 #include "BrachyDetectorConstruction.hh"
48 #include "BrachyPhysicsList.hh" 49 #include "BrachyPhysicsList.hh"
49 #include "BrachyPrimaryGeneratorAction.hh" 50 #include "BrachyPrimaryGeneratorAction.hh"
50 #include "G4SDManager.hh" 51 #include "G4SDManager.hh"
51 #include "Randomize.hh" 52 #include "Randomize.hh"
52 #include "G4RunManager.hh" 53 #include "G4RunManager.hh"
53 #include "G4SDManager.hh" 54 #include "G4SDManager.hh"
54 #include "G4UImanager.hh" 55 #include "G4UImanager.hh"
55 #include "G4UImessenger.hh" 56 #include "G4UImessenger.hh"
56 57
57 #include "G4ScoringManager.hh" 58 #include "G4ScoringManager.hh"
58 #include "G4UIExecutive.hh" 59 #include "G4UIExecutive.hh"
59 60
60 #include "G4ScoringManager.hh" 61 #include "G4ScoringManager.hh"
61 #include "BrachyUserScoreWriter.hh" 62 #include "BrachyUserScoreWriter.hh"
62 63
63 int main(int argc ,char ** argv) 64 int main(int argc ,char ** argv)
64 65
65 { 66 {
66 // Construct the default run manager 67 // Construct the default run manager
67 // 68 //
68 auto* pRunManager = G4RunManagerFactory::Cre 69 auto* pRunManager = G4RunManagerFactory::CreateRunManager();
69 G4int nThreads = 4; 70 G4int nThreads = 4;
70 pRunManager->SetNumberOfThreads(nThreads); 71 pRunManager->SetNumberOfThreads(nThreads);
71 72
72 G4cout << "***********************" << G4end 73 G4cout << "***********************" << G4endl;
73 G4cout << "*** Seed: " << G4Random::getTheSe 74 G4cout << "*** Seed: " << G4Random::getTheSeed() << " ***" << G4endl;
74 G4cout << "***********************" << G4end 75 G4cout << "***********************" << G4endl;
75 76
76 // Access to the Scoring Manager pointer 77 // Access to the Scoring Manager pointer
77 78
78 G4ScoringManager* scoringManager = G4Scoring 79 G4ScoringManager* scoringManager = G4ScoringManager::GetScoringManager();
79 80
80 // Overwrite the default output file with us 81 // Overwrite the default output file with user-defined one
81 scoringManager->SetScoreWriter(new BrachyUse 82 scoringManager->SetScoreWriter(new BrachyUserScoreWriter());
82 83
83 // Initialize the physics component 84 // Initialize the physics component
84 pRunManager -> SetUserInitialization(new Bra 85 pRunManager -> SetUserInitialization(new BrachyPhysicsList);
85 86
86 // Initialize the detector component 87 // Initialize the detector component
87 auto pDetectorConstruction = new BrachyDete << 88 BrachyDetectorConstruction *pDetectorConstruction = new BrachyDetectorConstruction();
88 pRunManager -> SetUserInitialization(pDetect 89 pRunManager -> SetUserInitialization(pDetectorConstruction);
89 90
90 // User action initialization 91 // User action initialization
91 92
92 auto actions = new BrachyActionInitialization << 93 BrachyActionInitialization* actions = new BrachyActionInitialization();
93 pRunManager->SetUserInitialization(actions); 94 pRunManager->SetUserInitialization(actions);
94 95
95 //Initialize G4 kernel 96 //Initialize G4 kernel
96 // pRunManager -> Initialize(); 97 // pRunManager -> Initialize();
97 98
98 // Visualization manager 99 // Visualization manager
99 G4VisManager* visManager = new G4VisExecutiv 100 G4VisManager* visManager = new G4VisExecutive;
100 visManager->Initialize(); 101 visManager->Initialize();
101 102
102 // get the pointer to the User Interface man 103 // get the pointer to the User Interface manager
103 G4UImanager* UImanager = G4UImanager::GetUIp 104 G4UImanager* UImanager = G4UImanager::GetUIpointer();
104 if (argc == 1) // Define UI session for in 105 if (argc == 1) // Define UI session for interactive mode.
105 { 106 {
106 auto ui = new G4UIExecutive(argc, argv); << 107 G4UIExecutive* ui = new G4UIExecutive(argc, argv);
107 G4cout << " UI session starts ..." << G4 108 G4cout << " UI session starts ..." << G4endl;
108 UImanager -> ApplyCommand("/control/exec 109 UImanager -> ApplyCommand("/control/execute VisualisationMacro.mac");
109 ui -> SessionStart(); 110 ui -> SessionStart();
110 delete ui; 111 delete ui;
111 } 112 }
112 else // Batch mode 113 else // Batch mode
113 { 114 {
114 G4String command = "/control/execute "; 115 G4String command = "/control/execute ";
115 G4String fileName = argv[1]; 116 G4String fileName = argv[1];
116 UImanager -> ApplyCommand(command+fileNa 117 UImanager -> ApplyCommand(command+fileName);
117 } 118 }
118 119
119 // Job termination 120 // Job termination
120 // Close the root file 121 // Close the root file
>> 122
>> 123 delete G4AnalysisManager::Instance();
121 124
122 delete visManager; 125 delete visManager;
123 delete pRunManager; 126 delete pRunManager;
124 127
125 return 0; 128 return 0;
126 } 129 }
127 130