Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 /// \file parallel/ThreadsafeScorers/ts_scorers.cc 27 /// \brief Main of the ThreadsafeScorers example 28 // 29 // 30 // 31 // 32 /// ts_scorers example shows how to use global scorers. The benefit of using 33 /// global scorers in memory-savings for problems with very large amounts 34 /// of scoring volumes. Additionally, the global scorers are more precise 35 /// w.r.t. the serial solution because of the lack of compounding 36 /// round-off error from multiple threads 37 /// 38 /// In this example, the global scorers are implemented as static member 39 /// variables in TSRun because TSRun is thread-local. The G4atomic 40 /// class is the core of the thread-safe scorers and can be uses 41 // 42 // 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 #include "G4RunManagerFactory.hh" 47 #include "G4Threading.hh" 48 #include "G4Types.hh" 49 #include "Randomize.hh" 50 51 // User Defined Classes 52 #include "TSActionInitialization.hh" 53 #include "TSDetectorConstruction.hh" 54 #include "TSPhysicsList.hh" 55 56 #include "G4Step.hh" 57 #include "G4Track.hh" 58 #include "G4UIExecutive.hh" 59 #include "G4UImanager.hh" 60 #include "G4VisExecutive.hh" 61 62 // for std::system(const char*) 63 #include <cstdlib> 64 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 66 67 void message(G4RunManager* runmanager) 68 { 69 G4MTRunManager* man = dynamic_cast<G4MTRunManager*>(runmanager); 70 if (man) { 71 man->SetNumberOfThreads(G4Threading::G4GetNumberOfCores()); 72 G4cout << "\n\n\t--> Running in multithreaded mode with " << man->GetNumberOfThreads() 73 << " threads\n\n" 74 << G4endl; 75 } 76 else { 77 G4cout << "\n\n\t--> Running in serial mode\n\n" << G4endl; 78 } 79 } 80 81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 82 83 int main(int argc, char** argv) 84 { 85 G4String macro; 86 if (argc > 1) macro = argv[argc - 1]; 87 88 // Detect interactive mode (if no arguments) and define UI session 89 // 90 G4UIExecutive* ui = 0; 91 if (macro.empty()) ui = new G4UIExecutive(argc, argv); 92 93 // Set the random seed 94 CLHEP::HepRandom::setTheSeed(1245214UL); 95 96 G4RunManager* runmanager = G4RunManagerFactory::CreateRunManager(G4RunManagerType::Tasking); 97 98 message(runmanager); 99 100 runmanager->SetUserInitialization(new TSDetectorConstruction); 101 102 runmanager->SetUserInitialization(new TSPhysicsList); 103 104 runmanager->SetUserInitialization(new TSActionInitialization); 105 106 runmanager->Initialize(); 107 108 // Initialize visualization 109 // 110 G4VisManager* visManager = new G4VisExecutive; 111 // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance. 112 // G4VisManager* visManager = new G4VisExecutive("Quiet"); 113 visManager->Initialize(); 114 115 // Get the pointer to the User Interface manager 116 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 117 118 // Process macro or start UI session 119 // 120 if (!ui) { 121 // batch mode 122 G4String command = "/control/execute "; 123 UImanager->ApplyCommand(command + macro); 124 } 125 else { 126 ui->SessionStart(); 127 } 128 129 // Job termination 130 // Free the store: user actions, physics_list and detector_description are 131 // owned and deleted by the run manager, so they should not be deleted 132 // in the main() program ! 133 delete visManager; 134 delete runmanager; 135 136 return 0; 137 } 138