| 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 #include "ActionInitialization.hh"
27 #include "DetectorConstruction.hh"
28 #include "PhysicsList.hh"
29
30 #include "G4AnalysisManager.hh"
31 #include "G4RunManagerFactory.hh"
32 #include "G4UIExecutive.hh"
33 #include "G4UImanager.hh"
34 #include "G4VisExecutive.hh"
35
36 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
37
38 int main(int argc, char** argv)
39 {
40 G4int seed(0);
41 // Detect interactive mode (if no arguments) and define UI session
42 G4UIExecutive* ui = nullptr;
43 if (argc == 1) {
44 ui = new G4UIExecutive(argc, argv);
45 }
46 if (argc > 2) {
47 seed = std::stoi(argv[2]);
48 long enterseed = seed;
49 G4Random::setTheSeed(enterseed);
50 G4Random::showEngineStatus();
51 G4cout << "seed number : " << seed << G4endl;
52 }
53 auto runManager = G4RunManagerFactory::CreateRunManager();
54 // Set mandatory initialization classes
55
56 auto pDetector = new DetectorConstruction();
57 runManager->SetUserInitialization(pDetector);
58 runManager->SetUserInitialization(new PhysicsList(pDetector));
59 runManager->SetUserInitialization(new ActionInitialization(pDetector));
60
61 ///
62 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
63 analysisManager->SetDefaultFileType("root");
64 analysisManager->SetNtupleDirectoryName("ntuple");
65 // open output file
66 //
67 std::string name = "Dose_" + std::to_string(seed);
68 G4bool fileOpen = analysisManager->OpenFile(name);
69 if (!fileOpen) {
70 G4cout << "\n---> HistoManager::book(): cannot open " << name << G4endl;
71 return 0;
72 }
73 analysisManager->SetFirstNtupleId(1);
74 ///////
75 // Visualization
76 G4VisExecutive* visManager = nullptr;
77 // Get the pointer to the User Interface manager
78 G4UImanager* UImanager = G4UImanager::GetUIpointer();
79 if (nullptr == ui) {
80 // Batch mode
81 G4String command = "/control/execute ";
82 G4String fileName = argv[1];
83 UImanager->ApplyCommand(command + fileName);
84 }
85 else {
86 visManager = new G4VisExecutive;
87 visManager->Initialize();
88 UImanager->ApplyCommand("/control/execute vis.mac");
89 ui->SessionStart();
90 delete ui;
91 delete visManager;
92 }
93 analysisManager->Write();
94 analysisManager->CloseFile();
95 delete runManager;
96 return 0;
97 }
98
99 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
100