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