G4empCrossSection Class Reference

#include <G4empCrossSection.hh>

Inheritance diagram for G4empCrossSection:

Public Member Functions
G4double	CrossSection (G4int Z, G4AtomicShellEnumerator shell, G4double incidentEnergy, G4double mass, const G4Material *mat) override
	G4empCrossSection (const G4empCrossSection &)=delete
	G4empCrossSection (const G4String &nam="")
std::vector< G4double >	GetCrossSection (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat) override
const G4String &	GetName () const
G4empCrossSection &	operator= (const G4empCrossSection &right)=delete
std::vector< G4double >	Probabilities (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat) override
G4int	SelectRandomShell (G4int Z, G4double incidentEnergy, G4double mass, G4double deltaEnergy, const G4Material *mat)
void	SetTotalCS (G4double) override
virtual	~G4empCrossSection ()

Private Attributes
G4int	flag
G4String	name
G4OrlicLiXsModel *	orlicShellLi
G4PaulKxsModel *	paulShellK
G4double	totalCS

Detailed Description

Definition at line 49 of file G4empCrossSection.hh.

Constructor & Destructor Documentation

G4empCrossSection() [1/2]

G4empCrossSection::G4empCrossSection ( const G4String &  nam = "" )

explicit

Definition at line 40 of file G4empCrossSection.cc.

  :G4VhShellCrossSection(nam),totalCS(0.0)
{ 
 if (nam == "Empirical") 
  {
    paulShellK = new G4PaulKxsModel();
    orlicShellLi = new G4OrlicLiXsModel();
    flag=0;
  }
  else 
  { 
    G4cout << "G4empCrossSection::G4empCrossSection: " 
       << "ERROR in G4empCrossSection name; Paul+Orlic is selected." 
       << G4endl; 
    paulShellK = new G4PaulKxsModel();
    orlicShellLi = new G4OrlicLiXsModel();
    flag=0;
  }
}

References flag, G4cout, G4endl, orlicShellLi, and paulShellK.

~G4empCrossSection()

G4empCrossSection::~G4empCrossSection ( )

virtual

Definition at line 62 of file G4empCrossSection.cc.

{ 
  delete paulShellK;
  delete orlicShellLi;
}

References orlicShellLi, and paulShellK.

G4empCrossSection() [2/2]

G4empCrossSection::G4empCrossSection ( const G4empCrossSection &  )

delete

Member Function Documentation

CrossSection()

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

overridevirtual

Implements G4VhShellCrossSection.

Definition at line 107 of file G4empCrossSection.cc.

{
  G4double res = 0.0;
  G4ParticleDefinition* aProton = G4Proton::Proton();
  
  if(fKShell == shell) { 
    res = paulShellK->CalculateKCrossSection(Z, mass, incidentEnergy);
  } 
  // this check should be done in the Orlic class, that can handle only protons;
  // however this would lead up tp three checks of the mass, while here we have only one
  // moreover, at the present time,this class handles explicitly Paul and Orlic models,
  // so it can hadle the responsibility of this check too
 
  else if (mass == aProton->GetPDGMass()) {
    
    if(fL1Shell == shell) { 
      if (flag==0) res =   orlicShellLi->CalculateL1CrossSection(Z, incidentEnergy);
    } 
    else if(fL2Shell == shell) { 
      if (flag==0) res =   orlicShellLi->CalculateL2CrossSection(Z, incidentEnergy);
    } 
    else if(fL3Shell == shell) { 
      if (flag==0) res =   orlicShellLi->CalculateL3CrossSection(Z, incidentEnergy);
    } 
  }
  return res;
}

References G4PaulKxsModel::CalculateKCrossSection(), G4OrlicLiXsModel::CalculateL1CrossSection(), G4OrlicLiXsModel::CalculateL2CrossSection(), G4OrlicLiXsModel::CalculateL3CrossSection(), fKShell, fL1Shell, fL2Shell, fL3Shell, flag, G4ParticleDefinition::GetPDGMass(), orlicShellLi, paulShellK, G4Proton::Proton(), and Z.

GetCrossSection()

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

overridevirtual

Implements G4VhShellCrossSection.

Definition at line 70 of file G4empCrossSection.cc.

{
  std::vector<G4double> crossSections;
  G4ParticleDefinition* aProton = G4Proton::Proton();
 
  crossSections.push_back( paulShellK->CalculateKCrossSection(Z, mass, incidentEnergy) );
 
  // this check should be done in the Orlic class, that can handle only protons;
  // however this would lead up tp three checks of the mass, while here we have only one
  // moreover, at the present time,this class handles explicitly Paul and Orlic models,
  // so it can hadle the responsibility of this check too
 
  if (mass == aProton->GetPDGMass()) {
 
    if (flag==0)
    {
      crossSections.push_back( orlicShellLi->CalculateL1CrossSection(Z, incidentEnergy) );
      crossSections.push_back( orlicShellLi->CalculateL2CrossSection(Z, incidentEnergy) );
      crossSections.push_back( orlicShellLi->CalculateL3CrossSection(Z, incidentEnergy) );
    }
 
  }
  else {
    crossSections.push_back( 0. );
    crossSections.push_back( 0. );
    crossSections.push_back( 0. );
  }  
  return crossSections;
 
}

References G4PaulKxsModel::CalculateKCrossSection(), G4OrlicLiXsModel::CalculateL1CrossSection(), G4OrlicLiXsModel::CalculateL2CrossSection(), G4OrlicLiXsModel::CalculateL3CrossSection(), flag, G4ParticleDefinition::GetPDGMass(), orlicShellLi, paulShellK, G4Proton::Proton(), 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().

operator=()

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

delete

Probabilities()

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

overridevirtual

Implements G4VhShellCrossSection.

Definition at line 140 of file G4empCrossSection.cc.

{  
  std::vector<G4double> crossSections = GetCrossSection(Z, incidentEnergy, mass, deltaEnergy,mat);
 
  for (size_t i=0; i<crossSections.size(); i++ ) {
    
    if (totalCS) {
      crossSections[i] = crossSections[i]/totalCS;
    }
    
  }
  return crossSections;
}

References GetCrossSection(), totalCS, 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;
}

References G4UniformRand, G4VhShellCrossSection::Probabilities(), and Z.

SetTotalCS()

void G4empCrossSection::SetTotalCS ( G4double  val )

overridevirtual

Reimplemented from G4VhShellCrossSection.

Definition at line 160 of file G4empCrossSection.cc.

                                              {
 
  totalCS = val;
 
}

References totalCS.

Field Documentation

flag

G4int G4empCrossSection::flag

private

Definition at line 82 of file G4empCrossSection.hh.

Referenced by CrossSection(), G4empCrossSection(), and GetCrossSection().

name

G4String G4VhShellCrossSection::name

privateinherited

Definition at line 98 of file G4VhShellCrossSection.hh.

Referenced by source.g4viscp.G4Scene::create_scene(), G4VhShellCrossSection::GetName(), mcscore.MCParticle::printout(), and source.g4viscp.G4Scene::update_scene().

orlicShellLi

G4OrlicLiXsModel* G4empCrossSection::orlicShellLi

private

Definition at line 79 of file G4empCrossSection.hh.

Referenced by CrossSection(), G4empCrossSection(), GetCrossSection(), and ~G4empCrossSection().

paulShellK

G4PaulKxsModel* G4empCrossSection::paulShellK

private

Definition at line 78 of file G4empCrossSection.hh.

Referenced by CrossSection(), G4empCrossSection(), GetCrossSection(), and ~G4empCrossSection().

totalCS

G4double G4empCrossSection::totalCS

private

Definition at line 81 of file G4empCrossSection.hh.

Referenced by Probabilities(), and SetTotalCS().

---

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

• source/processes/electromagnetic/lowenergy/include/G4empCrossSection.hh
• source/processes/electromagnetic/lowenergy/src/G4empCrossSection.cc