G4VModularPhysicsList* G4FastSimPhysicsList::getPhysicsList() {
// ownership passed to G4SimSvc which will register it in G4RunManager. To be deleted in ~G4RunManager()
G4VModularPhysicsList* physicsList = m_physicsListTool->getPhysicsList();
// Coupled transportation enables calculation of particle trajectory in envelopes with fast sim models attached
G4PhysicsListHelper::GetPhysicsListHelper()->UseCoupledTransportation();
// Attach Fast Simulation Process (will take over normal transportation if FastSimModel triggered
physicsList->RegisterPhysics(new sim::FastSimPhysics);
return physicsList;
}
G4VModularPhysicsList* GeantFastSimConfig::getPhysicsList() {
// we pass the ownership to Geant4Simulation which will register it in G4RunManager and will be deleted in ~G4RunManager()
G4VModularPhysicsList* physicsList = new FTFP_BERT;
// Coupled transportation enables calculation of particle trajectory in envelopes with fast sim models attached
G4PhysicsListHelper::GetPhysicsListHelper()->UseCoupledTransportation();
// Attach Fast Simulation Process (will take over normal transportation if FastSimModel triggered
physicsList->RegisterPhysics(new FastSimPhysics);
return physicsList;
}
/// Extend physics list from factory:
G4VUserPhysicsList* Geant4PhysicsListActionSequence::extensionList() {
G4VModularPhysicsList* physics = ( m_extends.empty() )
? new EmptyPhysics()
: G4PhysListFactory().GetReferencePhysList(m_extends);
// Register all physics constructors with the physics list
constructPhysics(physics);
// Ensure the particles and processes declared are also invoked.
// Hence: Create a special physics constructor doing so.
// Ownership is transferred to the physics list.
// Do not delete this pointer afterwards....
physics->RegisterPhysics(new ParticlePhysics(this,physics));
return physics;
}
int main(int argc, char **argv)
{
G4String macro = "";
G4bool showAll = false;
for (G4int i = 1; i < argc; i++) {
if (G4String(argv[i]) == "-m") macro = argv[++i];
else if (G4String(argv[i]) == "-a") showAll = true;
else {
PrintUsage();
return 1;
}
}
// Remove?
// Choose the Random engine
CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);
G4int seed = time(0);
G4cout << "\nseed = " << seed << G4endl;
CLHEP::HepRandom::setTheSeed(seed);
G4Random::setTheSeed(seed);
// Construct the default run manager
//
#ifdef G4MULTITHREADED
G4MTRunManager *runManager = new G4MTRunManager();
runManager->SetNumberOfThreads(G4Threading::G4GetNumberOfCores());
//runManager->SetNumberOfThreads(1);
#else
G4RunManager *runManager = new G4RunManager();
#endif
// Set mandatory initialization classes
//
// Detector construction
runManager->SetUserInitialization(new BIDetectorConstruction());
G4cout << "Here" << G4endl;
// Physics list
//G4VModularPhysicsList *physicsList = new FTFP_BERT;
//G4VModularPhysicsList *physicsList = new QGSP_BERT_HP;
G4VModularPhysicsList *physicsList = new Shielding;
//G4VModularPhysicsList *physicsList = new BIPhysicsList;
physicsList->SetVerboseLevel(0);
//physicsList->SetCutValue(1.*um, "proton");
//physicsList->SetCuts();
//physicsList->SetDefaultCutValue(100.*um);
runManager->SetUserInitialization(physicsList);
// Primary generator action and User action intialization
runManager->SetUserInitialization(new BIActionInitialization());
// Initialize G4 kernel
//
runManager->Initialize();
#ifdef G4VIS_USE
// Initialize visualization
G4VisManager *visManager = new G4VisExecutive;
visManager->Initialize();
if (!showAll) { //Show only proton
G4TrajectoryParticleFilter *filterp = new G4TrajectoryParticleFilter;
filterp->Add("proton");
visManager->RegisterModel(filterp);
}
#endif
// Get the pointer to the User Interface manager
G4UImanager *UImanager = G4UImanager::GetUIpointer();
if (macro != "") {
// batch mode
G4String command = "/control/execute ";
UImanager->ApplyCommand(command + macro);
} else {
// interactive mode : define UI session
#ifdef G4UI_USE
G4UIExecutive *ui = new G4UIExecutive(argc, argv);
#ifdef G4VIS_USE
UImanager->ApplyCommand("/control/execute init_vis.mac");
#else
UImanager->ApplyCommand("/control/execute init.mac");
#endif
if (ui->IsGUI()) {
UImanager->ApplyCommand("/control/execute icons.mac");
}
ui->SessionStart();
delete ui;
#endif
}
// Job termination
// Free the store: user actions, physics_list and detector_description are
// owned and deleted by the run manager, so they should not be deleted
// in the main() program !
#ifdef G4VIS_USE
delete visManager;
#endif
delete runManager;
return 0;
}