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 // This example is provided by the Geant4-DNA 27 // Any report or published results obtained us 28 // shall cite the following Geant4-DNA collabo 29 // Med. Phys. 37 (2010) 4692-4708 30 // J. Comput. Phys. 274 (2014) 841-882 31 // The Geant4-DNA web site is available at htt 32 // 33 // 34 /// \file chem1.cc 35 /// \brief Chem1 example 36 37 #include "ActionInitialization.hh" 38 #include "CommandLineParser.hh" 39 #include "DetectorConstruction.hh" 40 #include "PhysicsList.hh" 41 42 #include "G4DNAChemistryManager.hh" 43 #include "G4RunManagerFactory.hh" 44 #include "G4UIExecutive.hh" 45 #include "G4UImanager.hh" 46 #include "G4VisExecutive.hh" 47 48 //....oooOO0OOooo........oooOO0OOooo........oo 49 50 /* 51 * WARNING : Geant4 was initially not intended 52 * This code is delivered as a prototype 53 * We will be happy to hear from you, do not h 54 * communicate on the difficulties you may enc 55 * The user interface may change in the next r 56 * code has started 57 */ 58 59 using namespace G4DNAPARSER; 60 CommandLineParser* parser(0); 61 62 void Parse(int& argc, char** argv); 63 64 int main(int argc, char** argv) 65 { 66 ////////// 67 // Parse options given in commandLine 68 // 69 Parse(argc, argv); 70 71 ////////// 72 // Construct the run manager according to wh 73 // 74 Command* commandLine(nullptr); 75 76 auto* runManager = G4RunManagerFactory::Crea 77 if ((commandLine = parser->GetCommandIfActiv 78 int nThreads = 2; 79 if (commandLine->GetOption() == "NMAX") { 80 nThreads = G4Threading::G4GetNumberOfCor 81 } 82 else { 83 nThreads = G4UIcommand::ConvertToInt(com 84 } 85 runManager->SetNumberOfThreads(nThreads); 86 G4cout << "===== Chem1 is started with " < 87 << " threads =====" << G4endl; 88 } 89 90 ////////// 91 // Set mandatory user initialization classes 92 // 93 DetectorConstruction* detector = new Detecto 94 runManager->SetUserInitialization(new Physic 95 runManager->SetUserInitialization(detector); 96 runManager->SetUserInitialization(new Action 97 98 // Initialize visualization 99 G4VisManager* visManager = nullptr; 100 101 // Get the pointer to the User Interface man 102 G4UImanager* UImanager = G4UImanager::GetUIp 103 G4UIExecutive* ui(nullptr); 104 105 // interactive mode : define UI session 106 if ((commandLine = parser->GetCommandIfActiv 107 visManager = new G4VisExecutive; 108 visManager->Initialize(); 109 110 ui = new G4UIExecutive(argc, argv, command 111 112 if (parser->GetCommandIfActive("-novis") = 113 // visualization is used by default 114 if ((commandLine = parser->GetCommandIfA 115 // select a visualization driver if ne 116 UImanager->ApplyCommand(G4String("/vis 117 } 118 else { 119 // by default OGL is used 120 UImanager->ApplyCommand("/vis/open OGL 121 } 122 UImanager->ApplyCommand("/control/execut 123 } 124 125 if (ui->IsGUI()) { 126 UImanager->ApplyCommand("/control/execut 127 } 128 } 129 else if ((commandLine = parser->GetCommandIf 130 // to be use visualization file (= store t 131 // an external file: 132 // ASCIITree ; DAWNFILE ; HepRepFile ; VR 133 visManager = new G4VisExecutive; 134 visManager->Initialize(); 135 136 ui = new G4UIExecutive(argc, argv, command 137 UImanager->ApplyCommand(G4String("/vis/ope 138 UImanager->ApplyCommand("/control/execute 139 } 140 141 if ((commandLine = parser->GetCommandIfActiv 142 G4String command = "/control/execute "; 143 UImanager->ApplyCommand(command + commandL 144 } 145 else { 146 UImanager->ApplyCommand("/control/execute 147 } 148 149 if (ui) { 150 ui->SessionStart(); 151 delete ui; 152 } 153 154 // Job termination 155 // Free the store: user actions, physics_lis 156 // owned and deleted by the run manager, so 157 // in the main() program ! 158 159 delete visManager; 160 delete runManager; 161 162 CommandLineParser::DeleteInstance(); 163 164 return 0; 165 } 166 167 void Parse(int& argc, char** argv) 168 { 169 ////////// 170 // Parse options given in commandLine 171 // 172 parser = CommandLineParser::GetParser(); 173 174 parser->AddCommand("-gui", Command::OptionNo 175 "Select geant4 UI or just 176 177 parser->AddCommand("-mac", Command::WithOpti 178 179 // You cann your own command, as for instanc 180 // parser->AddCommand("-seed", 181 // Command::WithOption, 182 // "Give a seed value in 183 // it is then up to you to manage this optio 184 185 parser->AddCommand("-mt", Command::WithOptio 186 "Launch in MT mode (event 187 " NOT RECOMMANDED WITH CH 188 "2"); 189 190 parser->AddCommand("-chemOFF", Command::With 191 192 parser->AddCommand("-vis", Command::WithOpti 193 "OGL 600x600-0+0"); 194 195 parser->AddCommand("-novis", Command::Withou 196 197 ////////// 198 // If -h or --help is given in option : prin 199 // 200 if (parser->Parse(argc, argv) != 0) // help 201 { 202 // if you are using ROOT, create a TApplic 203 // to print the help from ROOT as well 204 CommandLineParser::DeleteInstance(); 205 std::exit(0); 206 } 207 208 /////////// 209 // Kill application if wrong argument in com 210 // 211 if (parser->CheckIfNotHandledOptionsExists(a 212 // if you are using ROOT, you should initi 213 // before this condition 214 abort(); 215 } 216 } 217