| 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 //
>> 27 // $Id: Brachy.cc
>> 28 // GEANT4 tag $Name: geant4-09-02-patch-03 $
>> 29 //
26 // ------------------------------------------- 30 // --------------------------------------------------------------
27 // GEANT 4 - Brachytherapy exa 31 // GEANT 4 - Brachytherapy example
28 // ------------------------------------------- 32 // --------------------------------------------------------------
29 // 33 //
30 // Code developed currently by: << 34 // Code developed by:
31 // S.Guatelli & D. Cutajar << 35 // S. Agostinelli, F. Foppiano, S. Garelli , M. Tropeano, S.Guatelli
32 36
33 // 37 //
34 // ******************************* 38 // *******************************
35 // * * 39 // * *
36 // * Brachy.cc * 40 // * Brachy.cc *
37 // * * 41 // * *
38 // ******************************* 42 // *******************************
39 // 43 //
40 // 44 //
41 #include "G4Types.hh" << 45
42 #include "G4RunManagerFactory.hh" << 46 #include "G4RunManager.hh"
43 #include "G4UImanager.hh" 47 #include "G4UImanager.hh"
44 #include "G4UIExecutive.hh" << 48 #include "G4UIterminal.hh"
45 #include "BrachyActionInitialization.hh" << 49 #include "G4UItcsh.hh"
>> 50 #include "BrachyFactoryIr.hh"
>> 51 #ifdef G4UI_USE_XM
>> 52 #include "G4UIXm.hh"
>> 53 #endif
>> 54
>> 55 #ifdef G4VIS_USE
46 #include "G4VisExecutive.hh" 56 #include "G4VisExecutive.hh"
>> 57 #endif
>> 58
>> 59 #include "BrachyEventAction.hh"
47 #include "BrachyDetectorConstruction.hh" 60 #include "BrachyDetectorConstruction.hh"
48 #include "BrachyPhysicsList.hh" 61 #include "BrachyPhysicsList.hh"
>> 62 #include "BrachyPhantomSD.hh"
49 #include "BrachyPrimaryGeneratorAction.hh" 63 #include "BrachyPrimaryGeneratorAction.hh"
50 #include "G4SDManager.hh" 64 #include "G4SDManager.hh"
51 #include "Randomize.hh" << 65 #include "BrachyRunAction.hh"
>> 66 #include "Randomize.hh"
52 #include "G4RunManager.hh" 67 #include "G4RunManager.hh"
53 #include "G4SDManager.hh" 68 #include "G4SDManager.hh"
54 #include "G4UImanager.hh" 69 #include "G4UImanager.hh"
55 #include "G4UImessenger.hh" 70 #include "G4UImessenger.hh"
56 71
57 #include "G4ScoringManager.hh" << 72 #ifdef G4ANALYSIS_USE
58 #include "G4UIExecutive.hh" << 73 #include "BrachyAnalysisManager.hh"
59 << 74 #endif
60 #include "G4ScoringManager.hh" <<
61 #include "BrachyUserScoreWriter.hh" <<
62 75
63 int main(int argc ,char ** argv) 76 int main(int argc ,char ** argv)
64 77
65 { 78 {
66 // Construct the default run manager << 79 G4RunManager* pRunManager = new G4RunManager;
67 // <<
68 auto* pRunManager = G4RunManagerFactory::Cre <<
69 G4int nThreads = 4; <<
70 pRunManager->SetNumberOfThreads(nThreads); <<
71 <<
72 G4cout << "***********************" << G4end <<
73 G4cout << "*** Seed: " << G4Random::getTheSe <<
74 G4cout << "***********************" << G4end <<
75 <<
76 // Access to the Scoring Manager pointer <<
77 <<
78 G4ScoringManager* scoringManager = G4Scoring <<
79 80
80 // Overwrite the default output file with us << 81 G4String sensitiveDetectorName = "Phantom";
81 scoringManager->SetScoreWriter(new BrachyUse <<
82 82
83 // Initialize the physics component 83 // Initialize the physics component
84 pRunManager -> SetUserInitialization(new Bra 84 pRunManager -> SetUserInitialization(new BrachyPhysicsList);
85 85
86 // Initialize the detector component 86 // Initialize the detector component
87 auto pDetectorConstruction = new BrachyDete << 87 BrachyDetectorConstruction *pDetectorConstruction = new BrachyDetectorConstruction(sensitiveDetectorName);
88 pRunManager -> SetUserInitialization(pDetect 88 pRunManager -> SetUserInitialization(pDetectorConstruction);
89 89
90 // User action initialization << 90 // Initialize the primary particles
>> 91 BrachyPrimaryGeneratorAction* primary = new BrachyPrimaryGeneratorAction();
>> 92 pRunManager -> SetUserAction(primary);
>> 93
>> 94 // Initialize Optional User Action
>> 95 BrachyEventAction *pEventAction = new BrachyEventAction();
>> 96 pRunManager -> SetUserAction(pEventAction );
91 97
92 auto actions = new BrachyActionInitialization << 98 BrachyRunAction *pRunAction = new BrachyRunAction();
93 pRunManager->SetUserInitialization(actions); << 99 pRunManager -> SetUserAction(pRunAction);
94 <<
95 //Initialize G4 kernel <<
96 // pRunManager -> Initialize(); <<
97 100
>> 101 //// Initialize the Visualization component
>> 102 #ifdef G4VIS_USE
98 // Visualization manager 103 // Visualization manager
99 G4VisManager* visManager = new G4VisExecutiv 104 G4VisManager* visManager = new G4VisExecutive;
100 visManager->Initialize(); 105 visManager->Initialize();
101 << 106 #endif
102 // get the pointer to the User Interface man << 107
103 G4UImanager* UImanager = G4UImanager::GetUIp << 108 // Output environment variables concerning the analysis:
>> 109 #ifdef G4ANALYSIS_USE
>> 110 G4cout << G4endl << G4endl << G4endl
>> 111 << " User Environment " << G4endl
>> 112 << " Using AIDA 3.2.1 analysis " << G4endl;
>> 113 # else
>> 114 G4cout << G4endl << G4endl << G4endl
>> 115 << " User Environment " << G4endl
>> 116 << " G4ANALYSIS_USE environment variable not set, NO ANALYSIS "
>> 117 << G4endl;
>> 118 #endif
>> 119
>> 120 // Initialize the interactive session
>> 121 G4UIsession* session = 0;
104 if (argc == 1) // Define UI session for in 122 if (argc == 1) // Define UI session for interactive mode.
105 { 123 {
106 auto ui = new G4UIExecutive(argc, argv); << 124 session = new G4UIterminal();
>> 125 }
>> 126
>> 127 //Initialize G4 kernel
>> 128 pRunManager -> Initialize();
>> 129
>> 130 // get the pointer to the User Interface manager
>> 131 G4UImanager* UI = G4UImanager::GetUIpointer();
>> 132 if (session) // Define UI session for interactive mode.
>> 133 {
107 G4cout << " UI session starts ..." << G4 134 G4cout << " UI session starts ..." << G4endl;
108 UImanager -> ApplyCommand("/control/exec << 135 UI -> ApplyCommand("/control/execute VisualisationMacro.mac");
109 ui -> SessionStart(); << 136 session -> SessionStart();
110 delete ui; << 137 delete session;
111 } 138 }
112 else // Batch mode 139 else // Batch mode
113 { << 140 {
114 G4String command = "/control/execute "; 141 G4String command = "/control/execute ";
115 G4String fileName = argv[1]; 142 G4String fileName = argv[1];
116 UImanager -> ApplyCommand(command+fileNa << 143 UI -> ApplyCommand(command+fileName);
117 } << 144 }
118 145
>> 146 // Close the output file containing histograms and n-tuple with the
>> 147 // results of the simulation.
>> 148 // Finish the analysis.
>> 149 #ifdef G4ANALYSIS_USE
>> 150 BrachyAnalysisManager* analysis = BrachyAnalysisManager::getInstance();
>> 151 analysis -> finish();
>> 152 #endif
>> 153
119 // Job termination 154 // Job termination
120 // Close the root file << 155 #ifdef G4VIS_USE
121 <<
122 delete visManager; 156 delete visManager;
>> 157 #endif
>> 158
123 delete pRunManager; 159 delete pRunManager;
124 160
125 return 0; 161 return 0;
126 } 162 }
127 163