G4HadronElasticPhysics Class Reference

#include <G4HadronElasticPhysics.hh>

Inheritance diagram for G4HadronElasticPhysics:

Collaboration diagram for G4HadronElasticPhysics:

Public Member Functions
	G4HadronElasticPhysics (G4int ver=0)
virtual	~G4HadronElasticPhysics ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
G4HadronElastic *	GetNeutronModel ()
G4HadronicProcess *	GetNeutronProcess ()
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
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()
Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess *process, G4ParticleDefinition *particle)
G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const
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 52 of file G4HadronElasticPhysics.hh.

Constructor & Destructor Documentation

G4HadronElasticPhysics::G4HadronElasticPhysics

(

G4int

ver = 0

)

Definition at line 89 of file G4HadronElasticPhysics.cc.

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

Here is the call graph for this function:

G4HadronElasticPhysics::~G4HadronElasticPhysics ( )

virtual

Definition at line 98 of file G4HadronElasticPhysics.cc.

{}

Member Function Documentation

void G4HadronElasticPhysics::ConstructParticle ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 101 of file G4HadronElasticPhysics.cc.

{
  // G4cout << "G4HadronElasticPhysics::ConstructParticle" << G4endl;
  G4MesonConstructor pMesonConstructor;
  pMesonConstructor.ConstructParticle();

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();

  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4HadronElasticPhysics::ConstructProcess ( void  )

virtual

Implements G4VPhysicsConstructor.

Reimplemented in HadronElasticPhysicsHP, HadronElasticPhysicsHP, and HadronElasticPhysicsHP.

Definition at line 114 of file G4HadronElasticPhysics.cc.

{
  if(wasActivated) { return; }
  wasActivated = true;

  const G4double elimitPi = 1.0*GeV;
  const G4double elimitAntiNuc = 100.*MeV;
  const G4double delta = 0.1*MeV;
  if(verbose > 1) {
    G4cout << "### HadronElasticPhysics::ConstructProcess: Elimit for pi " 
       << elimitPi/GeV << " GeV" << G4endl;
    G4cout << "                                         for anti-neuclei " 
       << elimitAntiNuc/GeV << " GeV"      << G4endl;
  }

  G4AntiNuclElastic* anuc = new G4AntiNuclElastic();
  anuc->SetMinEnergy(elimitAntiNuc);
  G4CrossSectionElastic* anucxs = 
    new G4CrossSectionElastic(anuc->GetComponentCrossSection());

  G4HadronElastic* lhep0 = new G4HadronElastic();
  G4HadronElastic* lhep1 = new G4HadronElastic();
  G4HadronElastic* lhep2 = new G4HadronElastic();
  lhep1->SetMaxEnergy(elimitPi+delta);
  lhep2->SetMaxEnergy(elimitAntiNuc+delta);

  G4ChipsElasticModel* chipsp = new G4ChipsElasticModel();
  neutronModel = new G4ChipsElasticModel();

  G4ElasticHadrNucleusHE* he = new G4ElasticHadrNucleusHE(); 
  he->SetMinEnergy(elimitPi);

  G4VCrossSectionDataSet* theComponentGGHadronNucleusData = 
    new G4CrossSectionElastic( new G4ComponentGGHadronNucleusXsc );

  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_neutron" ||
       pname == "anti_omega-"  || 
       pname == "anti_sigma-"  || 
       pname == "anti_sigma+"  || 
       pname == "anti_xi-"  || 
       pname == "anti_xi0"  || 
       pname == "lambda"    || 
       pname == "omega-"    || 
       pname == "sigma-"    || 
       pname == "sigma+"    || 
       pname == "xi-"       || 
       pname == "alpha"     ||
       pname == "deuteron"  ||
       pname == "triton"   
       ) {
      
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->RegisterMe(lhep0);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if(pname == "He3") {
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      G4VCrossSectionDataSet* theComponentGGNuclNuclData = 
        new G4CrossSectionElastic(new G4ComponentGGNuclNuclXsc());
      hel->AddDataSet(theComponentGGNuclNuclData);
      hel->RegisterMe(lhep0);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if(pname == "proton") {   

      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      //hel->AddDataSet(new G4BGGNucleonElasticXS(particle));

      //      hel->AddDataSet(new G4ChipsProtonElasticXS());
      hel->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsProtonElasticXS::Default_Name()));
     
      hel->RegisterMe(chipsp);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if(pname == "neutron") {   

      neutronProcess = new G4HadronElasticProcess();
      neutronProcess->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NeutronElasticXS::Default_Name()));
      neutronProcess->RegisterMe(neutronModel);
      pmanager->AddDiscreteProcess(neutronProcess);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " 
           << neutronProcess->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if (pname == "pi+" || pname == "pi-") { 

      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet(new G4BGGPionElasticXS(particle));
      hel->RegisterMe(lhep1);
      hel->RegisterMe(he);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if(pname == "kaon-"     || 
          pname == "kaon+"     || 
          pname == "kaon0S"    || 
          pname == "kaon0L" 
          ) {
      
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      //AR-14Aug2017 : Replaced Gheisha elastic kaon cross sections with
      //               Grichine's Glauber-Gribov ones. In this way, the
      //               total (elastic + inelastic) kaon cross sections
      //               are consistent with the PDG ones.
      //               For the time being, kept Gheisha elastic as
      //               final-state model.
      hel->AddDataSet( theComponentGGHadronNucleusData );
      hel->RegisterMe(lhep0);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if(
       pname == "anti_proton"    || 
       pname == "anti_alpha"     ||
       pname == "anti_deuteron"  ||
       pname == "anti_triton"    ||
       pname == "anti_He3"       ) {

      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet(anucxs);
      hel->RegisterMe(lhep2);
      hel->RegisterMe(anuc);
      pmanager->AddDiscreteProcess(hel);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4HadronElastic * G4HadronElasticPhysics::GetNeutronModel

(

)

Definition at line 270 of file G4HadronElasticPhysics.cc.

{
  return neutronModel;
}

Here is the caller graph for this function:

G4HadronicProcess * G4HadronElasticPhysics::GetNeutronProcess

(

)

Definition at line 275 of file G4HadronElasticPhysics.cc.

{
  return neutronProcess;
}

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/source/physics_lists/constructors/hadron_elastic/include/G4HadronElasticPhysics.hh
• source/geant4.10.03.p03/source/physics_lists/constructors/hadron_elastic/src/G4HadronElasticPhysics.cc