hadronic/FissionFragment/FissionFragment.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 // -------------------------------------------------------------
0027 //      GEANT4 FF_Neutron_HP
0028 //
0029 //  Command line options:
0030 //      -i ARG      : run in batch mode from script file ARG
0031 //      -o ARG      : write output to file ARG
0032 //                    (defaults to FF_Neutron_HP.out)
0033 //      -n ARG      : multithreading with ARG number of threads
0034 //                    (only works if Geant4 was compiled with
0035 //                    multithreading enables)
0036 //
0037 //  =============== Begin Documentation Comments ===============
0038 //!
0039 //! \file       FissionFragment.cc
0040 //! \author     B. Wendt (brycen.linn.wendt@cern.ch)
0041 //! \date       June 06, 2014
0042 //!
0043 //! \brief      Main program of the FissionFragment example
0044 //!
0045 //! \details    Application demonstrating the Fission Fragment model as used
0046 //!                 within the neutron_hp model. It demostrates the capability
0047 //!                 for fission product containment by the cladding in a water
0048 //!                 moderated sub-critical assembly.
0049 //!             It could also be further extended to calculate the effective
0050 //!                 multiplication factor of the subcritical assembly for
0051 //!                 various loading schemes.
0052 //!
0053 //  ================ End Documentation Comments ================
0054 //
0055 //  Modified:
0056 //
0057 //  05-08-20                                              ARibon
0058 //  Replaced deprecated HP environmental variables with UI commands
0059 //  23-06-14                                              BWendt
0060 //  Added check for NeutronHP fission generator environment variable
0061 //
0062 // -------------------------------------------------------------
0063
0064 #include "FFActionInitialization.hh"
0065 #include "FFDetectorConstruction.hh"
0066 #include "QGSP_BIC_HP.hh"
0067
0068 #include "G4ParticleHPManager.hh"
0069 #include "G4RunManagerFactory.hh"
0070 #include "G4UIExecutive.hh"
0071 #include "G4UImanager.hh"
0072 #include "G4VisExecutive.hh"
0073 #include "Randomize.hh"
0074 #include "globals.hh"
0075
0076 // Entry point
0077 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0078 int main(int argc, char* argv[])
0079 {
0080   int result;
0081   unsigned int numberOfThreads = 1;
0082
0083   G4String scriptFileName = "";
0084   G4String outputFileName = "FF_Neutron_HP.out";
0085   G4UImanager* UIManager = NULL;
0086
0087   // Activate production of fission fragments in neutronHP
0088   G4ParticleHPManager::GetInstance()->SetProduceFissionFragments(true);
0089
0090   char Force[] = "G4FORCENUMBEROFTHREADS";
0091   if (std::getenv(Force) != NULL) {
0092     char doNotForce[] = "G4FORCENUMBEROFTHREADS=1";
0093     putenv(doNotForce);
0094   }
0095
0096   // Indicate the example is starting
0097   G4cout << "####   Starting: " << argv[0] << "    ####" << G4endl;
0098
0099   //  Parse the command line arguments, if any
0100   for (int i = 1; i < argc; i += 2) {
0101     // Ensure that this is actually a command
0102     if (argv[i][0] != '-') {
0103       G4cerr << G4endl << "!!!!" << G4endl;
0104       G4cerr << "!!!! Error in argument " << i + 1 << G4endl;
0105       G4cerr << "!!!! A command-line option was expected, but \"" << argv[i] << "\" was found"
0106              << G4endl;
0107       G4cerr << "!!!! " << argv[0] << " will now terminate" << G4endl;
0108       G4cerr << "!!!!" << G4endl << G4endl;
0109
0110       return EXIT_FAILURE;
0111     }
0112
0113     // Ensure that the command-line option has an associated argument
0114     if (!(i + 1 < argc)) {
0115       G4cerr << G4endl << "!!!!" << G4endl;
0116       G4cerr << "!!!! Error in argument " << i + 2 << G4endl;
0117       G4cerr << "!!!! An argument was expected, but \"" << argv[i + 1] << "\" was found" << G4endl;
0118       G4cerr << "!!!! Ensure that a space is used to separate the "
0119                 "option and argument"
0120              << G4endl;
0121       G4cerr << "!!!! " << argv[0] << " will now terminate" << G4endl;
0122       G4cerr << "!!!!" << G4endl << G4endl;
0123
0124       return EXIT_FAILURE;
0125     }
0126
0127     switch (argv[i][1]) {
0128       case 'i':
0129         scriptFileName = "/control/execute ";
0130         scriptFileName.append(argv[i + 1]);
0131         break;
0132
0133       case 'o':
0134         outputFileName = argv[i + 1];
0135         break;
0136
0137       case 'n':
0138         result = sscanf(argv[i + 1], "%u", &numberOfThreads);
0139         if (result != 1) {
0140           G4cerr << G4endl << "!!!!" << G4endl;
0141           G4cerr << "!!!! Error in argument " << i + 2 << G4endl;
0142           G4cerr << "!!!! An positive number was expected, but \"" << argv[i + 1] << "\" was found"
0143                  << G4endl;
0144           G4cerr << "!!!! " << argv[0] << " will now terminate" << G4endl;
0145           G4cerr << "!!!!" << G4endl << G4endl;
0146
0147           return EXIT_FAILURE;
0148         }
0149         break;
0150
0151       default:
0152         G4cout << G4endl << "!!!!" << G4endl;
0153         G4cout << "!!!! Warning for command " << i + 1 << G4endl;
0154         G4cout << "!!!! \"" << argv[i] << "\" is not a valid command" << G4endl;
0155         G4cout << "!!!! " << argv[0] << " will ignore \"" << argv[i] << "\" and \"" << argv[i + 1]
0156                << "\"" << G4endl;
0157         G4cout << "!!!!" << G4endl << G4endl;
0158     }
0159   }
0160
0161   // Instantiate G4UIExecutive if interactive mode
0162   G4UIExecutive* ui = nullptr;
0163   if (scriptFileName.length() == 0) {
0164     ui = new G4UIExecutive(argc, argv);
0165   }
0166
0167   // Set the Random engine
0168   // A seed of 62737819 produced a maximum number of 67 events on the
0169   // author's system before timing out the nightly test
0170   const G4long seed = 62737819;
0171 #ifndef NDEBUG
0172   G4cout << "MT RNG Seed: " << seed << G4endl;
0173 #endif  // NDEBUG
0174   G4Random::setTheEngine(new CLHEP::MTwistEngine(seed));
0175
0176   // Initialize the multithreaded run manager
0177   auto* runManager = G4RunManagerFactory::CreateRunManager();
0178   runManager->SetNumberOfThreads(numberOfThreads);
0179   G4cout << "    Threads requested:    " << numberOfThreads << G4endl;
0180   G4cout << "    Threads started:      " << runManager->GetNumberOfThreads() << G4endl;
0181
0182   // Set mandatory initialization classes
0183   runManager->SetUserInitialization(new FFDetectorConstruction());
0184   runManager->SetUserInitialization(new QGSP_BIC_HP());
0185   runManager->SetUserInitialization(new FFActionInitialization());
0186
0187   // Initialize the Geant4 kernel
0188   runManager->Initialize();
0189
0190   // Initialize visualization
0191   G4VisManager* visManager = new G4VisExecutive();
0192   visManager->Initialize();
0193
0194   // Get the pointer to the User Interface manager
0195   UIManager = G4UImanager::GetUIpointer();
0196
0197   if (!ui) {
0198     // Batch mode
0199     UIManager->ApplyCommand(scriptFileName);
0200   }
0201   else {
0202     // Interactive mode
0203     ui->SessionStart();
0204     delete ui;
0205   }
0206
0207   // Job termination
0208   // Free the store: user actions, physics_list and detector_description are
0209   // owned and deleted by the run manager, so they should not be deleted
0210   // in the main() program !
0211   delete visManager;
0212   delete runManager;
0213
0214   return 0;
0215 }