#include <G4HadronPhysicsFTFP_BERT.hh>

Inheritance diagram for G4HadronPhysicsFTFP_BERT:

Public Member Functions
	G4HadronPhysicsFTFP_BERT (G4int verbose=1)
	G4HadronPhysicsFTFP_BERT (const G4String &name, G4bool quasiElastic=false)
	~G4HadronPhysicsFTFP_BERT () override=default
void	ConstructParticle () override
void	ConstructProcess () override
	G4HadronPhysicsFTFP_BERT (G4HadronPhysicsFTFP_BERT &)=delete
G4HadronPhysicsFTFP_BERT &	operator= (const G4HadronPhysicsFTFP_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	Pion ()
virtual void	Kaon ()
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	minFTFP_pion
G4double	maxBERT_pion
G4double	minFTFP_kaon
G4double	maxBERT_kaon
G4double	minFTFP_proton
G4double	maxBERT_proton
G4double	minFTFP_neutron
G4double	maxBERT_neutron
G4double	minBERT_proton
G4double	minBERT_neutron
G4bool	QuasiElastic
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 48 of file G4HadronPhysicsFTFP_BERT.hh.

Constructor & Destructor Documentation

◆ G4HadronPhysicsFTFP_BERT() [1/3]

G4HadronPhysicsFTFP_BERT::G4HadronPhysicsFTFP_BERT ( G4int verbose = 1 )

Definition at line 84 of file G4HadronPhysicsFTFP_BERT.cc.

                                                             :
  G4HadronPhysicsFTFP_BERT("hInelastic FTFP_BERT",false) 
{
  G4HadronicParameters::Instance()->SetVerboseLevel(verb);
}

◆ G4HadronPhysicsFTFP_BERT() [2/3]

G4HadronPhysicsFTFP_BERT::G4HadronPhysicsFTFP_BERT	(	const G4String &	name,
		G4bool	quasiElastic = false )

Definition at line 90 of file G4HadronPhysicsFTFP_BERT.cc.

  : G4VPhysicsConstructor(name), QuasiElastic(qe)
{
  SetPhysicsType(bHadronInelastic);
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  minFTFP_pion = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_pion = param->GetMaxEnergyTransitionFTF_Cascade();
  minFTFP_kaon = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_kaon = param->GetMaxEnergyTransitionFTF_Cascade();
  minFTFP_proton = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_proton = param->GetMaxEnergyTransitionFTF_Cascade();
  minFTFP_neutron = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_neutron = param->GetMaxEnergyTransitionFTF_Cascade();
  minBERT_proton = minBERT_neutron = 0.0;
  param->SetEnableBCParticles(true);
}

◆ ~G4HadronPhysicsFTFP_BERT()

G4HadronPhysicsFTFP_BERT::~G4HadronPhysicsFTFP_BERT ( )

overridedefault

◆ G4HadronPhysicsFTFP_BERT() [3/3]

G4HadronPhysicsFTFP_BERT::G4HadronPhysicsFTFP_BERT ( G4HadronPhysicsFTFP_BERT & )

delete

Member Function Documentation

◆ ConstructParticle()

void G4HadronPhysicsFTFP_BERT::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 107 of file G4HadronPhysicsFTFP_BERT.cc.

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

◆ ConstructProcess()

void G4HadronPhysicsFTFP_BERT::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Reimplemented in G4HadronPhysicsFTFP_BERT_ATL, G4HadronPhysicsFTFQGSP_BERT, G4HadronPhysicsINCLXX, and G4HadronPhysicsNuBeam.

Definition at line 290 of file G4HadronPhysicsFTFP_BERT.cc.

{
  // allow changing of parameters at PreInit
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  minFTFP_pion = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_pion = param->GetMaxEnergyTransitionFTF_Cascade();
  minFTFP_kaon = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_kaon = param->GetMaxEnergyTransitionFTF_Cascade();
  minFTFP_proton = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_proton = param->GetMaxEnergyTransitionFTF_Cascade();
  minFTFP_neutron = param->GetMinEnergyTransitionFTF_Cascade();
  maxBERT_neutron = param->GetMaxEnergyTransitionFTF_Cascade();
 
  if (G4Threading::IsMasterThread() && param->GetVerboseLevel() > 0) {
    DumpBanner();
  }
  CreateModels();
}

◆ CreateModels()

void G4HadronPhysicsFTFP_BERT::CreateModels ( )

protected

Definition at line 133 of file G4HadronPhysicsFTFP_BERT.cc.

{
  Neutron();
  Proton();
  Pion();
  Kaon();
  Others();
} 

Referenced by ConstructProcess(), G4HadronPhysicsFTFP_BERT_ATL::ConstructProcess(), G4HadronPhysicsINCLXX::ConstructProcess(), and G4HadronPhysicsNuBeam::ConstructProcess().

◆ DumpBanner()

void G4HadronPhysicsFTFP_BERT::DumpBanner ( )

protectedvirtual

Reimplemented in G4HadronPhysicsFTFQGSP_BERT.

Definition at line 122 of file G4HadronPhysicsFTFP_BERT.cc.

{
  G4cout << G4endl
       << " " << GetPhysicsName() << " : threshold between BERT and FTFP is over the interval " << G4endl
       << " for pions :   " << minFTFP_pion/GeV << " to " << maxBERT_pion/GeV  << " GeV" << G4endl
       << " for kaons :   " << minFTFP_kaon/GeV << " to " << maxBERT_kaon/GeV  << " GeV" << G4endl
       << " for proton :  " << minFTFP_proton/GeV << " to " << maxBERT_proton/GeV  << " GeV" << G4endl
       << " for neutron : " << minFTFP_neutron/GeV << " to " << maxBERT_neutron/GeV  << " GeV" << G4endl
       << G4endl;
}

Referenced by ConstructProcess(), G4HadronPhysicsFTFP_BERT_ATL::ConstructProcess(), G4HadronPhysicsINCLXX::ConstructProcess(), G4HadronPhysicsNuBeam::ConstructProcess(), and G4HadronPhysicsFTFQGSP_BERT::DumpBanner().

◆ Kaon()

void G4HadronPhysicsFTFP_BERT::Kaon ( )

protectedvirtual

Reimplemented in G4HadronPhysicsFTF_BIC, and G4HadronPhysicsINCLXX.

Definition at line 233 of file G4HadronPhysicsFTFP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto k = new G4KaonBuilder;
  AddBuilder(k);
  auto ftfpk = new G4FTFPKaonBuilder(QuasiElastic);
  AddBuilder(ftfpk);
  k->RegisterMe(ftfpk);
  ftfpk->SetMinEnergy(minFTFP_kaon);
  auto bertk  = new G4BertiniKaonBuilder;
  AddBuilder(bertk);
  k->RegisterMe(bertk);
  bertk->SetMaxEnergy(maxBERT_kaon);
  k->Build();
 
  // add cross section factor
  if( useFactorXS ) {
    G4ParticleTable* table = G4ParticleTable::GetParticleTable();
    for( auto & pdg : G4HadParticles::GetKaons() ) {
      auto part = table->FindParticle( pdg );
      if ( part == nullptr ) { continue; }
      G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(part);
      if(nullptr != inel) { 
        inel->MultiplyCrossSectionBy( param->XSFactorHadronInelastic() );
      }
    }
  }
}

Referenced by CreateModels().

◆ Neutron()

void G4HadronPhysicsFTFP_BERT::Neutron ( )

protectedvirtual

Reimplemented in G4HadronPhysicsFTF_BIC, G4HadronPhysicsFTFP_BERT_HP, G4HadronPhysicsINCLXX, and G4HadronPhysicsShielding.

Definition at line 142 of file G4HadronPhysicsFTFP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
  //General schema:
  // 1) Create a builder
  // 2) Call AddBuilder
  // 3) Configure the builder, possibly with sub-builders
  // 4) Call builder->Build()
  auto neu = new G4NeutronBuilder;
  AddBuilder(neu);
  auto ftfpn = new G4FTFPNeutronBuilder(QuasiElastic);
  AddBuilder( ftfpn );
  neu->RegisterMe(ftfpn);
  ftfpn->SetMinEnergy(minFTFP_neutron);
  auto bertn = new G4BertiniNeutronBuilder;
  AddBuilder(bertn);
  neu->RegisterMe(bertn);
  bertn->SetMinEnergy(minBERT_neutron);
  bertn->SetMaxEnergy(maxBERT_neutron);
  neu->Build();
 
  const G4ParticleDefinition* neutron = G4Neutron::Neutron();
  G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(neutron);
  if(nullptr != inel) { 
    inel->AddDataSet(new G4NeutronInelasticXS()); 
    if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
  G4HadronicProcess* capture = G4PhysListUtil::FindCaptureProcess(neutron);
  if (nullptr != capture) {
    capture->RegisterMe(new G4NeutronRadCapture());
  }
}

Referenced by CreateModels().

◆ operator=()

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

delete

◆ Others()

void G4HadronPhysicsFTFP_BERT::Others ( )

protectedvirtual

Reimplemented in G4HadronPhysicsINCLXX.

Definition at line 264 of file G4HadronPhysicsFTFP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
 
  // high energy particles
  if( param->GetMaxEnergy() > param->EnergyThresholdForHeavyHadrons() ) {
 
    // anti light ions
    G4HadronicBuilder::BuildAntiLightIonsFTFP();
 
    // hyperons
    G4HadronicBuilder::BuildHyperonsFTFP_BERT();
 
    // b-, c- baryons and mesons
    if( param->EnableBCParticles() ) {
      G4HadronicBuilder::BuildBCHadronsFTFP_BERT();
    }
 
    // light hypernuclei and anti-hypernuclei
    if ( param->EnableHyperNuclei() ) {
      G4HadronicBuilder::BuildHyperNucleiFTFP_BERT();
      G4HadronicBuilder::BuildHyperAntiNucleiFTFP_BERT();
    }
  }
}

Referenced by CreateModels().

◆ Pion()

void G4HadronPhysicsFTFP_BERT::Pion ( )

protectedvirtual

Reimplemented in G4HadronPhysicsFTF_BIC, and G4HadronPhysicsINCLXX.

Definition at line 201 of file G4HadronPhysicsFTFP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto pi = new G4PionBuilder;
  AddBuilder(pi);
  auto ftfppi = new G4FTFPPionBuilder(QuasiElastic);
  AddBuilder(ftfppi);
  pi->RegisterMe(ftfppi);
  ftfppi->SetMinEnergy(minFTFP_pion);
  auto bertpi = new G4BertiniPionBuilder;
  AddBuilder(bertpi);
  pi->RegisterMe(bertpi);
  bertpi->SetMaxEnergy(maxBERT_pion);
  pi->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() );
    }
  }
}

Referenced by CreateModels().

◆ Proton()

void G4HadronPhysicsFTFP_BERT::Proton ( )

protectedvirtual

Reimplemented in G4HadronPhysicsFTF_BIC, G4HadronPhysicsINCLXX, and G4HadronPhysicsNuBeam.

Definition at line 176 of file G4HadronPhysicsFTFP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto pro = new G4ProtonBuilder;
  AddBuilder(pro);
  auto ftfpp = new G4FTFPProtonBuilder(QuasiElastic);
  AddBuilder(ftfpp);
  pro->RegisterMe(ftfpp);
  ftfpp->SetMinEnergy(minFTFP_proton);
  auto bertp = new G4BertiniProtonBuilder;
  AddBuilder(bertp);
  pro->RegisterMe(bertp);
  bertp->SetMinEnergy(minBERT_proton);
  bertp->SetMaxEnergy(maxBERT_proton);
  pro->Build();
 
  const G4ParticleDefinition* proton = G4Proton::Proton();
  G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(proton);
  if(nullptr != inel) { 
    if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
}