G4EnergyRangeManager Class Reference

#include <G4EnergyRangeManager.hh>

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

Detailed Description

Definition at line 39 of file G4EnergyRangeManager.hh.

Constructor & Destructor Documentation

G4EnergyRangeManager()

G4EnergyRangeManager::G4EnergyRangeManager ( )

explicit

Definition at line 39 of file G4EnergyRangeManager.cc.

 : theHadronicInteractionCounter(0)
{}

~G4EnergyRangeManager()

G4EnergyRangeManager::~G4EnergyRangeManager

(

)

Definition at line 43 of file G4EnergyRangeManager.cc.

{}

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 );
  }
}

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;
    }
  }
}

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;
} 

Referenced by G4HadronicProcess::ChooseHadronicInteraction().

GetHadronicInteractionList()

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

(

)

Definition at line 156 of file G4EnergyRangeManager.cc.

{
  return theHadronicInteraction;
}

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

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;
}

Referenced by G4HadronicProcess::RegisterMe().

---

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

• G4EnergyRangeManager.hh
• G4EnergyRangeManager.cc