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 std::vector< G4double >	GetCrossSection (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)=0
virtual G4double	CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat)=0
virtual std::vector< G4double >	Probabilities (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)=0
virtual void	SetTotalCS (G4double)
const G4String &	GetName () const

Detailed Description

Definition at line 55 of file G4LivermoreIonisationCrossSection.hh.

Constructor & Destructor Documentation

G4LivermoreIonisationCrossSection()

G4LivermoreIonisationCrossSection::G4LivermoreIonisationCrossSection

(

const G4String &

nam = "LivermorePIXE"

)

Definition at line 47 of file G4LivermoreIonisationCrossSection.cc.

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

~G4LivermoreIonisationCrossSection()

G4LivermoreIonisationCrossSection::~G4LivermoreIonisationCrossSection ( )

virtual

Definition at line 62 of file G4LivermoreIonisationCrossSection.cc.

{
  delete crossSectionHandler;
}

Member Function Documentation

CrossSection()

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.

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

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  )

virtual

Implements G4VhShellCrossSection.

Definition at line 108 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;
}

Referenced by Probabilities().

Initialise()

void G4LivermoreIonisationCrossSection::Initialise

(

)

Definition at line 69 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;
}

Referenced by G4LivermoreIonisationCrossSection().

Probabilities()

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.

{
  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