#include <G4LivermoreIonisationCrossSection.hh>

Inheritance diagram for G4LivermoreIonisationCrossSection:

Public Member Functions
G4double	CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass=0.0, const G4Material *mat=0) override

	G4LivermoreIonisationCrossSection (const G4LivermoreIonisationCrossSection &)=delete

	G4LivermoreIonisationCrossSection (const G4String &nam="LivermorePIXE")

std::vector< G4double >	GetCrossSection (G4int Z, G4double incidentEnergy, G4double mass=0.0, G4double deltaEnergy=0.0, const G4Material *mat=0) override

const G4String &	GetName () const

void	Initialise ()

G4LivermoreIonisationCrossSection &	operator= (const G4LivermoreIonisationCrossSection &right)=delete

std::vector< G4double >	Probabilities (G4int Z, G4double incidentEnergy, G4double mass=0.0, G4double deltaEnergy=0, const G4Material *mat=0) override

G4int	SelectRandomShell (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)

virtual void	SetTotalCS (G4double)

virtual	~G4LivermoreIonisationCrossSection ()

Private Attributes
G4VCrossSectionHandler *	crossSectionHandler

G4double	fHighEnergyLimit

G4double	fLowEnergyLimit

G4String	name

const G4AtomicTransitionManager *	transitionManager

G4int	verboseLevel

Detailed Description

Definition at line 54 of file G4LivermoreIonisationCrossSection.hh.

Constructor & Destructor Documentation

◆ G4LivermoreIonisationCrossSection() [1/2]

G4LivermoreIonisationCrossSection::G4LivermoreIonisationCrossSection ( const G4String & nam = "LivermorePIXE" )

explicit

Definition at line 46 of file G4LivermoreIonisationCrossSection.cc.

                       : G4VhShellCrossSection(nam), crossSectionHandler(nullptr)
{
  fLowEnergyLimit  = 10.0*eV;
  fHighEnergyLimit = 100.0*GeV;
 
  transitionManager = G4AtomicTransitionManager::Instance();
  verboseLevel = 0;
  
  Initialise();
}

References eV, fHighEnergyLimit, fLowEnergyLimit, GeV, Initialise(), G4AtomicTransitionManager::Instance(), transitionManager, and verboseLevel.

◆ ~G4LivermoreIonisationCrossSection()

G4LivermoreIonisationCrossSection::~G4LivermoreIonisationCrossSection ( )

virtual

Definition at line 60 of file G4LivermoreIonisationCrossSection.cc.

{
  delete crossSectionHandler;
}

References crossSectionHandler.

◆ G4LivermoreIonisationCrossSection() [2/2]

G4LivermoreIonisationCrossSection::G4LivermoreIonisationCrossSection ( const G4LivermoreIonisationCrossSection & )

delete

Member Function Documentation

◆ CrossSection()

G4double G4LivermoreIonisationCrossSection::CrossSection	(	G4int	Z,
		G4AtomicShellEnumerator	shell,
		G4double	incidentEnergy,
		G4double	mass = `0.0`,
		const G4Material *	mat = `0`
	)

overridevirtual

Implements G4VhShellCrossSection.

Definition at line 88 of file G4LivermoreIonisationCrossSection.cc.

{
  G4double cross = 0.0;
  G4int n = G4int(shell);
  G4int nmax = std::min(9,transitionManager->NumberOfShells(Z));
  if(Z > 6 && Z < 93 && n < nmax && 
     kinEnergy >= fLowEnergyLimit && kinEnergy <= fHighEnergyLimit) {
    cross = crossSectionHandler->FindValue(Z, kinEnergy, n);
  }
  return cross;
}

References crossSectionHandler, fHighEnergyLimit, G4VCrossSectionHandler::FindValue(), fLowEnergyLimit, G4INCL::Math::min(), CLHEP::detail::n, G4AtomicTransitionManager::NumberOfShells(), transitionManager, and Z.

Referenced by GetCrossSection().

◆ GetCrossSection()

std::vector< G4double > G4LivermoreIonisationCrossSection::GetCrossSection	(	G4int	Z,
		G4double	incidentEnergy,
		G4double	mass = `0.0`,
		G4double	deltaEnergy = `0.0`,
		const G4Material *	mat = `0`
	)

overridevirtual

Implements G4VhShellCrossSection.

Definition at line 105 of file G4LivermoreIonisationCrossSection.cc.

{
  G4int nmax = std::min(9,transitionManager->NumberOfShells(Z));
  std::vector<G4double> vec(nmax,0.0); 
  for(G4int i=0; i<nmax; ++i) {
    vec[i] = CrossSection(Z, G4AtomicShellEnumerator(i), kinEnergy); 
  }
  return vec;
}

References CrossSection(), G4INCL::Math::min(), G4AtomicTransitionManager::NumberOfShells(), transitionManager, and Z.

Referenced by Probabilities().

◆ GetName()

const G4String & G4VhShellCrossSection::GetName ( ) const

inlineinherited

Definition at line 102 of file G4VhShellCrossSection.hh.

{
  return name;
} 

References G4VhShellCrossSection::name.

Referenced by G4UAtomicDeexcitation::InitialiseForNewRun().

◆ Initialise()

void G4LivermoreIonisationCrossSection::Initialise ( )

Definition at line 67 of file G4LivermoreIonisationCrossSection.cc.

{
  const G4int binForFluo = 20;
  G4int nbin = G4int(std::log10(fHighEnergyLimit/fLowEnergyLimit) + 0.5);
  if(nbin <= 0) { nbin = 1; }
  nbin *= binForFluo;
 
  // Data on shell ionisation x-sections
  if (crossSectionHandler) { 
    crossSectionHandler->Clear();
    delete crossSectionHandler; 
  }
 
  G4VDataSetAlgorithm* inter = new G4SemiLogInterpolation();
  crossSectionHandler = 
    new G4eCrossSectionHandler(inter,fLowEnergyLimit,fHighEnergyLimit,nbin);
  crossSectionHandler->LoadShellData("ioni/ion-ss-cs-");
  //G4cout << "!!!  G4LivermoreIonisationCrossSection::Initialise()" << G4endl;
}

References G4VCrossSectionHandler::Clear(), crossSectionHandler, fHighEnergyLimit, fLowEnergyLimit, and G4VCrossSectionHandler::LoadShellData().

Referenced by G4LivermoreIonisationCrossSection().

◆ operator=()

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

delete

◆ Probabilities()

std::vector< G4double > G4LivermoreIonisationCrossSection::Probabilities	(	G4int	Z,
		G4double	incidentEnergy,
		G4double	mass = `0.0`,
		G4double	deltaEnergy = `0`,
		const G4Material *	mat = `0`
	)

overridevirtual

Implements G4VhShellCrossSection.

Definition at line 121 of file G4LivermoreIonisationCrossSection.cc.

{
  std::vector<G4double> vec = GetCrossSection(Z, kinEnergy);
  size_t n = vec.size();
  size_t i;
  G4double sum = 0.0;
  for(i=0; i<n; ++i) { sum += vec[i]; }
  if(sum > 0.0) { 
    sum = 1.0/sum; 
    for(i=0; i<n; ++i) { vec[i] = vec[i]*sum; }
  }
  return vec;
}

References GetCrossSection(), CLHEP::detail::n, and Z.

◆ SelectRandomShell()

G4int G4VhShellCrossSection::SelectRandomShell	(	G4int	Z,
		G4double	incidentEnergy,
		G4double	mass,
		G4double	deltaEnergy,
		const G4Material *	mat
	)

inherited

Definition at line 63 of file G4VhShellCrossSection.cc.

{
  std::vector<G4double> p = 
    Probabilities(Z,incidentEnergy,mass,deltaEnergy,mat);
  G4int shell = -1;
  size_t nShells = p.size();
  G4double q = G4UniformRand();
  for (size_t i=0; i<nShells; i++) {
    
    if (p[i] >= q) {
      shell = i;
      break;
    }
    q -= p[i];
  }
  return shell;
}