Gabriele Cosmo
434 lines
16 KiB
C++
434 lines
16 KiB
C++
//
|
|
// ********************************************************************
|
|
// * 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. *
|
|
// ********************************************************************
|
|
//
|
|
//
|
|
//---------------------------------------------------------------------------
|
|
//
|
|
// ClassName: G4HadronicParameters
|
|
//
|
|
// Author: 2018 Alberto Ribon
|
|
//
|
|
// Description: Singleton to keep global hadronic parameters.
|
|
//
|
|
// Modified:
|
|
//
|
|
//----------------------------------------------------------------------------
|
|
//
|
|
#ifndef G4HadronicParameters_h
|
|
#define G4HadronicParameters_h 1
|
|
|
|
#include "globals.hh"
|
|
|
|
class G4HadronicParametersMessenger;
|
|
|
|
|
|
class G4HadronicParameters {
|
|
public:
|
|
|
|
static G4HadronicParameters* Instance();
|
|
~G4HadronicParameters();
|
|
|
|
// printing
|
|
void StreamInfo(std::ostream& os) const;
|
|
void Dump() const;
|
|
|
|
inline G4double GetMaxEnergy() const;
|
|
void SetMaxEnergy( const G4double val );
|
|
// Getter/Setter for the upper limit for Geant4 hadronic physics, for any application.
|
|
// Any hadronic model, physics list builder and constructor should use this method
|
|
// instead of putting an arbitrary value in the code.
|
|
// Any application which tries to use hadronic physics for an energy higher than this limit
|
|
// will get a run-time crash, because no model is found.
|
|
|
|
inline G4double GetMinEnergyTransitionFTF_Cascade() const;
|
|
inline G4double GetMaxEnergyTransitionFTF_Cascade() const;
|
|
void SetMinEnergyTransitionFTF_Cascade( const G4double val );
|
|
void SetMaxEnergyTransitionFTF_Cascade( const G4double val );
|
|
// Getter/Setter of the recommended energy limits, for physics lists, of the
|
|
// transition region between the Fritiof (FTF) string model and the
|
|
// intranuclear cascade model, either Bertini (BERT) or Binary (BIC).
|
|
|
|
inline G4double GetMinEnergyTransitionQGS_FTF() const;
|
|
inline G4double GetMaxEnergyTransitionQGS_FTF() const;
|
|
void SetMinEnergyTransitionQGS_FTF( const G4double val );
|
|
void SetMaxEnergyTransitionQGS_FTF( const G4double val );
|
|
// Getter/Setter of the recommended energy limits, for physics lists, of the
|
|
// transition region between the two strings models - the Quark Gluon String (QGS)
|
|
// model and the Fritiof (FTF) model.
|
|
|
|
inline G4double GetMinEnergyINCLXX_Pbar() const;
|
|
inline G4double GetMaxEnergyINCLXX_Pbar() const;
|
|
void SetMinEnergyINCLXX_Pbar( const G4double val );
|
|
void SetMaxEnergyINCLXX_Pbar( const G4double val );
|
|
// Getter/Setter of the recommended energy limits, for physics lists, of the
|
|
// intranuclear cascade model INCLXX, for pbar interaction.
|
|
|
|
inline G4double EnergyThresholdForHeavyHadrons() const;
|
|
void SetEnergyThresholdForHeavyHadrons( G4double val );
|
|
// If max kinetic energy is below this limit, then EM and hadronic physics are not
|
|
// instantiated for hyperons, anti-hyperons, anti light ions, b-, c- particles.
|
|
|
|
inline G4double XSFactorNucleonInelastic() const;
|
|
void SetXSFactorNucleonInelastic( G4double val );
|
|
inline G4double XSFactorNucleonElastic() const;
|
|
void SetXSFactorNucleonElastic( G4double val );
|
|
// Cross section factor for protons and neutrons.
|
|
|
|
inline G4double XSFactorPionInelastic() const;
|
|
void SetXSFactorPionInelastic( G4double val );
|
|
inline G4double XSFactorPionElastic() const;
|
|
void SetXSFactorPionElastic( G4double val );
|
|
// Cross section factor for pions.
|
|
|
|
inline G4double XSFactorHadronInelastic() const;
|
|
void SetXSFactorHadronInelastic( G4double val );
|
|
inline G4double XSFactorHadronElastic() const;
|
|
void SetXSFactorHadronElastic( G4double val );
|
|
// Cross section factor for other hadrons and ions.
|
|
|
|
inline G4double XSFactorEM() const;
|
|
void SetXSFactorEM( G4double val );
|
|
// Cross section factor for gamma and leptons.
|
|
|
|
inline G4bool EnableBCParticles() const;
|
|
void SetEnableBCParticles( G4bool val );
|
|
// Baryons and mesons with c- and b- quarks may be enabled/disabled.
|
|
// This flag is used both by EM and hadronic physics constructors.
|
|
|
|
inline G4bool EnableHyperNuclei() const;
|
|
void SetEnableHyperNuclei( G4bool val );
|
|
// Light hyper-nuclei may be enabled/disabled.
|
|
// This flag is used both by EM and hadronic physics constructors.
|
|
|
|
inline G4bool ApplyFactorXS() const;
|
|
void SetApplyFactorXS( G4bool val );
|
|
// Flag enabling cross section factor definition.
|
|
|
|
inline G4int GetVerboseLevel() const;
|
|
void SetVerboseLevel( const G4int val );
|
|
// Getter/Setter of the general verbosity level for hadronics.
|
|
|
|
inline G4bool EnableCRCoalescence() const;
|
|
void SetEnableCRCoalescence( G4bool val );
|
|
// Boolean switch that allows to apply the Cosmic Ray (CR) coalescence algorithm
|
|
// to the secondaries produced by a string model. By default it is disabled.
|
|
|
|
inline G4bool EnableIntegralInelasticXS() const;
|
|
inline G4bool EnableIntegralElasticXS() const;
|
|
void SetEnableIntegralInelasticXS( G4bool val );
|
|
void SetEnableIntegralElasticXS( G4bool val );
|
|
// Enable/disable integral method for main types of hadrons.
|
|
|
|
inline G4bool EnableDiffDissociationForBGreater10() const;
|
|
// For nucleon-hadron interactions, it's not decided what to do with diffraction
|
|
// dissociation. For the moment, they are turned off. This option allows it to
|
|
// be turned back on. Applies to Baryon Number > 10 or # target nucleons > 10.
|
|
void SetEnableDiffDissociationForBGreater10(G4bool val);
|
|
|
|
inline G4bool EnableCoherentChargeExchange() const;
|
|
void SetEnableCoherentChargeExchange( G4bool val );
|
|
// Coherent Charge exchange process may be enabled/disabled.
|
|
|
|
inline G4bool EnableNeutronGeneralProcess() const;
|
|
void SetEnableNeutronGeneralProcess( G4bool val );
|
|
// Neutron general process may be enabled/disabled.
|
|
|
|
inline G4bool EnableNUDEX() const;
|
|
void SetEnableNUDEX( G4bool val );
|
|
// NUDEX gamma de-excitation is enabled/disabled.
|
|
|
|
inline const G4String& GetTypeTablePT() const;
|
|
void SetTypeTablePT( const G4String& typeTablePT );
|
|
// Specify the type of PT table - between "calendf" and "njoy" for
|
|
// the URR (Unresolved Resonance Region) treatment of low-energy neutrons.
|
|
// ( Note that there is no default: an empty string "" is returned if
|
|
// it is not set explicitly. )
|
|
|
|
inline G4double GetEPRelativeLevel() const;
|
|
inline G4double GetEPAbsoluteLevel() const;
|
|
inline G4int GetEPReportLevel() const;
|
|
inline G4bool GetBinaryDebug() const;
|
|
inline const G4String& GetDirPARTICLEXS() const;
|
|
inline const G4String& GetPhysListDocDir() const;
|
|
inline const G4String& GetPhysListName() const;
|
|
// Access to environment variables.
|
|
|
|
inline G4double GetNeutronKineticEnergyThresholdForSVT() const;
|
|
void SetNeutronKineticEnergyThresholdForSVT( const G4double val );
|
|
// Getter/Setter for the neutron kinetic energy threshold for
|
|
// applying the SVT (Sampling of the Velocity of the Target) algorithm.
|
|
|
|
inline G4double GetTimeThresholdForRadioactiveDecay() const;
|
|
void SetTimeThresholdForRadioactiveDecay( const G4double val );
|
|
// Getter/Setter for the time threshold of radioactive decays
|
|
// (i.e. radioactive decays that happen later than this value are ignored).
|
|
|
|
inline G4bool IsBertiniAs11_2() const;
|
|
void SetBertiniAs11_2( G4bool val );
|
|
inline G4bool IsBertiniAngularEmissionsAs11_2() const;
|
|
void SetBertiniAngularEmissionsAs11_2( G4bool val );
|
|
inline G4bool IsBertiniNucleiModelAs11_2() const;
|
|
void SetBertiniNucleiModelAs11_2( G4bool val );
|
|
// Getter/Setter for the Bertini model behavior with respect to the
|
|
// Geant4 version 11.2.
|
|
// The first two methods "*BertiniAs11_2" refer to the overall behavior
|
|
// of the Bertini model - the one which should matter the most for the
|
|
// majority of applications.
|
|
// The methods "*BertiniAngularEmissionsAs11_2" refer to the angular
|
|
// distributions of the >= 4-body final state emissions.
|
|
// The last two methods "*BertiniNucleiModelAs11_2" refer to the
|
|
// modelling of nuclei.
|
|
// The overall behavior of Bertini with respect to the version Geant4
|
|
// version 11.2 depends on both the angular distributios of the
|
|
// >= 4-body final state emissions, and the modelling of nuclei.
|
|
// "True" in these methods means that the corresponding behavior of
|
|
// the Geant4 version 11.2 is retrieved;
|
|
// "False" means that the corresponding behavior of the Geant4 version
|
|
// 11.3 is retrieved.
|
|
|
|
inline G4bool UseRFilesForXS() const;
|
|
void SetUseRFilesForXS( G4bool val );
|
|
|
|
private:
|
|
|
|
G4HadronicParameters();
|
|
|
|
G4bool IsLocked() const;
|
|
|
|
static G4HadronicParameters* sInstance;
|
|
|
|
G4HadronicParametersMessenger* fMessenger;
|
|
|
|
G4double fMaxEnergy;
|
|
G4double fMinEnergyTransitionFTF_Cascade;
|
|
G4double fMaxEnergyTransitionFTF_Cascade;
|
|
G4double fMinEnergyTransitionQGS_FTF;
|
|
G4double fMaxEnergyTransitionQGS_FTF;
|
|
G4double fMinEnergyINCLXX_Pbar;
|
|
G4double fMaxEnergyINCLXX_Pbar;
|
|
G4double fEnergyThresholdForHeavyHadrons;
|
|
G4double fXSFactorNucleonInelastic = 1.0;
|
|
G4double fXSFactorPionInelastic = 1.0;
|
|
G4double fXSFactorHadronInelastic = 1.0;
|
|
G4double fXSFactorNucleonElastic = 1.0;
|
|
G4double fXSFactorPionElastic = 1.0;
|
|
G4double fXSFactorHadronElastic = 1.0;
|
|
G4double fXSFactorEM = 1.0;
|
|
G4double fXSFactorLimit = 0.2;
|
|
G4double fRelativeDiff = DBL_MAX;
|
|
G4double fAbsoluteDiff = DBL_MAX;
|
|
G4double fNeutronEkinThresholdForSVT = -1.0;
|
|
G4double fTimeThresholdForRadioactiveDecays = -1.0;
|
|
|
|
G4int fVerboseLevel = 1;
|
|
G4int fReportLevel = 0;
|
|
|
|
G4bool fEnableBC = false;
|
|
G4bool fEnableHyperNuclei = false;
|
|
G4bool fApplyFactorXS = false;
|
|
G4bool fEnableCRCoalescence = false;
|
|
G4bool fEnableIntegralInelasticXS = true;
|
|
G4bool fEnableIntegralElasticXS = true;
|
|
G4bool fEnableDiffDissociationForBGreater10 = false;
|
|
G4bool fEnableNUDEX = false;
|
|
G4bool fNeutronGeneral = false;
|
|
G4bool fChargeExchange = false;
|
|
G4bool fBinaryDebug = false;
|
|
G4bool fBertiniAngularEmissionsAs11_2 = false; // Keep the new G4 11.3 behavior
|
|
G4bool fBertiniNucleiModelAs11_2 = false; // Keep the new G4 11.3 behavior
|
|
G4bool fUseRFilesForXS = false;
|
|
|
|
G4String fTypeTablePT = "";
|
|
G4String fDirPARTICLEXS = "";
|
|
G4String fPhysListDocDir = "";
|
|
G4String fPhysListName = "";
|
|
};
|
|
|
|
inline G4double G4HadronicParameters::GetMaxEnergy() const {
|
|
return fMaxEnergy;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetMinEnergyTransitionFTF_Cascade() const {
|
|
return fMinEnergyTransitionFTF_Cascade;
|
|
}
|
|
inline G4double G4HadronicParameters::GetMaxEnergyTransitionFTF_Cascade() const {
|
|
return fMaxEnergyTransitionFTF_Cascade;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetMinEnergyTransitionQGS_FTF() const {
|
|
return fMinEnergyTransitionQGS_FTF;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetMaxEnergyTransitionQGS_FTF() const {
|
|
return fMaxEnergyTransitionQGS_FTF;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetMinEnergyINCLXX_Pbar() const {
|
|
return fMinEnergyINCLXX_Pbar;
|
|
}
|
|
inline G4double G4HadronicParameters::GetMaxEnergyINCLXX_Pbar() const {
|
|
return fMaxEnergyINCLXX_Pbar;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::EnergyThresholdForHeavyHadrons() const {
|
|
return fEnergyThresholdForHeavyHadrons;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorNucleonInelastic() const {
|
|
return fXSFactorNucleonInelastic;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorNucleonElastic() const {
|
|
return fXSFactorNucleonElastic;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorPionInelastic() const {
|
|
return fXSFactorPionInelastic;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorPionElastic() const {
|
|
return fXSFactorPionElastic;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorHadronInelastic() const {
|
|
return fXSFactorHadronInelastic;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorHadronElastic() const {
|
|
return fXSFactorHadronElastic;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::XSFactorEM() const {
|
|
return fXSFactorEM;
|
|
}
|
|
|
|
inline G4int G4HadronicParameters::GetVerboseLevel() const {
|
|
return fVerboseLevel;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableBCParticles() const {
|
|
return fEnableBC;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableHyperNuclei() const {
|
|
return fEnableHyperNuclei;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::ApplyFactorXS() const {
|
|
return fApplyFactorXS;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableCRCoalescence() const {
|
|
return fEnableCRCoalescence;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableIntegralInelasticXS() const {
|
|
return fEnableIntegralInelasticXS;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableIntegralElasticXS() const {
|
|
return fEnableIntegralElasticXS;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableDiffDissociationForBGreater10() const {
|
|
return fEnableDiffDissociationForBGreater10;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableNeutronGeneralProcess() const {
|
|
return fNeutronGeneral;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableNUDEX() const {
|
|
return fEnableNUDEX;
|
|
}
|
|
|
|
inline const G4String& G4HadronicParameters::GetTypeTablePT() const {
|
|
return fTypeTablePT;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::EnableCoherentChargeExchange() const {
|
|
return fChargeExchange;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::GetBinaryDebug() const {
|
|
return fBinaryDebug;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::UseRFilesForXS() const {
|
|
return fUseRFilesForXS;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetEPRelativeLevel() const {
|
|
return fRelativeDiff;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetEPAbsoluteLevel() const {
|
|
return fAbsoluteDiff;
|
|
}
|
|
|
|
inline G4int G4HadronicParameters::GetEPReportLevel() const {
|
|
return fReportLevel;
|
|
}
|
|
|
|
inline const G4String& G4HadronicParameters::GetDirPARTICLEXS() const {
|
|
return fDirPARTICLEXS;
|
|
}
|
|
|
|
inline const G4String& G4HadronicParameters::GetPhysListDocDir() const {
|
|
return fPhysListDocDir;
|
|
}
|
|
|
|
inline const G4String& G4HadronicParameters::GetPhysListName() const {
|
|
return fPhysListName;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetNeutronKineticEnergyThresholdForSVT() const {
|
|
return fNeutronEkinThresholdForSVT;
|
|
}
|
|
|
|
inline G4double G4HadronicParameters::GetTimeThresholdForRadioactiveDecay() const {
|
|
return fTimeThresholdForRadioactiveDecays;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::IsBertiniAs11_2() const {
|
|
return ( fBertiniAngularEmissionsAs11_2 && fBertiniNucleiModelAs11_2 );
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::IsBertiniAngularEmissionsAs11_2() const {
|
|
return fBertiniAngularEmissionsAs11_2;
|
|
}
|
|
|
|
inline G4bool G4HadronicParameters::IsBertiniNucleiModelAs11_2() const {
|
|
return fBertiniNucleiModelAs11_2;
|
|
}
|
|
|
|
inline std::ostream& operator<<(std::ostream& os, const G4HadronicParameters& p)
|
|
{
|
|
p.StreamInfo(os);
|
|
return os;
|
|
}
|
|
|
|
#endif
|