#include <G4HadronElasticPhysics.hh>

Inheritance diagram for G4HadronElasticPhysics:

Public Member Functions
	G4HadronElasticPhysics (G4int ver=0, const G4String &nam="hElasticWEL_CHIPS_XS")

virtual	~G4HadronElasticPhysics ()

void	ConstructParticle () override

void	ConstructProcess () override

G4HadronicProcess *	GetElasticProcess (const G4ParticleDefinition *part) const

G4HadronElastic *	GetElasticModel (const G4ParticleDefinition *part) const

G4HadronicProcess *	GetNeutronProcess () const

G4HadronElastic *	GetNeutronModel () const

void	AddXSection (const G4ParticleDefinition , G4VCrossSectionDataSet ) const

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 Attributes
G4int	verbose

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

Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess process, G4ParticleDefinition particle)

G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const

PhysicsBuilder_V	GetBuilders () const

void	AddBuilder (G4PhysicsBuilderInterface *bld)

Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 54 of file G4HadronElasticPhysics.hh.

Constructor & Destructor Documentation

◆ G4HadronElasticPhysics()

G4HadronElasticPhysics::G4HadronElasticPhysics	(	G4int	ver = `0`,
		const G4String &	nam = `"hElasticWEL_CHIPS_XS"`
	)

explicit

Definition at line 78 of file G4HadronElasticPhysics.cc.

  : G4VPhysicsConstructor(nam), verbose(ver)
{
  if(verbose > 1) { 
    G4cout << "### G4HadronElasticPhysics: " << GetPhysicsName() 
       << G4endl; 
  }
}

◆ ~G4HadronElasticPhysics()

G4HadronElasticPhysics::~G4HadronElasticPhysics ( )

virtual

Definition at line 87 of file G4HadronElasticPhysics.cc.

88{}

Member Function Documentation

◆ AddXSection()

void G4HadronElasticPhysics::AddXSection	(	const G4ParticleDefinition *	part,
		G4VCrossSectionDataSet *	cross
	)		const

Definition at line 237 of file G4HadronElasticPhysics.cc.

{
  G4HadronicProcess* hel = GetElasticProcess(part);
  if(hel) { hel->AddDataSet(cross); }
}

◆ ConstructParticle()

void G4HadronElasticPhysics::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 90 of file G4HadronElasticPhysics.cc.

{
  // G4cout << "G4HadronElasticPhysics::ConstructParticle" << G4endl;
  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();
 
  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();
 
  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  
}

◆ ConstructProcess()

void G4HadronElasticPhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Reimplemented in G4HadronElasticPhysicsHP, G4HadronElasticPhysicsLEND, G4HadronElasticPhysicsPHP, and G4HadronElasticPhysicsXS.

Definition at line 103 of file G4HadronElasticPhysics.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
  G4ParticleTable* table = G4ParticleTable::GetParticleTable();
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
  const G4double elimitAntiNuc = 100.*MeV;
  const G4double delta = 0.1*MeV;
  G4double emax = std::max(param->GetMaxEnergy(), elimitAntiNuc+delta);
  if(verbose > 1) {
    G4cout << "### HadronElasticPhysics::ConstructProcess: "
       << "Elimit for for anti-neuclei " << elimitAntiNuc/CLHEP::GeV << " GeV"
       << " for all hadrons Emax(GeV)= " << emax/CLHEP::GeV
           << G4endl;
  }
 
  G4HadronElastic* lhep0 = new G4HadronElastic();
  G4HadronElastic* lhep2 = new G4HadronElastic();
  lhep0->SetMaxEnergy(emax);
  lhep2->SetMaxEnergy(elimitAntiNuc+delta);
 
  G4ElasticHadrNucleusHE* he = new G4ElasticHadrNucleusHE(); 
  he->SetMaxEnergy(emax);
 
  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
  anuc->SetMinEnergy(elimitAntiNuc);
  anuc->SetMaxEnergy(emax);
 
  auto anucxs = G4HadProcesses::ElasticXS("AntiAGlauber");
  auto xsNN = G4HadProcesses::ElasticXS("Glauber-Gribov Nucl-nucl");
 
  // p
  G4ParticleDefinition* particle = G4Proton::Proton();
  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
  hel->RegisterMe(new G4ChipsElasticModel());
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorNucleonElastic() );
  ph->RegisterProcess(hel, particle);
 
  // n
  particle = G4Neutron::Neutron();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4NeutronElasticXS());
  hel->RegisterMe(new G4ChipsElasticModel());
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorNucleonElastic() );
  ph->RegisterProcess(hel, particle);
 
  // pi+
  particle = G4PionPlus::PionPlus();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGPionElasticXS(particle));
  hel->RegisterMe(he);
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorPionElastic() );
  ph->RegisterProcess(hel, particle);
 
  // pi-
  particle = G4PionMinus::PionMinus();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGPionElasticXS(particle));
  hel->RegisterMe(he);
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorPionElastic() );
  ph->RegisterProcess(hel, particle);
 
  // kaons
  G4HadronicBuilder::BuildElastic( G4HadParticles::GetKaons() );
 
  // d, t, He3, alpha
  for( auto & pdg : G4HadParticles::GetLightIons() ) {
    particle = table->FindParticle( pdg );
    if ( particle == nullptr ) { continue; }
 
    hel = new G4HadronElasticProcess();
    hel->AddDataSet(xsNN);
    hel->RegisterMe(lhep0);
    if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorHadronElastic() );
    ph->RegisterProcess(hel, particle);
  }
 
  // high energy particles
  if( emax > param->EnergyThresholdForHeavyHadrons() ) {
 
    // pbar, nbar, anti light ions
    for( auto & pdg : G4HadParticles::GetLightAntiIons() ) {
      particle = table->FindParticle( pdg );
      if ( particle == nullptr ) { continue; }
 
      hel = new G4HadronElasticProcess();
      hel->RegisterMe(lhep2);
      hel->RegisterMe(anuc);
      hel->AddDataSet(anucxs);
      if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorHadronElastic() );
      ph->RegisterProcess(hel, particle);
    }
 
    // hyperons
    G4HadronicBuilder::BuildElastic( G4HadParticles::GetHyperons() );
    G4HadronicBuilder::BuildElastic( G4HadParticles::GetAntiHyperons() );
 
    // b-, c- baryons and mesons
    if( G4HadronicParameters::Instance()->EnableBCParticles() ) {
      G4HadronicBuilder::BuildElastic( G4HadParticles::GetBCHadrons() );
    }
  }
}

Referenced by G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), and G4HadronElasticPhysicsXS::ConstructProcess().

◆ GetElasticModel()

G4HadronElastic * G4HadronElasticPhysics::GetElasticModel ( const G4ParticleDefinition * part ) const

Definition at line 216 of file G4HadronElasticPhysics.cc.

{
  G4HadronElastic* mod = nullptr;
  G4HadronicProcess* hel = GetElasticProcess(part);
  if(hel) {
    std::vector<G4HadronicInteraction*>& hi = hel->GetHadronicInteractionList();
    if( !hi.empty() ) { mod = static_cast<G4HadronElastic*>(hi[0]); }
  }
  return mod;
}

Referenced by G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), and GetNeutronModel().

◆ GetElasticProcess()

G4HadronicProcess * G4HadronElasticPhysics::GetElasticProcess ( const G4ParticleDefinition * part ) const

Definition at line 210 of file G4HadronElasticPhysics.cc.

{
  return G4PhysListUtil::FindElasticProcess(part);
}

Referenced by AddXSection(), G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), GetElasticModel(), and GetNeutronProcess().

◆ GetNeutronModel()

G4HadronElastic * G4HadronElasticPhysics::GetNeutronModel ( ) const

Definition at line 232 of file G4HadronElasticPhysics.cc.

{
  return GetElasticModel(G4Neutron::Neutron());
}

◆ GetNeutronProcess()

G4HadronicProcess * G4HadronElasticPhysics::GetNeutronProcess ( ) const

Definition at line 227 of file G4HadronElasticPhysics.cc.

{
  return GetElasticProcess(G4Neutron::Neutron());
}

Member Data Documentation

◆ verbose

G4int G4HadronElasticPhysics::verbose

protected

Definition at line 89 of file G4HadronElasticPhysics.hh.

Referenced by G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), ConstructProcess(), G4HadronElasticPhysics(), G4HadronElasticPhysicsHP::G4HadronElasticPhysicsHP(), G4HadronElasticPhysicsLEND::G4HadronElasticPhysicsLEND(), G4HadronElasticPhysicsPHP::G4HadronElasticPhysicsPHP(), and G4HadronElasticPhysicsXS::G4HadronElasticPhysicsXS().

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

Public Member Functions

Protected Attributes

Additional Inherited Members