From cd7daa4e293869e803217d1623b1ca306c951960 Mon Sep 17 00:00:00 2001 From: Yessica Dominguez <yessicadomin@gmail.com> Date: Tue, 6 Jun 2023 20:56:57 +0000 Subject: [PATCH] Upload New File --- CODE/sim.cc | 46 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 46 insertions(+) create mode 100644 CODE/sim.cc diff --git a/CODE/sim.cc b/CODE/sim.cc new file mode 100644 index 0000000..dae2a67 --- /dev/null +++ b/CODE/sim.cc @@ -0,0 +1,46 @@ +#include <iostream> + +#include "G4RunManager.hh" +#include "G4MTRunManager.hh" +#include "G4UImanager.hh" +#include "QGSP_BERT.hh" +#include "G4UImanager.hh" + +#include "construction.hh" +#include "physics.hh" +#include "action.hh" + +int main(int argc, char** argv) +{ + + G4UImanager* UImanager = G4UImanager::GetUIpointer(); + + #ifdef G4MULTITHREADED + G4MTRunManager* runManager = new G4MTRunManager; + #else + G4RunManager* runManager = new G4RunManager; + #endif + runManager->SetUserInitialization(new MyDetectorConstruction()); + runManager->SetUserInitialization(new MyPhysicsList()); + runManager->SetUserInitialization(new MyActionInitialization()); + + G4VModularPhysicsList* physics = new QGSP_BERT(); + physics->RegisterPhysics(new G4DecayPhysics()); + runManager->SetUserInitialization(physics); + + runManager->Initialize(); + + if (argc > 1) + { + G4String command = "/control/execute "; + G4String fileName = argv[1]; + G4UImanager::GetUIpointer()->ApplyCommand(command + fileName); + } + + G4String command = "/control/execute input.in"; + UImanager->ApplyCommand(command); + delete runManager; + + return 0; +} + -- GitLab