Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/parameterisations/Par03/examplePar03.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 //
 27 //
 28 //
 29 // --------------------------------------------------------------
 30 //      GEANT 4 - examplePar03
 31 // --------------------------------------------------------------
 32 // Comments
 33 //
 34 // Example of a main program demonstrating multiple energy deposition
 35 // from fast simulation model.
 36 //
 37 //-------------------------------------------------------------------
 38 
 39 #include "FTFP_BERT.hh"
 40 #include "Par03ActionInitialisation.hh"
 41 #include "Par03DetectorConstruction.hh"
 42 
 43 #include "G4EmParameters.hh"
 44 #include "G4FastSimulationPhysics.hh"
 45 #include "G4HadronicProcessStore.hh"
 46 #include "G4RunManagerFactory.hh"
 47 #include "G4Types.hh"
 48 #include "G4UIExecutive.hh"
 49 #include "G4UImanager.hh"
 50 #include "G4VisExecutive.hh"
 51 
 52 #include <sstream>
 53 
 54 int main(int argc, char** argv)
 55 {
 56   // Macro name from arguments
 57   G4String batchMacroName;
 58   G4bool useInteractiveMode = true;
 59   G4String helpMsg(
 60     "Usage: " + G4String(argv[0]) +
 61     " [option(s)] \n No additional arguments triggers an interactive mode "
 62     "executing vis.mac macro. \n Options:\n\t-h\t\tdisplay this help "
 63     "message\n\t-m MACRO\ttriggers a batch mode executing MACRO\n");
 64   for (G4int i = 1; i < argc; ++i) {
 65     G4String argument(argv[i]);
 66     if (argument == "-h" || argument == "--help") {
 67       G4cout << helpMsg << G4endl;
 68       return 0;
 69     }
 70     else if (argument == "-m") {
 71       batchMacroName = G4String(argv[i + 1]);
 72       useInteractiveMode = false;
 73       ++i;
 74     }
 75     else {
 76       G4Exception(
 77         "main", "Unknown argument", FatalErrorInArgument,
 78         ("Unknown argument passed to " + G4String(argv[0]) + " : " + argument + "\n" + helpMsg)
 79           .c_str());
 80     }
 81   }
 82 
 83   // Initialization of default Run manager
 84   auto* runManager = G4RunManagerFactory::CreateRunManager(G4RunManagerType::Default);
 85 
 86   // Detector geometry:
 87   auto detector = new Par03DetectorConstruction();
 88   runManager->SetUserInitialization(detector);
 89 
 90   // Physics list
 91   auto physicsList = new FTFP_BERT();
 92   // Add fast simulation physics
 93   auto fastSimulationPhysics = new G4FastSimulationPhysics();
 94   fastSimulationPhysics->BeVerbose();
 95   fastSimulationPhysics->ActivateFastSimulation("e-");
 96   fastSimulationPhysics->ActivateFastSimulation("e+");
 97   fastSimulationPhysics->ActivateFastSimulation("gamma");
 98   physicsList->RegisterPhysics(fastSimulationPhysics);
 99   // reduce verbosity of physics lists
100   G4EmParameters::Instance()->SetVerbose(0);
101   runManager->SetUserInitialization(physicsList);
102   G4HadronicProcessStore::Instance()->SetVerbose(0);
103 
104   //-------------------------------
105   // UserAction classes
106   //-------------------------------
107   runManager->SetUserInitialization(new Par03ActionInitialisation(detector));
108 
109   //----------------
110   // Visualization:
111   //----------------
112   G4cout << "Instantiating Visualization Manager......." << G4endl;
113   G4VisManager* visManager = new G4VisExecutive;
114   visManager->Initialize();
115   G4UImanager* UImanager = G4UImanager::GetUIpointer();
116 
117   if (useInteractiveMode) {
118     auto ui = new G4UIExecutive(argc, argv);
119     UImanager->ApplyCommand("/control/execute vis.mac");
120     ui->SessionStart();
121     delete ui;
122   }
123   else {
124     G4String command = "/control/execute ";
125     UImanager->ApplyCommand(command + batchMacroName);
126   }
127 
128   // Free the store: user actions, physics_list and detector_description are
129   //                 owned and deleted by the run manager, so they should not
130   //                 be deleted in the main() program !
131 
132   delete visManager;
133   delete runManager;
134 
135   return 0;
136 }
137