#include <G4HadronPhysicsQGSP_BERT.hh>

Inheritance diagram for G4HadronPhysicsQGSP_BERT:

Public Member Functions
	G4HadronPhysicsQGSP_BERT (G4int verbose=1)
	G4HadronPhysicsQGSP_BERT (const G4String &name, G4bool quasiElastic=true)
	~G4HadronPhysicsQGSP_BERT () override=default
void	ConstructParticle () override
void	ConstructProcess () override
	G4HadronPhysicsQGSP_BERT (G4HadronPhysicsQGSP_BERT &)=delete
G4HadronPhysicsQGSP_BERT &	operator= (const G4HadronPhysicsQGSP_BERT &right)=delete
Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")
	G4VPhysicsConstructor (const G4String &name, G4int physics_type)
virtual	~G4VPhysicsConstructor ()
void	SetPhysicsName (const G4String &="")
const G4String &	GetPhysicsName () const
void	SetPhysicsType (G4int)
G4int	GetPhysicsType () const
G4int	GetInstanceID () const
virtual void	TerminateWorker ()
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const

Protected Member Functions
void	CreateModels ()
virtual void	Neutron ()
virtual void	Proton ()
virtual void	PiK ()
virtual void	Others ()
virtual void	DumpBanner ()
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess process, G4ParticleDefinition particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
PhysicsBuilder_V	GetBuilders () const
void	AddBuilder (G4PhysicsBuilderInterface *bld)

Protected Attributes
G4double	minQGSP_proton
G4double	minQGSP_neutron
G4double	minQGSP_pik
G4double	minFTFP_proton
G4double	minFTFP_neutron
G4double	minFTFP_pik
G4double	maxFTFP_proton
G4double	maxFTFP_neutron
G4double	maxFTFP_pik
G4double	minBERT_proton {0.0}
G4double	minBERT_neutron {0.0}
G4double	minBERT_pik {0.0}
G4double	maxBERT_proton
G4double	maxBERT_neutron
G4double	maxBERT_pik
G4double	maxBIC_proton {0.0}
G4double	maxBIC_neutron {0.0}
G4double	minBIC_neutron {0.0}
G4bool	QuasiElasticFTF {false}
G4bool	QuasiElasticQGS {true}
Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel = 0
G4String	namePhysics = ""
G4int	typePhysics = 0
G4ParticleTable *	theParticleTable = nullptr
G4int	g4vpcInstanceID = 0

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Types inherited from G4VPhysicsConstructor
using	PhysicsBuilder_V = G4VPCData::PhysicsBuilders_V
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 52 of file G4HadronPhysicsQGSP_BERT.hh.

Constructor & Destructor Documentation

◆ G4HadronPhysicsQGSP_BERT() [1/3]

G4HadronPhysicsQGSP_BERT::G4HadronPhysicsQGSP_BERT ( G4int verbose = 1 )

Definition at line 92 of file G4HadronPhysicsQGSP_BERT.cc.

    :  G4HadronPhysicsQGSP_BERT("hInelastic QGSP_BERT", true)
{
  G4HadronicParameters::Instance()->SetVerboseLevel(verb);
}

Referenced by G4HadronPhysicsQGS_BIC::G4HadronPhysicsQGS_BIC(), G4HadronPhysicsQGSP_BERT(), G4HadronPhysicsQGSP_BERT(), G4HadronPhysicsQGSP_BERT_HP::G4HadronPhysicsQGSP_BERT_HP(), G4HadronPhysicsQGSP_BIC::G4HadronPhysicsQGSP_BIC(), G4HadronPhysicsQGSP_FTFP_BERT::G4HadronPhysicsQGSP_FTFP_BERT(), and operator=().

◆ G4HadronPhysicsQGSP_BERT() [2/3]

G4HadronPhysicsQGSP_BERT::G4HadronPhysicsQGSP_BERT	(	const G4String &	name,
		G4bool	quasiElastic = true )

Definition at line 98 of file G4HadronPhysicsQGSP_BERT.cc.

    :  G4VPhysicsConstructor(name)
{
  SetPhysicsType(bHadronInelastic);
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  minQGSP_proton = minQGSP_neutron = minQGSP_pik = 
    param->GetMinEnergyTransitionQGS_FTF();
  maxFTFP_proton = maxFTFP_neutron = maxFTFP_pik = 
    param->GetMaxEnergyTransitionQGS_FTF();
  minFTFP_proton = minFTFP_neutron = minFTFP_pik = 
    param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_proton = maxBERT_neutron = maxBERT_pik = 
    param->GetMaxEnergyTransitionFTF_Cascade();
  param->SetEnableBCParticles(true);
}

◆ ~G4HadronPhysicsQGSP_BERT()

G4HadronPhysicsQGSP_BERT::~G4HadronPhysicsQGSP_BERT ( )

overridedefault

◆ G4HadronPhysicsQGSP_BERT() [3/3]

G4HadronPhysicsQGSP_BERT::G4HadronPhysicsQGSP_BERT ( G4HadronPhysicsQGSP_BERT & )

delete

Member Function Documentation

◆ ConstructParticle()

void G4HadronPhysicsQGSP_BERT::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 270 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();
 
  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();
 
  G4ShortLivedConstructor pShortLivedConstructor;
  pShortLivedConstructor.ConstructParticle();  
 
  G4IonConstructor pIonConstructor;
  pIonConstructor.ConstructParticle();
}

◆ ConstructProcess()

void G4HadronPhysicsQGSP_BERT::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 285 of file G4HadronPhysicsQGSP_BERT.cc.

{
  // allow changing of parameters at PreInit
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  minQGSP_proton = minQGSP_neutron = minQGSP_pik = 
    param->GetMinEnergyTransitionQGS_FTF();
  maxFTFP_proton = maxFTFP_neutron = maxFTFP_pik = 
    param->GetMaxEnergyTransitionQGS_FTF();
  minFTFP_proton = minFTFP_neutron = minFTFP_pik = 
    param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_proton = maxBERT_neutron = maxBERT_pik =
    param->GetMaxEnergyTransitionFTF_Cascade();
 
  if ( G4Threading::IsMasterThread() && param->GetVerboseLevel() > 0 )
    DumpBanner();
 
  // build models and x-sections
  CreateModels();
}

◆ CreateModels()

void G4HadronPhysicsQGSP_BERT::CreateModels ( )

protected

Definition at line 114 of file G4HadronPhysicsQGSP_BERT.cc.

{
  Neutron();
  Proton();
  PiK();
  Others();
}

Referenced by ConstructProcess().

◆ DumpBanner()

void G4HadronPhysicsQGSP_BERT::DumpBanner ( )

protectedvirtual

Definition at line 305 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4cout << G4endl;
  G4cout << " " << GetPhysicsName() << "  Thresholds: " << G4endl;
  if ( maxBIC_proton > 0.0 ||  maxBIC_neutron > 0.0) {
    G4cout << "    0) between BIC and BERT for p, n over the interval " 
       << minBERT_proton/CLHEP::GeV << " to "
       << maxBIC_proton/CLHEP::GeV << " GeV. " << G4endl;
  }
  G4cout << "    1) between BERT  and FTF/P over the interval " 
     << minFTFP_proton/CLHEP::GeV << " to "
     << maxBERT_proton/CLHEP::GeV << " GeV. " << G4endl;
  G4cout << "    2) between FTF/P and QGS/P over the interval " 
     << minQGSP_proton/CLHEP::GeV << " to "
     << maxFTFP_proton/CLHEP::GeV << " GeV. " << G4endl;
  G4cout << "  -- quasiElastic: " << QuasiElasticQGS << " for QGS "
     << " and " << QuasiElasticFTF << " for FTF" << G4endl;
}

Referenced by ConstructProcess().

◆ Neutron()

void G4HadronPhysicsQGSP_BERT::Neutron ( )

protectedvirtual

Reimplemented in G4HadronPhysicsQGS_BIC, G4HadronPhysicsQGSP_BERT_HP, and G4HadronPhysicsQGSP_BIC_HP.

Definition at line 122 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  const G4ParticleDefinition* neutron = G4Neutron::Neutron();
  auto inel = new G4HadronInelasticProcess( "neutronInelastic", neutron );
  neutron->GetProcessManager()->AddDiscreteProcess( inel );
 
  G4QGSPNeutronBuilder qgs( QuasiElasticQGS );
  qgs.SetMinEnergy( minQGSP_neutron );
  qgs.Build( inel );
 
  G4FTFPNeutronBuilder ftf( QuasiElasticFTF );
  ftf.SetMinEnergy( minFTFP_neutron );
  ftf.SetMaxEnergy( maxFTFP_neutron );
  ftf.Build( inel );
 
  if ( maxBERT_neutron > minBERT_neutron ) {
    G4BertiniNeutronBuilder bert;
    bert.SetMinEnergy( minBERT_neutron );
    bert.SetMaxEnergy( maxBERT_neutron );
    bert.Build( inel );
  }
  
  if ( maxBIC_neutron > 0.0 ) {
    G4BinaryNeutronBuilder bic;
    bic.SetMinEnergy( minBIC_neutron );
    bic.SetMaxEnergy( maxBIC_neutron );
    bic.Build( inel );
  }
 
  inel->AddDataSet( new G4NeutronInelasticXS() ); 
  if ( useFactorXS ) {
    inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
 
  auto capture = new G4NeutronCaptureProcess( "nCaptureXS" );
  neutron->GetProcessManager()->AddDiscreteProcess(capture);
  capture->RegisterMe( new G4NeutronRadCapture() );
}

Referenced by CreateModels().

◆ operator=()

G4HadronPhysicsQGSP_BERT & G4HadronPhysicsQGSP_BERT::operator= ( const G4HadronPhysicsQGSP_BERT & right )

delete

◆ Others()

void G4HadronPhysicsQGSP_BERT::Others ( )

protectedvirtual

Definition at line 250 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
 
  // high energy particles
  if( param->GetMaxEnergy() > param->EnergyThresholdForHeavyHadrons() ) {
 
    // anti light ions
    G4HadronicBuilder::BuildAntiLightIonsFTFP();
 
    // hyperons
    G4HadronicBuilder::BuildHyperonsQGSP_FTFP_BERT(true);
 
    // b-, c- baryons and mesons
    if( param->EnableBCParticles() ) {
      G4HadronicBuilder::BuildBCHadronsQGSP_FTFP_BERT(true);
    }
  }
}

Referenced by CreateModels().

◆ PiK()

void G4HadronPhysicsQGSP_BERT::PiK ( )

protectedvirtual

Definition at line 203 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  G4PiKBuilder pik;
 
  G4QGSPPiKBuilder qgs( QuasiElasticQGS );
  qgs.SetMinEnergy( minQGSP_pik );
  pik.RegisterMe( &qgs );
 
  G4FTFPPiKBuilder ftf( QuasiElasticFTF );
  ftf.SetMaxEnergy( maxFTFP_pik );
  ftf.SetMinEnergy( minFTFP_pik );
  pik.RegisterMe( &ftf );
 
  G4BertiniPiKBuilder bert;
  bert.SetMaxEnergy( maxBERT_pik );
  pik.RegisterMe( &bert );
 
  // build all models
  pik.Build();
  
  // add cross section factor
  if ( useFactorXS ) {
    const G4ParticleDefinition* pion = G4PionPlus::PionPlus();
    G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess( pion );
    if ( nullptr != inel ) {
      inel->MultiplyCrossSectionBy( param->XSFactorPionInelastic() );
    }
    pion = G4PionMinus::PionMinus();
    inel = G4PhysListUtil::FindInelasticProcess( pion );
    if ( nullptr != inel ) {
      inel->MultiplyCrossSectionBy( param->XSFactorPionInelastic() );
    }
    G4ParticleTable* table = G4ParticleTable::GetParticleTable();
    for ( auto const & pdg : G4HadParticles::GetKaons() ) {
      auto part = table->FindParticle( pdg );
      if ( part == nullptr ) { continue; }
      inel = G4PhysListUtil::FindInelasticProcess(part);
      if ( nullptr != inel ) { 
        inel->MultiplyCrossSectionBy( param->XSFactorHadronInelastic() );
      }
    }
  }
}

Referenced by CreateModels().

◆ Proton()

void G4HadronPhysicsQGSP_BERT::Proton ( )

protectedvirtual

Reimplemented in G4HadronPhysicsQGS_BIC, and G4HadronPhysicsQGSP_BIC_AllHP.

Definition at line 164 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  const G4ParticleDefinition* proton = G4Proton::Proton();
  auto inel = new G4HadronInelasticProcess( "protonInelastic", proton );
  proton->GetProcessManager()->AddDiscreteProcess( inel );
 
  G4QGSPProtonBuilder qgs(QuasiElasticQGS);
  qgs.SetMinEnergy(minQGSP_proton);
  qgs.Build( inel );
 
  G4FTFPProtonBuilder ftf(QuasiElasticFTF);
  ftf.SetMinEnergy( minFTFP_proton );
  ftf.SetMaxEnergy( maxFTFP_proton );
  ftf.Build( inel );
 
  if ( maxBERT_proton >  minBERT_proton ) {
    G4BertiniProtonBuilder bert;
    bert.SetMinEnergy( minBERT_proton );
    bert.SetMaxEnergy( maxBERT_proton );
    bert.Build( inel );
  }
 
  if ( maxBIC_proton > 0.0 ) {
    G4BinaryProtonBuilder bic;
    bic.SetMaxEnergy( maxBIC_proton);
    bic.Build( inel );
  }
 
  auto xsinel = new G4BGGNucleonInelasticXS( proton );
  inel->AddDataSet( xsinel );
 
  if ( useFactorXS ) {
    inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
}