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 // 27 /// \file molecular.cc 28 /// \brief Molecular level simulation of DNA 29 30 #include "ActionInitialization.hh" 31 #include "DetectorConstruction.hh" 32 #include "PhysicsList.hh" 33 34 #include "G4DNAChemistryManager.hh" 35 #include "G4MoleculeGun.hh" 36 #include "G4RunManagerFactory.hh" 37 #include "G4UIExecutive.hh" 38 #include "G4UImanager.hh" 39 #include "G4VisExecutive.hh" 40 41 #include <ctime> 42 43 //....oooOO0OOooo........oooOO0OOooo........oo 44 45 namespace 46 { 47 void PrintUsage() 48 { 49 G4cout << " Usage: " << G4endl; 50 G4cout << " molecular [-m macro ] [-t nThrea 51 G4cout << " -p is the G4DNA Physics List o 52 << " G4EmDNAPhysics" << G4endl; 53 G4cout << " note: -t option is available o 54 } 55 } // namespace 56 57 int main(int argc, char** argv) 58 { 59 if (argc > 7) { 60 PrintUsage(); 61 return 1; 62 } 63 64 G4String macro; 65 G4int phys_option = 2; 66 67 G4int nThreads = 2; 68 for (G4int ii = 1; ii < argc; ii = ii + 2) { 69 if (G4String(argv[ii]) == "-m") { 70 macro = argv[ii + 1]; 71 } 72 else if (G4String(argv[ii]) == "-p") { 73 phys_option = G4UIcommand::ConvertToInt( 74 } 75 else if (G4String(argv[ii]) == "-t") { 76 nThreads = G4UIcommand::ConvertToInt(arg 77 } 78 else { 79 PrintUsage(); 80 return 1; 81 } 82 } 83 84 G4UIExecutive* ui = nullptr; 85 if (!macro.size()) { 86 ui = new G4UIExecutive(argc, argv); 87 } 88 89 /* 90 time_t timeStart; 91 time(&timeStart); 92 long seed = timeStart; 93 G4Random::setTheSeed(seed); 94 */ 95 96 auto* runManager = G4RunManagerFactory::Crea 97 if (nThreads > 0) { 98 runManager->SetNumberOfThreads(nThreads); 99 } 100 101 runManager->SetUserInitialization(new Detect 102 G4VModularPhysicsList* physicsList = new Phy 103 runManager->SetUserInitialization(physicsLis 104 runManager->SetUserInitialization(new Action 105 G4DNAChemistryManager::Instance()->Initializ 106 107 G4VisExecutive* visManager = nullptr; 108 109 // Get the pointer to the User Interface man 110 G4UImanager* UImanager = G4UImanager::GetUIp 111 if (!macro.empty()) { 112 // batch mode 113 G4String command = "/control/execute "; 114 UImanager->ApplyCommand(command + macro); 115 } 116 else { 117 visManager = new G4VisExecutive; 118 visManager->Initialize(); 119 UImanager->ApplyCommand("/control/execute 120 ui->SessionStart(); 121 delete ui; 122 delete visManager; 123 } 124 delete runManager; 125 return 0; 126 } 127 128 //....oooOO0OOooo........oooOO0OOooo........oo 129