Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // 27 /// \file medical/GammaTherapy/GammaTherapy.cc << 28 /// \brief Main program of the medical/GammaTh << 29 // << 30 // << 31 // ------------------------------------------- 24 // ------------------------------------------------------------- 32 // GEANT4 ibrem 25 // GEANT4 ibrem 33 // 26 // 34 // History: based on object model of 27 // History: based on object model of 35 // 2nd December 1995, G.Cosmo 28 // 2nd December 1995, G.Cosmo 36 // ---------- ibrem --------------------- 29 // ---------- ibrem ---------------------------- 37 // 30 // 38 // Modified: 05.04.01 Vladimir Ivanchenko new 31 // Modified: 05.04.01 Vladimir Ivanchenko new design of ibrem 39 // 32 // 40 // ------------------------------------------- 33 // ------------------------------------------------------------- 41 34 42 //....oooOO0OOooo........oooOO0OOooo........oo 35 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 43 //....oooOO0OOooo........oooOO0OOooo........oo 36 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 44 37 45 #include "ActionInitialization.hh" << 38 #include "G4RunManager.hh" 46 #include "DetectorConstruction.hh" << 47 #include "PhysicsList.hh" << 48 << 49 #include "G4RunManagerFactory.hh" << 50 #include "G4Types.hh" << 51 #include "G4UIExecutive.hh" << 52 #include "G4UImanager.hh" 39 #include "G4UImanager.hh" >> 40 #include "G4UIterminal.hh" >> 41 #include "G4UItcsh.hh" >> 42 #include "Randomize.hh" >> 43 53 #include "G4VisExecutive.hh" 44 #include "G4VisExecutive.hh" 54 45 >> 46 #include "DetectorConstruction.hh" >> 47 #include "PhysicsList.hh" >> 48 #include "PrimaryGeneratorAction.hh" >> 49 #include "EventAction.hh" >> 50 #include "TrackingAction.hh" >> 51 #include "RunAction.hh" >> 52 #include "Histo.hh" >> 53 >> 54 #include <string> >> 55 #include <fstream> >> 56 #include <iomanip> >> 57 55 //....oooOO0OOooo........oooOO0OOooo........oo 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 56 //....oooOO0OOooo........oooOO0OOooo........oo 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 57 60 58 int main(int argc, char** argv) << 61 int main(int argc,char** argv) { 59 { << 62 60 // detect interactive mode (if no arguments) << 63 G4int verbose = 1; 61 G4UIExecutive* ui = nullptr; << 64 //choose the Random engine 62 if (argc == 1) ui = new G4UIExecutive(argc, << 65 CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine); >> 66 >> 67 // G4String fName = argv[1]; >> 68 // (Histo::GetInstance)->LoadParameters(fName); 63 69 64 // Construct the default run manager 70 // Construct the default run manager 65 auto* runManager = G4RunManagerFactory::Crea << 71 G4RunManager * runManager = new G4RunManager(); 66 G4int nThreads = std::min(G4Threading::G4Get << 67 if (argc == 3) nThreads = G4UIcommand::Conve << 68 runManager->SetNumberOfThreads(nThreads); << 69 G4cout << "===== GammaTherapy is started wit << 70 << " threads =====" << G4endl; << 71 72 72 // set mandatory initialization classes 73 // set mandatory initialization classes 73 DetectorConstruction* det = new DetectorCons 74 DetectorConstruction* det = new DetectorConstruction(); 74 runManager->SetUserInitialization(det); 75 runManager->SetUserInitialization(det); 75 76 76 runManager->SetUserInitialization(new Physic 77 runManager->SetUserInitialization(new PhysicsList()); 77 78 78 // set user action classes << 79 // visualization manager 79 runManager->SetUserInitialization(new Action << 80 G4VisManager* visManager = new G4VisExecutive(); >> 81 visManager->Initialize(); 80 82 81 // initialize visualization << 83 // set user action classes 82 G4VisManager* visManager = nullptr; << 84 runManager->SetUserAction(new PrimaryGeneratorAction(det)); >> 85 runManager->SetUserAction(new RunAction()); >> 86 runManager->SetUserAction(new EventAction()); >> 87 runManager->SetUserAction(new TrackingAction()); 83 88 84 // get the pointer to the User Interface man 89 // get the pointer to the User Interface manager 85 G4UImanager* UImanager = G4UImanager::GetUIp << 90 G4UImanager* UI = G4UImanager::GetUIpointer(); >> 91 if(1 < verbose) UI->ListCommands("/testem/"); 86 92 87 if (ui) { << 93 if (argc==1) // Define UI terminal for interactive mode 88 // interactive mode << 94 { 89 visManager = new G4VisExecutive; << 95 G4UIsession * session = new G4UIterminal; 90 visManager->Initialize(); << 96 #ifdef G4UI_USE_TCSH 91 UImanager->ApplyCommand("/control/execute << 97 session = new G4UIterminal(new G4UItcsh); 92 ui->SessionStart(); << 98 #else 93 delete ui; << 99 session = new G4UIterminal(); 94 } << 100 #endif 95 else { << 101 session->SessionStart(); 96 // batch mode << 102 delete session; 97 G4String command = "/control/execute "; << 103 } 98 G4String fileName = argv[1]; << 104 else if (argc>1) // Batch mode with 1 or more files 99 UImanager->ApplyCommand(command + fileName << 105 { 100 } << 106 if(verbose >0) G4cout << "UI interface is started" << G4endl; >> 107 G4String command = "/control/execute "; >> 108 G4String fileName = argv[1]; >> 109 UI->ApplyCommand(command+fileName); >> 110 } 101 111 102 // job termination 112 // job termination >> 113 103 delete visManager; 114 delete visManager; 104 delete runManager; 115 delete runManager; >> 116 >> 117 return 0; 105 } 118 } 106 119 107 //....oooOO0OOooo........oooOO0OOooo........oo << 120 108 121