#include <G4HadronDElasticPhysics.hh>

Inheritance diagram for G4HadronDElasticPhysics:

Public Member Functions
	G4HadronDElasticPhysics (G4int ver=1)

virtual	~G4HadronDElasticPhysics ()

void	ConstructProcess () override

	G4HadronDElasticPhysics (G4HadronDElasticPhysics &)=delete

G4HadronDElasticPhysics &	operator= (const G4HadronDElasticPhysics &right)=delete

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

	~G4HadronElasticPhysics () override=default

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

	G4HadronElasticPhysics (G4HadronElasticPhysics &)=delete

G4HadronElasticPhysics &	operator= (const G4HadronElasticPhysics &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

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)

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

Definition at line 43 of file G4HadronDElasticPhysics.hh.

Constructor & Destructor Documentation

◆ G4HadronDElasticPhysics() [1/2]

G4HadronDElasticPhysics::G4HadronDElasticPhysics ( G4int ver = 1 )

explicit

Definition at line 91 of file G4HadronDElasticPhysics.cc.

  : G4HadronElasticPhysics(ver, "hElasticDIFFUSE")
{
  if(ver > 1) { 
    G4cout << "### G4HadronDElasticPhysics: " << GetPhysicsName() 
       << G4endl; 
  }
}

◆ ~G4HadronDElasticPhysics()

G4HadronDElasticPhysics::~G4HadronDElasticPhysics ( )

virtual

Definition at line 100 of file G4HadronDElasticPhysics.cc.

101{}

◆ G4HadronDElasticPhysics() [2/2]

G4HadronDElasticPhysics::G4HadronDElasticPhysics ( G4HadronDElasticPhysics & )

delete

Member Function Documentation

◆ ConstructProcess()

void G4HadronDElasticPhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 103 of file G4HadronDElasticPhysics.cc.

{
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
  G4ParticleTable* table = G4ParticleTable::GetParticleTable();
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
  const G4double eminDElastic  = 10.*MeV;
  const G4double elimitAntiNuc = 100.*MeV;
  const G4double delta = 0.1*MeV;
  G4double emax = std::max(param->GetMaxEnergy(), elimitAntiNuc+delta);
  if(param->GetVerboseLevel() > 1) {
    G4cout << "### HadronDElasticPhysics Construct Processes " 
       << " for anti-neuclei " 
       << elimitAntiNuc/GeV << " GeV"      << G4endl;
  }
 
  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
  anuc->SetMinEnergy(elimitAntiNuc);
  anuc->SetMaxEnergy(emax);
 
  auto anucxs = G4HadProcesses::ElasticXS("AntiAGlauber");
  auto xsNN = G4HadProcesses::ElasticXS("Glauber-Gribov Nucl-nucl");
 
  G4HadronElastic* lhep0 = new G4HadronElastic();
  G4HadronElastic* lhep1 = new G4HadronElastic();
  lhep1->SetMaxEnergy(eminDElastic + delta);
  G4HadronElastic* lhep2 = new G4HadronElastic();
  lhep2->SetMaxEnergy(elimitAntiNuc);
 
  G4DiffuseElastic* model = nullptr;
 
  // p
  G4ParticleDefinition* particle = G4Proton::Proton();
  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
  model = new G4DiffuseElastic();
  model->SetMinEnergy(eminDElastic);
  hel->RegisterMe(lhep1);
  hel->RegisterMe(model);  
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorNucleonElastic() );
  ph->RegisterProcess(hel, particle);
 
  // n
  particle = G4Neutron::Neutron();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4NeutronElasticXS());
  model = new G4DiffuseElastic();
  model->SetMinEnergy(eminDElastic);
  hel->RegisterMe(lhep1);
  hel->RegisterMe(model);
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorNucleonElastic() );
  ph->RegisterProcess(hel, particle);
 
  // pi+
  particle = G4PionPlus::PionPlus();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGPionElasticXS(particle));
  model = new G4DiffuseElastic();
  model->SetMinEnergy(eminDElastic);
  hel->RegisterMe(lhep1);
  hel->RegisterMe(model);  
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorPionElastic() );
  ph->RegisterProcess(hel, particle);
 
  // pi-
  particle = G4PionMinus::PionMinus();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGPionElasticXS(particle));
  model = new G4DiffuseElastic();
  model->SetMinEnergy(eminDElastic);
  hel->RegisterMe(lhep1);
  hel->RegisterMe(model);  
  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() );
    }
  }
}

◆ operator=()

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