61int main(
int argc,
char** argv)
72 if ( argc == 2 ) macro = argv[1];
75 for ( G4int i=1; i<argc; i=i+2 )
77 if (
G4String(argv[i]) ==
"-m" ) macro = argv[i+1];
78 else if (
G4String(argv[i]) ==
"-b" ) onOffBiasing = argv[i+1];
87 if ( onOffBiasing ==
"" ) onOffBiasing =
"on";
96 auto* runManager = G4RunManagerFactory::CreateRunManager();
97 runManager->SetNumberOfThreads(4);
101 runManager->SetUserInitialization(detector);
106 if ( onOffBiasing ==
"on" )
108 biasingPhysics->Bias(
"gamma");
109 biasingPhysics->Bias(
"neutron");
110 physicsList->RegisterPhysics(biasingPhysics);
111 G4cout <<
" ********************************************************* " << G4endl;
112 G4cout <<
" ********** processes are wrapped for biasing ************ " << G4endl;
113 G4cout <<
" ********************************************************* " << G4endl;
117 G4cout <<
" ************************************************* " << G4endl;
118 G4cout <<
" ********** processes are not wrapped ************ " << G4endl;
119 G4cout <<
" ************************************************* " << G4endl;
121 runManager->SetUserInitialization(physicsList);
126 runManager->Initialize();
131 visManager->Initialize();
134 G4UImanager* UImanager = G4UImanager::GetUIpointer();
138 G4String command =
"/control/execute ";
139 UImanager->ApplyCommand(command+macro);