81int main(
int argc,
char **argv)
84 if (argc>1) mode = atoi(argv[1]);
86 G4cout <<
" mode not used yet, refer to B01 to see WeightWindow technique "
90 G4int numberOfEvents = 100;
91 G4long myseed = 345354;
93 auto* runManager = G4RunManagerFactory::CreateRunManager();
94 runManager->SetNumberOfThreads(2);
96 G4Random::setTheSeed(myseed);
100 runManager->SetUserInitialization(detector);
102 G4String parallelName(
"ParallelBiasingWorld");
105 detector->RegisterParallelWorld(pdet);
109 pgs.SetParallel(
true);
115 runManager->SetUserInitialization(physicsList);
120 runManager->SetUserInitialization(actions);
122 runManager->Initialize();
127 G4UImanager* UImanager = G4UImanager::GetUIpointer();
128 G4String command1 =
"/control/cout/setCoutFile threadOut";
129 UImanager->ApplyCommand(command1);
130 G4String command2 =
"/run/beamOn " +
131 G4UIcommand::ConvertToString(numberOfEvents);
132 UImanager->ApplyCommand(command2);
136 G4GeometryManager::GetInstance()->OpenGeometry();