G4HadronHElasticPhysics Class Reference

#include <G4HadronHElasticPhysics.hh>

Inheritance diagram for G4HadronHElasticPhysics:

Public Member Functions
	G4HadronHElasticPhysics (G4int ver=1, G4bool diffraction=false)
virtual	~G4HadronHElasticPhysics ()
void	ConstructProcess () override
void	SetDiffraction (G4bool val)
	G4HadronHElasticPhysics (G4HadronHElasticPhysics &)=delete
G4HadronHElasticPhysics &	operator= (const G4HadronHElasticPhysics &right)=delete
Public Member Functions inherited from G4HadronElasticPhysics
	G4HadronElasticPhysics (G4int ver=1, const G4String &nam="hElasticWEL_CHIPS_XS")
	~G4HadronElasticPhysics () override=default
void	ConstructParticle () 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

Public Attributes
G4bool	fDiffraction

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 G4HadronHElasticPhysics.hh.

Constructor & Destructor Documentation

G4HadronHElasticPhysics() [1/2]

G4HadronHElasticPhysics::G4HadronHElasticPhysics	(	G4int	ver = 1,
		G4bool	diffraction = false )

Definition at line 84 of file G4HadronHElasticPhysics.cc.

  : G4HadronElasticPhysics(ver, "hElastic_BEST"), 
    fDiffraction(diffraction) 
{
  if (ver > 1) { 
    G4cout << "### G4HadronHElasticPhysics: " << GetPhysicsName() 
       << "  low-mass diffraction: " << fDiffraction << G4endl; 
  }
}

~G4HadronHElasticPhysics()

G4HadronHElasticPhysics::~G4HadronHElasticPhysics ( )

virtual

Definition at line 94 of file G4HadronHElasticPhysics.cc.

{}

G4HadronHElasticPhysics() [2/2]

G4HadronHElasticPhysics::G4HadronHElasticPhysics ( G4HadronHElasticPhysics & )

delete

Member Function Documentation

ConstructProcess()

void G4HadronHElasticPhysics::ConstructProcess ( )

overridevirtual

Reimplemented from G4HadronElasticPhysics.

Definition at line 96 of file G4HadronHElasticPhysics.cc.

                                               {
  G4HadronicParameters* param = G4HadronicParameters::Instance();
  G4bool useFactorXS = param->ApplyFactorXS();
  G4ParticleTable* table = G4ParticleTable::GetParticleTable();
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
 
  const G4double elimitDiffuse = 10.0;
  const G4double elimitAntiNuc = 100.0*MeV;
  const G4double delta = 0.1*MeV;
  G4double emax = std::max(param->GetMaxEnergy(), elimitAntiNuc+delta);
  if (param->GetVerboseLevel() > 1 ) {
    G4cout << "### HadronHElasticPhysics::ConstructProcess: lower energy limit for DiffuseElastic : " 
       << elimitDiffuse/GeV << " GeV" << G4endl
           << "                                             transition energy for anti-nuclei : " 
       << elimitAntiNuc/GeV << " GeV" << G4endl;
  }
  G4HadronElastic* lhep0 = new G4HadronElastic();
  G4HadronElastic* lhep1 = new G4HadronElastic();
  G4HadronElastic* lhep2 = new G4HadronElastic();
  lhep0->SetMaxEnergy(emax);
  lhep1->SetMaxEnergy(elimitDiffuse+delta);
  lhep2->SetMaxEnergy(elimitAntiNuc+delta);
 
  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
  anuc->SetMinEnergy( elimitAntiNuc );
  anuc->SetMaxEnergy(emax);
 
  auto anucxs = G4HadProcesses::ElasticXS("AntiAGlauber");
  auto xsNN = G4HadProcesses::ElasticXS("Glauber-Gribov Nucl-nucl");
 
  G4LMsdGenerator* diffGen = nullptr;
  G4DiffElasticRatio* diffRatio = nullptr;
  if( fDiffraction ) {
    diffGen = new G4LMsdGenerator("LMsdDiffraction");
    diffRatio = new G4DiffElasticRatio();
  }
 
  // Use Chips elastic model only for the hydrogen element and above an energy threshold
  G4ChipsElasticModel* chipsH = new G4ChipsElasticModel();
  chipsH->SetMinEnergy( elimitDiffuse );
  const G4ElementTable* theElementTable = G4Element::GetElementTable();
  for ( size_t i_ele = 0; i_ele < theElementTable->size(); i_ele++ ) {
    G4Element* element = (*theElementTable)[ i_ele ];
    if ( element->GetZ() > 1.0 ) chipsH->DeActivateFor( element );
  }
 
  // p
  G4ParticleDefinition* particle = G4Proton::Proton();
  G4HadronElasticProcess* hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
  G4DiffuseElastic* protonDiffuseElastic = new G4DiffuseElastic();
  protonDiffuseElastic->SetMinEnergy( elimitDiffuse );
  hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
  hel->RegisterMe( protonDiffuseElastic );
  hel->RegisterMe( lhep1 );
  if( fDiffraction) hel->SetDiffraction(diffGen, diffRatio);
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorNucleonElastic() );
  ph->RegisterProcess(hel, particle);
 
  // n
  particle = G4Neutron::Neutron();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4NeutronElasticXS());
  G4DiffuseElastic* neutronDiffuseElastic = new G4DiffuseElastic();
  neutronDiffuseElastic->SetMinEnergy( elimitDiffuse );
  hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
  hel->RegisterMe( neutronDiffuseElastic );
  hel->RegisterMe( lhep1 );
  if( fDiffraction) hel->SetDiffraction(diffGen, diffRatio);
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorNucleonElastic() );
  ph->RegisterProcess(hel, particle);
 
  // pi+
  particle = G4PionPlus::PionPlus();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGPionElasticXS(particle));
  G4DiffuseElastic* dElastic = new G4DiffuseElastic();
  dElastic->SetMinEnergy( elimitDiffuse );
  hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
  hel->RegisterMe( dElastic );
  hel->RegisterMe( lhep1 );
  if( fDiffraction) hel->SetDiffraction(diffGen, diffRatio);
  if( useFactorXS ) hel->MultiplyCrossSectionBy( param->XSFactorPionElastic() );
  ph->RegisterProcess(hel, particle);
 
  // pi-
  particle = G4PionMinus::PionMinus();
  hel = new G4HadronElasticProcess();
  hel->AddDataSet(new G4BGGPionElasticXS(particle));
  dElastic = new G4DiffuseElastic();
  dElastic->SetMinEnergy( elimitDiffuse );
  hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
  hel->RegisterMe( dElastic );
  hel->RegisterMe( lhep1 );
  if( fDiffraction) hel->SetDiffraction(diffGen, diffRatio);
  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=()

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

delete

SetDiffraction()

void G4HadronHElasticPhysics::SetDiffraction ( G4bool val )

inline

Definition at line 53 of file G4HadronHElasticPhysics.hh.

{ fDiffraction = val; };

Member Data Documentation

fDiffraction

G4bool G4HadronHElasticPhysics::fDiffraction

Definition at line 59 of file G4HadronHElasticPhysics.hh.

Referenced by ConstructProcess(), G4HadronHElasticPhysics(), and SetDiffraction().

---

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

• G4HadronHElasticPhysics.hh
• G4HadronHElasticPhysics.cc