G4HadronPhysicsINCLXX Class Reference

#include <G4HadronPhysicsINCLXX.hh>

Inheritance diagram for G4HadronPhysicsINCLXX:

Public Member Functions
	G4HadronPhysicsINCLXX (G4int verbose=1)
	G4HadronPhysicsINCLXX (const G4String &name, const G4bool quasiElastic=true, const G4bool neutronHP=false, const G4bool ftfp=false)
virtual	~G4HadronPhysicsINCLXX ()
void	ConstructProcess () override
void	SetQuasiElastic (G4bool value)
	G4HadronPhysicsINCLXX (G4HadronPhysicsINCLXX &)=delete
G4HadronPhysicsINCLXX &	operator= (const G4HadronPhysicsINCLXX &right)=delete
Public Member Functions inherited from G4HadronPhysicsFTFP_BERT
	G4HadronPhysicsFTFP_BERT (G4int verbose=1)
	G4HadronPhysicsFTFP_BERT (const G4String &name, G4bool quasiElastic=false)
virtual	~G4HadronPhysicsFTFP_BERT ()
void	ConstructParticle () 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	Neutron () override
void	Proton () override
void	Pion () override
void	Kaon () override
void	Others () override
Protected Member Functions inherited from G4HadronPhysicsFTFP_BERT
void	CreateModels ()
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)

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
Protected Attributes inherited from G4HadronPhysicsFTFP_BERT
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
Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Build hadronic physics using INCL++, high-energy models (QGSP or FTFP) and possibly NeutronHP.

See also

G4INCLXXProtonBuilder

G4INCLXXNeutronBuilder

G4INCLXXPionBuilder

G4IonINCLXXBuilder

Definition at line 62 of file G4HadronPhysicsINCLXX.hh.

Constructor & Destructor Documentation

G4HadronPhysicsINCLXX() [1/3]

G4HadronPhysicsINCLXX::G4HadronPhysicsINCLXX

(

G4int

verbose = 1

)

Definition at line 95 of file G4HadronPhysicsINCLXX.cc.

    : G4HadronPhysicsINCLXX("hInelastic INCLXX")
{
  G4HadronicParameters::Instance()->SetVerboseLevel(verb);
}

G4HadronPhysicsINCLXX() [2/3]

G4HadronPhysicsINCLXX::G4HadronPhysicsINCLXX	(	const G4String &	name,
		const G4bool	quasiElastic = true,
		const G4bool	neutronHP = false,
		const G4bool	ftfp = false )

Definition at line 101 of file G4HadronPhysicsINCLXX.cc.

  : G4HadronPhysicsFTFP_BERT(name, quasiElastic), 
    withNeutronHP(neutronHP),
    withFTFP(ftfp)
{
  QuasiElastic = withFTFP ? false : true;
  minBERT_neutron = withNeutronHP ? 19.9*MeV : 0.0;
}

~G4HadronPhysicsINCLXX()

G4HadronPhysicsINCLXX::~G4HadronPhysicsINCLXX ( )

virtual

Definition at line 320 of file G4HadronPhysicsINCLXX.cc.

{}

G4HadronPhysicsINCLXX() [3/3]

G4HadronPhysicsINCLXX::G4HadronPhysicsINCLXX ( G4HadronPhysicsINCLXX & )

delete

Member Function Documentation

ConstructProcess()

void G4HadronPhysicsINCLXX::ConstructProcess ( )

overridevirtual

Reimplemented from G4HadronPhysicsFTFP_BERT.

Definition at line 323 of file G4HadronPhysicsINCLXX.cc.

{
  if(G4Threading::IsMasterThread() &&
     G4HadronicParameters::Instance()->GetVerboseLevel() > 0) {
      DumpBanner();
  }
  CreateModels();
}

Kaon()

void G4HadronPhysicsINCLXX::Kaon ( )

overrideprotectedvirtual

Reimplemented from G4HadronPhysicsFTFP_BERT.

Definition at line 241 of file G4HadronPhysicsINCLXX.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto k = new G4KaonBuilder;
  AddBuilder(k);
  G4PhysicsBuilderInterface* string = nullptr;
  if(withFTFP) {
    string = new G4FTFPKaonBuilder(QuasiElastic);
  } else {
    string = new G4QGSPKaonBuilder(QuasiElastic);
  }
  string->SetMinEnergy(14.*GeV);
  AddBuilder(string);
  k->RegisterMe(string);
 
  auto bert = new G4BertiniKaonBuilder;
  bert->SetMinEnergy(0.0*GeV);
  bert->SetMaxEnergy(15.0*GeV);
  AddBuilder(bert);
  k->RegisterMe(bert);
 
  k->Build();
 
  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() );
      }
    }
  }
}

Neutron()

void G4HadronPhysicsINCLXX::Neutron ( )

overrideprotectedvirtual

Reimplemented from G4HadronPhysicsFTFP_BERT.

Definition at line 110 of file G4HadronPhysicsINCLXX.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( withNeutronHP );
  AddBuilder(neu);
  G4PhysicsBuilderInterface* string = nullptr;
  if(withFTFP) {
    string = new G4FTFPNeutronBuilder(QuasiElastic);
  } else {
    string = new G4QGSPNeutronBuilder(QuasiElastic);
  }
  string->SetMinEnergy(15.*GeV);
  AddBuilder(string);
  neu->RegisterMe(string);
 
  auto inclxxn = new G4INCLXXNeutronBuilder;
  inclxxn->SetMaxEnergy(20.*GeV);
  AddBuilder(inclxxn);
  neu->RegisterMe(inclxxn);
 
  if(withNeutronHP) {
      inclxxn->UsePreCompound(false);
      inclxxn->SetMinEnergy(minBERT_neutron);
      auto hpn = new G4NeutronPHPBuilder;
      AddBuilder(hpn);
      neu->RegisterMe(hpn);
  } else {
      inclxxn->UsePreCompound(true);
      inclxxn->SetMinPreCompoundEnergy(0.0*MeV);
      inclxxn->SetMaxPreCompoundEnergy(2.0*MeV);
      inclxxn->SetMinEnergy(1.0*MeV);
  }
 
  neu->Build();
 
  const G4ParticleDefinition* neutron = G4Neutron::Neutron();
  G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(neutron);
  if(nullptr != inel) { 
    if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
  G4HadronicProcess* capture = G4PhysListUtil::FindCaptureProcess(neutron);
  if (nullptr != capture) {
    G4NeutronRadCapture* theNeutronRadCapture = new G4NeutronRadCapture(); 
    theNeutronRadCapture->SetMinEnergy( minBERT_neutron ); 
    capture->RegisterMe( theNeutronRadCapture );
  }
  G4HadronicProcess* fission = G4PhysListUtil::FindFissionProcess(neutron);
  if (nullptr != fission) {
    G4LFission* theNeutronLEPFission = new G4LFission();
    theNeutronLEPFission->SetMinEnergy( minBERT_neutron );
    theNeutronLEPFission->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
    fission->RegisterMe( theNeutronLEPFission );
  }
}

operator=()

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

delete

Others()

void G4HadronPhysicsINCLXX::Others ( )

overrideprotectedvirtual

Reimplemented from G4HadronPhysicsFTFP_BERT.

Definition at line 279 of file G4HadronPhysicsINCLXX.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
 
  // high energy particles
  if( param->GetMaxEnergy() > param->EnergyThresholdForHeavyHadrons() ) {
 
    // anti light ions
    G4HadronicBuilder::BuildAntiLightIonsINCLXX();
 
    if ( param->EnableHyperNuclei() ) {
      // INCLXX is currently capable of handling light hypernuclei projectiles,
      // but not light anti-hypernuclei projectiles, therefore FTFP must be used
      // for the latter.
      // Note that the QGSP string model cannot currently handle nuclear projectiles
      // of any kind, so only the FTFP string model can be used together with INCLXX
      // for the simulation of nuclear interactions light hypernuclei.
      G4HadronicBuilder::BuildHyperAntiNucleiFTFP_BERT();
      G4HadronicBuilder::BuildHyperNucleiFTFP_INCLXX();
    }
 
    if(withFTFP) {
      // hyperons
      G4HadronicBuilder::BuildHyperonsFTFP_BERT();
 
      // b-, c- baryons and mesons
      if( param->EnableBCParticles() ) {
    G4HadronicBuilder::BuildBCHadronsFTFP_BERT();
      }
    } else {
      // hyperons
      G4HadronicBuilder::BuildHyperonsQGSP_FTFP_BERT(true);
 
      // b-, c- baryons and mesons
      if( param->EnableBCParticles() ) {
    G4HadronicBuilder::BuildBCHadronsQGSP_FTFP_BERT(true);
      }
    }
  }
}

Pion()

void G4HadronPhysicsINCLXX::Pion ( )

overrideprotectedvirtual

Reimplemented from G4HadronPhysicsFTFP_BERT.

Definition at line 202 of file G4HadronPhysicsINCLXX.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto pi = new G4PionBuilder;
  AddBuilder(pi);
  G4PhysicsBuilderInterface* string = nullptr;
  if(withFTFP) {
    string = new G4FTFPPionBuilder(QuasiElastic);
  } else {
    string = new G4QGSPPionBuilder(QuasiElastic);
  }
  string->SetMinEnergy(15.*GeV);
  AddBuilder(string);
  pi->RegisterMe(string);
 
  auto inclxx = new G4INCLXXPionBuilder;
  inclxx->SetMinEnergy(0.0*GeV);
  inclxx->SetMaxEnergy(20.*GeV);
  AddBuilder(inclxx);
  pi->RegisterMe(inclxx);
 
  pi->Build();
 
  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() );
    }
  }
}

Proton()

void G4HadronPhysicsINCLXX::Proton ( )

overrideprotectedvirtual

Reimplemented from G4HadronPhysicsFTFP_BERT.

Definition at line 171 of file G4HadronPhysicsINCLXX.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto pro =new G4ProtonBuilder;
  AddBuilder(pro);
  G4PhysicsBuilderInterface* string = nullptr;
  if(withFTFP) {
    string = new G4FTFPProtonBuilder(QuasiElastic);
  } else {
    string = new G4QGSPProtonBuilder(QuasiElastic);
  }
  string->SetMinEnergy(15.*GeV);
  AddBuilder(string);
  pro->RegisterMe(string);
 
  auto inclxxp = new G4INCLXXProtonBuilder;
  AddBuilder(inclxxp);
  inclxxp->SetMinEnergy(1.0*MeV);
  inclxxp->SetMaxEnergy(20.0*GeV);
  pro->RegisterMe(inclxxp);
  pro->Build();
 
  const G4ParticleDefinition* proton = G4Proton::Proton();
  G4HadronicProcess* inel = G4PhysListUtil::FindInelasticProcess(proton);
  if(nullptr != inel) { 
    if( useFactorXS ) inel->MultiplyCrossSectionBy( param->XSFactorNucleonInelastic() );
  }
}

SetQuasiElastic()

void G4HadronPhysicsINCLXX::SetQuasiElastic ( G4bool value )

inline

Definition at line 71 of file G4HadronPhysicsINCLXX.hh.

{QuasiElastic = value;}; 

---

The documentation for this class was generated from the following files:

• G4HadronPhysicsINCLXX.hh
• G4HadronPhysicsINCLXX.cc