LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
exampleGB05.cc File Reference

Main program of example GB05. More...

#include "G4Types.hh"
#include "G4RunManager.hh"
#include "GB05ActionInitialization.hh"
#include "G4UImanager.hh"
#include "GB05DetectorConstruction.hh"
#include "GB05PrimaryGeneratorAction.hh"
#include "FTFP_BERT.hh"
#include "G4GenericBiasingPhysics.hh"
#include "G4VisExecutive.hh"
#include "G4UIExecutive.hh"

Go to the source code of this file.

Functions

int main (int argc, char **argv)
 

Detailed Description

Main program of example GB05.

Definition in file exampleGB05.cc.

Function Documentation

int main ( int  argc,
char **  argv 
)

Definition at line 66 of file exampleGB05.cc.

References PrintUsage().

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
115  GB05DetectorConstruction* detector = new GB05DetectorConstruction();
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:
141  runManager->SetUserInitialization(new GB05ActionInitialization);
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 }
void PrintUsage()