G4INCLConfig.hh

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
// INCL++ intra-nuclear cascade model
// Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
// Davide Mancusi, CEA
// Alain Boudard, CEA
// Sylvie Leray, CEA
// Joseph Cugnon, University of Liege
//
#define INCLXX_IN_GEANT4_MODE 1

#include "globals.hh"

#ifndef G4INCLConfig_hh
#define G4INCLConfig_hh 1

#include "G4INCLParticleSpecies.hh"
#include "G4INCLConfigEnums.hh"
#include "G4INCLRandomSeedVector.hh"
#include <iostream>
#include <string>
#include <sstream>
// #include <cassert>

#if defined(HAS_BOOST_PROGRAM_OPTIONS) && !defined(INCLXX_IN_GEANT4_MODE)
#include <boost/program_options/options_description.hpp>
#include <boost/program_options/parsers.hpp>
#include <boost/program_options/variables_map.hpp>
#include <fstream>
#include <cstdlib>

namespace po = boost::program_options;
#endif

namespace G4INCL {

  class Config {
  public:
    Config();

    Config(G4int, G4int, ParticleSpecies, G4double);

    Config(G4int argc, char *argv[], G4bool isFullRun);

    ~Config();

    void init();

    std::string summary();

    G4int getVerbosity() const { return verbosity; }

    std::string const &getCalculationTitle() const { return title; }

    std::string const &getOutputFileRoot() const { return outputFileRoot; }

    G4int getNumberOfShots() const { return nShots; }

    G4bool isNaturalTarget() const { return naturalTarget; }

    G4int getTargetA() const { return targetSpecies.theA; }

    G4int getTargetZ() const { return targetSpecies.theZ; }

    void setTargetA(G4int A) { targetSpecies.theA = A; }

    void setTargetZ(G4int Z) { targetSpecies.theZ = Z; }

    ParticleType getProjectileType() const { return projectileSpecies.theType; }

    ParticleSpecies getProjectileSpecies() const { return projectileSpecies; }

    void setProjectileSpecies(ParticleSpecies const &ps) { projectileSpecies=ps; }

    G4double getProjectileKineticEnergy() const { return projectileKineticEnergy; }

    void setProjectileKineticEnergy(G4double const kinE) { projectileKineticEnergy=kinE; }

    G4int getVerboseEvent() const { return verboseEvent; }

    static std::string const getVersionID();

    static std::string const getVersionHash();

    static std::string const getVersionString() {
      std::stringstream ss;
      ss << getVersionID() << "-" << getVersionHash();
      return ss.str();
    }

    Random::SeedVector getRandomSeeds() const {
      Random::SeedVector s;
      s.push_back(randomSeed1);
      s.push_back(randomSeed2);
      return s;
    }

    PauliType getPauliType() const { return pauliType; }

    G4bool getCDPP() const { return CDPP; }

    G4bool getPionPotential() const { return pionPotential; }

    CoulombType getCoulombType() const { return coulombType; }

    PotentialType getPotentialType() const { return potentialType; }

    void setPotentialType(PotentialType type) { potentialType = type; }

    void setPionPotential(const G4bool pionPot) { pionPotential = pionPot; }

    LocalEnergyType getLocalEnergyBBType() const { return localEnergyBBType; }

    LocalEnergyType getLocalEnergyPiType() const { return localEnergyPiType; }

    void setLocalEnergyBBType(const LocalEnergyType t) { localEnergyBBType=t; }

    void setLocalEnergyPiType(const LocalEnergyType t) { localEnergyPiType=t; }

    std::string const &getLogFileName() const { return logFileName; }

    DeExcitationType getDeExcitationType() const { return deExcitationType; }

    std::string getDeExcitationString() const { return deExcitationString; }

    ClusterAlgorithmType getClusterAlgorithm() const { return clusterAlgorithmType; }

    G4int getClusterMaxMass() const { return clusterMaxMass; }

    void setClusterMaxMass(const G4int m){ clusterMaxMass=m; }

    G4bool getBackToSpectator() const { return backToSpectator; }

    G4bool getUseRealMasses() const { return useRealMasses; }

    void setUseRealMasses(G4bool use) { useRealMasses = use; }

    std::string const &getINCLXXDataFilePath() const {
      return INCLXXDataFilePath;
    }

#ifdef INCL_DEEXCITATION_ABLAXX
    std::string const &getABLAv3pCxxDataFilePath() const {
      return ablav3pCxxDataFilePath;
    }
#endif

#ifdef INCL_DEEXCITATION_ABLA07
    std::string const &getABLA07DataFilePath() const {
      return abla07DataFilePath;
    }
#endif
#ifdef INCL_DEEXCITATION_GEMINIXX
    std::string const &getGEMINIXXDataFilePath() const {
      return geminixxDataFilePath;
    }
#endif

    G4double getImpactParameter() const { return impactParameter; }

    SeparationEnergyType getSeparationEnergyType() const { return separationEnergyType; }

    FermiMomentumType getFermiMomentumType() const { return fermiMomentumType; }

    void setFermiMomentumType(FermiMomentumType const f) { fermiMomentumType=f; }

    G4double getCutNN() const { return cutNN; }

#ifdef INCL_ROOT_USE
    std::string const &getROOTSelectionString() const {
      return rootSelectionString;
    }
#endif

#ifdef INCL_DEEXCITATION_FERMI_BREAKUP
    G4int getMaxMassFermiBreakUp() const {
      return maxMassFermiBreakUp;
    }
#endif

    G4double getRPCorrelationCoefficient(const ParticleType t) const {
// assert(t==Proton || t==Neutron);
      return ((t==Proton) ? rpCorrelationCoefficientProton : rpCorrelationCoefficientNeutron);
    }

    void setRPCorrelationCoefficient(const ParticleType t, const G4double corrCoeff) {
// assert(t==Proton || t==Neutron);
      if(t==Proton)
        rpCorrelationCoefficientProton=corrCoeff;
      else
        rpCorrelationCoefficientNeutron=corrCoeff;
    }

    void setRPCorrelationCoefficient(const G4double corrCoeff) {
      setRPCorrelationCoefficient(Proton,corrCoeff);
      setRPCorrelationCoefficient(Neutron,corrCoeff);
    }

    G4double getNeutronSkinThickness() const { return neutronSkinThickness; }

    void setNeutronSkinThickness(const G4double d) { neutronSkinThickness=d; }

    G4double getNeutronSkinAdditionalDiffuseness() const { return neutronSkinAdditionalDiffuseness; }

    void setNeutronSkinAdditionalDiffuseness(const G4double d) { neutronSkinAdditionalDiffuseness=d; }

    G4bool getRefraction() const { return refraction; }

    void setRefraction(const G4bool r) { refraction = r; }

#if defined(HAS_BOOST_PROGRAM_OPTIONS) && !defined(INCLXX_IN_GEANT4_MODE)
    std::string const echo() const;
#endif

  private:

#if defined(HAS_BOOST_PROGRAM_OPTIONS) && !defined(INCLXX_IN_GEANT4_MODE)
    std::string echoOptionsDescription(const po::options_description &aDesc) const;

    po::options_description runOptDesc;
    po::options_description hiddenOptDesc;
    po::options_description genericOptDesc;
    po::options_description physicsOptDesc;
    po::variables_map variablesMap;
#endif

    G4int verbosity;
    std::string inputFileName;
    std::string title;
    std::string outputFileRoot;
    std::string fileSuffix;
    std::string logFileName;

    G4int nShots;

    std::string targetString;
    ParticleSpecies targetSpecies;
    G4bool naturalTarget;

    std::string projectileString;
    ParticleSpecies projectileSpecies;
    G4double projectileKineticEnergy;

    G4int verboseEvent;

    G4int randomSeed1, randomSeed2;
    static const G4int randomSeedMin, randomSeedMax;

    std::string pauliString;
    PauliType pauliType;
    G4bool CDPP;

    std::string coulombString;
    CoulombType coulombType;

    std::string potentialString;
    PotentialType potentialType;
    G4bool pionPotential;

    std::string localEnergyBBString;
    LocalEnergyType localEnergyBBType;

    std::string localEnergyPiString;
    LocalEnergyType localEnergyPiType;

    std::string deExcitationModelList;
    std::string deExcitationOptionDescription;
    std::string deExcitationString;
    DeExcitationType deExcitationType;
#ifdef INCL_DEEXCITATION_ABLAXX
    std::string ablav3pCxxDataFilePath;
#endif
#ifdef INCL_DEEXCITATION_ABLA07
    std::string abla07DataFilePath;
#endif
#ifdef INCL_DEEXCITATION_GEMINIXX
    std::string geminixxDataFilePath;
#endif
    std::string INCLXXDataFilePath;

    std::string clusterAlgorithmString;
    ClusterAlgorithmType clusterAlgorithmType;

    G4int clusterMaxMass;

    G4bool backToSpectator;

    G4bool useRealMasses;

    G4double impactParameter;

    std::string separationEnergyString;
    SeparationEnergyType separationEnergyType;

    std::string fermiMomentumString;
    FermiMomentumType fermiMomentumType;

    G4double cutNN;

#ifdef INCL_ROOT_USE
    std::string rootSelectionString; 
#endif

#ifdef INCL_DEEXCITATION_FERMI_BREAKUP
    G4int maxMassFermiBreakUp;
#endif

    G4double rpCorrelationCoefficient;
    G4double rpCorrelationCoefficientProton;
    G4double rpCorrelationCoefficientNeutron;

    G4double neutronSkinThickness;
    G4double neutronSkinAdditionalDiffuseness;

    G4bool refraction;
  };

}

#endif