ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
exampleGB04.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file exampleGB04.cc
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
//
29
//
30
31
#include "
G4Types.hh
"
32
33
#ifdef G4MULTITHREADED
34
#include "
G4MTRunManager.hh
"
35
#else
36
#include "
G4RunManager.hh
"
37
#endif
38
#include "
GB04ActionInitialization.hh
"
39
40
#include "
G4UImanager.hh
"
41
42
#include "
GB04DetectorConstruction.hh
"
43
#include "
GB04PrimaryGeneratorAction.hh
"
44
45
#include "
FTFP_BERT.hh
"
46
#include "
G4GenericBiasingPhysics.hh
"
47
48
#include "
G4VisExecutive.hh
"
49
#include "
G4UIExecutive.hh
"
50
51
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52
53
namespace
{
54
void
PrintUsage() {
55
G4cerr
<<
" Usage: "
<<
G4endl
;
56
G4cerr
<<
" ./exampleGB04 [-m macro ] "
57
<<
" [-b biasing {'on','off'}]"
58
<<
"\n or\n ./exampleGB04 [macro.mac]"
59
<<
G4endl
;
60
}
61
}
62
63
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
64
65
int
main
(
int
argc,
char
** argv)
66
{
67
// Evaluate arguments
68
//
69
if
( argc > 5 ) {
70
PrintUsage();
71
return
1;
72
}
73
74
G4String
macro(
""
);
75
G4String
onOffBiasing(
""
);
76
if
( argc == 2 ) macro = argv[1];
77
else
78
{
79
for
(
G4int
i=1; i<argc; i=i+2 )
80
{
81
if
(
G4String
(argv[i]) ==
"-m"
) macro = argv[i+1];
82
else
if
(
G4String
(argv[i]) ==
"-b"
) onOffBiasing = argv[i+1];
83
else
84
{
85
PrintUsage();
86
return
1;
87
}
88
}
89
}
90
91
if
( onOffBiasing ==
""
) onOffBiasing =
"on"
;
92
93
// Instantiate G4UIExecutive if interactive mode
94
G4UIExecutive
* ui =
nullptr
;
95
if
( macro ==
""
) {
96
ui =
new
G4UIExecutive
(argc, argv);
97
}
98
99
// -- Construct the run manager : MT or sequential one
100
#ifdef G4MULTITHREADED
101
G4MTRunManager
* runManager =
new
G4MTRunManager
;
102
G4cout
<<
" ********** Run Manager constructed in MT mode ************ "
103
<<
G4endl
;
104
// -- Choose 4 threads:
105
runManager->
SetNumberOfThreads
(4);
106
#else
107
G4RunManager
* runManager =
new
G4RunManager
;
108
G4cout
<<
" ********** Run Manager constructed in sequential mode ************ "
109
<<
G4endl
;
110
#endif
111
112
// -- Set mandatory initialization classes
113
GB04DetectorConstruction
* detector =
new
GB04DetectorConstruction
();
114
runManager->
SetUserInitialization
(detector);
115
// -- Select a physics list
116
FTFP_BERT
* physicsList =
new
FTFP_BERT
;
117
// -- And augment it with biasing facilities:
118
G4GenericBiasingPhysics
* biasingPhysics =
new
G4GenericBiasingPhysics
();
119
if
( onOffBiasing ==
"on"
)
120
{
121
// -- Create list of physics processes to be biased: only brem. in this case:
122
std::vector< G4String > processToBias;
123
processToBias.push_back(
"eBrem"
);
124
// -- Pass the list to the G4GenericBiasingPhysics, which will wrap the eBrem
125
// -- process of e- and e+ to activate the biasing of it:
126
biasingPhysics->
PhysicsBias
(
"e-"
, processToBias);
127
biasingPhysics->
PhysicsBias
(
"e+"
, processToBias);
128
physicsList->
RegisterPhysics
(biasingPhysics);
129
G4cout
<<
" ********************************************************* "
130
<<
G4endl
;
131
G4cout
<<
" ********** processes are wrapped for biasing ************ "
132
<<
G4endl
;
133
G4cout
<<
" ********************************************************* "
134
<<
G4endl
;
135
}
136
else
137
{
138
G4cout
<<
" ************************************************* "
<<
G4endl
;
139
G4cout
<<
" ********** processes are not wrapped ************ "
<<
G4endl
;
140
G4cout
<<
" ************************************************* "
<<
G4endl
;
141
}
142
runManager->
SetUserInitialization
(physicsList);
143
// -- Action initialization:
144
runManager->
SetUserInitialization
(
new
GB04ActionInitialization
);
145
146
// Initialize G4 kernel
147
runManager->
Initialize
();
148
149
// Initialize visualization
150
G4VisManager
* visManager =
new
G4VisExecutive
;
151
// G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
152
visManager->
Initialize
();
153
154
// Get the pointer to the User Interface manager
155
G4UImanager
* UImanager =
G4UImanager::GetUIpointer
();
156
157
if
( !ui )
// batch mode
158
{
159
G4String
command =
"/control/execute "
;
160
UImanager->
ApplyCommand
(command+macro);
161
}
162
else
163
{
// interactive mode : define UI session
164
// UImanager->ApplyCommand("/control/execute vis.mac");
165
// if (ui->IsGUI())
166
// UImanager->ApplyCommand("/control/execute gui.mac");
167
ui->
SessionStart
();
168
delete
ui;
169
}
170
171
delete
visManager;
172
delete
runManager;
173
174
175
return
0;
176
}
177
178
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
geant4
tree
geant4-10.6-release
examples
extended
biasing
GB04
exampleGB04.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:02
using
1.8.2 with
ECCE GitHub integration