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 // 27 /// \file ExGflasha.cc 28 /// \brief Main program of the parameterisations/gflash/gflasha example 29 // 30 31 // G4 includes 32 #include "G4PhysListFactory.hh" 33 #include "G4RunManagerFactory.hh" 34 #include "G4Timer.hh" 35 #include "G4Types.hh" 36 #include "G4UImanager.hh" 37 #include "G4ios.hh" 38 39 // my project 40 #include "ExGflashActionInitialization.hh" 41 #include "ExGflashDetectorConstruction.hh" 42 43 #include "G4FastSimulationPhysics.hh" 44 #include "G4UIExecutive.hh" 45 #include "G4VisExecutive.hh" 46 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 48 49 int main(int argc, char** argv) 50 { 51 // Instantiate G4UIExecutive if interactive mode 52 G4UIExecutive* ui = nullptr; 53 if (argc == 1) { 54 ui = new G4UIExecutive(argc, argv); 55 } 56 57 // timer to see GFlash performance 58 G4Timer timer; 59 timer.Start(); 60 61 G4cout << "+-------------------------------------------------------+" << G4endl; 62 G4cout << "| |" << G4endl; 63 G4cout << "| This is an example of Shower |" << G4endl; 64 G4cout << "| Parameterization with GFLASH |" << G4endl; 65 G4cout << "+-------------------------------------------------------+" << G4endl; 66 67 auto* runManager = G4RunManagerFactory::CreateRunManager(); 68 69 // UserInitialization classes (mandatory) 70 auto detector = new ExGflashDetectorConstruction; 71 runManager->SetUserInitialization(detector); 72 73 // -- Select a physics list: 74 G4PhysListFactory listFactory; 75 G4String name; 76 auto list_name = std::getenv("PHYSLIST"); 77 if (list_name == nullptr || std::strlen(list_name) == 0) { 78 name = "FTFP_BERT"; 79 } 80 else { 81 name = list_name; 82 } 83 G4VModularPhysicsList* physicsList = listFactory.GetReferencePhysList(name); 84 85 // -- Create a fast simulation physics constructor, used to augment 86 // -- the above physics list to allow for fast simulation: 87 auto fastSimulationPhysics = new G4FastSimulationPhysics(); 88 89 // -- We now configure the fastSimulationPhysics object. 90 // -- The gflash model (GFlashShowerModel, see 91 // ExGflashDetectorConstruction.cc) 92 // -- is applicable to e+ and e- : we augment the physics list for these 93 // -- particles (by adding a G4FastSimulationManagerProcess with below's 94 // -- calls), this will make the fast simulation to be activated: 95 fastSimulationPhysics->ActivateFastSimulation("e-"); 96 fastSimulationPhysics->ActivateFastSimulation("e+"); 97 // -- Register this fastSimulationPhysics to the physicsList, 98 // -- when the physics list will be called by the run manager 99 // -- (will happen at initialization of the run manager) 100 // -- for physics process construction, the fast simulation 101 // -- configuration will be applied as well. 102 physicsList->RegisterPhysics(fastSimulationPhysics); 103 runManager->SetUserInitialization(physicsList); 104 105 // Action initialization: 106 runManager->SetUserInitialization(new ExGflashActionInitialization(detector)); 107 108 G4VisManager* visManager = new G4VisExecutive; 109 visManager->Initialize(); 110 111 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 112 113 if (ui != nullptr) // Define UI terminal for interactive mode 114 { 115 UImanager->ApplyCommand("/control/execute vis.mac"); 116 ui->SessionStart(); 117 delete ui; 118 } 119 else // Batch mode 120 { 121 G4String s = *(argv + 1); 122 UImanager->ApplyCommand("/control/execute " + s); 123 } 124 125 delete visManager; 126 delete runManager; 127 128 timer.Stop(); 129 G4cout << G4endl; 130 G4cout << "******************************************"; 131 G4cout << G4endl; 132 G4cout << "Total Real Elapsed Time is: " << timer.GetRealElapsed(); 133 G4cout << G4endl; 134 G4cout << "Total System Elapsed Time: " << timer.GetSystemElapsed(); 135 G4cout << G4endl; 136 G4cout << "Total GetUserElapsed Time: " << timer.GetUserElapsed(); 137 G4cout << G4endl; 138 G4cout << "******************************************"; 139 G4cout << G4endl; 140 141 return 0; 142 } 143 144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 145