G4LivermoreIonisationCrossSection Class Reference

#include <G4LivermoreIonisationCrossSection.hh>

Inheritance diagram for G4LivermoreIonisationCrossSection:

Public Member Functions
	G4LivermoreIonisationCrossSection (const G4String &nam="LivermorePIXE")
virtual	~G4LivermoreIonisationCrossSection ()
void	Initialise ()
G4double	CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass=0.0, const G4Material *mat=0)
std::vector< G4double >	GetCrossSection (G4int Z, G4double incidentEnergy, G4double mass=0.0, G4double deltaEnergy=0.0, const G4Material *mat=0)
std::vector< G4double >	Probabilities (G4int Z, G4double incidentEnergy, G4double mass=0.0, G4double deltaEnergy=0, const G4Material *mat=0)
Public Member Functions inherited from G4VhShellCrossSection
	G4VhShellCrossSection (const G4String &xname="")
virtual	~G4VhShellCrossSection ()
G4int	SelectRandomShell (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
virtual void	SetTotalCS (G4double)
const G4String &	GetName () const

Detailed Description

Definition at line 55 of file G4LivermoreIonisationCrossSection.hh.

Constructor & Destructor Documentation

G4LivermoreIonisationCrossSection::G4LivermoreIonisationCrossSection

(

const G4String &

nam = "LivermorePIXE"

)

Definition at line 47 of file G4LivermoreIonisationCrossSection.cc.

References python.hepunit::eV, python.hepunit::GeV, Initialise(), and G4AtomicTransitionManager::Instance().

                       : G4VhShellCrossSection(nam), crossSectionHandler(0)
{
  fLowEnergyLimit  = 10.0*eV;
  fHighEnergyLimit = 100.0*GeV;

  transitionManager = G4AtomicTransitionManager::Instance();

  verboseLevel = 0;
  
  Initialise();
}

G4LivermoreIonisationCrossSection::~G4LivermoreIonisationCrossSection ( )

virtual

Definition at line 62 of file G4LivermoreIonisationCrossSection.cc.

{
  delete crossSectionHandler;
}

Member Function Documentation

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

virtual

Implements G4VhShellCrossSection.

Definition at line 90 of file G4LivermoreIonisationCrossSection.cc.

References G4VCrossSectionHandler::FindValue(), G4INCL::Math::min(), n, nmax, and G4AtomicTransitionManager::NumberOfShells().

Referenced by GetCrossSection().

{
  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) {
    //G4cout << "Z= " << Z << "  n= " << n << " E(MeV)= " << kinEnergy/MeV << G4endl;
    cross = crossSectionHandler->FindValue(Z, kinEnergy, n);
  }
  return cross;
}

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

virtual

Implements G4VhShellCrossSection.

Definition at line 108 of file G4LivermoreIonisationCrossSection.cc.

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

Referenced by Probabilities().

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

void G4LivermoreIonisationCrossSection::Initialise

(

)

Definition at line 69 of file G4LivermoreIonisationCrossSection.cc.

References G4VCrossSectionHandler::Clear(), and G4VCrossSectionHandler::LoadShellData().

Referenced by G4LivermoreIonisationCrossSection().

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

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

virtual

Implements G4VhShellCrossSection.

Definition at line 124 of file G4LivermoreIonisationCrossSection.cc.

References GetCrossSection(), and n.

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

---

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

• G4LivermoreIonisationCrossSection.hh
• G4LivermoreIonisationCrossSection.cc