G4HadronHElasticPhysics Class Reference

#include <G4HadronHElasticPhysics.hh>

Inheritance diagram for G4HadronHElasticPhysics:

Public Member Functions
	G4HadronHElasticPhysics (G4int ver=0)
virtual	~G4HadronHElasticPhysics ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
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)
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 45 of file G4HadronHElasticPhysics.hh.

Constructor & Destructor Documentation

G4HadronHElasticPhysics::G4HadronHElasticPhysics

(

G4int

ver = 0

)

Definition at line 81 of file G4HadronHElasticPhysics.cc.

References G4cout, G4endl, and G4VPhysicsConstructor::GetPhysicsName().

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

G4HadronHElasticPhysics::~G4HadronHElasticPhysics ( )

virtual

Definition at line 90 of file G4HadronHElasticPhysics.cc.

{}

Member Function Documentation

void G4HadronHElasticPhysics::ConstructParticle ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 93 of file G4HadronHElasticPhysics.cc.

References G4MesonConstructor::ConstructParticle(), G4IonConstructor::ConstructParticle(), and G4BaryonConstructor::ConstructParticle().

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

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();

  //  Construct light ions
  G4IonConstructor pConstructor;
  pConstructor.ConstructParticle();  
}

void G4HadronHElasticPhysics::ConstructProcess ( void  )

virtual

Implements G4VPhysicsConstructor.

Definition at line 107 of file G4HadronHElasticPhysics.cc.

References G4HadronicProcess::AddDataSet(), G4ProcessManager::AddDiscreteProcess(), aParticleIterator, G4cout, G4endl, G4AntiNuclElastic::GetComponentCrossSection(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetProcessManager(), G4VProcess::GetProcessName(), python.hepunit::GeV, python.hepunit::MeV, eplot::pname, G4HadronicProcess::RegisterMe(), G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().

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

  const G4double elimit = 1.0*GeV;
  const G4double elimitAntiNuc = 100*MeV;
  const G4double delta = 0.1*MeV;
  if(verbose > 1) {
    G4cout << "### HadronElasticPhysics::ConstructProcess: Elimit for HE " 
       << elimit/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(elimit+delta);
  lhep2->SetMaxEnergy(elimitAntiNuc+delta);

  G4ChipsElasticModel* chips = new G4ChipsElasticModel();
  chips->SetMaxEnergy(elimit+delta);

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

  aParticleIterator->reset();
  while( (*aParticleIterator)() )
  {
    G4ParticleDefinition* particle = aParticleIterator->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 == "proton") {   

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

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

      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet(new G4NeutronElasticXS());
      hel->RegisterMe(chips);
      hel->RegisterMe(he);
      pmanager->AddDiscreteProcess(hel);
      if(verbose > 1) {
    G4cout << "### HadronElasticPhysics: " 
           << hel->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();
      hel->RegisterMe(lhep0);
      //hel->AddDataSet(new G4CHIPSElasticXS());
      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);
    }
  }
}

---

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

• G4HadronHElasticPhysics.hh
• G4HadronHElasticPhysics.cc