G4HadronPhysicsQGS_BIC Class Reference

#include <G4HadronPhysicsQGS_BIC.hh>

Inheritance diagram for G4HadronPhysicsQGS_BIC:

Public Member Functions
	G4HadronPhysicsQGS_BIC (G4int verbose=1)
	G4HadronPhysicsQGS_BIC (const G4String &name, G4bool quasiElastic=true)
virtual	~G4HadronPhysicsQGS_BIC ()
	G4HadronPhysicsQGS_BIC (G4HadronPhysicsQGS_BIC &)=delete
G4HadronPhysicsQGS_BIC &	operator= (const G4HadronPhysicsQGS_BIC &right)=delete
Public Member Functions inherited from G4HadronPhysicsQGSP_BERT
	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
virtual void	Neutron ()
virtual void	Proton ()
virtual void	Pion ()
Protected Member Functions inherited from G4HadronPhysicsQGSP_BERT
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)

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

Detailed Description

Definition at line 44 of file G4HadronPhysicsQGS_BIC.hh.

Constructor & Destructor Documentation

G4HadronPhysicsQGS_BIC() [1/3]

G4HadronPhysicsQGS_BIC::G4HadronPhysicsQGS_BIC

(

G4int

verbose = 1

)

Definition at line 81 of file G4HadronPhysicsQGS_BIC.cc.

    : G4HadronPhysicsQGS_BIC("hInelastic QGS_BIC",true) {}

G4HadronPhysicsQGS_BIC() [2/3]

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

Definition at line 84 of file G4HadronPhysicsQGS_BIC.cc.

  : G4HadronPhysicsQGSP_BERT(name, qe) 
{
  minBERT_pion = 1.0*GeV;
  maxBIC_pion  = 1.5*GeV;
}

~G4HadronPhysicsQGS_BIC()

G4HadronPhysicsQGS_BIC::~G4HadronPhysicsQGS_BIC ( )

virtual

Definition at line 91 of file G4HadronPhysicsQGS_BIC.cc.

{}

G4HadronPhysicsQGS_BIC() [3/3]

G4HadronPhysicsQGS_BIC::G4HadronPhysicsQGS_BIC ( G4HadronPhysicsQGS_BIC &  )

delete

Member Function Documentation

Neutron()

void G4HadronPhysicsQGS_BIC::Neutron ( )

protectedvirtual

Reimplemented from G4HadronPhysicsQGSP_BERT.

Definition at line 94 of file G4HadronPhysicsQGS_BIC.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 G4QGSBinaryNeutronBuilder(QuasiElasticQGS);
  AddBuilder(qgs);
  qgs->SetMinEnergy(minQGSP_neutron);
  neu->RegisterMe(qgs);
  auto ftf = new G4FTFBinaryNeutronBuilder(QuasiElasticFTF);
  AddBuilder(ftf);
  ftf->SetMinEnergy(minFTFP_neutron);
  ftf->SetMaxEnergy(maxFTFP_neutron);
  neu->RegisterMe(ftf);
  auto bicn = new G4BinaryNeutronBuilder;
  AddBuilder(bicn);
  bicn->SetMaxEnergy(maxBERT_neutron);
  neu->RegisterMe(bicn);
  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());
  }
}

operator=()

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

delete

Pion()

void G4HadronPhysicsQGS_BIC::Pion ( )

protectedvirtual

Reimplemented from G4HadronPhysicsQGSP_BERT.

Definition at line 161 of file G4HadronPhysicsQGS_BIC.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
 
  auto pi = new G4PionBuilder;
  AddBuilder(pi);
  auto qgs = new G4QGSBinaryPionBuilder(QuasiElasticQGS);
  AddBuilder(qgs);
  qgs->SetMinEnergy(minQGSP_pik);
  pi->RegisterMe(qgs);
  auto ftf = new G4FTFBinaryPionBuilder(QuasiElasticFTF);
  AddBuilder(ftf);
  ftf->SetMinEnergy(minFTFP_pik);
  ftf->SetMaxEnergy(maxFTFP_pik);
  pi->RegisterMe(ftf);
  auto bert = new G4BertiniPionBuilder;
  AddBuilder(bert);
  bert->SetMinEnergy(minBERT_pion);
  bert->SetMaxEnergy(maxBERT_pik);
  pi->RegisterMe(bert);
  auto bic = new G4BinaryPionBuilder;
  AddBuilder(bic);
  bic->SetMaxEnergy(maxBIC_pion);
  pi->RegisterMe(bic);
  pi->Build();
 
  auto k = new G4KaonBuilder;
  AddBuilder(k);
  auto qgsk = new G4QGSBinaryKaonBuilder(QuasiElasticQGS);
  AddBuilder(qgsk);
  qgsk->SetMinEnergy(minQGSP_pik);
  k->RegisterMe(qgsk);
  auto ftfk = new G4FTFBinaryKaonBuilder(QuasiElasticFTF);
  AddBuilder(ftfk);
  ftfk->SetMaxEnergy(maxFTFP_pik);
  k->RegisterMe(ftfk);
  auto bertk = new G4BertiniKaonBuilder;
  AddBuilder(bertk);
  bertk->SetMaxEnergy(maxBERT_pik);
  k->RegisterMe(bertk);
  k->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() );
      }
    }
  }
}

Proton()

void G4HadronPhysicsQGS_BIC::Proton ( )

protectedvirtual

Reimplemented from G4HadronPhysicsQGSP_BERT.

Definition at line 132 of file G4HadronPhysicsQGS_BIC.cc.

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

---

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

• G4HadronPhysicsQGS_BIC.hh
• G4HadronPhysicsQGS_BIC.cc