#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

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 ()

virtual void	ConstructParticle ()=0

virtual void	ConstructProcess ()=0

void	SetPhysicsName (const G4String &="")

const G4String &	GetPhysicsName () const

void	SetPhysicsType (G4int)

G4int	GetPhysicsType () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

G4int	GetInstanceID () const

virtual void	TerminateWorker ()

Protected Member Functions
void	CreateModels ()

virtual void	Neutron ()

virtual void	Proton ()

virtual void	Pion ()

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

G4double	minBERT_neutron

G4double	minBERT_pik

G4double	maxBERT_proton

G4double	maxBERT_neutron

G4double	maxBERT_pik

G4bool	QuasiElasticFTF

G4bool	QuasiElasticQGS

Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel

G4String	namePhysics

G4int	typePhysics

G4ParticleTable *	theParticleTable

G4int	g4vpcInstanceID

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 51 of file G4HadronPhysicsQGSP_BERT.hh.

Constructor & Destructor Documentation

◆ G4HadronPhysicsQGSP_BERT() [1/3]

G4HadronPhysicsQGSP_BERT::G4HadronPhysicsQGSP_BERT ( G4int verbose = 1 )

Definition at line 90 of file G4HadronPhysicsQGSP_BERT.cc.

91 : G4HadronPhysicsQGSP_BERT("hInelastic QGSP_BERT",true) {}

G4HadronPhysicsQGSP_BERT

Definition: G4HadronPhysicsQGSP_BERT.hh:52

◆ G4HadronPhysicsQGSP_BERT() [2/3]

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

Definition at line 93 of file G4HadronPhysicsQGSP_BERT.cc.

    :  G4VPhysicsConstructor(name)
{
  QuasiElasticFTF= false;   // Use built-in quasi-elastic (not add-on)
  QuasiElasticQGS= true;    // For QGS, it must use it.
  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();
  minBERT_proton = minBERT_neutron = minBERT_pik = 0.0;    
  param->SetEnableBCParticles(true);
}

◆ ~G4HadronPhysicsQGSP_BERT()

G4HadronPhysicsQGSP_BERT::~G4HadronPhysicsQGSP_BERT ( )

override

Definition at line 255 of file G4HadronPhysicsQGSP_BERT.cc.

256{}

◆ 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 258 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 273 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_pik = param->GetMaxEnergyTransitionFTF_Cascade();
 
  if(G4Threading::IsMasterThread()) {
      DumpBanner();
  }
  CreateModels();
}

◆ CreateModels()

void G4HadronPhysicsQGSP_BERT::CreateModels ( )

protected

Definition at line 111 of file G4HadronPhysicsQGSP_BERT.cc.

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

Referenced by ConstructProcess().

◆ DumpBanner()

void G4HadronPhysicsQGSP_BERT::DumpBanner ( )

protectedvirtual

Definition at line 291 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4cout << G4endl;
  G4cout << " " << GetPhysicsName() << "  Thresholds: " << G4endl;
  G4cout << "    1) between BERT  and FTF/P over the interval " 
     << minFTFP_proton/GeV << " to " << maxBERT_proton/GeV << " GeV. " << G4endl;
  G4cout << "    2) between FTF/P and QGS/P over the interval " 
     << minQGSP_proton/GeV << " to " << maxFTFP_proton/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, and G4HadronPhysicsQGSP_BERT_HP.

Definition at line 119 of file G4HadronPhysicsQGSP_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 qgs = new G4QGSPNeutronBuilder(QuasiElasticQGS);
  AddBuilder(qgs);
  qgs->SetMinEnergy(minQGSP_neutron);
  neu->RegisterMe(qgs);
  auto ftf = new G4FTFPNeutronBuilder(QuasiElasticFTF);
  AddBuilder(ftf);
  ftf->SetMinEnergy(minFTFP_neutron);
  ftf->SetMaxEnergy(maxFTFP_neutron);
  neu->RegisterMe(ftf);
  auto bert = new G4BertiniNeutronBuilder;
  AddBuilder(bert);
  bert->SetMinEnergy(minBERT_neutron);
  bert->SetMaxEnergy(maxBERT_neutron);
  neu->RegisterMe(bert);
  neu->Build();
 
  const G4ParticleDefinition* neutron = G4Neutron::Neutron();
  G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(neutron);
  if(inel) { 
    inel->AddDataSet(new G4NeutronInelasticXS()); 
    if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
  G4HadronicProcess* capture = G4PhysListUtil::FindCaptureProcess(neutron);
  if (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 235 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().

◆ Pion()

void G4HadronPhysicsQGSP_BERT::Pion ( )

protectedvirtual

Reimplemented in G4HadronPhysicsQGS_BIC.

Definition at line 188 of file G4HadronPhysicsQGSP_BERT.cc.

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

Referenced by CreateModels().

◆ Proton()

void G4HadronPhysicsQGSP_BERT::Proton ( )

protectedvirtual

Reimplemented in G4HadronPhysicsQGS_BIC.

Definition at line 158 of file G4HadronPhysicsQGSP_BERT.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto pro = new G4ProtonBuilder;
  AddBuilder(pro);
  auto qgs = new G4QGSPProtonBuilder(QuasiElasticQGS);
  AddBuilder(qgs);
  qgs->SetMinEnergy(minQGSP_proton);
  pro->RegisterMe(qgs);
  auto ftf = new G4FTFPProtonBuilder(QuasiElasticFTF);
  AddBuilder(ftf);
  ftf->SetMinEnergy(minFTFP_proton);
  ftf->SetMaxEnergy(maxFTFP_proton);
  pro->RegisterMe(ftf);
  auto bert = new G4BertiniProtonBuilder;
  AddBuilder(bert);
  bert->SetMinEnergy(minBERT_proton);
  bert->SetMaxEnergy(maxBERT_proton);
  pro->RegisterMe(bert);
  pro->Build();
 
  const G4ParticleDefinition* proton = G4Proton::Proton();
  G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(proton);
  if(inel) { 
    if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
}