Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/physicslists/factory/factory.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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 factory.cc
 27 /// \brief Main program of the physicslists/factory example
 28 //
 29 //
 30 //
 31 // -------------------------------------------------------------
 32 //      factory
 33 //
 34 //  Application demonstrating the usage of G4PhysListFactory
 35 //  to build the concrete physics list
 36 //
 37 //  Modified from hadronic/Hadr00/Hadr00.cc
 38 //     Author: V.Ivanchenko 20 June 2008
 39 // -------------------------------------------------------------
 40 //
 41 //
 42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 44 
 45 #include "ActionInitialization.hh"
 46 #include "DetectorConstruction.hh"
 47 #include "PrimaryGeneratorAction.hh"
 48 
 49 #include "G4PhysListFactory.hh"
 50 #include "G4RunManagerFactory.hh"
 51 #include "G4UIExecutive.hh"
 52 #include "G4UImanager.hh"
 53 #include "G4VModularPhysicsList.hh"
 54 #include "G4VisExecutive.hh"
 55 #include "Randomize.hh"
 56 
 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 58 
 59 namespace
 60 {
 61 void PrintUsage()
 62 {
 63   G4cerr << " Usage: " << G4endl;
 64   G4cerr << " factory [-m macro ] [-p physList ] [-u UIsession] [-t nThreads]" << G4endl;
 65   G4cerr << "   note: -t option is available only for multi-threaded mode." << G4endl;
 66   G4cerr << G4endl;
 67 }
 68 }  // namespace
 69 
 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 71 
 72 int main(int argc, char** argv)
 73 {
 74   // Evaluate arguments
 75   //
 76   if (argc > 9) {
 77     PrintUsage();
 78     return 1;
 79   }
 80 
 81   G4String macro;
 82   G4String session;
 83   G4String physListName;
 84   G4String gdmlFileName;
 85 #ifdef G4MULTITHREADED
 86   G4int nofThreads = 0;
 87 #endif
 88   for (G4int i = 1; i < argc; i = i + 2) {
 89     if (G4String(argv[i]) == "-m")
 90       macro = argv[i + 1];
 91     else if (G4String(argv[i]) == "-u")
 92       session = argv[i + 1];
 93     else if (G4String(argv[i]) == "-p")
 94       physListName = argv[i + 1];
 95 #ifdef G4MULTITHREADED
 96     else if (G4String(argv[i]) == "-t") {
 97       nofThreads = G4UIcommand::ConvertToInt(argv[i + 1]);
 98     }
 99 #endif
100     else {
101       PrintUsage();
102       return 1;
103     }
104   }
105 
106   // Detect interactive mode (if no arguments) and define UI session
107   //
108   G4UIExecutive* ui = nullptr;
109   if (!macro.size()) {
110     ui = new G4UIExecutive(argc, argv, session);
111   }
112 
113   // Choose the Random engine  //choose the Random engine
114   G4Random::setTheEngine(new CLHEP::RanecuEngine());
115 
116   // Construct the run manager
117   auto* runManager = G4RunManagerFactory::CreateRunManager();
118 #ifdef G4MULTITHREADED
119   if (nofThreads > 0) {
120     runManager->SetNumberOfThreads(nofThreads);
121   }
122 #endif
123 
124   // Physics list factory
125   G4PhysListFactory factory;
126   G4VModularPhysicsList* physList = nullptr;
127 
128   // Get physics list name
129   if (!physListName.size()) {
130     // Physics List is defined via environment variable PHYSLIST
131     char* physListNameEnv = std::getenv("PHYSLIST");
132     if (physListNameEnv) {
133       physListName = G4String(physListNameEnv);
134     }
135   }
136 
137   // Check if the name is known to the factory
138   if (physListName.size() && (!factory.IsReferencePhysList(physListName))) {
139     G4cerr << "Physics list " << physListName << " is not available in PhysListFactory." << G4endl;
140     physListName.clear();
141   }
142 
143   // If name is not defined use FTFP_BERT
144   if (!physListName.size()) {
145     physListName = "FTFP_BERT";
146   }
147 
148   // Reference PhysicsList via its name
149   physList = factory.GetReferencePhysList(physListName);
150 
151   // Set mandatory initialization classes
152   runManager->SetUserInitialization(new DetectorConstruction());
153   runManager->SetUserInitialization(physList);
154 
155   // set user action classes
156   runManager->SetUserInitialization(new ActionInitialization("factory"));
157 
158   // Initialize visualization
159   G4VisManager* visManager = new G4VisExecutive;
160   // G4VisExecutive can take a verbosity argument - see /vis/verbose guidance.
161   // G4VisManager* visManager = new G4VisExecutive("Quiet");
162   visManager->Initialize();
163 
164   // Get the pointer to the User Interface manager
165   G4UImanager* UImanager = G4UImanager::GetUIpointer();
166 
167   if (macro.size()) {
168     // batch mode
169     G4String command = "/control/execute ";
170     UImanager->ApplyCommand(command + macro);
171   }
172   else {
173     // interactive mode : define UI session
174     UImanager->ApplyCommand("/control/execute init_vis.mac");
175     ui->SessionStart();
176     delete ui;
177   }
178 
179   // Job termination
180   // Free the store: user actions, physics_list and detector_description are
181   // owned and deleted by the run manager, so they should not be deleted
182   // in the main() program !
183 
184   delete visManager;
185   delete runManager;
186 }
187 
188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
189