hadronic/Hadr01/Hadr01.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /// \file hadronic/Hadr01/Hadr01.cc
0027 /// \brief Main program of the hadronic/Hadr01 example
0028 //
0029 //
0030 //
0031 // -------------------------------------------------------------
0032 //      GEANT4 Hadr01
0033 //
0034 //  Application demonstrating Geant4 hadronic physics:
0035 //  beam interaction with a target
0036 //
0037 //  Authors: A.Bagulya, I.Gudowska, V.Ivanchenko, N.Starkov
0038 //
0039 //  Modified:
0040 //  29.12.2009 V.Ivanchenko introduced access to reference PhysLists
0041 //
0042 // -------------------------------------------------------------
0043 //
0044 //
0045 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0046 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0047
0048 #include "DetectorConstruction.hh"
0049 #include "EventAction.hh"
0050 #include "PhysicsList.hh"
0051 #include "PrimaryGeneratorAction.hh"
0052 #include "RunAction.hh"
0053 #include "StackingAction.hh"
0054
0055 #include "G4HadronicParameters.hh"
0056 #include "G4PhysListFactory.hh"
0057 #include "G4RunManagerFactory.hh"
0058 #include "G4UIExecutive.hh"
0059 #include "G4UImanager.hh"
0060 #include "G4VModularPhysicsList.hh"
0061 #include "G4VisExecutive.hh"
0062
0063 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0064
0065 int main(int argc, char** argv)
0066 {
0067   // detect interactive mode (if no arguments) and define UI session
0068   G4UIExecutive* ui = nullptr;
0069   if (argc == 1) {
0070     ui = new G4UIExecutive(argc, argv);
0071   }
0072
0073   // Construct a serial run manager
0074   auto* runManager = G4RunManagerFactory::CreateRunManager(G4RunManagerType::SerialOnly);
0075
0076   // set mandatory initialization classes
0077   runManager->SetUserInitialization(new DetectorConstruction());
0078
0079   G4PhysListFactory factory;
0080   G4VModularPhysicsList* phys = nullptr;
0081   G4String physName = "";
0082
0083   // Physics List name defined via 3nd argument
0084   if (argc >= 3) {
0085     physName = argv[2];
0086   }
0087
0088   // Physics List name defined via environment variable
0089   if ("" == physName) {
0090     char* path = std::getenv("PHYSLIST");
0091     if (nullptr != path) {
0092       physName = G4String(path);
0093     }
0094   }
0095
0096   G4cout << "PhysicsList: " << physName << G4endl;
0097
0098   // reference PhysicsList via its name
0099   if ("" != physName && factory.IsReferencePhysList(physName)) {
0100     phys = factory.GetReferencePhysList(physName);
0101   }
0102
0103   // local Physics List
0104   if (nullptr == phys) {
0105     phys = new PhysicsList();
0106   }
0107
0108   // optional change of overlap cascade/string
0109   if (argc >= 5) {
0110     auto param = G4HadronicParameters::Instance();
0111     G4double e1 = CLHEP::GeV * std::strtod(argv[3], 0);
0112     G4double e2 = CLHEP::GeV * std::strtod(argv[4], 0);
0113     G4cout << "### Bertini/FTFP limits: e1(GeV)=" << e1 / CLHEP::GeV
0114            << "  e2(GeV)=" << e2 / CLHEP::GeV << G4endl;
0115     param->SetMinEnergyTransitionFTF_Cascade(e1);
0116     param->SetMaxEnergyTransitionFTF_Cascade(e2);
0117   }
0118
0119   if (argc >= 6) {
0120     auto param = G4HadronicParameters::Instance();
0121     param->SetEnableNUDEX(true);
0122   }
0123
0124   // define physics
0125   runManager->SetUserInitialization(phys);
0126   runManager->SetUserAction(new PrimaryGeneratorAction());
0127
0128   // set user action classes
0129   runManager->SetUserAction(new RunAction());
0130   runManager->SetUserAction(new EventAction());
0131   runManager->SetUserAction(new StackingAction());
0132
0133   // initialize visualization
0134   G4VisManager* visManager = nullptr;
0135
0136   // get the pointer to the User Interface manager
0137   G4UImanager* UImanager = G4UImanager::GetUIpointer();
0138
0139   if (ui) {
0140     // interactive mode
0141     visManager = new G4VisExecutive;
0142     visManager->Initialize();
0143     ui->SessionStart();
0144     delete ui;
0145   }
0146   else {
0147     // batch mode
0148     G4String command = "/control/execute ";
0149     G4String fileName = argv[1];
0150     UImanager->ApplyCommand(command + fileName);
0151   }
0152
0153   // job termination
0154   delete visManager;
0155   delete runManager;
0156 }
0157
0158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......