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" )
109 std::vector< G4String > processToBias;
110 processToBias.push_back(
"eBrem");
113 biasingPhysics->PhysicsBias(
"e-", processToBias);
114 biasingPhysics->PhysicsBias(
"e+", processToBias);
115 physicsList->RegisterPhysics(biasingPhysics);
116 G4cout <<
" ********************************************************* "
118 G4cout <<
" ********** processes are wrapped for biasing ************ "
120 G4cout <<
" ********************************************************* "
125 G4cout <<
" ************************************************* " << G4endl;
126 G4cout <<
" ********** processes are not wrapped ************ " << G4endl;
127 G4cout <<
" ************************************************* " << G4endl;
129 runManager->SetUserInitialization(physicsList);
134 runManager->Initialize();
139 visManager->Initialize();
142 G4UImanager* UImanager = G4UImanager::GetUIpointer();
146 G4String command =
"/control/execute ";
147 UImanager->ApplyCommand(command+macro);