G4EnergyRangeManager Class Reference

#include <G4EnergyRangeManager.hh>

Public Member Functions
void	BuildPhysicsTable (const G4ParticleDefinition &)
void	Dump (G4int verbose=0)
	G4EnergyRangeManager ()
G4HadronicInteraction *	GetHadronicInteraction (const G4HadProjectile &aHadProjectile, G4Nucleus &aTargetNucleus, const G4Material *aMaterial, const G4Element *anElement) const
std::vector< G4HadronicInteraction * > &	GetHadronicInteractionList ()
void	RegisterMe (G4HadronicInteraction *a)
	~G4EnergyRangeManager ()

Private Member Functions
	G4EnergyRangeManager (const G4EnergyRangeManager &right)
G4bool	operator!= (const G4EnergyRangeManager &right) const
G4EnergyRangeManager &	operator= (const G4EnergyRangeManager &right)
G4bool	operator== (const G4EnergyRangeManager &right) const

Private Attributes
std::vector< G4HadronicInteraction * >	theHadronicInteraction
G4int	theHadronicInteractionCounter

Detailed Description

Definition at line 39 of file G4EnergyRangeManager.hh.

Constructor & Destructor Documentation

G4EnergyRangeManager() [1/2]

G4EnergyRangeManager::G4EnergyRangeManager ( )

explicit

Definition at line 39 of file G4EnergyRangeManager.cc.

 : theHadronicInteractionCounter(0)
{}

~G4EnergyRangeManager()

G4EnergyRangeManager::~G4EnergyRangeManager

(

)

Definition at line 43 of file G4EnergyRangeManager.cc.

{}

G4EnergyRangeManager() [2/2]

G4EnergyRangeManager::G4EnergyRangeManager ( const G4EnergyRangeManager &  right )

private

Member Function Documentation

BuildPhysicsTable()

void G4EnergyRangeManager::BuildPhysicsTable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 178 of file G4EnergyRangeManager.cc.

{
  for (auto & hadi : theHadronicInteraction) {
    hadi->BuildPhysicsTable( aParticleType );
  }
}

References theHadronicInteraction.

Referenced by G4HadronicProcess::BuildPhysicsTable().

Dump()

void G4EnergyRangeManager::Dump

(

G4int

verbose = 0

)

Definition at line 161 of file G4EnergyRangeManager.cc.

{
  G4cout << "G4EnergyRangeManager " << this << G4endl;
  for (G4int i = 0 ; i < theHadronicInteractionCounter; i++) {
    G4cout << "   HadronicModel " << i <<":"
           << theHadronicInteraction[i]->GetModelName() << G4endl;
    if (verbose > 0) {
      G4cout << "      Minimum Energy " 
         << theHadronicInteraction[i]->GetMinEnergy()/GeV << " [GeV], "
             << "Maximum Energy " 
         << theHadronicInteraction[i]->GetMaxEnergy()/GeV << " [GeV]"
             << G4endl;
    }
  }
}

References G4cout, G4endl, GeV, theHadronicInteraction, and theHadronicInteractionCounter.

GetHadronicInteraction()

G4HadronicInteraction * G4EnergyRangeManager::GetHadronicInteraction	(	const G4HadProjectile &	aHadProjectile,
		G4Nucleus &	aTargetNucleus,
		const G4Material *	aMaterial,
		const G4Element *	anElement
	)		const

Definition at line 59 of file G4EnergyRangeManager.cc.

{
  // VI shortcut: if only one interaction is registered skip all checks
  if(1 == theHadronicInteractionCounter) { return theHadronicInteraction[0]; }
  else if(0 == theHadronicInteractionCounter) {
    G4cout << "G4EnergyRangeManager::GetHadronicInteraction: "
       << "no models defined for a process" << G4endl;
    return nullptr;
  }
 
  G4double kineticEnergy = aHadProjectile.GetKineticEnergy();
  // For ions, get kinetic energy per nucleon
  if ( std::abs( aHadProjectile.GetDefinition()->GetBaryonNumber() ) > 1 ) {
    kineticEnergy /= static_cast< G4double >( std::abs( aHadProjectile.GetDefinition()->GetBaryonNumber() ) );
  }
 
  G4int cou = 0, memory = 0, memor2 = 0;
  G4double emi1 = 0.0, ema1 = 0.0, emi2 = 0.0, ema2 = 0.0;
 
  for (G4int i = 0; i<theHadronicInteractionCounter; ++i) {
    if ( theHadronicInteraction[i]->IsApplicable( aHadProjectile, aTargetNucleus ) ) {
      G4double low  = theHadronicInteraction[i]->GetMinEnergy( aMaterial, anElement );
      G4double high = theHadronicInteraction[i]->GetMaxEnergy( aMaterial, anElement );
      if (low <= kineticEnergy && high >= kineticEnergy) {
        ++cou;
        emi2 = emi1;
        ema2 = ema1;
        emi1 = low;
        ema1 = high;
        memor2 = memory;
        memory = i;
      }
    }
  }
 
  G4HadronicInteraction* hi = nullptr;
  switch (cou) {
    case 0:
      G4cout << "No model found out of " << theHadronicInteractionCounter << G4endl;
      for( G4int j=0; j<theHadronicInteractionCounter; ++j) {
    G4HadronicInteraction* hint=theHadronicInteraction[j];
    G4cout << "   "<< j << ".  Elow= " << hint->GetMinEnergy(aMaterial,anElement)
           <<", Ehigh= " << hint->GetMaxEnergy(aMaterial,anElement)
           <<"   " << hint->GetModelName() << G4endl;
      }
      break;
 
    case 1:
      hi = theHadronicInteraction[memory];
      break;
 
    case 2:
      if( (emi2<=emi1 && ema2>=ema1) || (emi2>=emi1 && ema2<=ema1) ) {
    G4cout << "Energy ranges of two models fully overlapping " << G4endl;
    for( G4int j=0; j<theHadronicInteractionCounter; ++j) {
      G4HadronicInteraction* hint=theHadronicInteraction[j];
      G4cout << "   "<< j << ".  Elow= " << hint->GetMinEnergy(aMaterial,anElement)
         <<", Ehigh= " << hint->GetMaxEnergy(aMaterial,anElement)
         <<"   " << hint->GetModelName() << G4endl;
    }
      } else {
    G4double rand = G4UniformRand();
    G4int mem;
    if( emi1 < emi2 ) {
      if( (ema1-kineticEnergy) < rand*(ema1-emi2) ) {
        mem = memor2;
      } else {
        mem = memory;
      }
    } else {
      if( (ema2-kineticEnergy) < rand*(ema2-emi1) ) {
        mem = memory;
      } else {
        mem = memor2;
      }
    }
    hi = theHadronicInteraction[mem];
      }
      break;
 
    default:
      G4cout << "More than two competing models for this energy" << G4endl;
      for( G4int j=0; j<theHadronicInteractionCounter; ++j) {
    G4HadronicInteraction* hint=theHadronicInteraction[j];
    G4cout << "   "<< j << ".  Elow= " << hint->GetMinEnergy(aMaterial,anElement)
           <<", Ehigh= " << hint->GetMaxEnergy(aMaterial,anElement)
           <<"   " << hint->GetModelName() << G4endl;
      }
      break;
  }
  return hi;
} 

References G4cout, G4endl, G4UniformRand, G4ParticleDefinition::GetBaryonNumber(), G4HadProjectile::GetDefinition(), G4HadProjectile::GetKineticEnergy(), G4HadronicInteraction::GetMaxEnergy(), G4HadronicInteraction::GetMinEnergy(), G4HadronicInteraction::GetModelName(), theHadronicInteraction, and theHadronicInteractionCounter.

Referenced by G4HadronicProcess::ChooseHadronicInteraction().

GetHadronicInteractionList()

std::vector< G4HadronicInteraction * > & G4EnergyRangeManager::GetHadronicInteractionList

(

)

Definition at line 156 of file G4EnergyRangeManager.cc.

{
  return theHadronicInteraction;
}

References theHadronicInteraction.

Referenced by G4HadronicProcess::GetHadronicInteractionList(), and G4HadronicProcess::GetHadronicModel().

operator!=()

G4bool G4EnergyRangeManager::operator!= ( const G4EnergyRangeManager &  right ) const

private

operator=()

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

private

operator==()

G4bool G4EnergyRangeManager::operator== ( const G4EnergyRangeManager &  right ) const

private

RegisterMe()

void G4EnergyRangeManager::RegisterMe

(

G4HadronicInteraction *

a

)

Definition at line 46 of file G4EnergyRangeManager.cc.

{
  if(!a) { return; }
  if(0 < theHadronicInteractionCounter) {
    for(G4int i=0; i<theHadronicInteractionCounter; ++i) {
      if(a == theHadronicInteraction[i]) { return; }
    }
  }
  theHadronicInteraction.push_back(a);
  ++theHadronicInteractionCounter;
}

References theHadronicInteraction, and theHadronicInteractionCounter.

Referenced by G4HadronicProcess::RegisterMe().

Field Documentation

theHadronicInteraction

std::vector<G4HadronicInteraction*> G4EnergyRangeManager::theHadronicInteraction

private

Definition at line 68 of file G4EnergyRangeManager.hh.

Referenced by BuildPhysicsTable(), Dump(), GetHadronicInteraction(), GetHadronicInteractionList(), and RegisterMe().

theHadronicInteractionCounter

G4int G4EnergyRangeManager::theHadronicInteractionCounter

private

Definition at line 67 of file G4EnergyRangeManager.hh.

Referenced by Dump(), GetHadronicInteraction(), and RegisterMe().

---

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

• source/processes/hadronic/management/include/G4EnergyRangeManager.hh
• source/processes/hadronic/management/src/G4EnergyRangeManager.cc