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 // Phys. Med. Biol. 63(10) (2018) 105014-12pp 32 // The Geant4-DNA web site is available at htt 33 // 34 // 35 /// \file chem5.cc 36 /// \brief Chem5 example 37 38 #include "ActionInitialization.hh" 39 #include "CommandLineParser.hh" 40 #include "DetectorConstruction.hh" 41 #include "PhysicsList.hh" 42 43 #include "G4DNAChemistryManager.hh" 44 #include "G4RunManagerFactory.hh" 45 #include "G4UIExecutive.hh" 46 #include "G4UImanager.hh" 47 #include "G4VisExecutive.hh" 48 49 /* 50 * WARNING : Geant4 was initially not intended 51 * This code is delivered as a prototype 52 * We will be happy to hear from you, do not h 53 * communicate on the difficulties you may enc 54 * The user interface may change in the next r 55 * code has started 56 */ 57 58 //....oooOO0OOooo........oooOO0OOooo........oo 59 60 using namespace G4DNAPARSER; 61 CommandLineParser* parser(0); 62 long seed = 0; 63 64 unsigned int noise(); 65 void SetSeed(); 66 67 void Parse(int& argc, char** argv); 68 69 //....oooOO0OOooo........oooOO0OOooo........oo 70 71 int main(int argc, char** argv) 72 { 73 // Parse options given in commandLine 74 Parse(argc, argv); 75 Command* commandLine(0); 76 SetSeed(); 77 78 // Construct the run manager according to wh 79 // 80 auto* runManager = G4RunManagerFactory::Crea 81 82 if ((commandLine = parser->GetCommandIfActiv 83 int nThreads = 2; 84 if (commandLine->GetOption() == "NMAX") { 85 nThreads = G4Threading::G4GetNumberOfCor 86 } 87 else { 88 nThreads = G4UIcommand::ConvertToInt(com 89 } 90 91 runManager->SetNumberOfThreads(nThreads); 92 } 93 94 G4cout << "********************************* 95 << "******\n===== Chem5 is started wi 96 << G4Threading::G4GetNumberOfCores() 97 << " available threads =====\n\n***** 98 << "*******************************" 99 100 // Set mandatory initialization classes 101 runManager->SetUserInitialization(new Physic 102 runManager->SetUserInitialization(new Detect 103 runManager->SetUserInitialization(new Action 104 105 // Initialize G4 kernel 106 runManager->Initialize(); 107 108 // Initialize visualization 109 G4VisManager* visManager = new G4VisExecutiv 110 // G4VisExecutive can take a verbosity argum 111 // G4VisManager* visManager = new G4VisExecu 112 visManager->Initialize(); 113 114 // Get the pointer to the User Interface man 115 G4UImanager* UImanager = G4UImanager::GetUIp 116 G4UIExecutive* ui(0); 117 118 // interactive mode : define UI session 119 if ((commandLine = parser->GetCommandIfActiv 120 ui = new G4UIExecutive(argc, argv, command 121 122 if (ui->IsGUI()) UImanager->ApplyCommand(" 123 124 if (parser->GetCommandIfActive("-novis") = 125 // visualization is used by default 126 if ((commandLine = parser->GetCommandIfA 127 // select a visualization driver if ne 128 UImanager->ApplyCommand(G4String("/vis 129 } 130 else { 131 // by default OGL is used 132 UImanager->ApplyCommand("/vis/open OGL 133 } 134 UImanager->ApplyCommand("/control/execut 135 } 136 } 137 else { 138 // to be use visualization file (= store t 139 // an external file: 140 // ASCIITree ; DAWNFILE ; HepRepFile ; VR 141 if ((commandLine = parser->GetCommandIfAct 142 UImanager->ApplyCommand(G4String("/vis/o 143 UImanager->ApplyCommand("/control/execut 144 } 145 } 146 147 if ((commandLine = parser->GetCommandIfActiv 148 G4String command = "/control/execute "; 149 UImanager->ApplyCommand(command + commandL 150 } 151 else { 152 UImanager->ApplyCommand("/control/execute 153 } 154 155 if ((commandLine = parser->GetCommandIfActiv 156 ui->SessionStart(); 157 delete ui; 158 } 159 160 // Job termination 161 // Free the store: user actions, physics_lis 162 // owned and deleted by the run manager, so 163 // in the main() program ! 164 delete visManager; 165 delete runManager; 166 CommandLineParser::DeleteInstance(); 167 return 0; 168 } 169 170 //....oooOO0OOooo........oooOO0OOooo........oo 171 172 bool IsBracket(char c) 173 { 174 switch (c) { 175 case '[': 176 case ']': 177 return true; 178 default: 179 return false; 180 } 181 } 182 183 //....oooOO0OOooo........oooOO0OOooo........oo 184 185 void SetSeed() 186 { 187 Command* commandLine(0); 188 189 if ((commandLine = parser->GetCommandIfActiv 190 seed = atoi(commandLine->GetOption().c_str 191 } 192 193 if (seed == 0) // If no seed given in argum 194 { 195 long jobID_int = 0; 196 long noice = 0; 197 198 //____________________________________ 199 // In case on cluster 200 if ((commandLine = parser->GetCommandIfAct 201 noice = labs((long)noise()); 202 203 const char* env = std::getenv("PBS_JOBID 204 205 if (env) { 206 G4String buffer(env); 207 G4String jobID_string = buffer.substr( 208 jobID_string.erase(std::remove_if(jobI 209 jobID_string.end()) 210 jobID_int = atoi(jobID_string.c_str()) 211 } 212 else { 213 env = std::getenv("SGE_TASK_ID"); 214 if (env) jobID_int = atoi(env); 215 } 216 } // end cluster 217 218 //____________________________________ 219 seed = ((long)time(NULL)) + jobID_int + no 220 } 221 222 G4cout << "Seed used : " << seed << G4endl; 223 G4Random::setTheEngine(new CLHEP::MixMaxRng( 224 G4Random::setTheSeed(seed); 225 226 // Choose the Random engine 227 // CLHEP::HepRandom::setTheEngine(new CLHEP: 228 } 229 230 //....oooOO0OOooo........oooOO0OOooo........oo 231 232 unsigned int noise() 233 { 234 #if defined(WIN32) || defined(_WIN32) || defin 235 // TODO: MS Win 236 return std::time(0); 237 #else 238 unsigned int random_seed, random_seed_a, ran 239 std::ifstream file("/dev/urandom", std::ios: 240 if (file.is_open()) { 241 char* memblock; 242 int size = sizeof(int); 243 memblock = new char[size]; 244 file.read(memblock, size); 245 file.close(); 246 random_seed_a = *reinterpret_cast<int*>(me 247 delete[] memblock; 248 } // end if 249 else { 250 random_seed_a = 0; 251 } 252 random_seed_b = std::time(0); 253 random_seed = random_seed_a xor random_seed_ 254 return random_seed; 255 #endif 256 } 257 258 //....oooOO0OOooo........oooOO0OOooo........oo 259 260 void Parse(int& argc, char** argv) 261 { 262 ////////// 263 // Parse options given in commandLine 264 // 265 parser = CommandLineParser::GetParser(); 266 267 parser->AddCommand("-gui", Command::OptionNo 268 "Select geant4 UI or just 269 "session", 270 "qt"); 271 272 parser->AddCommand("-mac", Command::WithOpti 273 274 parser->AddCommand("-seed", Command::WithOpt 275 "seed"); 276 277 parser->AddCommand("-mt", Command::WithOptio 278 "Launch in MT mode (event 279 " NOT RECOMMANDED WITH CH 280 "2"); 281 282 parser->AddCommand("-chemOFF", Command::With 283 284 parser->AddCommand("-vis", Command::WithOpti 285 "OGL 600x600-0+0"); 286 287 parser->AddCommand("-novis", Command::Withou 288 289 parser->AddCommand("-cluster", Command::With 290 "Launch the code on a clu 291 292 ////////// 293 // If -h or --help is given in option : prin 294 // 295 if (parser->Parse(argc, argv) != 0) // help 296 { 297 // if you are using ROOT, create a TApplic 298 // to print the help from ROOT as well 299 CommandLineParser::DeleteInstance(); 300 std::exit(0); 301 } 302 303 /////////// 304 // Kill application if wrong argument in com 305 // 306 if (parser->CheckIfNotHandledOptionsExists(a 307 // if you are using ROOT, you should initi 308 // before this condition 309 std::exit(0); 310 } 311 } 312