G4HadronDElasticPhysics Class Reference

#include <G4HadronDElasticPhysics.hh>

Inheritance diagram for G4HadronDElasticPhysics:

Public Member Functions
	G4HadronDElasticPhysics (G4int ver=1)
virtual	~G4HadronDElasticPhysics ()
void	ConstructParticle () override
void	ConstructProcess () override
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 ()

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

Constructor & Destructor Documentation

G4HadronDElasticPhysics()

G4HadronDElasticPhysics::G4HadronDElasticPhysics ( G4int  ver = 1 )

explicit

Definition at line 88 of file G4HadronDElasticPhysics.cc.

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

~G4HadronDElasticPhysics()

G4HadronDElasticPhysics::~G4HadronDElasticPhysics ( )

virtual

Definition at line 97 of file G4HadronDElasticPhysics.cc.

{}

Member Function Documentation

ConstructParticle()

void G4HadronDElasticPhysics::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 100 of file G4HadronDElasticPhysics.cc.

{
  // G4cout << "G4HadronDElasticPhysics::ConstructParticle" << G4endl;
  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();
 
  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();
 
  //  Construct light ions
  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  
}

ConstructProcess()

void G4HadronDElasticPhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 114 of file G4HadronDElasticPhysics.cc.

{
  const G4double elimitAntiNuc = 100.1*MeV;
  if(verbose > 1) {
    G4cout << "### HadronDElasticPhysics Construct Processes " 
       << " for anti-neuclei " 
       << elimitAntiNuc/GeV << " GeV"      << G4endl;
  }
 
  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
  anuc->SetMinEnergy(elimitAntiNuc);
  G4CrossSectionElastic* anucxs = 
    new G4CrossSectionElastic(anuc->GetComponentCrossSection());
 
  G4VCrossSectionDataSet* theComponentGGHadronNucleusData = 
    new G4CrossSectionElastic( new G4ComponentGGHadronNucleusXsc );
 
  G4VCrossSectionDataSet* theComponentGGNuclNuclData = 
    new G4CrossSectionElastic( new G4ComponentGGNuclNuclXsc );
 
  G4HadronElastic* lhep0 = new G4HadronElastic();
  G4HadronElastic* lhep1 = new G4HadronElastic();
  lhep1->SetMaxEnergy(10.1*MeV);
  G4HadronElastic* lhep2 = new G4HadronElastic();
  lhep2->SetMaxEnergy(elimitAntiNuc);
 
  G4DiffuseElastic* model = nullptr;
  G4HadronElasticProcess* hel = nullptr;
 
  auto myParticleIterator=GetParticleIterator();
  myParticleIterator->reset();
  while( (*myParticleIterator)() )
  {
    G4ParticleDefinition* particle = myParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String pname = particle->GetParticleName();
    if(pname == "anti_lambda"  ||
       pname == "anti_omega-"  || 
       pname == "anti_sigma-"  || 
       pname == "anti_sigma0"  || 
       pname == "anti_sigma+"  || 
       pname == "anti_xi-"     || 
       pname == "anti_xi0"     || 
       pname == "lambda"       || 
       pname == "omega-"       || 
       pname == "sigma-"       || 
       pname == "sigma0"       ||
       pname == "sigma+"       || 
       pname == "xi-"          ||       
       pname == "xi0"
      ) {
      
      hel = new G4HadronElasticProcess();
      hel->RegisterMe(lhep0);
      hel->AddDataSet( theComponentGGHadronNucleusData );
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronDElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
 
    } else if(pname == "proton") {   
 
      hel = new G4HadronElasticProcess();
      hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
      model = new G4DiffuseElastic();
      hel->RegisterMe(lhep1);
      hel->RegisterMe(model);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronDElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
    } else if(pname == "neutron") {   
 
      hel = new G4HadronElasticProcess();
      hel->AddDataSet(new G4NeutronElasticXS());
      model = new G4DiffuseElastic();
      hel->RegisterMe(lhep1);
      hel->RegisterMe(model);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronDElasticPhysics: " 
           << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
    } else if (pname == "pi+" || pname == "pi-") { 
 
      hel = new G4HadronElasticProcess();
      hel->AddDataSet(new G4BGGPionElasticXS(particle));
      model = new G4DiffuseElastic();
      hel->RegisterMe(lhep1);
      hel->RegisterMe(model);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronDElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
    } else if(pname == "kaon-"     || 
          pname == "kaon+"     || 
          pname == "kaon0S"    || 
          pname == "kaon0L" 
          ) {
      
      hel = new G4HadronElasticProcess();
              model = new G4DiffuseElastic();
      hel->RegisterMe(lhep1);
      hel->RegisterMe(model);
      hel->AddDataSet( theComponentGGHadronNucleusData );
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
    } else if(pname == "alpha"     ||
              pname == "deuteron"  ||
              pname == "triton"    ||
              pname == "He3"
             ) {
      hel = new G4HadronElasticProcess();
      hel->AddDataSet(theComponentGGNuclNuclData);
      hel->RegisterMe(lhep0);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
 
    } else if(
       pname == "anti_neutron"   ||
       pname == "anti_proton"    || 
       pname == "anti_alpha"     ||
       pname == "anti_deuteron"  ||
       pname == "anti_triton"    ||
       pname == "anti_He3"       ) {
 
      hel = new G4HadronElasticProcess();
      hel->AddDataSet(anucxs);
      hel->RegisterMe(lhep2);
      hel->RegisterMe(anuc);
      pmanager->AddDiscreteProcess(hel);
    }
  }
  if(verbose > 1) {
    G4cout << "### HadronDElasticPhysics Construct Processes " << G4endl;
  }
}

---

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

• G4HadronDElasticPhysics.hh
• G4HadronDElasticPhysics.cc