#include <G4PenelopeIonisationCrossSection.hh>

Inheritance diagram for G4PenelopeIonisationCrossSection:

Public Member Functions
	G4PenelopeIonisationCrossSection ()
	Constructor. More...

	~G4PenelopeIonisationCrossSection ()
	Destructor. Clean all tables. More...

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

G4double	CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat)

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

void	SetVerbosityLevel (G4int vl)
	Getter/setter for the verbosity level. More...

G4int	GetVerbosityLevel ()

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 58 of file G4PenelopeIonisationCrossSection.hh.

Constructor & Destructor Documentation

G4PenelopeIonisationCrossSection::G4PenelopeIonisationCrossSection ( )

Constructor.

NOTICE: working only for e- at the moment (no interface available for e+)

Definition at line 46 of file G4PenelopeIonisationCrossSection.cc.

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

                                                                    : 
   G4VhShellCrossSection("Penelope"),shellIDTable(0),
   theCrossSectionHandler(0)
 { 
   oscManager = G4PenelopeOscillatorManager::GetOscillatorManager();
   nMaxLevels = 9;
 
   // Verbosity scale:
   // 0 = nothing
   // 1 = calculation of cross sections, file openings, sampling of atoms
   // 2 = entering in methods
   verboseLevel = 0;
 
   fLowEnergyLimit  = 10.0*eV;
   fHighEnergyLimit = 100.0*GeV;
 
   transitionManager = G4AtomicTransitionManager::Instance();
 }

G4PenelopeIonisationCrossSection::~G4PenelopeIonisationCrossSection ( )

Destructor. Clean all tables.

Definition at line 67 of file G4PenelopeIonisationCrossSection.cc.

 {
   if (theCrossSectionHandler)
     delete theCrossSectionHandler;
 }

Member Function Documentation

G4double G4PenelopeIonisationCrossSection::CrossSection	(	G4int	Z,
		G4AtomicShellEnumerator	shell,
		G4double	incidentEnergy,
		G4double	mass,
		const G4Material *	mat
	)

virtual

Purely virtual method from the base interface. Returns the cross section for the given shell in the element Z of material mat at the specified energy

Implements G4VhShellCrossSection.

Definition at line 75 of file G4PenelopeIonisationCrossSection.cc.

References python.hepunit::barn, G4PenelopeIonisationXSHandler::BuildXSTable(), G4Electron::Electron(), python.hepunit::eV, FatalException, G4cout, G4endl, G4Exception(), G4PenelopeIonisationXSHandler::GetCrossSectionTableForCouple(), G4PenelopeOscillator::GetIonisationEnergy(), G4Material::GetName(), G4PenelopeOscillatorManager::GetNumberOfZAtomsPerMolecule(), G4PenelopeOscillatorManager::GetOscillatorIonisation(), G4PenelopeOscillator::GetParentZ(), G4PenelopeOscillator::GetResonanceEnergy(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopeOscillator::GetShellFlag(), JustWarning, python.hepunit::keV, G4INCL::Math::min(), nmax, and G4AtomicTransitionManager::NumberOfShells().

Referenced by GetCrossSection().

 {
   if (verboseLevel > 1)    
     G4cout << "Entering in method G4PenelopeIonisationCrossSection::CrossSection()" << G4endl;
         
   G4double cross = 0.;
 
   //Material pointer is not available
   if (!material)
     {
       //CRASH!
       G4ExceptionDescription ed;
       ed << "The method has been called with a NULL G4Material pointer" << G4endl;
       G4Exception("G4PenelopeIonisationCrossSection::CrossSection()","em2042",
           FatalException,ed);
       return cross;
     }
 
   if (!theCrossSectionHandler)
     theCrossSectionHandler = new G4PenelopeIonisationXSHandler();
 
   theCrossSectionHandler->BuildXSTable(material,0.,G4Electron::Electron());
 
   G4int nmax = std::min(nMaxLevels,transitionManager->NumberOfShells(Z));
 
   if(G4int(shell) < nmax && 
      incidentEnergy >= fLowEnergyLimit && incidentEnergy <= fHighEnergyLimit) 
     {
       //The shells in Penelope are organized per *material*, rather than per 
       //element, so given a material one has to find the proper index for the 
       //given Z and shellID. An appropriate lookup table is used to avoid 
       //recalculation.
       G4int index = FindShellIDIndex(material,Z,shell);
 
       //Index is not available!
       if (index < 0)
     return cross;
 
       const G4PenelopeCrossSection*  theXS = 
     theCrossSectionHandler->GetCrossSectionTableForCouple(G4Electron::Electron(),
                                   material,
                                   0.);
 
       //Cross check that everything is fine:
       G4PenelopeOscillator* theOsc = oscManager->GetOscillatorIonisation(material,index);
       if (theOsc->GetParentZ() != Z || theOsc->GetShellFlag()-1 != G4int(shell))
     {
       //something went wrong!
       G4ExceptionDescription ed;
       ed << "There is something wrong here: it looks like the index is wrong" << G4endl;
       ed << "Requested: shell " << G4int(shell) << " and Z = " << Z << G4endl;
       ed << "Retrieved: " << theOsc->GetShellFlag()-1 << " and Z = " << theOsc->GetParentZ() << G4endl;   
       G4Exception("G4PenelopeIonisationCrossSection::CrossSection()","em2043",
               JustWarning,ed);
       return cross;
     }
 
 
       G4double crossPerMolecule = (theXS) ? theXS->GetShellCrossSection(index,incidentEnergy) : 0.;
 
       //Now it must be converted in cross section per atom. I need the number of 
       //atoms of the given Z per molecule.
       G4double atomsPerMolec = oscManager->GetNumberOfZAtomsPerMolecule(material,Z);
       if (atomsPerMolec)
     cross = crossPerMolecule/atomsPerMolec;
      
 
       if (verboseLevel > 0)
     {
       G4cout << "Cross section of shell " << G4int(shell) << " and Z= " << Z;
       G4cout << " of material: " << material->GetName() << " and energy = " << incidentEnergy/keV << " keV" << G4endl;
       G4cout << "--> " << cross/barn << " barn" << G4endl;
       G4cout << "Shell binding energy: " << theOsc->GetIonisationEnergy()/eV << " eV;" ;
       G4cout << " resonance energy: " << theOsc->GetResonanceEnergy()/eV << "eV" << G4endl;
        if (verboseLevel > 2)
          {
            G4cout << "Cross section per molecule: " << crossPerMolecule/barn << " barn" << G4endl;
            G4cout << "Atoms " << Z << " per molecule: " << atomsPerMolec << G4endl;
          }
     }     
     }  
   
   return cross;
 }

std::vector< G4double > G4PenelopeIonisationCrossSection::GetCrossSection	(	G4int	Z,
		G4double	incidentEnergy,
		G4double	mass,
		G4double	deltaEnergy,
		const G4Material *	mat
	)

virtual

Purely virtual method from the base interface. Returns the cross section for all levels of element Z in material mat at the given energy

Implements G4VhShellCrossSection.

Definition at line 166 of file G4PenelopeIonisationCrossSection.cc.

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

Referenced by Probabilities().

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

G4int G4PenelopeIonisationCrossSection::GetVerbosityLevel ( )

inline

Definition at line 95 of file G4PenelopeIonisationCrossSection.hh.

95 {return verboseLevel;};

std::vector< G4double > G4PenelopeIonisationCrossSection::Probabilities	(	G4int	Z,
		G4double	incidentEnergy,
		G4double	mass,
		G4double	deltaEnergy,
		const G4Material *	mat
	)

virtual

Purely virtual method from the base interface. Returns the shell ionisation probabilities for the given Z in the material mat at the specified energy.

Implements G4VhShellCrossSection.

Definition at line 182 of file G4PenelopeIonisationCrossSection.cc.

References GetCrossSection(), and n.

 {
   std::vector<G4double> vec = GetCrossSection(Z, kinEnergy,0,0,mat);
   size_t n = vec.size();
   size_t i=0;
   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;
 }

void G4PenelopeIonisationCrossSection::SetVerbosityLevel ( G4int vl )

inline

Getter/setter for the verbosity level.

Definition at line 94 of file G4PenelopeIonisationCrossSection.hh.

94 {verboseLevel = vl;};

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation