Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 /// \file Hadr07.cc 27 /// \brief Main program of the hadronic/Hadr07 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 #include "ActionInitialization.hh" 34 #include "DetectorConstruction.hh" 35 #include "PhysicsList.hh" 36 37 #include "G4ParticleHPManager.hh" 38 #include "G4RunManagerFactory.hh" 39 #include "G4SteppingVerbose.hh" 40 #include "G4Types.hh" 41 #include "G4UIExecutive.hh" 42 #include "G4UImanager.hh" 43 #include "G4VisExecutive.hh" 44 #include "Randomize.hh" 45 46 //....oooOO0OOooo........oooOO0OOooo........oo 47 48 int main(int argc, char** argv) 49 { 50 // detect interactive mode (if no arguments) 51 G4UIExecutive* ui = nullptr; 52 if (argc == 1) ui = new G4UIExecutive(argc, 53 54 // choose the Random engine 55 G4Random::setTheEngine(new CLHEP::RanecuEngi 56 57 // use SteppingVerbose with units 58 G4int precision = 4; 59 G4SteppingVerbose::UseBestUnit(precision); 60 61 // construct the run manager 62 auto runManager = G4RunManagerFactory::Creat 63 if (argc == 3) { 64 G4int nThreads = G4UIcommand::ConvertToInt 65 runManager->SetNumberOfThreads(nThreads); 66 } 67 68 // set mandatory initialization classes 69 DetectorConstruction* det = new DetectorCons 70 runManager->SetUserInitialization(det); 71 72 PhysicsList* phys = new PhysicsList; 73 runManager->SetUserInitialization(phys); 74 runManager->SetUserInitialization(new Action 75 76 // Replaced HP environmental variables with 77 G4ParticleHPManager::GetInstance()->SetSkipM 78 G4ParticleHPManager::GetInstance()->SetDoNot 79 G4ParticleHPManager::GetInstance()->SetUseOn 80 G4ParticleHPManager::GetInstance()->SetNegle 81 G4ParticleHPManager::GetInstance()->SetProdu 82 G4ParticleHPManager::GetInstance()->SetUseWe 83 G4ParticleHPManager::GetInstance()->SetUseNR 84 85 // initialize visualization 86 G4VisManager* visManager = nullptr; 87 88 // get the pointer to the User Interface man 89 G4UImanager* UImanager = G4UImanager::GetUIp 90 91 if (ui) { 92 // interactive mode 93 visManager = new G4VisExecutive; 94 visManager->Initialize(); 95 ui->SessionStart(); 96 delete ui; 97 } 98 else { 99 // batch mode 100 G4String command = "/control/execute "; 101 G4String fileName = argv[1]; 102 UImanager->ApplyCommand(command + fileName 103 } 104 105 // job termination 106 delete visManager; 107 delete runManager; 108 } 109 110 //....oooOO0OOooo........oooOO0OOooo........oo 111