G4HadronHElasticPhysics Class Reference

#include <G4HadronHElasticPhysics.hh>

Inheritance diagram for G4HadronHElasticPhysics:

Collaboration diagram for G4HadronHElasticPhysics:

Public Member Functions
	G4HadronHElasticPhysics (G4int ver=0, G4bool diffraction=false)
virtual	~G4HadronHElasticPhysics ()
virtual void	ConstructParticle ()
virtual void	ConstructProcess ()
void	SetDiffraction (G4bool val)
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 47 of file G4HadronHElasticPhysics.hh.

Constructor & Destructor Documentation

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

Definition at line 89 of file G4HadronHElasticPhysics.cc.

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

Here is the call graph for this function:

G4HadronHElasticPhysics::~G4HadronHElasticPhysics ( )

virtual

Definition at line 99 of file G4HadronHElasticPhysics.cc.

{}

Member Function Documentation

void G4HadronHElasticPhysics::ConstructParticle ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 102 of file G4HadronHElasticPhysics.cc.

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

  G4BaryonConstructor pBaryonConstructor;
  pBaryonConstructor.ConstructParticle();

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

Here is the call graph for this function:

void G4HadronHElasticPhysics::ConstructProcess ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 115 of file G4HadronHElasticPhysics.cc.

                                               {

  const G4double elimitDiffuse = 0.0;
  const G4double elimitAntiNuc = 100.0*MeV;
  const G4double delta = 0.1*MeV;

  if ( verbose > 1 ) {
    G4cout << "### HadronHElasticPhysics::ConstructProcess: lower energy limit for DiffuseElastic : " 
       << elimitDiffuse/GeV << " GeV" << G4endl
           << "                                             transition energy for anti-nuclei : " 
       << elimitAntiNuc/GeV << " GeV" << G4endl;
  }

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

  G4HadronElastic* lhep = new G4HadronElastic();
  lhep->SetMaxEnergy( elimitAntiNuc + delta );

  // Three instances of Chips elastic model: one used everywhere,
  // one used below a energy threshold, and one used only for the
  // hydrogen element.
  G4ChipsElasticModel* chips1 = new G4ChipsElasticModel();
  G4ChipsElasticModel* chips2 = new G4ChipsElasticModel();
  chips2->SetMaxEnergy( elimitAntiNuc + delta );
  G4ChipsElasticModel* chipsH = new G4ChipsElasticModel();
  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 );
  }

  G4NuclNuclDiffuseElastic* diffuseNuclNuclElastic = new G4NuclNuclDiffuseElastic();
  diffuseNuclNuclElastic->SetMinEnergy( elimitDiffuse );

  G4VCrossSectionDataSet* theComponentGGNuclNuclData = 
    new G4CrossSectionElastic( new G4ComponentGGNuclNuclXsc() );

  G4LMsdGenerator* diffGen = 0;
  if(fDiffraction) {
    diffGen = new G4LMsdGenerator("LMsdDiffraction");
    diffRatio = new G4DiffElasticRatio();
    G4AutoDelete::Register(diffRatio);    
  }

  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_sigma-"  ||
         pname == "anti_sigma0"  || 
         pname == "anti_sigma+"  || 
         pname == "anti_xi-"     || 
         pname == "anti_xi0"     ||
         pname == "anti_omega-"
       ) {
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4ChipsAntiBaryonElasticXS::Default_Name() ) );
      hel->RegisterMe( chips1 );
      pmanager->AddDiscreteProcess( hel );
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }
      
    } else if ( pname == "lambda"  || 
                pname == "sigma-"  ||
                pname == "sigma0"  || 
                pname == "sigma+"  || 
                pname == "xi-"     || 
                pname == "xi0"     ||
                pname == "omega-"
              ) {
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4ChipsHyperonElasticXS::Default_Name() ) );
      hel->RegisterMe( chips1 );
      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 ) );
      // To preserve reproducibility, a different instance of
      // G4DiffuseElastic must be used for each particle type.
      G4DiffuseElastic* protonDiffuseElastic = new G4DiffuseElastic();
      protonDiffuseElastic->SetMinEnergy( elimitDiffuse );
      hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
      hel->RegisterMe( protonDiffuseElastic );
      pmanager->AddDiscreteProcess( hel );
      if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( pname == "neutron" ) {   
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet(G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4NeutronElasticXS::Default_Name()) );
      // To preserve reproducibility, a different instance of
      // G4DiffuseElastic must be used for each particle type.
      G4DiffuseElastic* neutronDiffuseElastic = new G4DiffuseElastic();
      neutronDiffuseElastic->SetMinEnergy( elimitDiffuse );
      hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
      hel->RegisterMe( neutronDiffuseElastic );
      pmanager->AddDiscreteProcess( hel );
      if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " 
           << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( pname == "pi-" ) { 
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( new G4BGGPionElasticXS( particle ) );
      // To preserve reproducibility, a different instance of
      // G4DiffuseElastic must be used for each particle type.
      G4DiffuseElastic* pionMinusDiffuseElastic = new G4DiffuseElastic();
      pionMinusDiffuseElastic->SetMinEnergy( elimitDiffuse );
      hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
      hel->RegisterMe( pionMinusDiffuseElastic );
      pmanager->AddDiscreteProcess( hel );
      if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( pname == "pi+" ) { 
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( new G4BGGPionElasticXS( particle ) );
      // To preserve reproducibility, a different instance of
      // G4DiffuseElastic must be used for each particle type.
      G4DiffuseElastic* pionPlusDiffuseElastic = new G4DiffuseElastic();
      hel->RegisterMe( chipsH );  // Use Chips only for Hydrogen element
      pionPlusDiffuseElastic->SetMinEnergy( elimitDiffuse );
      hel->RegisterMe( pionPlusDiffuseElastic );
      pmanager->AddDiscreteProcess( hel );
      if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
      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();
      if ( pname == "kaon-" ) {
        hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsKaonMinusElasticXS::Default_Name() ) );
      } else if ( pname == "kaon+" ) {
        hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet(G4ChipsKaonPlusElasticXS::Default_Name() ) );
      } else {
        hel->AddDataSet( G4CrossSectionDataSetRegistry::Instance()->GetCrossSectionDataSet( G4ChipsKaonZeroElasticXS::Default_Name() ) );
      }
      hel->RegisterMe( chips1 );
      pmanager->AddDiscreteProcess( hel );
      if(fDiffraction) { hel->SetDiffraction(diffGen, diffRatio); }
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( 
                pname == "deuteron"  ||
                pname == "triton"    ||
                pname == "He3"       ||
                pname == "alpha"
              ) {
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( theComponentGGNuclNuclData );
      // To preserve reproducibility, replace temporarily 
      // G4NuclNuclDiffuseElastic with the Gheisha elastic model.
      //hel->RegisterMe( diffuseNuclNuclElastic );
      G4HadronElastic* lhepLightIon = new G4HadronElastic();     
      hel->RegisterMe( lhepLightIon );
      pmanager->AddDiscreteProcess( hel );
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( pname == "anti_proton"  ||  pname == "anti_neutron" ) {
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( anucxs );
      hel->RegisterMe( chips2 );
      hel->RegisterMe( anuc );
      pmanager->AddDiscreteProcess( hel );
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( pname == "anti_deuteron"  ||
                pname == "anti_triton"    ||
                pname == "anti_He3"       ||
                pname == "anti_alpha"
              ) {
      G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      hel->AddDataSet( anucxs );
      hel->RegisterMe( lhep );
      hel->RegisterMe( anuc );
      pmanager->AddDiscreteProcess( hel );
      if ( verbose > 1 ) {
    G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
           << " added for " << particle->GetParticleName() << G4endl;
      }

    } else if ( pname == "GenericIon" ) {
      // To preserve reproducibility, disable temporarily 
      // G4NuclNuclDiffuseElastic.
      //G4HadronElasticProcess* hel = new G4HadronElasticProcess();
      //hel->AddDataSet( theComponentGGNuclNuclData );
      //hel->RegisterMe( diffuseNuclNuclElastic );
      //pmanager->AddDiscreteProcess( hel );
      //if ( verbose > 1 ) {
      //  G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
      //         << " added for " << particle->GetParticleName() << G4endl;
      //}

    }

  }

}

Here is the call graph for this function:

void G4HadronHElasticPhysics::SetDiffraction ( G4bool  val )

inline

Definition at line 64 of file G4HadronHElasticPhysics.hh.

{fDiffraction = val;};

---

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

• geant4.10.03.p01/source/physics_lists/constructors/hadron_elastic/include/G4HadronHElasticPhysics.hh
• geant4.10.03.p01/source/physics_lists/constructors/hadron_elastic/src/G4HadronHElasticPhysics.cc