#include <G4PenelopeIonisationModel.hh>

Inheritance diagram for G4PenelopeIonisationModel:

Public Member Functions
virtual G4double	ChargeSquareRatio (const G4Track &)

G4double	ComputeCrossSectionPerAtom (const G4ParticleDefinition , const G4Element , G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

G4double	ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double, G4double, G4double, G4double, G4double) override

virtual G4double	ComputeCrossSectionPerShell (const G4ParticleDefinition *, G4int Z, G4int shellIdx, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

G4double	ComputeDEDX (const G4MaterialCutsCouple , const G4ParticleDefinition , G4double kineticEnergy, G4double cutEnergy=DBL_MAX)

G4double	ComputeDEDXPerVolume (const G4Material , const G4ParticleDefinition , G4double kineticEnergy, G4double cutEnergy) override

G4double	ComputeMeanFreePath (const G4ParticleDefinition , G4double kineticEnergy, const G4Material , G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

virtual void	CorrectionsAlongStep (const G4MaterialCutsCouple , const G4DynamicParticle , const G4double &length, G4double &eloss)

G4double	CrossSection (const G4MaterialCutsCouple , const G4ParticleDefinition , G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

G4double	CrossSectionPerVolume (const G4Material material, const G4ParticleDefinition theParticle, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy=DBL_MAX) override

G4bool	DeexcitationFlag () const

virtual void	DefineForRegion (const G4Region *)

virtual void	FillNumberOfSecondaries (G4int &numberOfTriplets, G4int &numberOfRecoil)

G4bool	ForceBuildTableFlag () const

	G4PenelopeIonisationModel (const G4ParticleDefinition *p=nullptr, const G4String &processName="PenIoni")

	G4PenelopeIonisationModel (const G4PenelopeIonisationModel &)=delete

G4VEmAngularDistribution *	GetAngularDistribution ()

virtual G4double	GetChargeSquareRatio (const G4ParticleDefinition , const G4Material , G4double kineticEnergy)

G4PhysicsTable *	GetCrossSectionTable ()

const G4Element *	GetCurrentElement () const

const G4Isotope *	GetCurrentIsotope () const

G4ElementData *	GetElementData ()

std::vector< G4EmElementSelector * > *	GetElementSelectors ()

G4VEmFluctuationModel *	GetModelOfFluctuations ()

const G4String &	GetName () const

virtual G4double	GetPartialCrossSection (const G4Material , G4int level, const G4ParticleDefinition , G4double kineticEnergy)

virtual G4double	GetParticleCharge (const G4ParticleDefinition , const G4Material , G4double kineticEnergy)

G4VEmModel *	GetTripletModel ()

G4int	GetVerbosityLevel ()

G4double	HighEnergyActivationLimit () const

G4double	HighEnergyLimit () const

void	Initialise (const G4ParticleDefinition *, const G4DataVector &) override

void	InitialiseElementSelectors (const G4ParticleDefinition *, const G4DataVector &)

virtual void	InitialiseForElement (const G4ParticleDefinition *, G4int Z)

virtual void	InitialiseForMaterial (const G4ParticleDefinition , const G4Material )

void	InitialiseLocal (const G4ParticleDefinition , G4VEmModel ) override

G4bool	IsActive (G4double kinEnergy) const

G4bool	IsLocked () const

G4bool	IsMaster () const

G4double	LowEnergyActivationLimit () const

G4double	LowEnergyLimit () const

G4bool	LPMFlag () const

G4double	MaxSecondaryKinEnergy (const G4DynamicParticle *dynParticle)

G4double	MinEnergyCut (const G4ParticleDefinition , const G4MaterialCutsCouple ) override

virtual G4double	MinPrimaryEnergy (const G4Material , const G4ParticleDefinition , G4double cut=0.0)

virtual void	ModelDescription (std::ostream &outFile) const

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

G4double	PolarAngleLimit () const

void	SampleSecondaries (std::vector< G4DynamicParticle * > , const G4MaterialCutsCouple , const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override

G4double	SecondaryThreshold () const

G4int	SelectIsotopeNumber (const G4Element *)

const G4Element *	SelectRandomAtom (const G4Material , const G4ParticleDefinition , G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

const G4Element *	SelectRandomAtom (const G4MaterialCutsCouple , const G4ParticleDefinition , G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

G4int	SelectRandomAtomNumber (const G4Material *)

const G4Element *	SelectTargetAtom (const G4MaterialCutsCouple , const G4ParticleDefinition , G4double kineticEnergy, G4double logKineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)

void	SetActivationHighEnergyLimit (G4double)

void	SetActivationLowEnergyLimit (G4double)

void	SetAngularDistribution (G4VEmAngularDistribution *)

void	SetAngularGeneratorFlag (G4bool)

void	SetCrossSectionTable (G4PhysicsTable *, G4bool isLocal)

void	SetCurrentCouple (const G4MaterialCutsCouple *)

void	SetDeexcitationFlag (G4bool val)

void	SetElementSelectors (std::vector< G4EmElementSelector * > *)

void	SetFluctuationFlag (G4bool val)

void	SetForceBuildTable (G4bool val)

void	SetHighEnergyLimit (G4double)

void	SetLocked (G4bool)

void	SetLowEnergyLimit (G4double)

void	SetLPMFlag (G4bool val)

void	SetMasterThread (G4bool val)

void	SetParticleChange (G4VParticleChange , G4VEmFluctuationModel f=nullptr)

void	SetPolarAngleLimit (G4double)

void	SetSecondaryThreshold (G4double)

void	SetTripletModel (G4VEmModel *)

virtual void	SetupForMaterial (const G4ParticleDefinition , const G4Material , G4double kineticEnergy)

void	SetUseBaseMaterials (G4bool val)

void	SetVerbosityLevel (G4int lev)

virtual void	StartTracking (G4Track *)

G4bool	UseAngularGeneratorFlag () const

G4bool	UseBaseMaterials () const

virtual G4double	Value (const G4MaterialCutsCouple , const G4ParticleDefinition , G4double kineticEnergy)

virtual	~G4PenelopeIonisationModel ()

Protected Member Functions
const G4MaterialCutsCouple *	CurrentCouple () const

G4ParticleChangeForGamma *	GetParticleChangeForGamma ()

G4ParticleChangeForLoss *	GetParticleChangeForLoss ()

virtual G4double	MaxSecondaryEnergy (const G4ParticleDefinition *, G4double kineticEnergy)

void	SetCurrentElement (const G4Element *)

Protected Attributes
size_t	basedCoupleIndex = 0

size_t	currentCoupleIndex = 0

G4ElementData *	fElementData = nullptr

const G4ParticleDefinition *	fParticle

G4ParticleChangeForLoss *	fParticleChange

G4double	inveplus

G4bool	lossFlucFlag = true

const G4Material *	pBaseMaterial = nullptr

G4double	pFactor = 1.0

G4VParticleChange *	pParticleChange = nullptr

const std::vector< G4double > *	theDensityFactor = nullptr

const std::vector< G4int > *	theDensityIdx = nullptr

G4PhysicsTable *	xSectionTable = nullptr

Private Member Functions
void	SampleFinalStateElectron (const G4Material *, G4double cutEnergy, G4double kineticEnergy)

void	SampleFinalStatePositron (const G4Material *, G4double cutEnergy, G4double kineticEnergy)

void	SetParticle (const G4ParticleDefinition *)

Private Attributes
G4VEmAngularDistribution *	anglModel = nullptr

std::vector< G4EmElementSelector * > *	elmSelectors = nullptr

G4double	eMaxActive = DBL_MAX

G4double	eMinActive = 0.0

G4VAtomDeexcitation *	fAtomDeexcitation

G4double	fCosThetaPrimary

G4double	fCosThetaSecondary

G4PenelopeIonisationXSHandler *	fCrossSectionHandler

const G4MaterialCutsCouple *	fCurrentCouple = nullptr

const G4Element *	fCurrentElement = nullptr

const G4Isotope *	fCurrentIsotope = nullptr

G4LossTableManager *	fEmManager

G4double	fEnergySecondary

G4double	fIntrinsicHighEnergyLimit

G4double	fIntrinsicLowEnergyLimit

G4bool	fIsInitialised

G4double	fKineticEnergy1

G4bool	flagDeexcitation = false

G4bool	flagForceBuildTable = false

G4bool	fLocalTable

G4VEmFluctuationModel *	flucModel = nullptr

size_t	fNBins

G4PenelopeOscillatorManager *	fOscManager

G4bool	fPIXEflag

G4int	fTargetOscillator

G4VEmModel *	fTripletModel = nullptr

G4int	fVerboseLevel

G4double	highLimit

G4bool	isLocked = false

G4bool	isMaster = true

G4bool	localElmSelectors = true

G4bool	localTable = true

G4double	lowLimit

const G4String	name

G4int	nsec = 5

G4int	nSelectors = 0

G4double	polarAngleLimit

G4double	secondaryThreshold = DBL_MAX

G4bool	theLPMflag = false

G4bool	useAngularGenerator = false

G4bool	useBaseMaterials = false

std::vector< G4double >	xsec

Detailed Description

Definition at line 65 of file G4PenelopeIonisationModel.hh.

Constructor & Destructor Documentation

◆ G4PenelopeIonisationModel() [1/2]

G4PenelopeIonisationModel::G4PenelopeIonisationModel	(	const G4ParticleDefinition *	p = `nullptr`,
		const G4String &	processName = `"PenIoni"`
	)

explicit

Definition at line 73 of file G4PenelopeIonisationModel.cc.

  :G4VEmModel(nam),fParticleChange(nullptr),fParticle(nullptr),
   fCrossSectionHandler(nullptr), 
   fAtomDeexcitation(nullptr), fKineticEnergy1(0.*eV),
   fCosThetaPrimary(1.0),fEnergySecondary(0.*eV),
   fCosThetaSecondary(0.0),fTargetOscillator(-1),
   fIsInitialised(false),fPIXEflag(false),fLocalTable(false)
{
  fIntrinsicLowEnergyLimit = 100.0*eV;
  fIntrinsicHighEnergyLimit = 100.0*GeV;
  //  SetLowEnergyLimit(fIntrinsicLowEnergyLimit);
  SetHighEnergyLimit(fIntrinsicHighEnergyLimit);
  fNBins = 200;
 
  if (part)
    SetParticle(part);
 
  //
  fOscManager = G4PenelopeOscillatorManager::GetOscillatorManager();
  //
  fVerboseLevel= 0;   
  // Verbosity scale:
  // 0 = nothing
  // 1 = warning for energy non-conservation
  // 2 = details of energy budget
  // 3 = calculation of cross sections, file openings, sampling of atoms
  // 4 = entering in methods
 
  // Atomic deexcitation model activated by default
  SetDeexcitationFlag(true);
}

References eV, fIntrinsicHighEnergyLimit, fIntrinsicLowEnergyLimit, fNBins, fOscManager, fVerboseLevel, G4PenelopeOscillatorManager::GetOscillatorManager(), GeV, G4VEmModel::SetDeexcitationFlag(), G4VEmModel::SetHighEnergyLimit(), and SetParticle().

◆ ~G4PenelopeIonisationModel()

G4PenelopeIonisationModel::~G4PenelopeIonisationModel ( )

virtual

Definition at line 108 of file G4PenelopeIonisationModel.cc.

{
  if (IsMaster() || fLocalTable)
    {
      if (fCrossSectionHandler)
    delete fCrossSectionHandler;
    }
}

References fCrossSectionHandler, fLocalTable, and G4VEmModel::IsMaster().

◆ G4PenelopeIonisationModel() [2/2]

G4PenelopeIonisationModel::G4PenelopeIonisationModel ( const G4PenelopeIonisationModel & )

delete

Member Function Documentation

◆ ChargeSquareRatio()

G4double G4VEmModel::ChargeSquareRatio ( const G4Track & track )

virtualinherited

Reimplemented in G4BraggIonGasModel, and G4BetheBlochIonGasModel.

Definition at line 374 of file G4VEmModel.cc.

{
  return GetChargeSquareRatio(track.GetParticleDefinition(), 
                              track.GetMaterial(), track.GetKineticEnergy());
}

References G4VEmModel::GetChargeSquareRatio(), G4Track::GetKineticEnergy(), G4Track::GetMaterial(), and G4Track::GetParticleDefinition().

Referenced by G4VEnergyLossProcess::PostStepGetPhysicalInteractionLength().

◆ ComputeCrossSectionPerAtom() [1/2]

G4double G4VEmModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	part,
		const G4Element *	elm,
		G4double	kinEnergy,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

inlineinherited

Definition at line 566 of file G4VEmModel.hh.

{
  SetCurrentElement(elm);
  return ComputeCrossSectionPerAtom(part,kinEnergy,elm->GetZ(),elm->GetN(),
                                    cutEnergy,maxEnergy);
}

References G4VEmModel::ComputeCrossSectionPerAtom(), G4Element::GetN(), G4Element::GetZ(), and G4VEmModel::SetCurrentElement().

◆ ComputeCrossSectionPerAtom() [2/2]

G4double G4PenelopeIonisationModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	,
		G4double	,
		G4double	,
		G4double	,
		G4double	,
		G4double
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 340 of file G4PenelopeIonisationModel.cc.

{
  G4cout << "*** G4PenelopeIonisationModel -- WARNING ***" << G4endl;
  G4cout << "Penelope Ionisation model v2008 does not calculate cross section _per atom_ " << G4endl;
  G4cout << "so the result is always zero. For physics values, please invoke " << G4endl;
  G4cout << "GetCrossSectionPerVolume() or GetMeanFreePath() via the G4EmCalculator" << G4endl;
  return 0;
}

References G4cout, and G4endl.

◆ ComputeCrossSectionPerShell()

G4double G4VEmModel::ComputeCrossSectionPerShell	(	const G4ParticleDefinition *	,
		G4int	Z,
		G4int	shellIdx,
		G4double	kinEnergy,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

virtualinherited

Definition at line 351 of file G4VEmModel.cc.

{
  return 0.0;
}

Referenced by G4EmCalculator::ComputeCrossSectionPerShell().

◆ ComputeDEDX()

G4double G4VEmModel::ComputeDEDX	(	const G4MaterialCutsCouple *	couple,
		const G4ParticleDefinition *	part,
		G4double	kineticEnergy,
		G4double	cutEnergy = `DBL_MAX`
	)

inlineinherited

Definition at line 528 of file G4VEmModel.hh.

{
  SetCurrentCouple(couple);
  return pFactor*ComputeDEDXPerVolume(pBaseMaterial,part,kinEnergy,cutEnergy);
}

References G4VEmModel::ComputeDEDXPerVolume(), G4VEmModel::pBaseMaterial, G4VEmModel::pFactor, and G4VEmModel::SetCurrentCouple().

◆ ComputeDEDXPerVolume()

G4double G4PenelopeIonisationModel::ComputeDEDXPerVolume	(	const G4Material *	material,
		const G4ParticleDefinition *	theParticle,
		G4double	kineticEnergy,
		G4double	cutEnergy
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 356 of file G4PenelopeIonisationModel.cc.

{
  // Penelope model v2008 to calculate the stopping power for soft inelastic collisions
  // below the threshold. It makes use of the Generalised Oscillator Strength (GOS)
  // model from
  //  D. Liljequist, J. Phys. D: Appl. Phys. 16 (1983) 1567
  //
  // The stopping power is calculated analytically using the dsigma/dW cross
  // section from the GOS models, which includes separate contributions from
  // distant longitudinal collisions, distant transverse collisions and
  // close collisions. Only the atomic oscillators that come in the play
  // (according to the threshold) are considered for the calculation.
  // Differential cross sections have a different form for e+ and e-.
  //
  // Fermi density correction is calculated analytically according to
  //  U. Fano, Ann. Rev. Nucl. Sci. 13 (1963),1
 
  if (fVerboseLevel > 3)
    G4cout << "Calling ComputeDEDX() of G4PenelopeIonisationModel" << G4endl;
 
  //Either Initialize() was not called, or we are in a slave and InitializeLocal() was 
  //not invoked
  if (!fCrossSectionHandler)
    {
      //create a **thread-local** version of the table. Used only for G4EmCalculator and 
      //Unit Tests
      fLocalTable = true;
      fCrossSectionHandler = new G4PenelopeIonisationXSHandler(fNBins); 
    }
 
  const G4PenelopeCrossSection* theXS = 
    fCrossSectionHandler->GetCrossSectionTableForCouple(theParticle,material,
                              cutEnergy);
  if (!theXS)
    {
      //If we are here, it means that Initialize() was inkoved, but the MaterialTable was 
      //not filled up. This can happen in a UnitTest or via G4EmCalculator
      if (fVerboseLevel > 0)
    {   
      //Issue a G4Exception (warning) only in verbose mode
      G4ExceptionDescription ed;
      ed << "Unable to retrieve the cross section table for " << 
        theParticle->GetParticleName() <<
        " in " << material->GetName() << ", cut = " << cutEnergy/keV << " keV " << G4endl;
      ed << "This can happen only in Unit Tests or via G4EmCalculator" << G4endl;
      G4Exception("G4PenelopeIonisationModel::ComputeDEDXPerVolume()",
              "em2038",JustWarning,ed);
    }
      //protect file reading via autolock
      G4AutoLock lock(&PenelopeIonisationModelMutex);
      fCrossSectionHandler->BuildXSTable(material,cutEnergy,theParticle);
      lock.unlock();
      //now it should be ok
      theXS = 
    fCrossSectionHandler->GetCrossSectionTableForCouple(theParticle,
                              material,
                              cutEnergy);
    }
 
  G4double sPowerPerMolecule = 0.0;
  if (theXS)
    sPowerPerMolecule = theXS->GetSoftStoppingPower(kineticEnergy);
 
  G4double atomDensity = material->GetTotNbOfAtomsPerVolume();
  G4double atPerMol =  fOscManager->GetAtomsPerMolecule(material);
                                                                                        
  G4double moleculeDensity = 0.; 
  if (atPerMol)
    moleculeDensity = atomDensity/atPerMol;
  G4double sPowerPerVolume = sPowerPerMolecule*moleculeDensity;
 
  if (fVerboseLevel > 2)
    {
      G4cout << "G4PenelopeIonisationModel " << G4endl;
      G4cout << "Stopping power < " << cutEnergy/keV << " keV at " <<
        kineticEnergy/keV << " keV = " << 
    sPowerPerVolume/(keV/mm) << " keV/mm" << G4endl;
    }
  return sPowerPerVolume;
}

References G4PenelopeIonisationXSHandler::BuildXSTable(), fCrossSectionHandler, fLocalTable, fNBins, fOscManager, fVerboseLevel, G4cout, G4endl, G4Exception(), G4PenelopeOscillatorManager::GetAtomsPerMolecule(), G4PenelopeIonisationXSHandler::GetCrossSectionTableForCouple(), G4ParticleDefinition::GetParticleName(), G4PenelopeCrossSection::GetSoftStoppingPower(), JustWarning, keV, eplot::material, mm, anonymous_namespace{G4PenelopeIonisationModel.cc}::PenelopeIonisationModelMutex, and G4TemplateAutoLock< _Mutex_t >::unlock().

◆ ComputeMeanFreePath()

G4double G4VEmModel::ComputeMeanFreePath	(	const G4ParticleDefinition *	part,
		G4double	kineticEnergy,
		const G4Material *	material,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

inlineinherited

Definition at line 553 of file G4VEmModel.hh.

{
  G4double cross = CrossSectionPerVolume(material,part,ekin,emin,emax);
  return (cross > 0.0) ? 1./cross : DBL_MAX;
}

References G4VEmModel::CrossSectionPerVolume(), DBL_MAX, emax, and eplot::material.

◆ CorrectionsAlongStep()

void G4VEmModel::CorrectionsAlongStep	(	const G4MaterialCutsCouple *	,
		const G4DynamicParticle *	,
		const G4double &	length,
		G4double &	eloss
	)

virtualinherited

Reimplemented in G4IonParametrisedLossModel, G4AtimaEnergyLossModel, G4BraggIonModel, G4ICRU73QOModel, G4BetheBlochModel, and G4LindhardSorensenIonModel.

Definition at line 399 of file G4VEmModel.cc.

402{}

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4VEnergyLossProcess::AlongStepDoIt(), G4EmCalculator::ComputeDEDX(), and G4EmCalculator::GetDEDX().

◆ CrossSection()

G4double G4VEmModel::CrossSection	(	const G4MaterialCutsCouple *	couple,
		const G4ParticleDefinition *	part,
		G4double	kineticEnergy,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

inlineinherited

Definition at line 539 of file G4VEmModel.hh.

{
  SetCurrentCouple(couple);
  return pFactor*CrossSectionPerVolume(pBaseMaterial,part,kinEnergy,
                                       cutEnergy,maxEnergy);
}

References G4VEmModel::CrossSectionPerVolume(), G4VEmModel::pBaseMaterial, G4VEmModel::pFactor, and G4VEmModel::SetCurrentCouple().

Referenced by G4PolarizedAnnihilation::ComputeAsymmetry(), G4PolarizedIonisation::ComputeAsymmetry(), G4PolarizedCompton::ComputeAsymmetry(), G4EmModelManager::FillLambdaVector(), and G4EmMultiModel::SampleSecondaries().

◆ CrossSectionPerVolume()

G4double G4PenelopeIonisationModel::CrossSectionPerVolume	(	const G4Material *	material,
		const G4ParticleDefinition *	theParticle,
		G4double	kineticEnergy,
		G4double	cutEnergy,
		G4double	maxEnergy = `DBL_MAX`
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 236 of file G4PenelopeIonisationModel.cc.

{
  // Penelope model v2008 to calculate the cross section for inelastic collisions above the
  // threshold. It makes use of the Generalised Oscillator Strength (GOS) model from
  //  D. Liljequist, J. Phys. D: Appl. Phys. 16 (1983) 1567
  //
  // The total cross section is calculated analytically by taking
  // into account the atomic oscillators coming into the play for a given threshold.
  //
  // For incident e- the maximum energy allowed for the delta-rays is energy/2.
  // because particles are undistinghishable.
  //
  // The contribution is splitted in three parts: distant longitudinal collisions,
  // distant transverse collisions and close collisions. Each term is described by
  // its own analytical function.
  // Fermi density correction is calculated analytically according to
  //  U. Fano, Ann. Rev. Nucl. Sci. 13 (1963),1
  //
  if (fVerboseLevel > 3)
    G4cout << "Calling CrossSectionPerVolume() of G4PenelopeIonisationModel" << G4endl;
 
  SetupForMaterial(theParticle, material, energy);
   
  G4double totalCross = 0.0;
  G4double crossPerMolecule = 0.;
 
  //Either Initialize() was not called, or we are in a slave and InitializeLocal() was 
  //not invoked
  if (!fCrossSectionHandler)
    {
      //create a **thread-local** version of the table. Used only for G4EmCalculator and 
      //Unit Tests
      fLocalTable = true;
      fCrossSectionHandler = new G4PenelopeIonisationXSHandler(fNBins); 
    }
  
  const G4PenelopeCrossSection* theXS = 
    fCrossSectionHandler->GetCrossSectionTableForCouple(theParticle,
                              material,
                              cutEnergy);
  if (!theXS)
    {
      //If we are here, it means that Initialize() was inkoved, but the MaterialTable was 
      //not filled up. This can happen in a UnitTest or via G4EmCalculator
      if (fVerboseLevel > 0)
    {
      //Issue a G4Exception (warning) only in verbose mode
      G4ExceptionDescription ed;
      ed << "Unable to retrieve the cross section table for " << 
        theParticle->GetParticleName() << 
        " in " << material->GetName() << ", cut = " << cutEnergy/keV << " keV " << G4endl;
      ed << "This can happen only in Unit Tests or via G4EmCalculator" << G4endl;
      G4Exception("G4PenelopeIonisationModel::CrossSectionPerVolume()",
              "em2038",JustWarning,ed);
    }
      //protect file reading via autolock
      G4AutoLock lock(&PenelopeIonisationModelMutex);
      fCrossSectionHandler->BuildXSTable(material,cutEnergy,theParticle);
      lock.unlock();
      //now it should be ok
      theXS = 
    fCrossSectionHandler->GetCrossSectionTableForCouple(theParticle,
                              material,
                              cutEnergy);
    }
 
  if (theXS)
    crossPerMolecule = theXS->GetHardCrossSection(energy);
 
  G4double atomDensity = material->GetTotNbOfAtomsPerVolume();
  G4double atPerMol =  fOscManager->GetAtomsPerMolecule(material);
 
  if (fVerboseLevel > 3)
    G4cout << "Material " << material->GetName() << " has " << atPerMol <<
      "atoms per molecule" << G4endl;
 
  G4double moleculeDensity = 0.; 
  if (atPerMol)
    moleculeDensity = atomDensity/atPerMol;
  G4double crossPerVolume = crossPerMolecule*moleculeDensity;
 
  if (fVerboseLevel > 2)
  {
    G4cout << "G4PenelopeIonisationModel " << G4endl;
    G4cout << "Mean free path for delta emission > " << cutEnergy/keV << " keV at " <<
      energy/keV << " keV = " << (1./crossPerVolume)/mm << " mm" << G4endl;
    if (theXS)
      totalCross = (theXS->GetTotalCrossSection(energy))*moleculeDensity;
    G4cout << "Total free path for ionisation (no threshold) at " <<
      energy/keV << " keV = " << (1./totalCross)/mm << " mm" << G4endl;
  }
  return crossPerVolume;
}

References G4PenelopeIonisationXSHandler::BuildXSTable(), G4INCL::KinematicsUtils::energy(), fCrossSectionHandler, fLocalTable, fNBins, fOscManager, fVerboseLevel, G4cout, G4endl, G4Exception(), G4PenelopeOscillatorManager::GetAtomsPerMolecule(), G4PenelopeIonisationXSHandler::GetCrossSectionTableForCouple(), G4PenelopeCrossSection::GetHardCrossSection(), G4ParticleDefinition::GetParticleName(), G4PenelopeCrossSection::GetTotalCrossSection(), JustWarning, keV, eplot::material, mm, anonymous_namespace{G4PenelopeIonisationModel.cc}::PenelopeIonisationModelMutex, G4VEmModel::SetupForMaterial(), and G4TemplateAutoLock< _Mutex_t >::unlock().

◆ CurrentCouple()

const G4MaterialCutsCouple * G4VEmModel::CurrentCouple ( ) const

inlineprotectedinherited

Definition at line 490 of file G4VEmModel.hh.

{
  return fCurrentCouple;
}

References G4VEmModel::fCurrentCouple.

◆ DeexcitationFlag()

G4bool G4VEmModel::DeexcitationFlag ( ) const

inlineinherited

Definition at line 704 of file G4VEmModel.hh.

{
  return flagDeexcitation;
}

References G4VEmModel::flagDeexcitation.

Referenced by G4EmModelManager::DumpModelList().

◆ DefineForRegion()

void G4VEmModel::DefineForRegion ( const G4Region * )

virtualinherited

Reimplemented in G4PAIModel, and G4PAIPhotModel.

Definition at line 360 of file G4VEmModel.cc.

361{}

Referenced by G4EmModelManager::AddEmModel().

◆ FillNumberOfSecondaries()

void G4VEmModel::FillNumberOfSecondaries	(	G4int &	numberOfTriplets,
		G4int &	numberOfRecoil
	)

virtualinherited

Definition at line 365 of file G4VEmModel.cc.

{
  numberOfTriplets = 0;
  numberOfRecoil = 0;
}

Referenced by G4VEmProcess::PostStepDoIt(), and G4VEnergyLossProcess::PostStepDoIt().

◆ ForceBuildTableFlag()

G4bool G4VEmModel::ForceBuildTableFlag ( ) const

inlineinherited

Definition at line 711 of file G4VEmModel.hh.

{
  return flagForceBuildTable;
}

References G4VEmModel::flagForceBuildTable.

Referenced by G4VMscModel::GetParticleChangeForMSC().

◆ GetAngularDistribution()

G4VEmAngularDistribution * G4VEmModel::GetAngularDistribution ( )

inlineinherited

Definition at line 621 of file G4VEmModel.hh.

{
  return anglModel;
}

References G4VEmModel::anglModel.

◆ GetChargeSquareRatio()

G4double G4VEmModel::GetChargeSquareRatio	(	const G4ParticleDefinition *	p,
		const G4Material *	,
		G4double	kineticEnergy
	)

virtualinherited

Reimplemented in G4BraggIonModel, G4BraggModel, G4IonParametrisedLossModel, G4AtimaEnergyLossModel, G4BetheBlochModel, and G4LindhardSorensenIonModel.

Definition at line 382 of file G4VEmModel.cc.

{
  const G4double q = p->GetPDGCharge()*inveplus;
  return q*q;
}

References G4ParticleDefinition::GetPDGCharge(), and G4VEmModel::inveplus.

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4AdjointCSManager::BuildTotalSigmaTables(), G4VEmModel::ChargeSquareRatio(), G4EmCalculator::ComputeDEDX(), G4AdjointIonIonisationModel::CorrectPostStepWeight(), G4ContinuousGainOfEnergy::GetContinuousStepLimit(), and G4EmCalculator::GetDEDX().

◆ GetCrossSectionTable()

G4PhysicsTable * G4VEmModel::GetCrossSectionTable ( )

inlineinherited

Definition at line 870 of file G4VEmModel.hh.

{
  return xSectionTable;
}

References G4VEmModel::xSectionTable.

Referenced by G4VMultipleScattering::BuildPhysicsTable(), G4EmModelManager::DumpModelList(), G4EmCalculator::FindLambdaTable(), G4WentzelVIModel::Initialise(), and G4VMultipleScattering::StorePhysicsTable().

◆ GetCurrentElement()

const G4Element * G4VEmModel::GetCurrentElement ( ) const

inlineinherited

Definition at line 505 of file G4VEmModel.hh.

{
  return fCurrentElement;
}

References G4VEmModel::fCurrentElement.

Referenced by G4VEmProcess::GetCurrentElement(), G4VEnergyLossProcess::GetCurrentElement(), G4VEmProcess::GetTargetElement(), G4PolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4PolarizedPhotoElectricModel::SampleSecondaries(), G4PolarizedBremsstrahlungModel::SelectedAtom(), G4PolarizedGammaConversionModel::SelectedAtom(), and G4eBremParametrizedModel::SetCurrentElement().

◆ GetCurrentIsotope()

const G4Isotope * G4VEmModel::GetCurrentIsotope ( ) const

inlineinherited

Definition at line 512 of file G4VEmModel.hh.

{
  return fCurrentIsotope;
}

References G4VEmModel::fCurrentIsotope.

Referenced by G4VEmProcess::GetTargetIsotope().

◆ GetElementData()

G4ElementData * G4VEmModel::GetElementData ( )

inlineinherited

Definition at line 863 of file G4VEmModel.hh.

{
  return fElementData;
}

References G4VEmModel::fElementData.

Referenced by G4MuPairProductionModel::InitialiseLocal(), G4ePairProduction::StreamProcessInfo(), and G4MuPairProduction::StreamProcessInfo().

◆ GetElementSelectors()

std::vector< G4EmElementSelector * > * G4VEmModel::GetElementSelectors ( )

inlineinherited

Definition at line 844 of file G4VEmModel.hh.

{
  return elmSelectors;
}

References G4VEmModel::elmSelectors.

◆ GetModelOfFluctuations()

G4VEmFluctuationModel * G4VEmModel::GetModelOfFluctuations ( )

inlineinherited

Definition at line 614 of file G4VEmModel.hh.

{
  return flucModel;
}

References G4VEmModel::flucModel.

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4VEnergyLossProcess::AlongStepDoIt(), G4IonParametrisedLossModel::CorrectionsAlongStep(), G4AtimaEnergyLossModel::CorrectionsAlongStep(), G4BraggIonModel::CorrectionsAlongStep(), G4BetheBlochModel::CorrectionsAlongStep(), G4LindhardSorensenIonModel::CorrectionsAlongStep(), G4BraggModel::GetChargeSquareRatio(), G4VEnergyLossProcess::GetDEDXDispersion(), and G4EmMultiModel::Initialise().

◆ GetName()

const G4String & G4VEmModel::GetName ( ) const

inlineinherited

Definition at line 837 of file G4VEmModel.hh.

{
  return name;
}

References G4VEmModel::name.

Referenced by G4NIELCalculator::AddEmModel(), G4VLEPTSModel::BuildMeanFreePathTable(), G4VLEPTSModel::BuildPhysicsTable(), G4EmCalculator::ComputeDEDX(), G4MuPairProductionModel::DataCorrupted(), G4EmElementSelector::Dump(), G4EmModelManager::DumpModelList(), G4EmCalculator::FindEmModel(), G4NIELCalculator::G4NIELCalculator(), G4EmMultiModel::Initialise(), G4IonParametrisedLossModel::Initialise(), G4EmModelManager::Initialise(), G4IonParametrisedLossModel::PrintDEDXTable(), G4DNAAttachment::PrintInfo(), G4DNAChargeDecrease::PrintInfo(), G4DNAChargeIncrease::PrintInfo(), G4DNADissociation::PrintInfo(), G4DNAElastic::PrintInfo(), G4DNAExcitation::PrintInfo(), G4DNAIonisation::PrintInfo(), G4DNAPlasmonExcitation::PrintInfo(), G4DNAPositronium::PrintInfo(), G4DNARotExcitation::PrintInfo(), G4DNAVibExcitation::PrintInfo(), G4VLEPTSModel::ReadIXS(), G4LossTableManager::Register(), G4DNAModelInterface::RegisterModel(), G4IonParametrisedLossModel::RemoveDEDXTable(), G4VLEPTSModel::SampleEnergyLoss(), G4EmConfigurator::SetExtraEmModel(), G4EmConfigurator::SetModelForRegion(), and G4EmConfigurator::UpdateModelEnergyRange().

◆ GetPartialCrossSection()

G4double G4VEmModel::GetPartialCrossSection	(	const G4Material *	,
		G4int	level,
		const G4ParticleDefinition *	,
		G4double	kineticEnergy
	)

virtualinherited

Reimplemented in G4DNABornExcitationModel1, G4DNABornExcitationModel2, G4DNABornIonisationModel1, G4DNABornIonisationModel2, G4DNAMillerGreenExcitationModel, and G4DNARelativisticIonisationModel.

Definition at line 261 of file G4VEmModel.cc.

{ 
  return 0.0;
}

◆ GetParticleChangeForGamma()

G4ParticleChangeForGamma * G4VEmModel::GetParticleChangeForGamma ( )

protectedinherited

Definition at line 123 of file G4VEmModel.cc.

{
  G4ParticleChangeForGamma* p = nullptr;
  if (pParticleChange != nullptr) {
    p = static_cast<G4ParticleChangeForGamma*>(pParticleChange);
  } else {
    p = new G4ParticleChangeForGamma();
    pParticleChange = p;
  }
  if(fTripletModel != nullptr) { fTripletModel->SetParticleChange(p); }
  return p;
}

References G4VEmModel::fTripletModel, G4VEmModel::pParticleChange, and G4VEmModel::SetParticleChange().

◆ GetParticleChangeForLoss()

G4ParticleChangeForLoss * G4VEmModel::GetParticleChangeForLoss ( )

protectedinherited

Definition at line 108 of file G4VEmModel.cc.

{
  G4ParticleChangeForLoss* p = nullptr;
  if (pParticleChange != nullptr) {
    p = static_cast<G4ParticleChangeForLoss*>(pParticleChange);
  } else {
    p = new G4ParticleChangeForLoss();
    pParticleChange = p;
  }
  if(fTripletModel != nullptr) { fTripletModel->SetParticleChange(p); }
  return p;
}

References G4VEmModel::fTripletModel, G4VEmModel::pParticleChange, and G4VEmModel::SetParticleChange().

◆ GetParticleCharge()

G4double G4VEmModel::GetParticleCharge	(	const G4ParticleDefinition *	p,
		const G4Material *	,
		G4double	kineticEnergy
	)

virtualinherited

Reimplemented in G4IonParametrisedLossModel, G4BraggIonGasModel, G4BetheBlochIonGasModel, G4AtimaEnergyLossModel, G4BetheBlochModel, G4BraggIonModel, G4BraggModel, and G4LindhardSorensenIonModel.

Definition at line 391 of file G4VEmModel.cc.

{
  return p->GetPDGCharge();
}

References G4ParticleDefinition::GetPDGCharge().

Referenced by G4VEnergyLossProcess::AlongStepDoIt().

◆ GetTripletModel()

G4VEmModel * G4VEmModel::GetTripletModel ( )

inlineinherited

Definition at line 638 of file G4VEmModel.hh.

{
  return fTripletModel;
}

References G4VEmModel::fTripletModel.

Referenced by G4eBremsstrahlungRelModel::Initialise(), G4SeltzerBergerModel::Initialise(), and G4eBremsstrahlungRelModel::SampleSecondaries().

◆ GetVerbosityLevel()

G4int G4PenelopeIonisationModel::GetVerbosityLevel ( )

inline

Definition at line 109 of file G4PenelopeIonisationModel.hh.

109{return fVerboseLevel;};

References fVerboseLevel.

◆ HighEnergyActivationLimit()

G4double G4VEmModel::HighEnergyActivationLimit ( ) const

inlineinherited

Definition at line 669 of file G4VEmModel.hh.

{
  return eMaxActive;
}

References G4VEmModel::eMaxActive.

Referenced by G4EmModelManager::DumpModelList(), G4VMscModel::GetParticleChangeForMSC(), G4WentzelVIModel::Initialise(), and G4EmModelManager::Initialise().

◆ HighEnergyLimit()

G4double G4VEmModel::HighEnergyLimit ( ) const

inlineinherited

Definition at line 655 of file G4VEmModel.hh.

{
  return highLimit;
}

References G4VEmModel::highLimit.

◆ Initialise()

void G4PenelopeIonisationModel::Initialise	(	const G4ParticleDefinition *	particle,
		const G4DataVector &	theCuts
	)

overridevirtual

Implements G4VEmModel.

Definition at line 119 of file G4PenelopeIonisationModel.cc.

{
  if (fVerboseLevel > 3)
    G4cout << "Calling G4PenelopeIonisationModel::Initialise()" << G4endl;
 
  fAtomDeexcitation = G4LossTableManager::Instance()->AtomDeexcitation();
  //Issue warning if the AtomicDeexcitation has not been declared
  if (!fAtomDeexcitation)
    {
      G4cout << G4endl;
      G4cout << "WARNING from G4PenelopeIonisationModel " << G4endl;
      G4cout << "Atomic de-excitation module is not instantiated, so there will not be ";
      G4cout << "any fluorescence/Auger emission." << G4endl;
      G4cout << "Please make sure this is intended" << G4endl;
    }
 
  if (fAtomDeexcitation)
    fPIXEflag = fAtomDeexcitation->IsPIXEActive();
 
  //If the PIXE flag is active, the PIXE interface will take care of the 
  //atomic de-excitation. The model does not need to do that.
  //Issue warnings here
  if (fPIXEflag && IsMaster() && particle==G4Electron::Electron())
    {
      G4String theModel = G4EmParameters::Instance()->PIXEElectronCrossSectionModel();
      G4cout << "======================================================================" << G4endl;
      G4cout << "The G4PenelopeIonisationModel is being used with the PIXE flag ON." << G4endl;
      G4cout << "Atomic de-excitation will be produced statistically by the PIXE " << G4endl;
      G4cout << "interface by using the shell cross section --> " << theModel << G4endl;
      G4cout << "The built-in model procedure for atomic de-excitation is disabled. " << G4endl;
      G4cout << "*Please be sure this is intended*, or disable PIXE by" << G4endl;
      G4cout << "/process/em/pixe false" << G4endl;    
      G4cout << "======================================================================" << G4endl;
    }
 
  SetParticle(particle);
 
  //Only the master model creates/manages the tables. All workers get the 
  //pointer to the table, and use it as readonly
  if (IsMaster() && particle == fParticle)
    {
      //Set the number of bins for the tables. 20 points per decade
      fNBins = (size_t) (20*std::log10(HighEnergyLimit()/LowEnergyLimit()));
      fNBins = std::max(fNBins,(size_t)100);
      
      //Clear and re-build the tables
      if (fCrossSectionHandler)
    {
      delete fCrossSectionHandler;
      fCrossSectionHandler = 0;
    }
      fCrossSectionHandler = new G4PenelopeIonisationXSHandler(fNBins);
      fCrossSectionHandler->SetVerboseLevel(fVerboseLevel);
 
      //Build tables for all materials
      G4ProductionCutsTable* theCoupleTable = 
    G4ProductionCutsTable::GetProductionCutsTable();      
      for (size_t i=0;i<theCoupleTable->GetTableSize();i++)
    {
      const G4Material* theMat = 
        theCoupleTable->GetMaterialCutsCouple(i)->GetMaterial();
      //Forces the building of the cross section tables
      fCrossSectionHandler->BuildXSTable(theMat,theCuts.at(i),particle,
                           IsMaster());  
    }
 
      if (fVerboseLevel > 2) {
    G4cout << "Penelope Ionisation model v2008 is initialized " << G4endl
           << "Energy range: "
           << LowEnergyLimit() / keV << " keV - "
           << HighEnergyLimit() / GeV << " GeV. Using " 
           << fNBins << " bins." 
           << G4endl;
      }
    }
 
  if(fIsInitialised) 
    return;
  fParticleChange = GetParticleChangeForLoss();
  fIsInitialised = true;
 
  return;
}

References G4LossTableManager::AtomDeexcitation(), G4PenelopeIonisationXSHandler::BuildXSTable(), G4Electron::Electron(), fAtomDeexcitation, fCrossSectionHandler, fIsInitialised, fNBins, fParticle, fParticleChange, fPIXEflag, fVerboseLevel, G4cout, G4endl, G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4VEmModel::GetParticleChangeForLoss(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), GeV, G4VEmModel::HighEnergyLimit(), G4EmParameters::Instance(), G4LossTableManager::Instance(), G4VEmModel::IsMaster(), G4VAtomDeexcitation::IsPIXEActive(), keV, G4VEmModel::LowEnergyLimit(), G4INCL::Math::max(), G4EmParameters::PIXEElectronCrossSectionModel(), SetParticle(), and G4PenelopeIonisationXSHandler::SetVerboseLevel().

◆ InitialiseElementSelectors()

void G4VEmModel::InitialiseElementSelectors	(	const G4ParticleDefinition *	part,
		const G4DataVector &	cuts
	)

inherited

Definition at line 138 of file G4VEmModel.cc.

{
  // using spline for element selectors should be investigated in details
  // because small number of points may provide biased results
  // large number of points requires significant increase of memory
  G4bool spline = false;
  
  //G4cout << "IES: for " << GetName() << " Emin(MeV)= " << lowLimit/MeV 
  //         << " Emax(MeV)= " << highLimit/MeV << G4endl;
  
  // two times less bins because probability functon is normalized 
  // so correspondingly is more smooth
  if(highLimit <= lowLimit) { return; }
 
  G4int nbinsPerDec = G4EmParameters::Instance()->NumberOfBinsPerDecade();
 
  G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  G4int numOfCouples = theCoupleTable->GetTableSize();
 
  // prepare vector
  if(!elmSelectors) {
    elmSelectors = new std::vector<G4EmElementSelector*>;
  }
  if(numOfCouples > nSelectors) { 
    for(G4int i=nSelectors; i<numOfCouples; ++i) { 
      elmSelectors->push_back(nullptr); 
    }
    nSelectors = numOfCouples;
  }
 
  // initialise vector
  for(G4int i=0; i<numOfCouples; ++i) {
 
    // no need in element selectors for infinite cuts
    if(cuts[i] == DBL_MAX) { continue; }
   
    auto couple = theCoupleTable->GetMaterialCutsCouple(i); 
    auto material = couple->GetMaterial();
    SetCurrentCouple(couple);
 
    // selector already exist then delete
    delete (*elmSelectors)[i];
 
    G4double emin = std::max(lowLimit, MinPrimaryEnergy(material, part, cuts[i]));
    G4double emax = std::max(highLimit, 10*emin);
    static const G4double invlog106 = 1.0/(6*G4Log(10.));
    G4int nbins = (G4int)(nbinsPerDec*G4Log(emax/emin)*invlog106);
    nbins = std::max(nbins, 3);
 
    (*elmSelectors)[i] = new G4EmElementSelector(this,material,nbins,
                         emin,emax,spline);
    ((*elmSelectors)[i])->Initialise(part, cuts[i]);
    /*      
      G4cout << "G4VEmModel::InitialiseElmSelectors i= " << i 
             << "  "  << part->GetParticleName() 
             << " for " << GetName() << "  cut= " << cuts[i] 
             << "  " << (*elmSelectors)[i] << G4endl;      
      ((*elmSelectors)[i])->Dump(part);
    */
  } 
}

References DBL_MAX, G4VEmModel::elmSelectors, emax, G4Log(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), G4VEmModel::highLimit, G4VEmModel::Initialise(), G4EmParameters::Instance(), G4VEmModel::lowLimit, eplot::material, G4INCL::Math::max(), G4VEmModel::MinPrimaryEnergy(), G4VEmModel::nSelectors, G4EmParameters::NumberOfBinsPerDecade(), and G4VEmModel::SetCurrentCouple().

◆ InitialiseForElement()

void G4VEmModel::InitialiseForElement	(	const G4ParticleDefinition *	,
		G4int	Z
	)

virtualinherited

◆ InitialiseForMaterial()

void G4VEmModel::InitialiseForMaterial	(	const G4ParticleDefinition *	part,
		const G4Material *	material
	)

virtualinherited

Definition at line 209 of file G4VEmModel.cc.

{
  if(material != nullptr) {
    size_t n = material->GetNumberOfElements();
    for(size_t i=0; i<n; ++i) {
      G4int Z = material->GetElement(i)->GetZasInt();
      InitialiseForElement(part, Z);
    }
  }
}

References G4VEmModel::InitialiseForElement(), eplot::material, CLHEP::detail::n, and Z.

Referenced by G4EmCalculator::FindEmModel().

◆ InitialiseLocal()

void G4PenelopeIonisationModel::InitialiseLocal	(	const G4ParticleDefinition *	part,
		G4VEmModel *	masterModel
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 206 of file G4PenelopeIonisationModel.cc.

{
  if (fVerboseLevel > 3)
    G4cout << "Calling  G4PenelopeIonisationModel::InitialiseLocal()" << G4endl;
  //
  //Check that particle matches: one might have multiple master models (e.g. 
  //for e+ and e-).
  //
  if (part == fParticle)
    {
      //Get the const table pointers from the master to the workers
      const G4PenelopeIonisationModel* theModel = 
    static_cast<G4PenelopeIonisationModel*> (masterModel);
      
      //Copy pointers to the data tables
      fCrossSectionHandler = theModel->fCrossSectionHandler;
      
      //copy data
      fNBins = theModel->fNBins;
      
      //Same verbosity for all workers, as the master
      fVerboseLevel = theModel->fVerboseLevel;
    }
  return;
}

References fCrossSectionHandler, fNBins, fParticle, fVerboseLevel, G4cout, and G4endl.

◆ IsActive()

G4bool G4VEmModel::IsActive ( G4double kinEnergy ) const

inlineinherited

Definition at line 795 of file G4VEmModel.hh.

{
  return (kinEnergy >= eMinActive && kinEnergy <= eMaxActive);
}

References G4VEmModel::eMaxActive, and G4VEmModel::eMinActive.

Referenced by G4VEnergyLossProcess::AlongStepDoIt(), G4NuclearStopping::AlongStepDoIt(), G4VMultipleScattering::AlongStepGetPhysicalInteractionLength(), G4VEnergyLossProcess::AlongStepGetPhysicalInteractionLength(), G4VEmProcess::PostStepDoIt(), G4VEnergyLossProcess::PostStepDoIt(), G4VEmProcess::PostStepGetPhysicalInteractionLength(), and G4VEnergyLossProcess::PostStepGetPhysicalInteractionLength().

◆ IsLocked()

G4bool G4VEmModel::IsLocked ( ) const

inlineinherited

Definition at line 877 of file G4VEmModel.hh.

{
  return isLocked;
}

References G4VEmModel::isLocked.

Referenced by G4VMscModel::InitialiseParameters(), G4VMscModel::SetGeomFactor(), G4VMscModel::SetLambdaLimit(), G4VMscModel::SetLateralDisplasmentFlag(), G4VMscModel::SetRangeFactor(), G4VMscModel::SetSafetyFactor(), G4VMscModel::SetSampleZ(), G4VMscModel::SetSkin(), and G4VMscModel::SetStepLimitType().

◆ IsMaster()

G4bool G4VEmModel::IsMaster ( ) const

inlineinherited

Definition at line 746 of file G4VEmModel.hh.

{
  return isMaster;
}

References G4VEmModel::isMaster.

◆ LowEnergyActivationLimit()

G4double G4VEmModel::LowEnergyActivationLimit ( ) const

inlineinherited

Definition at line 676 of file G4VEmModel.hh.

{
  return eMinActive;
}

References G4VEmModel::eMinActive.

Referenced by G4EmModelManager::DumpModelList(), G4VMscModel::GetParticleChangeForMSC(), G4WentzelVIModel::Initialise(), and G4EmModelManager::Initialise().

◆ LowEnergyLimit()

G4double G4VEmModel::LowEnergyLimit ( ) const

inlineinherited

Definition at line 662 of file G4VEmModel.hh.

{
  return lowLimit;
}

References G4VEmModel::lowLimit.

◆ LPMFlag()

G4bool G4VEmModel::LPMFlag ( ) const

inlineinherited

Definition at line 697 of file G4VEmModel.hh.

{
  return theLPMflag;
}

References G4VEmModel::theLPMflag.

Referenced by G4eBremsstrahlungRelModel::Initialise(), G4eBremsstrahlungRelModel::SetupForMaterial(), G4SeltzerBergerModel::SetupForMaterial(), and G4eBremsstrahlungRelModel::~G4eBremsstrahlungRelModel().

◆ MaxSecondaryEnergy()

G4double G4VEmModel::MaxSecondaryEnergy	(	const G4ParticleDefinition *	,
		G4double	kineticEnergy
	)

protectedvirtualinherited

Reimplemented in G4BraggIonModel, G4BraggModel, G4ICRU73QOModel, G4MollerBhabhaModel, G4PAIModel, G4PAIPhotModel, G4mplIonisationWithDeltaModel, G4MuBetheBlochModel, G4AtimaEnergyLossModel, G4BetheBlochModel, G4LindhardSorensenIonModel, and G4IonParametrisedLossModel.

Definition at line 432 of file G4VEmModel.cc.

{
  return kineticEnergy;
}

Referenced by G4VEmModel::MaxSecondaryKinEnergy().

◆ MaxSecondaryKinEnergy()

G4double G4VEmModel::MaxSecondaryKinEnergy ( const G4DynamicParticle * dynParticle )

inlineinherited

Definition at line 520 of file G4VEmModel.hh.

{
  return MaxSecondaryEnergy(dynPart->GetParticleDefinition(),
                            dynPart->GetKineticEnergy());
}

References G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetParticleDefinition(), and G4VEmModel::MaxSecondaryEnergy().

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4VEnergyLossProcess::AlongStepDoIt(), G4VEnergyLossProcess::GetDEDXDispersion(), G4MuBetheBlochModel::SampleSecondaries(), G4PolarizedIonisationModel::SampleSecondaries(), G4BraggIonModel::SampleSecondaries(), G4BraggModel::SampleSecondaries(), G4ICRU73QOModel::SampleSecondaries(), and G4IonParametrisedLossModel::SampleSecondaries().

◆ MinEnergyCut()

G4double G4PenelopeIonisationModel::MinEnergyCut	(	const G4ParticleDefinition *	,
		const G4MaterialCutsCouple *
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 442 of file G4PenelopeIonisationModel.cc.

{
  return fIntrinsicLowEnergyLimit;
}

References fIntrinsicLowEnergyLimit.

◆ MinPrimaryEnergy()

G4double G4VEmModel::MinPrimaryEnergy	(	const G4Material *	,
		const G4ParticleDefinition *	,
		G4double	cut = `0.0`
	)

virtualinherited

Reimplemented in G4eBremsstrahlungRelModel, G4BoldyshevTripletModel, G4eCoulombScatteringModel, G4hCoulombScatteringModel, G4LivermoreNuclearGammaConversionModel, G4LivermorePolarizedGammaConversionModel, G4MuBremsstrahlungModel, G4MuPairProductionModel, and G4eDPWACoulombScatteringModel.

Definition at line 415 of file G4VEmModel.cc.

{
  return 0.0;
}

Referenced by G4LossTableBuilder::BuildTableForModel(), and G4VEmModel::InitialiseElementSelectors().

◆ ModelDescription()

void G4VEmModel::ModelDescription ( std::ostream & outFile ) const

virtualinherited

Reimplemented in G4eeToHadronsMultiModel.

Definition at line 469 of file G4VEmModel.cc.

{
  outFile << "The description for this model has not been written yet.\n";
}

◆ operator=()

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

delete

◆ PolarAngleLimit()

G4double G4VEmModel::PolarAngleLimit ( ) const

inlineinherited

Definition at line 683 of file G4VEmModel.hh.

{
  return polarAngleLimit;
}

References G4VEmModel::polarAngleLimit.

Referenced by G4eCoulombScatteringModel::Initialise(), G4eDPWACoulombScatteringModel::Initialise(), G4hCoulombScatteringModel::Initialise(), and G4WentzelVIModel::Initialise().

◆ SampleFinalStateElectron()

void G4PenelopeIonisationModel::SampleFinalStateElectron	(	const G4Material *	mat,
		G4double	cutEnergy,
		G4double	kineticEnergy
	)

private

Definition at line 697 of file G4PenelopeIonisationModel.cc.

{
  // This method sets the final ionisation parameters
  // fKineticEnergy1, fCosThetaPrimary (= updates of the primary e-)
  // fEnergySecondary, fCosThetaSecondary (= info of delta-ray)
  // fTargetOscillator (= ionised oscillator)
  //
  // The method implements SUBROUTINE EINa of Penelope
  //
 
  G4PenelopeOscillatorTable* theTable = fOscManager->GetOscillatorTableIonisation(mat);
  size_t numberOfOscillators = theTable->size();
  const G4PenelopeCrossSection* theXS = 
    fCrossSectionHandler->GetCrossSectionTableForCouple(G4Electron::Electron(),mat,
                              cutEnergy);
  G4double delta = fCrossSectionHandler->GetDensityCorrection(mat,kineticEnergy);
 
  // Selection of the active oscillator
  G4double TST = G4UniformRand();
  fTargetOscillator = numberOfOscillators-1; //initialization, last oscillator
  G4double XSsum = 0.;
 
  for (size_t i=0;i<numberOfOscillators-1;i++)
    {     
      XSsum += theXS->GetNormalizedShellCrossSection(i,kineticEnergy); 
        
      if (XSsum > TST)
    {
      fTargetOscillator = (G4int) i;
      break;
    }
    }
 
  if (fVerboseLevel > 3)
    {
      G4cout << "SampleFinalStateElectron: sampled oscillator #" << 
    fTargetOscillator << "." << G4endl;
      G4cout << "Ionisation energy: " << 
    (*theTable)[fTargetOscillator]->GetIonisationEnergy()/eV << 
    " eV " << G4endl;
      G4cout << "Resonance energy: : " << 
    (*theTable)[fTargetOscillator]->GetResonanceEnergy()/eV << " eV "
         << G4endl;
    }
  //Constants
  G4double rb = kineticEnergy + 2.0*electron_mass_c2;
  G4double gam = 1.0+kineticEnergy/electron_mass_c2;
  G4double gam2 = gam*gam;
  G4double beta2 = (gam2-1.0)/gam2;
  G4double amol = ((gam-1.0)/gam)*((gam-1.0)/gam);
 
  //Partial cross section of the active oscillator
  G4double resEne = (*theTable)[fTargetOscillator]->GetResonanceEnergy();
  G4double invResEne = 1.0/resEne;
  G4double ionEne = (*theTable)[fTargetOscillator]->GetIonisationEnergy();
  G4double cutoffEne = (*theTable)[fTargetOscillator]->GetCutoffRecoilResonantEnergy();
  G4double XHDL = 0.;
  G4double XHDT = 0.;
  G4double QM = 0.;
  G4double cps = 0.;
  G4double cp = 0.;
 
  //Distant excitations
  if (resEne > cutEnergy && resEne < kineticEnergy)
    {
      cps = kineticEnergy*rb;
      cp = std::sqrt(cps);
      G4double XHDT0 = std::max(G4Log(gam2)-beta2-delta,0.);
      if (resEne > 1.0e-6*kineticEnergy)
    {
      G4double cpp = std::sqrt((kineticEnergy-resEne)*(kineticEnergy-resEne+2.0*electron_mass_c2));
      QM = std::sqrt((cp-cpp)*(cp-cpp)+electron_mass_c2*electron_mass_c2)-electron_mass_c2;
    }
      else
    {
      QM = resEne*resEne/(beta2*2.0*electron_mass_c2);
      QM *= (1.0-QM*0.5/electron_mass_c2);
    }
      if (QM < cutoffEne)
    {
      XHDL = G4Log(cutoffEne*(QM+2.0*electron_mass_c2)/(QM*(cutoffEne+2.0*electron_mass_c2)))
        *invResEne;
      XHDT = XHDT0*invResEne;     
    }
      else
    {
      QM = cutoffEne;
      XHDL = 0.;
      XHDT = 0.;
    }
    }
  else
    {
      QM = cutoffEne;
      cps = 0.;
      cp = 0.;
      XHDL = 0.;
      XHDT = 0.;
    }
 
  //Close collisions
  G4double EE = kineticEnergy + ionEne;
  G4double wmaxc = 0.5*EE;
  G4double wcl = std::max(cutEnergy,cutoffEne);
  G4double rcl = wcl/EE;
  G4double XHC = 0.;
  if (wcl < wmaxc)
    {
      G4double rl1 = 1.0-rcl;
      G4double rrl1 = 1.0/rl1;
      XHC = (amol*(0.5-rcl)+1.0/rcl-rrl1+
         (1.0-amol)*G4Log(rcl*rrl1))/EE;
    }
 
  //Total cross section per molecule for the active shell, in cm2
  G4double XHTOT = XHC + XHDL + XHDT;
 
  //very small cross section, do nothing
  if (XHTOT < 1.e-14*barn)
    {
      fKineticEnergy1=kineticEnergy;
      fCosThetaPrimary=1.0;
      fEnergySecondary=0.0;
      fCosThetaSecondary=1.0;
      fTargetOscillator = numberOfOscillators-1;
      return;
    }
  
  //decide which kind of interaction we'll have
  TST = XHTOT*G4UniformRand();
  
  // Hard close collision
  G4double TS1 = XHC;
  
  if (TST < TS1)
    {
      G4double A = 5.0*amol;
      G4double ARCL = A*0.5*rcl;
      G4double rk=0.;
      G4bool loopAgain = false;
      do
    {
      loopAgain = false;
      G4double fb = (1.0+ARCL)*G4UniformRand();
      if (fb < 1)
        rk = rcl/(1.0-fb*(1.0-(rcl+rcl)));       
      else
        rk = rcl + (fb-1.0)*(0.5-rcl)/ARCL;
      G4double rk2 = rk*rk;
      G4double rkf = rk/(1.0-rk);
      G4double phi = 1.0+rkf*rkf-rkf+amol*(rk2+rkf);
      if (G4UniformRand()*(1.0+A*rk2) > phi)
        loopAgain = true;
    }while(loopAgain);
      //energy and scattering angle (primary electron)
      G4double deltaE = rk*EE;
      
      fKineticEnergy1 = kineticEnergy - deltaE;
      fCosThetaPrimary = std::sqrt(fKineticEnergy1*rb/(kineticEnergy*(rb-deltaE)));
      //energy and scattering angle of the delta ray
      fEnergySecondary = deltaE - ionEne; //subtract ionisation energy
      fCosThetaSecondary= std::sqrt(deltaE*rb/(kineticEnergy*(deltaE+2.0*electron_mass_c2)));
      if (fVerboseLevel > 3)
    G4cout << "SampleFinalStateElectron: sampled close collision " << G4endl;
      return;                    
    }
 
  //Hard distant longitudinal collisions
  TS1 += XHDL;
  G4double deltaE = resEne;
  fKineticEnergy1 = kineticEnergy - deltaE;
 
  if (TST < TS1)
    {
      G4double QS = QM/(1.0+QM*0.5/electron_mass_c2);
      G4double Q = QS/(std::pow((QS/cutoffEne)*(1.0+cutoffEne*0.5/electron_mass_c2),G4UniformRand())
               - (QS*0.5/electron_mass_c2));
      G4double QTREV = Q*(Q+2.0*electron_mass_c2);
      G4double cpps = fKineticEnergy1*(fKineticEnergy1+2.0*electron_mass_c2);
      fCosThetaPrimary = (cpps+cps-QTREV)/(2.0*cp*std::sqrt(cpps));
      if (fCosThetaPrimary > 1.) 
    fCosThetaPrimary = 1.0;
      //energy and emission angle of the delta ray
      fEnergySecondary = deltaE - ionEne;
      fCosThetaSecondary = 0.5*(deltaE*(kineticEnergy+rb-deltaE)+QTREV)/std::sqrt(cps*QTREV);
      if (fCosThetaSecondary > 1.0)
    fCosThetaSecondary = 1.0;
      if (fVerboseLevel > 3)
    G4cout << "SampleFinalStateElectron: sampled distant longitudinal collision " << G4endl;
      return;      
    }
 
  //Hard distant transverse collisions
  fCosThetaPrimary = 1.0;
  //energy and emission angle of the delta ray
  fEnergySecondary = deltaE - ionEne;
  fCosThetaSecondary = 0.5;
  if (fVerboseLevel > 3)
    G4cout << "SampleFinalStateElectron: sampled distant transverse collision " << G4endl;
 
  return;
}

References A, barn, cp, G4Electron::Electron(), source.hepunit::electron_mass_c2, eV, fCosThetaPrimary, fCosThetaSecondary, fCrossSectionHandler, fEnergySecondary, fKineticEnergy1, fOscManager, fTargetOscillator, fVerboseLevel, G4cout, G4endl, G4Log(), G4UniformRand, G4InuclParticleNames::gam, G4PenelopeIonisationXSHandler::GetCrossSectionTableForCouple(), G4PenelopeIonisationXSHandler::GetDensityCorrection(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4PenelopeOscillatorManager::GetOscillatorTableIonisation(), G4INCL::Math::max(), and Q.

Referenced by SampleSecondaries().

◆ SampleFinalStatePositron()

void G4PenelopeIonisationModel::SampleFinalStatePositron	(	const G4Material *	mat,
		G4double	cutEnergy,
		G4double	kineticEnergy
	)

private

Definition at line 904 of file G4PenelopeIonisationModel.cc.

{
  // This method sets the final ionisation parameters
  // fKineticEnergy1, fCosThetaPrimary (= updates of the primary e-)
  // fEnergySecondary, fCosThetaSecondary (= info of delta-ray)
  // fTargetOscillator (= ionised oscillator)
  //
  // The method implements SUBROUTINE PINa of Penelope
  // 
 
  G4PenelopeOscillatorTable* theTable = fOscManager->GetOscillatorTableIonisation(mat);
  size_t numberOfOscillators = theTable->size();
  const G4PenelopeCrossSection* theXS = 
    fCrossSectionHandler->GetCrossSectionTableForCouple(G4Positron::Positron(),mat,
                            cutEnergy);
  G4double delta = fCrossSectionHandler->GetDensityCorrection(mat,kineticEnergy);
 
  // Selection of the active oscillator
  G4double TST = G4UniformRand();
  fTargetOscillator = numberOfOscillators-1; //initialization, last oscillator
  G4double XSsum = 0.;
  for (size_t i=0;i<numberOfOscillators-1;i++)
    {
      XSsum += theXS->GetNormalizedShellCrossSection(i,kineticEnergy);      
      if (XSsum > TST)
    {
      fTargetOscillator = (G4int) i;
      break;
    }
    }
 
  if (fVerboseLevel > 3)
    {
      G4cout << "SampleFinalStatePositron: sampled oscillator #" << 
    fTargetOscillator << "." << G4endl;
      G4cout << "Ionisation energy: " << (*theTable)[fTargetOscillator]->GetIonisationEnergy()/eV 
         << " eV " << G4endl; 
      G4cout << "Resonance energy: : " << (*theTable)[fTargetOscillator]->GetResonanceEnergy()/eV 
         << " eV " << G4endl;
    }
 
  //Constants
  G4double rb = kineticEnergy + 2.0*electron_mass_c2;
  G4double gam = 1.0+kineticEnergy/electron_mass_c2;
  G4double gam2 = gam*gam;
  G4double beta2 = (gam2-1.0)/gam2;
  G4double g12 = (gam+1.0)*(gam+1.0);
  G4double amol = ((gam-1.0)/gam)*((gam-1.0)/gam);
  //Bhabha coefficients
  G4double bha1 = amol*(2.0*g12-1.0)/(gam2-1.0);
  G4double bha2 = amol*(3.0+1.0/g12);
  G4double bha3 = amol*2.0*gam*(gam-1.0)/g12;
  G4double bha4 = amol*(gam-1.0)*(gam-1.0)/g12;
 
  //
  //Partial cross section of the active oscillator
  //
  G4double resEne = (*theTable)[fTargetOscillator]->GetResonanceEnergy();
  G4double invResEne = 1.0/resEne;
  G4double ionEne = (*theTable)[fTargetOscillator]->GetIonisationEnergy();
  G4double cutoffEne = (*theTable)[fTargetOscillator]->GetCutoffRecoilResonantEnergy();
 
  G4double XHDL = 0.;
  G4double XHDT = 0.;
  G4double QM = 0.;
  G4double cps = 0.;
  G4double cp = 0.;
 
  //Distant excitations XS (same as for electrons)
  if (resEne > cutEnergy && resEne < kineticEnergy)
    {
      cps = kineticEnergy*rb;
      cp = std::sqrt(cps);
      G4double XHDT0 = std::max(G4Log(gam2)-beta2-delta,0.);
      if (resEne > 1.0e-6*kineticEnergy)
    {
      G4double cpp = std::sqrt((kineticEnergy-resEne)*(kineticEnergy-resEne+2.0*electron_mass_c2));
      QM = std::sqrt((cp-cpp)*(cp-cpp)+electron_mass_c2*electron_mass_c2)-electron_mass_c2;
    }
      else
    {
      QM = resEne*resEne/(beta2*2.0*electron_mass_c2);
      QM *= (1.0-QM*0.5/electron_mass_c2);
    }
      if (QM < cutoffEne)
    {
      XHDL = G4Log(cutoffEne*(QM+2.0*electron_mass_c2)/(QM*(cutoffEne+2.0*electron_mass_c2)))
        *invResEne;
      XHDT = XHDT0*invResEne;     
    }
      else
    {
      QM = cutoffEne;
      XHDL = 0.;
      XHDT = 0.;
    }
    }
  else
    {
      QM = cutoffEne;
      cps = 0.;
      cp = 0.;
      XHDL = 0.;
      XHDT = 0.;
    }
  //Close collisions (Bhabha)
  G4double wmaxc = kineticEnergy;
  G4double wcl = std::max(cutEnergy,cutoffEne);
  G4double rcl = wcl/kineticEnergy;
  G4double XHC = 0.;
  if (wcl < wmaxc)
    {
      G4double rl1 = 1.0-rcl;
      XHC = ((1.0/rcl-1.0)+bha1*G4Log(rcl)+bha2*rl1
         + (bha3/2.0)*(rcl*rcl-1.0) 
         + (bha4/3.0)*(1.0-rcl*rcl*rcl))/kineticEnergy;
    }
 
  //Total cross section per molecule for the active shell, in cm2
  G4double XHTOT = XHC + XHDL + XHDT;
 
  //very small cross section, do nothing
  if (XHTOT < 1.e-14*barn)
    {
      fKineticEnergy1=kineticEnergy;
      fCosThetaPrimary=1.0;
      fEnergySecondary=0.0;
      fCosThetaSecondary=1.0;
      fTargetOscillator = numberOfOscillators-1;
      return;
    }
 
  //decide which kind of interaction we'll have
  TST = XHTOT*G4UniformRand();
  
  // Hard close collision
  G4double TS1 = XHC;
  if (TST < TS1)
    {
      G4double rl1 = 1.0-rcl;
      G4double rk=0.;
      G4bool loopAgain = false;
      do
    {
      loopAgain = false;
      rk = rcl/(1.0-G4UniformRand()*rl1);
      G4double phi = 1.0-rk*(bha1-rk*(bha2-rk*(bha3-bha4*rk)));
      if (G4UniformRand() > phi)
        loopAgain = true;
    }while(loopAgain);
      //energy and scattering angle (primary electron)
      G4double deltaE = rk*kineticEnergy;      
      fKineticEnergy1 = kineticEnergy - deltaE;
      fCosThetaPrimary = std::sqrt(fKineticEnergy1*rb/(kineticEnergy*(rb-deltaE)));
      //energy and scattering angle of the delta ray
      fEnergySecondary = deltaE - ionEne; //subtract ionisation energy
      fCosThetaSecondary= std::sqrt(deltaE*rb/(kineticEnergy*(deltaE+2.0*electron_mass_c2)));
      if (fVerboseLevel > 3)
    G4cout << "SampleFinalStatePositron: sampled close collision " << G4endl;
      return;                    
    }
 
  //Hard distant longitudinal collisions
  TS1 += XHDL;
  G4double deltaE = resEne;
  fKineticEnergy1 = kineticEnergy - deltaE;
  if (TST < TS1)
    {
      G4double QS = QM/(1.0+QM*0.5/electron_mass_c2);
      G4double Q = QS/(std::pow((QS/cutoffEne)*(1.0+cutoffEne*0.5/electron_mass_c2),G4UniformRand())
               - (QS*0.5/electron_mass_c2));
      G4double QTREV = Q*(Q+2.0*electron_mass_c2);
      G4double cpps = fKineticEnergy1*(fKineticEnergy1+2.0*electron_mass_c2);
      fCosThetaPrimary = (cpps+cps-QTREV)/(2.0*cp*std::sqrt(cpps));
      if (fCosThetaPrimary > 1.) 
    fCosThetaPrimary = 1.0;
      //energy and emission angle of the delta ray
      fEnergySecondary = deltaE - ionEne;
      fCosThetaSecondary = 0.5*(deltaE*(kineticEnergy+rb-deltaE)+QTREV)/std::sqrt(cps*QTREV);
      if (fCosThetaSecondary > 1.0)
    fCosThetaSecondary = 1.0;
      if (fVerboseLevel > 3)
    G4cout << "SampleFinalStatePositron: sampled distant longitudinal collision " << G4endl;
      return;      
    }
 
  //Hard distant transverse collisions
  fCosThetaPrimary = 1.0;
  //energy and emission angle of the delta ray
  fEnergySecondary = deltaE - ionEne;
  fCosThetaSecondary = 0.5;
 
  if (fVerboseLevel > 3)    
    G4cout << "SampleFinalStatePositron: sampled distant transverse collision " << G4endl;
    
  return;
}

References barn, cp, source.hepunit::electron_mass_c2, eV, fCosThetaPrimary, fCosThetaSecondary, fCrossSectionHandler, fEnergySecondary, fKineticEnergy1, fOscManager, fTargetOscillator, fVerboseLevel, G4cout, G4endl, G4Log(), G4UniformRand, G4InuclParticleNames::gam, G4PenelopeIonisationXSHandler::GetCrossSectionTableForCouple(), G4PenelopeIonisationXSHandler::GetDensityCorrection(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4PenelopeOscillatorManager::GetOscillatorTableIonisation(), G4INCL::Math::max(), G4Positron::Positron(), and Q.

Referenced by SampleSecondaries().

◆ SampleSecondaries()

void G4PenelopeIonisationModel::SampleSecondaries	(	std::vector< G4DynamicParticle * > *	fvect,
		const G4MaterialCutsCouple *	couple,
		const G4DynamicParticle *	aDynamicParticle,
		G4double	tmin,
		G4double	maxEnergy
	)

overridevirtual

Implements G4VEmModel.

Definition at line 450 of file G4PenelopeIonisationModel.cc.

{
  // Penelope model v2008 to sample the final state following an hard inelastic interaction.
  // It makes use of the Generalised Oscillator Strength (GOS) model from
  //  D. Liljequist, J. Phys. D: Appl. Phys. 16 (1983) 1567
  //
  // The GOS model is used to calculate the individual cross sections for all
  // the atomic oscillators coming in the play, taking into account the three
  // contributions (distant longitudinal collisions, distant transverse collisions and
  // close collisions). Then the target shell and the interaction channel are
  // sampled. Final state of the delta-ray (energy, angle) are generated according
  // to the analytical distributions (dSigma/dW) for the selected interaction
  // channels.
  // For e-, the maximum energy for the delta-ray is initialEnergy/2. (because
  // particles are indistinghusbable), while it is the full initialEnergy for
  // e+.
  // The efficiency of the random sampling algorithm (e.g. for close collisions)
  // decreases when initial and cutoff energy increase (e.g. 87% for 10-keV primary
  // and 1 keV threshold, 99% for 10-MeV primary and 10-keV threshold).
  // Differential cross sections have a different form for e+ and e-.
  //
  // WARNING: The model provides an _average_ description of inelastic collisions.
  // Anyway, the energy spectrum associated to distant excitations of a given
  // atomic shell is approximated as a single resonance. The simulated energy spectra
  // show _unphysical_ narrow peaks at energies that are multiple of the shell
  // resonance energies. The spurious speaks are automatically smoothed out after
  // multiple inelastic collisions.
  //
  // The model determines also the original shell from which the delta-ray is expelled,
  // in order to produce fluorescence de-excitation (from G4DeexcitationManager)
  //
  // Fermi density correction is calculated analytically according to
  //  U. Fano, Ann. Rev. Nucl. Sci. 13 (1963),1
 
  if (fVerboseLevel > 3)
    G4cout << "Calling SamplingSecondaries() of G4PenelopeIonisationModel" << G4endl;
 
  G4double kineticEnergy0 = aDynamicParticle->GetKineticEnergy();
  const G4ParticleDefinition* theParticle = aDynamicParticle->GetDefinition();
 
  if (kineticEnergy0 <= fIntrinsicLowEnergyLimit)
  {
    fParticleChange->SetProposedKineticEnergy(0.);
    fParticleChange->ProposeLocalEnergyDeposit(kineticEnergy0);
    return ;
  }
 
  const G4Material* material = couple->GetMaterial();
  const G4PenelopeOscillatorTable* theTable = fOscManager->GetOscillatorTableIonisation(material); 
 
  G4ParticleMomentum particleDirection0 = aDynamicParticle->GetMomentumDirection();
 
  //Initialise final-state variables. The proper values will be set by the methods
  // SampleFinalStateElectron() and SampleFinalStatePositron()
  fKineticEnergy1=kineticEnergy0;
  fCosThetaPrimary=1.0;
  fEnergySecondary=0.0;
  fCosThetaSecondary=1.0;
  fTargetOscillator = -1;
     
  if (theParticle == G4Electron::Electron())
    SampleFinalStateElectron(material,cutE,kineticEnergy0);
  else if (theParticle == G4Positron::Positron())
    SampleFinalStatePositron(material,cutE,kineticEnergy0);
  else
    {
      G4ExceptionDescription ed;
      ed << "Invalid particle " << theParticle->GetParticleName() << G4endl;
      G4Exception("G4PenelopeIonisationModel::SamplingSecondaries()",
          "em0001",FatalException,ed);
      
    }
  if (fEnergySecondary == 0) return;
 
  if (fVerboseLevel > 3)
  {
     G4cout << "G4PenelopeIonisationModel::SamplingSecondaries() for " << 
    theParticle->GetParticleName() << G4endl;
      G4cout << "Final eKin = " << fKineticEnergy1 << " keV" << G4endl;
      G4cout << "Final cosTheta = " << fCosThetaPrimary << G4endl;
      G4cout << "Delta-ray eKin = " << fEnergySecondary << " keV" << G4endl;
      G4cout << "Delta-ray cosTheta = " << fCosThetaSecondary << G4endl;
      G4cout << "Oscillator: " << fTargetOscillator << G4endl;
   }
   
  //Update the primary particle
  G4double sint = std::sqrt(1. - fCosThetaPrimary*fCosThetaPrimary);
  G4double phiPrimary  = twopi * G4UniformRand();
  G4double dirx = sint * std::cos(phiPrimary);
  G4double diry = sint * std::sin(phiPrimary);
  G4double dirz = fCosThetaPrimary;
 
  G4ThreeVector electronDirection1(dirx,diry,dirz);
  electronDirection1.rotateUz(particleDirection0);
   
  if (fKineticEnergy1 > 0)
    {
      fParticleChange->ProposeMomentumDirection(electronDirection1);
      fParticleChange->SetProposedKineticEnergy(fKineticEnergy1);
    }
  else    
    fParticleChange->SetProposedKineticEnergy(0.);
    
  //Generate the delta ray
  G4double ionEnergyInPenelopeDatabase = 
    (*theTable)[fTargetOscillator]->GetIonisationEnergy();
  
  //Now, try to handle fluorescence
  //Notice: merged levels are indicated with Z=0 and flag=30
  G4int shFlag = (*theTable)[fTargetOscillator]->GetShellFlag(); 
  G4int Z = (G4int) (*theTable)[fTargetOscillator]->GetParentZ();
 
  //initialize here, then check photons created by Atomic-Deexcitation, and the final state e-
  const G4AtomicTransitionManager* transitionManager = G4AtomicTransitionManager::Instance();
  G4double bindingEnergy = 0.*eV;  
 
  const G4AtomicShell* shell = nullptr;
  //Real level
  if (Z > 0 && shFlag<30)
    {
      shell = transitionManager->Shell(Z,shFlag-1);
      bindingEnergy = shell->BindingEnergy();
      //shellId = shell->ShellId();
    }
 
  //correct the fEnergySecondary to account for the fact that the Penelope 
  //database of ionisation energies is in general (slightly) different 
  //from the fluorescence database used in Geant4.
  fEnergySecondary += ionEnergyInPenelopeDatabase-bindingEnergy;
  
  G4double localEnergyDeposit = bindingEnergy;
  //testing purposes only
  G4double energyInFluorescence = 0;
  G4double energyInAuger = 0; 
 
  if (fEnergySecondary < 0)
    {
      //It means that there was some problem/mismatch between the two databases. 
      //In this case, the available energy is ok to excite the level according 
      //to the Penelope database, but not according to the Geant4 database
      //Full residual energy is deposited locally
      localEnergyDeposit += fEnergySecondary;
      fEnergySecondary = 0.0;
    }
 
  //Notice: shell might be nullptr (invalid!) if shFlag=30. Must be protected
  //Disable the built-in de-excitation of the PIXE flag is active. In this 
  //case, the PIXE interface takes care (statistically) of producing the 
  //de-excitation.
  //Now, take care of fluorescence, if required
  if (fAtomDeexcitation && !fPIXEflag && shell)
    {
      G4int index = couple->GetIndex();
      if (fAtomDeexcitation->CheckDeexcitationActiveRegion(index))
    {
      size_t nBefore = fvect->size();
      fAtomDeexcitation->GenerateParticles(fvect,shell,Z,index);
      size_t nAfter = fvect->size(); 
      
      if (nAfter>nBefore) //actual production of fluorescence
        {
          for (size_t j=nBefore;j<nAfter;j++) //loop on products
        {
          G4double itsEnergy = ((*fvect)[j])->GetKineticEnergy();
                  if (itsEnergy < localEnergyDeposit) // valid secondary, generate it
                    {
              localEnergyDeposit -= itsEnergy;
              if (((*fvect)[j])->GetParticleDefinition() == G4Gamma::Definition())
                energyInFluorescence += itsEnergy;
              else if (((*fvect)[j])->GetParticleDefinition() == G4Electron::Definition())
                energyInAuger += itsEnergy;
                    }
                  else //invalid secondary: takes more than the available energy: delete it
            {
              delete (*fvect)[j];
              (*fvect)[j] = nullptr;
            }           
        }
        }
    }
    }
 
  // Generate the delta ray --> to be done only if above cut
  if (fEnergySecondary > cutE)
    {
      G4DynamicParticle* electron = nullptr;
      G4double sinThetaE = std::sqrt(1.-fCosThetaSecondary*fCosThetaSecondary);
      G4double phiEl = phiPrimary+pi; //pi with respect to the primary electron/positron
      G4double xEl = sinThetaE * std::cos(phiEl);
      G4double yEl = sinThetaE * std::sin(phiEl);
      G4double zEl = fCosThetaSecondary;
      G4ThreeVector eDirection(xEl,yEl,zEl); //electron direction
      eDirection.rotateUz(particleDirection0);
      electron = new G4DynamicParticle (G4Electron::Electron(),
                    eDirection,fEnergySecondary) ;
      fvect->push_back(electron);
    }
  else
    {
      localEnergyDeposit += fEnergySecondary;
      fEnergySecondary = 0;
    }
 
  if (localEnergyDeposit < 0) //Should not be: issue a G4Exception (warning)
    {
      G4Exception("G4PenelopeIonisationModel::SampleSecondaries()",
          "em2099",JustWarning,"WARNING: Negative local energy deposit");
      localEnergyDeposit=0.;
    }
  fParticleChange->ProposeLocalEnergyDeposit(localEnergyDeposit);
 
  if (fVerboseLevel > 1)
    {
      G4cout << "-----------------------------------------------------------" << G4endl;
      G4cout << "Energy balance from G4PenelopeIonisation" << G4endl;
      G4cout << "Incoming primary energy: " << kineticEnergy0/keV << " keV" << G4endl;
      G4cout << "-----------------------------------------------------------" << G4endl;
      G4cout << "Outgoing primary energy: " << fKineticEnergy1/keV << " keV" << G4endl;
      G4cout << "Delta ray " << fEnergySecondary/keV << " keV" << G4endl;
      if (energyInFluorescence)
    G4cout << "Fluorescence x-rays: " << energyInFluorescence/keV << " keV" << G4endl;
      if (energyInAuger)
    G4cout << "Auger electrons: " << energyInAuger/keV << " keV" << G4endl;     
      G4cout << "Local energy deposit " << localEnergyDeposit/keV << " keV" << G4endl;
      G4cout << "Total final state: " << (fEnergySecondary+energyInFluorescence+fKineticEnergy1+
                      localEnergyDeposit+energyInAuger)/keV <<
    " keV" << G4endl;
      G4cout << "-----------------------------------------------------------" << G4endl;
    }
      
  if (fVerboseLevel > 0)
    {
      G4double energyDiff = std::fabs(fEnergySecondary+energyInFluorescence+fKineticEnergy1+
                      localEnergyDeposit+energyInAuger-kineticEnergy0);
      if (energyDiff > 0.05*keV)
    G4cout << "Warning from G4PenelopeIonisation: problem with energy conservation: " <<  
      (fEnergySecondary+energyInFluorescence+fKineticEnergy1+localEnergyDeposit+energyInAuger)/keV <<
      " keV (final) vs. " <<
      kineticEnergy0/keV << " keV (initial)" << G4endl;      
    }
    
}

References G4AtomicShell::BindingEnergy(), G4InuclSpecialFunctions::bindingEnergy(), G4VAtomDeexcitation::CheckDeexcitationActiveRegion(), G4Gamma::Definition(), G4Electron::Definition(), G4Electron::Electron(), G4InuclParticleNames::electron, eV, FatalException, fAtomDeexcitation, fCosThetaPrimary, fCosThetaSecondary, fEnergySecondary, fIntrinsicLowEnergyLimit, fKineticEnergy1, fOscManager, fParticleChange, fPIXEflag, fTargetOscillator, fVerboseLevel, G4cout, G4endl, G4Exception(), G4UniformRand, G4VAtomDeexcitation::GenerateParticles(), G4DynamicParticle::GetDefinition(), G4MaterialCutsCouple::GetIndex(), G4DynamicParticle::GetKineticEnergy(), G4MaterialCutsCouple::GetMaterial(), G4DynamicParticle::GetMomentumDirection(), G4PenelopeOscillatorManager::GetOscillatorTableIonisation(), G4ParticleDefinition::GetParticleName(), G4AtomicTransitionManager::Instance(), JustWarning, keV, eplot::material, pi, G4Positron::Positron(), G4VParticleChange::ProposeLocalEnergyDeposit(), G4ParticleChangeForLoss::ProposeMomentumDirection(), CLHEP::Hep3Vector::rotateUz(), SampleFinalStateElectron(), SampleFinalStatePositron(), G4ParticleChangeForLoss::SetProposedKineticEnergy(), G4AtomicTransitionManager::Shell(), twopi, and Z.

◆ SecondaryThreshold()

G4double G4VEmModel::SecondaryThreshold ( ) const

inlineinherited

Definition at line 690 of file G4VEmModel.hh.

{
  return secondaryThreshold;
}

References G4VEmModel::secondaryThreshold.

Referenced by G4LivermoreBremsstrahlungModel::SampleSecondaries(), G4eBremParametrizedModel::SampleSecondaries(), G4eBremsstrahlungRelModel::SampleSecondaries(), G4SeltzerBergerModel::SampleSecondaries(), G4MuBremsstrahlungModel::SampleSecondaries(), and G4MuPairProductionModel::SampleSecondaries().

◆ SelectIsotopeNumber()

G4int G4VEmModel::SelectIsotopeNumber ( const G4Element * elm )

inherited

Definition at line 319 of file G4VEmModel.cc.

{
  SetCurrentElement(elm);
  const size_t ni = elm->GetNumberOfIsotopes();
  fCurrentIsotope = elm->GetIsotope(0);
  size_t idx = 0;
  if(ni > 1) {
    const G4double* ab = elm->GetRelativeAbundanceVector();
    G4double x = G4UniformRand();
    for(; idx<ni; ++idx) {
      x -= ab[idx];
      if (x <= 0.0) { 
    fCurrentIsotope = elm->GetIsotope(idx);
    break; 
      }
    }
  }
  return fCurrentIsotope->GetN();
}

References ab, G4VEmModel::fCurrentIsotope, G4UniformRand, G4Element::GetIsotope(), G4Isotope::GetN(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), and G4VEmModel::SetCurrentElement().

Referenced by G4eSingleCoulombScatteringModel::SampleSecondaries(), G4IonCoulombScatteringModel::SampleSecondaries(), G4eCoulombScatteringModel::SampleSecondaries(), G4hCoulombScatteringModel::SampleSecondaries(), and G4BetheHeitler5DModel::SampleSecondaries().

◆ SelectRandomAtom() [1/2]

const G4Element * G4VEmModel::SelectRandomAtom	(	const G4Material *	mat,
		const G4ParticleDefinition *	pd,
		G4double	kineticEnergy,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

inherited

Definition at line 275 of file G4VEmModel.cc.

{
  size_t n = mat->GetNumberOfElements();
  fCurrentElement = mat->GetElement(0);
  if (n > 1) {
    const G4double x = G4UniformRand()*
      G4VEmModel::CrossSectionPerVolume(mat,pd,kinEnergy,tcut,tmax);
    for(size_t i=0; i<n; ++i) {
      if (x <= xsec[i]) {
        fCurrentElement = mat->GetElement(i);
        break;
      }
    }
  }
  return fCurrentElement;
}

References G4VEmModel::CrossSectionPerVolume(), G4VEmModel::fCurrentElement, G4UniformRand, G4Material::GetElement(), G4Material::GetNumberOfElements(), CLHEP::detail::n, and G4VEmModel::xsec.

◆ SelectRandomAtom() [2/2]

const G4Element * G4VEmModel::SelectRandomAtom	(	const G4MaterialCutsCouple *	couple,
		const G4ParticleDefinition *	part,
		G4double	kineticEnergy,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

inlineinherited

Definition at line 580 of file G4VEmModel.hh.

{
  SetCurrentCouple(couple);
  fCurrentElement = (nSelectors > 0) ?
    ((*elmSelectors)[couple->GetIndex()])->SelectRandomAtom(kinEnergy) :
    SelectRandomAtom(pBaseMaterial,part,kinEnergy,cutEnergy,maxEnergy);
  fCurrentIsotope = nullptr;
  return fCurrentElement;
}

References G4VEmModel::elmSelectors, G4VEmModel::fCurrentElement, G4VEmModel::fCurrentIsotope, G4MaterialCutsCouple::GetIndex(), G4VEmModel::nSelectors, G4VEmModel::pBaseMaterial, G4VEmModel::SelectRandomAtom(), and G4VEmModel::SetCurrentCouple().

◆ SelectRandomAtomNumber()

G4int G4VEmModel::SelectRandomAtomNumber ( const G4Material * mat )

inherited

Definition at line 297 of file G4VEmModel.cc.

{
  // this algorith assumes that cross section is proportional to
  // number electrons multiplied by number of atoms
  const size_t nn = mat->GetNumberOfElements();
  fCurrentElement = mat->GetElement(0);
  if(1 < nn) {
    const G4double* at = mat->GetVecNbOfAtomsPerVolume();
    G4double tot = mat->GetTotNbOfAtomsPerVolume()*G4UniformRand();
    for(size_t i=0; i<nn; ++i) {
      tot -= at[i];
      if(tot <= 0.0) { 
    fCurrentElement = mat->GetElement(i);
    break; 
      }
    }
  }
  return fCurrentElement->GetZasInt();
}

References G4VEmModel::fCurrentElement, G4UniformRand, G4Material::GetElement(), G4Material::GetNumberOfElements(), G4Material::GetTotNbOfAtomsPerVolume(), G4Material::GetVecNbOfAtomsPerVolume(), G4Element::GetZasInt(), and G4InuclParticleNames::nn.

Referenced by G4AtimaEnergyLossModel::SampleSecondaries(), G4BetheBlochModel::SampleSecondaries(), G4BraggIonModel::SampleSecondaries(), G4BraggModel::SampleSecondaries(), G4ICRU73QOModel::SampleSecondaries(), G4LindhardSorensenIonModel::SampleSecondaries(), G4MollerBhabhaModel::SampleSecondaries(), and G4IonParametrisedLossModel::SampleSecondaries().

◆ SelectTargetAtom()

const G4Element * G4VEmModel::SelectTargetAtom	(	const G4MaterialCutsCouple *	couple,
		const G4ParticleDefinition *	part,
		G4double	kineticEnergy,
		G4double	logKineticEnergy,
		G4double	cutEnergy = `0.0`,
		G4double	maxEnergy = `DBL_MAX`
	)

inlineinherited

Definition at line 597 of file G4VEmModel.hh.

{
  SetCurrentCouple(couple);
  fCurrentElement = (nSelectors > 0)
   ? ((*elmSelectors)[couple->GetIndex()])->SelectRandomAtom(kinEnergy,logKinE)
   : SelectRandomAtom(pBaseMaterial,part,kinEnergy,cutEnergy,maxEnergy);
  fCurrentIsotope = nullptr;
  return fCurrentElement;
}

References G4VEmModel::elmSelectors, G4VEmModel::fCurrentElement, G4VEmModel::fCurrentIsotope, G4MaterialCutsCouple::GetIndex(), G4VEmModel::nSelectors, G4VEmModel::pBaseMaterial, G4VEmModel::SelectRandomAtom(), and G4VEmModel::SetCurrentCouple().

◆ SetActivationHighEnergyLimit()

void G4VEmModel::SetActivationHighEnergyLimit ( G4double val )

inlineinherited

Definition at line 781 of file G4VEmModel.hh.

{
  eMaxActive = val;
}

References G4VEmModel::eMaxActive.

Referenced by G4EmModelActivator::ActivateMicroElec(), G4EmDNAPhysics_option7::ConstructProcess(), G4EmDNAPhysics_option8::ConstructProcess(), G4NuclearStopping::InitialiseProcess(), and G4EmConfigurator::SetExtraEmModel().

◆ SetActivationLowEnergyLimit()

void G4VEmModel::SetActivationLowEnergyLimit ( G4double val )

inlineinherited

Definition at line 788 of file G4VEmModel.hh.

{
  eMinActive = val;
}

References G4VEmModel::eMinActive.

◆ SetAngularDistribution()

void G4VEmModel::SetAngularDistribution ( G4VEmAngularDistribution * p )

inlineinherited

Definition at line 628 of file G4VEmModel.hh.

{
  if(p != anglModel) {
    delete anglModel;
    anglModel = p;
  }
}

References G4VEmModel::anglModel.

Referenced by G4EmLivermorePhysics::ConstructProcess(), G4EmLowEPPhysics::ConstructProcess(), G4EmStandardPhysics::ConstructProcess(), G4EmStandardPhysics_option3::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), G4EmStandardPhysicsSS::ConstructProcess(), G4BetheHeitlerModel::G4BetheHeitlerModel(), G4DNABornIonisationModel1::G4DNABornIonisationModel1(), G4DNABornIonisationModel2::G4DNABornIonisationModel2(), G4DNAEmfietzoglouIonisationModel::G4DNAEmfietzoglouIonisationModel(), G4DNARuddIonisationExtendedModel::G4DNARuddIonisationExtendedModel(), G4DNARuddIonisationModel::G4DNARuddIonisationModel(), G4eBremParametrizedModel::G4eBremParametrizedModel(), G4eBremsstrahlungRelModel::G4eBremsstrahlungRelModel(), G4IonParametrisedLossModel::G4IonParametrisedLossModel(), G4LivermoreBremsstrahlungModel::G4LivermoreBremsstrahlungModel(), G4LivermoreIonisationModel::G4LivermoreIonisationModel(), G4LivermorePhotoElectricModel::G4LivermorePhotoElectricModel(), G4LivermoreRayleighModel::G4LivermoreRayleighModel(), G4MicroElecInelasticModel::G4MicroElecInelasticModel(), G4MicroElecInelasticModel_new::G4MicroElecInelasticModel_new(), G4MuBremsstrahlungModel::G4MuBremsstrahlungModel(), G4MuPairProductionModel::G4MuPairProductionModel(), G4PAIModel::G4PAIModel(), G4PAIPhotModel::G4PAIPhotModel(), G4PairProductionRelModel::G4PairProductionRelModel(), G4PEEffectFluoModel::G4PEEffectFluoModel(), G4SeltzerBergerModel::G4SeltzerBergerModel(), G4AtimaEnergyLossModel::Initialise(), G4BetheBlochModel::Initialise(), G4BraggIonModel::Initialise(), G4BraggModel::Initialise(), G4ICRU73QOModel::Initialise(), G4LindhardSorensenIonModel::Initialise(), and G4MollerBhabhaModel::Initialise().

◆ SetAngularGeneratorFlag()

void G4VEmModel::SetAngularGeneratorFlag ( G4bool val )

inlineinherited

Definition at line 725 of file G4VEmModel.hh.

{
  useAngularGenerator = val;
}

References G4VEmModel::useAngularGenerator.

Referenced by G4VEnergyLossProcess::PreparePhysicsTable().

◆ SetCrossSectionTable()

void G4VEmModel::SetCrossSectionTable	(	G4PhysicsTable *	p,
		G4bool	isLocal
	)

inherited

Definition at line 455 of file G4VEmModel.cc.

{
  if(p != xSectionTable) {
    if(xSectionTable != nullptr && localTable) { 
      xSectionTable->clearAndDestroy(); 
      delete xSectionTable;
    }
    xSectionTable = p;
  }
  localTable = isLocal;
}

References G4PhysicsTable::clearAndDestroy(), G4VEmModel::localTable, and G4VEmModel::xSectionTable.

Referenced by G4VMultipleScattering::BuildPhysicsTable().

◆ SetCurrentCouple()

void G4VEmModel::SetCurrentCouple ( const G4MaterialCutsCouple * ptr )

inlineinherited

Definition at line 472 of file G4VEmModel.hh.

{
  if(fCurrentCouple != ptr) {
    fCurrentCouple = ptr;
    basedCoupleIndex = currentCoupleIndex = ptr->GetIndex();
    pBaseMaterial = ptr->GetMaterial();
    pFactor = 1.0;
    if(useBaseMaterials) {
      basedCoupleIndex = (*theDensityIdx)[currentCoupleIndex];
      if(nullptr != pBaseMaterial->GetBaseMaterial()) 
    pBaseMaterial = pBaseMaterial->GetBaseMaterial();
      pFactor = (*theDensityFactor)[currentCoupleIndex];
    }
  }
}

References G4VEmModel::basedCoupleIndex, G4VEmModel::currentCoupleIndex, G4VEmModel::fCurrentCouple, G4Material::GetBaseMaterial(), G4MaterialCutsCouple::GetIndex(), G4MaterialCutsCouple::GetMaterial(), G4VEmModel::pBaseMaterial, G4VEmModel::pFactor, and G4VEmModel::useBaseMaterials.

Referenced by G4VMultipleScattering::AlongStepGetPhysicalInteractionLength(), G4EmMultiModel::ComputeCrossSectionPerAtom(), G4VEmModel::ComputeDEDX(), G4TablesForExtrapolator::ComputeTrasportXS(), G4UrbanAdjointMscModel::ComputeTruePathLengthLimit(), G4GoudsmitSaundersonMscModel::ComputeTruePathLengthLimit(), G4UrbanMscModel::ComputeTruePathLengthLimit(), G4VEmModel::CrossSection(), G4AdjointPhotoElectricModel::DefineCurrentMaterialAndElectronEnergy(), G4WentzelVIModel::DefineMaterial(), G4WentzelVIRelModel::DefineMaterial(), G4EmCalculator::GetCrossSectionPerVolume(), G4LivermoreGammaConversion5DModel::Initialise(), G4VEmModel::InitialiseElementSelectors(), G4PEEffectFluoModel::SampleSecondaries(), G4EmMultiModel::SampleSecondaries(), G4VEnergyLossProcess::SelectModel(), G4VEmProcess::SelectModel(), G4VEmModel::SelectRandomAtom(), G4VEmModel::SelectTargetAtom(), and G4VEmModel::Value().

◆ SetCurrentElement()

void G4VEmModel::SetCurrentElement ( const G4Element * elm )

inlineprotectedinherited

Definition at line 497 of file G4VEmModel.hh.

{
  fCurrentElement = elm;
  fCurrentIsotope = nullptr;
}

References G4VEmModel::fCurrentElement, and G4VEmModel::fCurrentIsotope.

Referenced by G4VEmModel::ComputeCrossSectionPerAtom(), G4eBremsstrahlungRelModel::ComputeDEDXPerVolume(), G4eBremParametrizedModel::ComputeDEDXPerVolume(), and G4VEmModel::SelectIsotopeNumber().

◆ SetDeexcitationFlag()

void G4VEmModel::SetDeexcitationFlag ( G4bool val )

inlineinherited

Definition at line 823 of file G4VEmModel.hh.

{
  flagDeexcitation = val;
}

References G4VEmModel::flagDeexcitation.

Referenced by G4DNABornIonisationModel1::G4DNABornIonisationModel1(), G4DNABornIonisationModel2::G4DNABornIonisationModel2(), G4DNACPA100IonisationModel::G4DNACPA100IonisationModel(), G4DNAEmfietzoglouIonisationModel::G4DNAEmfietzoglouIonisationModel(), G4DNARelativisticIonisationModel::G4DNARelativisticIonisationModel(), G4DNARuddIonisationExtendedModel::G4DNARuddIonisationExtendedModel(), G4DNARuddIonisationModel::G4DNARuddIonisationModel(), G4KleinNishinaModel::G4KleinNishinaModel(), G4LEPTSIonisationModel::G4LEPTSIonisationModel(), G4LivermoreComptonModel::G4LivermoreComptonModel(), G4LivermorePhotoElectricModel::G4LivermorePhotoElectricModel(), G4LivermorePolarizedComptonModel::G4LivermorePolarizedComptonModel(), G4LowEPComptonModel::G4LowEPComptonModel(), G4LowEPPolarizedComptonModel::G4LowEPPolarizedComptonModel(), G4MicroElecInelasticModel::G4MicroElecInelasticModel(), G4MicroElecInelasticModel_new::G4MicroElecInelasticModel_new(), G4PEEffectFluoModel::G4PEEffectFluoModel(), G4PenelopeBremsstrahlungModel::G4PenelopeBremsstrahlungModel(), G4PenelopeComptonModel::G4PenelopeComptonModel(), G4PenelopeIonisationModel(), G4PenelopePhotoElectricModel::G4PenelopePhotoElectricModel(), G4AtimaEnergyLossModel::Initialise(), G4BetheBlochModel::Initialise(), G4BraggIonModel::Initialise(), G4BraggModel::Initialise(), G4ICRU73QOModel::Initialise(), and G4LindhardSorensenIonModel::Initialise().

◆ SetElementSelectors()

void G4VEmModel::SetElementSelectors ( std::vector< G4EmElementSelector * > * p )

inlineinherited

Definition at line 852 of file G4VEmModel.hh.

{
  if(p != elmSelectors) {
    elmSelectors = p;
    nSelectors = (nullptr != elmSelectors) ? G4int(elmSelectors->size()) : 0;
    localElmSelectors = false;
  }
}

References G4VEmModel::elmSelectors, G4VEmModel::localElmSelectors, and G4VEmModel::nSelectors.

◆ SetFluctuationFlag()

void G4VEmModel::SetFluctuationFlag ( G4bool val )

inlineinherited

Definition at line 732 of file G4VEmModel.hh.

{
  lossFlucFlag = val;
}

References G4VEmModel::lossFlucFlag.

Referenced by G4EmCalculator::ComputeNuclearDEDX().

◆ SetForceBuildTable()

void G4VEmModel::SetForceBuildTable ( G4bool val )

inlineinherited

Definition at line 830 of file G4VEmModel.hh.

{
  flagForceBuildTable = val;
}

References G4VEmModel::flagForceBuildTable.

◆ SetHighEnergyLimit()

void G4VEmModel::SetHighEnergyLimit ( G4double val )

inlineinherited

Definition at line 767 of file G4VEmModel.hh.

{
  highLimit = val;
}

References G4VEmModel::highLimit.

◆ SetLocked()

void G4VEmModel::SetLocked ( G4bool val )

inlineinherited

Definition at line 884 of file G4VEmModel.hh.

{
  isLocked = val;
}

References G4VEmModel::isLocked.

Referenced by G4EmModelActivator::ActivateEmOptions(), G4EmModelActivator::AddStandardScattering(), and G4EmModelActivator::SetMscParameters().

◆ SetLowEnergyLimit()

void G4VEmModel::SetLowEnergyLimit ( G4double val )

inlineinherited

Definition at line 774 of file G4VEmModel.hh.

{
  lowLimit = val;
}

References G4VEmModel::lowLimit.

◆ SetLPMFlag()

void G4VEmModel::SetLPMFlag ( G4bool val )

inlineinherited

Definition at line 816 of file G4VEmModel.hh.

{
  theLPMflag = val;
}

References G4VEmModel::theLPMflag.

Referenced by G4eBremsstrahlungRelModel::G4eBremsstrahlungRelModel(), G4LivermoreBremsstrahlungModel::G4LivermoreBremsstrahlungModel(), G4SeltzerBergerModel::G4SeltzerBergerModel(), and G4eBremsstrahlung::InitialiseEnergyLossProcess().

◆ SetMasterThread()

void G4VEmModel::SetMasterThread ( G4bool val )

inlineinherited

Definition at line 739 of file G4VEmModel.hh.

{
  isMaster = val;
}

References G4VEmModel::isMaster.

Referenced by G4VEmProcess::PreparePhysicsTable(), G4VEnergyLossProcess::PreparePhysicsTable(), and G4VMultipleScattering::PreparePhysicsTable().

◆ SetParticle()

void G4PenelopeIonisationModel::SetParticle ( const G4ParticleDefinition * p )

private

Definition at line 1106 of file G4PenelopeIonisationModel.cc.

{
  if(!fParticle) {
    fParticle = p;  
  }
}

References fParticle.

Referenced by G4PenelopeIonisationModel(), and Initialise().

◆ SetParticleChange()

void G4VEmModel::SetParticleChange	(	G4VParticleChange *	p,
		G4VEmFluctuationModel *	f = `nullptr`
	)

inherited

Definition at line 447 of file G4VEmModel.cc.

{
  if(p != nullptr && pParticleChange != p) { pParticleChange = p; }
  if(flucModel != f) { flucModel = f; }
}

References G4VEmModel::flucModel, and G4VEmModel::pParticleChange.

Referenced by G4VMultipleScattering::AddEmModel(), G4VEmProcess::AddEmModel(), G4VEnergyLossProcess::AddEmModel(), G4VEmModel::GetParticleChangeForGamma(), G4VEmModel::GetParticleChangeForLoss(), G4EmMultiModel::Initialise(), G4IonParametrisedLossModel::Initialise(), G4DNADummyModel::Initialise(), and G4NuclearStopping::InitialiseProcess().

◆ SetPolarAngleLimit()

void G4VEmModel::SetPolarAngleLimit ( G4double val )

inlineinherited

Definition at line 802 of file G4VEmModel.hh.

{
  if(!isLocked) { polarAngleLimit = val; }
}

References G4VEmModel::isLocked, and G4VEmModel::polarAngleLimit.

Referenced by G4EmModelActivator::ActivateEmOptions(), G4TablesForExtrapolator::ComputeTrasportXS(), G4CoulombScattering::InitialiseProcess(), G4VEmProcess::PreparePhysicsTable(), and G4VMultipleScattering::PreparePhysicsTable().

◆ SetSecondaryThreshold()

void G4VEmModel::SetSecondaryThreshold ( G4double val )

inlineinherited

Definition at line 809 of file G4VEmModel.hh.

{
  secondaryThreshold = val;
}

References G4VEmModel::secondaryThreshold.

Referenced by G4MuBremsstrahlung::InitialiseEnergyLossProcess(), G4MuPairProduction::InitialiseEnergyLossProcess(), and G4eBremsstrahlung::InitialiseEnergyLossProcess().

◆ SetTripletModel()

void G4VEmModel::SetTripletModel ( G4VEmModel * p )

inlineinherited

Definition at line 645 of file G4VEmModel.hh.

{
  if(p != fTripletModel) {
    delete fTripletModel;
    fTripletModel = p;
  }
}

References G4VEmModel::fTripletModel.

Referenced by G4eplusTo2GammaOKVIModel::G4eplusTo2GammaOKVIModel().

◆ SetupForMaterial()

void G4VEmModel::SetupForMaterial	(	const G4ParticleDefinition *	,
		const G4Material *	,
		G4double	kineticEnergy
	)

virtualinherited

Reimplemented in G4eBremParametrizedModel, G4eBremsstrahlungRelModel, G4PairProductionRelModel, and G4SeltzerBergerModel.

Definition at line 440 of file G4VEmModel.cc.

442{}

Referenced by G4VEmModel::CrossSectionPerVolume(), G4PenelopeComptonModel::CrossSectionPerVolume(), G4PenelopeBremsstrahlungModel::CrossSectionPerVolume(), CrossSectionPerVolume(), G4EmCalculator::FindEmModel(), and G4EmElementSelector::Initialise().

◆ SetUseBaseMaterials()

void G4VEmModel::SetUseBaseMaterials ( G4bool val )

inlineinherited

Definition at line 753 of file G4VEmModel.hh.

{
  useBaseMaterials = val;
}

References G4VEmModel::useBaseMaterials.

Referenced by G4VEmProcess::BuildPhysicsTable(), G4VEnergyLossProcess::BuildPhysicsTable(), G4VMultipleScattering::BuildPhysicsTable(), G4TablesForExtrapolator::ComputeElectronDEDX(), G4TablesForExtrapolator::ComputeMuonDEDX(), G4TablesForExtrapolator::ComputeProtonDEDX(), G4TablesForExtrapolator::ComputeTrasportXS(), G4VEmProcess::PreparePhysicsTable(), G4VEnergyLossProcess::PreparePhysicsTable(), and G4VMultipleScattering::PreparePhysicsTable().

◆ SetVerbosityLevel()

void G4PenelopeIonisationModel::SetVerbosityLevel ( G4int lev )

inline

Definition at line 108 of file G4PenelopeIonisationModel.hh.

108{fVerboseLevel = lev;};

References fVerboseLevel.

◆ StartTracking()

void G4VEmModel::StartTracking ( G4Track * )

virtualinherited

Reimplemented in G4UrbanAdjointMscModel, G4GoudsmitSaundersonMscModel, G4UrbanMscModel, and G4WentzelVIModel.

Definition at line 270 of file G4VEmModel.cc.

271{}

Referenced by G4VMultipleScattering::StartTracking().

◆ UseAngularGeneratorFlag()

G4bool G4VEmModel::UseAngularGeneratorFlag ( ) const

inlineinherited

Definition at line 718 of file G4VEmModel.hh.

{
  return useAngularGenerator;
}

References G4VEmModel::useAngularGenerator.

Referenced by G4AtimaEnergyLossModel::Initialise(), G4BetheBlochModel::Initialise(), G4BraggIonModel::Initialise(), G4BraggModel::Initialise(), G4ICRU73QOModel::Initialise(), G4LindhardSorensenIonModel::Initialise(), G4MollerBhabhaModel::Initialise(), G4AtimaEnergyLossModel::SampleSecondaries(), G4BetheBlochModel::SampleSecondaries(), G4BraggIonModel::SampleSecondaries(), G4BraggModel::SampleSecondaries(), G4ICRU73QOModel::SampleSecondaries(), G4LindhardSorensenIonModel::SampleSecondaries(), and G4MollerBhabhaModel::SampleSecondaries().

◆ UseBaseMaterials()

G4bool G4VEmModel::UseBaseMaterials ( ) const

inlineinherited

Definition at line 760 of file G4VEmModel.hh.

{
  return useBaseMaterials;
}

References G4VEmModel::useBaseMaterials.

◆ Value()

G4double G4VEmModel::Value	(	const G4MaterialCutsCouple *	couple,
		const G4ParticleDefinition *	p,
		G4double	kineticEnergy
	)

virtualinherited

Definition at line 406 of file G4VEmModel.cc.

{
  SetCurrentCouple(couple);
  return pFactor*e*e*CrossSectionPerVolume(pBaseMaterial,p,e,0.0,DBL_MAX);
}

References G4VEmModel::CrossSectionPerVolume(), DBL_MAX, G4VEmModel::pBaseMaterial, G4VEmModel::pFactor, and G4VEmModel::SetCurrentCouple().

Referenced by G4IonParametrisedLossModel::BuildRangeVector(), G4LossTableBuilder::BuildTableForModel(), G4LivermorePhotoElectricModel::ComputeCrossSectionPerAtom(), G4IonParametrisedLossModel::ComputeLossForStep(), G4VLEPTSModel::GetMeanFreePath(), G4DNABornExcitationModel2::GetPartialCrossSection(), G4DNABornExcitationModel2::Initialise(), G4DNABornExcitationModel2::RandomSelect(), G4SeltzerBergerModel::SampleEnergyTransfer(), G4LivermoreBremsstrahlungModel::SampleSecondaries(), and G4LivermorePhotoElectricModel::SampleSecondaries().

Field Documentation

◆ anglModel

G4VEmAngularDistribution* G4VEmModel::anglModel = nullptr

privateinherited

Definition at line 414 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetAngularDistribution(), G4VEmModel::SetAngularDistribution(), and G4VEmModel::~G4VEmModel().

◆ basedCoupleIndex

size_t G4VEmModel::basedCoupleIndex = 0

protectedinherited

Definition at line 449 of file G4VEmModel.hh.

Referenced by G4VMscModel::GetTransportMeanFreePath(), and G4VEmModel::SetCurrentCouple().

◆ currentCoupleIndex

size_t G4VEmModel::currentCoupleIndex = 0

protectedinherited

Definition at line 448 of file G4VEmModel.hh.

Referenced by G4VEmModel::SetCurrentCouple().

◆ elmSelectors

std::vector<G4EmElementSelector*>* G4VEmModel::elmSelectors = nullptr

privateinherited

Definition at line 419 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetElementSelectors(), G4VEmModel::InitialiseElementSelectors(), G4VEmModel::SelectRandomAtom(), G4VEmModel::SelectTargetAtom(), G4VEmModel::SetElementSelectors(), and G4VEmModel::~G4VEmModel().

◆ eMaxActive

G4double G4VEmModel::eMaxActive = DBL_MAX

privateinherited

Definition at line 439 of file G4VEmModel.hh.

Referenced by G4VEmModel::HighEnergyActivationLimit(), G4VEmModel::IsActive(), and G4VEmModel::SetActivationHighEnergyLimit().

◆ eMinActive

G4double G4VEmModel::eMinActive = 0.0

privateinherited

Definition at line 438 of file G4VEmModel.hh.

Referenced by G4VEmModel::IsActive(), G4VEmModel::LowEnergyActivationLimit(), and G4VEmModel::SetActivationLowEnergyLimit().

◆ fAtomDeexcitation

G4VAtomDeexcitation* G4PenelopeIonisationModel::fAtomDeexcitation

private

Definition at line 129 of file G4PenelopeIonisationModel.hh.

Referenced by Initialise(), and SampleSecondaries().

◆ fCosThetaPrimary

G4double G4PenelopeIonisationModel::fCosThetaPrimary

private

Definition at line 132 of file G4PenelopeIonisationModel.hh.

Referenced by SampleFinalStateElectron(), SampleFinalStatePositron(), and SampleSecondaries().

◆ fCosThetaSecondary

G4double G4PenelopeIonisationModel::fCosThetaSecondary

private

Definition at line 134 of file G4PenelopeIonisationModel.hh.

Referenced by SampleFinalStateElectron(), SampleFinalStatePositron(), and SampleSecondaries().

◆ fCrossSectionHandler

G4PenelopeIonisationXSHandler* G4PenelopeIonisationModel::fCrossSectionHandler

private

Definition at line 128 of file G4PenelopeIonisationModel.hh.

Referenced by ComputeDEDXPerVolume(), CrossSectionPerVolume(), Initialise(), InitialiseLocal(), SampleFinalStateElectron(), SampleFinalStatePositron(), and ~G4PenelopeIonisationModel().

◆ fCurrentCouple

const G4MaterialCutsCouple* G4VEmModel::fCurrentCouple = nullptr

privateinherited

Definition at line 416 of file G4VEmModel.hh.

Referenced by G4VEmModel::CurrentCouple(), and G4VEmModel::SetCurrentCouple().

◆ fCurrentElement

const G4Element* G4VEmModel::fCurrentElement = nullptr

privateinherited

Definition at line 417 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetCurrentElement(), G4VEmModel::SelectRandomAtom(), G4VEmModel::SelectRandomAtomNumber(), G4VEmModel::SelectTargetAtom(), and G4VEmModel::SetCurrentElement().

◆ fCurrentIsotope

const G4Isotope* G4VEmModel::fCurrentIsotope = nullptr

privateinherited

Definition at line 418 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetCurrentIsotope(), G4VEmModel::SelectIsotopeNumber(), G4VEmModel::SelectRandomAtom(), G4VEmModel::SelectTargetAtom(), and G4VEmModel::SetCurrentElement().

◆ fElementData

G4ElementData* G4VEmModel::fElementData = nullptr

protectedinherited

Definition at line 424 of file G4VEmModel.hh.

Referenced by G4MuPairProductionModel::FindScaledEnergy(), G4VEmModel::GetElementData(), G4MuPairProductionModel::Initialise(), G4MuPairProductionModel::InitialiseLocal(), G4MuPairProductionModel::MakeSamplingTables(), G4MuPairProductionModel::RetrieveTables(), G4MuPairProductionModel::StoreTables(), and G4VEmModel::~G4VEmModel().

◆ fEmManager

G4LossTableManager* G4VEmModel::fEmManager

privateinherited

Definition at line 420 of file G4VEmModel.hh.

Referenced by G4VEmModel::G4VEmModel(), and G4VEmModel::~G4VEmModel().

◆ fEnergySecondary

G4double G4PenelopeIonisationModel::fEnergySecondary

private

Definition at line 133 of file G4PenelopeIonisationModel.hh.

Referenced by SampleFinalStateElectron(), SampleFinalStatePositron(), and SampleSecondaries().

◆ fIntrinsicHighEnergyLimit

G4double G4PenelopeIonisationModel::fIntrinsicHighEnergyLimit

private

Definition at line 138 of file G4PenelopeIonisationModel.hh.

Referenced by G4PenelopeIonisationModel().

◆ fIntrinsicLowEnergyLimit

G4double G4PenelopeIonisationModel::fIntrinsicLowEnergyLimit

private

Definition at line 137 of file G4PenelopeIonisationModel.hh.

Referenced by G4PenelopeIonisationModel(), MinEnergyCut(), and SampleSecondaries().

◆ fIsInitialised

G4bool G4PenelopeIonisationModel::fIsInitialised

private

Definition at line 143 of file G4PenelopeIonisationModel.hh.

Referenced by Initialise().

◆ fKineticEnergy1

G4double G4PenelopeIonisationModel::fKineticEnergy1

private

Definition at line 131 of file G4PenelopeIonisationModel.hh.

Referenced by SampleFinalStateElectron(), SampleFinalStatePositron(), and SampleSecondaries().

◆ flagDeexcitation

G4bool G4VEmModel::flagDeexcitation = false

privateinherited

Definition at line 455 of file G4VEmModel.hh.

Referenced by G4VEmModel::DeexcitationFlag(), and G4VEmModel::SetDeexcitationFlag().

◆ flagForceBuildTable

G4bool G4VEmModel::flagForceBuildTable = false

privateinherited

Definition at line 456 of file G4VEmModel.hh.

Referenced by G4VEmModel::ForceBuildTableFlag(), and G4VEmModel::SetForceBuildTable().

◆ fLocalTable

G4bool G4PenelopeIonisationModel::fLocalTable

private

Definition at line 146 of file G4PenelopeIonisationModel.hh.

Referenced by ComputeDEDXPerVolume(), CrossSectionPerVolume(), and ~G4PenelopeIonisationModel().

◆ flucModel

G4VEmFluctuationModel* G4VEmModel::flucModel = nullptr

privateinherited

Definition at line 413 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetModelOfFluctuations(), and G4VEmModel::SetParticleChange().

◆ fNBins

size_t G4PenelopeIonisationModel::fNBins

private

Definition at line 142 of file G4PenelopeIonisationModel.hh.

Referenced by ComputeDEDXPerVolume(), CrossSectionPerVolume(), G4PenelopeIonisationModel(), Initialise(), and InitialiseLocal().

◆ fOscManager

G4PenelopeOscillatorManager* G4PenelopeIonisationModel::fOscManager

private

Definition at line 127 of file G4PenelopeIonisationModel.hh.

Referenced by ComputeDEDXPerVolume(), CrossSectionPerVolume(), G4PenelopeIonisationModel(), SampleFinalStateElectron(), SampleFinalStatePositron(), and SampleSecondaries().

◆ fParticle

const G4ParticleDefinition* G4PenelopeIonisationModel::fParticle

protected

Definition at line 116 of file G4PenelopeIonisationModel.hh.

Referenced by Initialise(), InitialiseLocal(), and SetParticle().

◆ fParticleChange

G4ParticleChangeForLoss* G4PenelopeIonisationModel::fParticleChange

protected

Definition at line 115 of file G4PenelopeIonisationModel.hh.

Referenced by Initialise(), and SampleSecondaries().

◆ fPIXEflag

G4bool G4PenelopeIonisationModel::fPIXEflag

private

Definition at line 144 of file G4PenelopeIonisationModel.hh.

Referenced by Initialise(), and SampleSecondaries().

◆ fTargetOscillator

G4int G4PenelopeIonisationModel::fTargetOscillator

private

Definition at line 141 of file G4PenelopeIonisationModel.hh.

Referenced by SampleFinalStateElectron(), SampleFinalStatePositron(), and SampleSecondaries().

◆ fTripletModel

G4VEmModel* G4VEmModel::fTripletModel = nullptr

privateinherited

Definition at line 415 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetParticleChangeForGamma(), G4VEmModel::GetParticleChangeForLoss(), G4VEmModel::GetTripletModel(), and G4VEmModel::SetTripletModel().

◆ fVerboseLevel

G4int G4PenelopeIonisationModel::fVerboseLevel

private

Definition at line 140 of file G4PenelopeIonisationModel.hh.

Referenced by ComputeDEDXPerVolume(), CrossSectionPerVolume(), G4PenelopeIonisationModel(), GetVerbosityLevel(), Initialise(), InitialiseLocal(), SampleFinalStateElectron(), SampleFinalStatePositron(), SampleSecondaries(), and SetVerbosityLevel().

◆ highLimit

G4double G4VEmModel::highLimit

privateinherited

Definition at line 437 of file G4VEmModel.hh.

Referenced by G4VEmModel::HighEnergyLimit(), G4VEmModel::InitialiseElementSelectors(), and G4VEmModel::SetHighEnergyLimit().

◆ inveplus

G4double G4VEmModel::inveplus

protectedinherited

Definition at line 431 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetChargeSquareRatio(), G4VMscModel::GetDEDX(), G4VMscModel::GetEnergy(), G4VMscModel::GetRange(), G4LindhardSorensenIonModel::SetupParameters(), and G4BetheBlochModel::SetupParameters().

◆ isLocked

G4bool G4VEmModel::isLocked = false

privateinherited

Definition at line 463 of file G4VEmModel.hh.

Referenced by G4VEmModel::IsLocked(), G4VEmModel::SetLocked(), and G4VEmModel::SetPolarAngleLimit().

◆ isMaster

G4bool G4VEmModel::isMaster = true

privateinherited

Definition at line 457 of file G4VEmModel.hh.

Referenced by G4VEmModel::IsMaster(), G4VEmModel::SetMasterThread(), and G4VEmModel::~G4VEmModel().

◆ localElmSelectors

G4bool G4VEmModel::localElmSelectors = true

privateinherited

Definition at line 460 of file G4VEmModel.hh.

Referenced by G4VEmModel::SetElementSelectors(), and G4VEmModel::~G4VEmModel().

◆ localTable

G4bool G4VEmModel::localTable = true

privateinherited

Definition at line 459 of file G4VEmModel.hh.

Referenced by G4VEmModel::SetCrossSectionTable(), and G4VEmModel::~G4VEmModel().

◆ lossFlucFlag

G4bool G4VEmModel::lossFlucFlag = true

protectedinherited

Definition at line 450 of file G4VEmModel.hh.

Referenced by G4ICRU49NuclearStoppingModel::NuclearStoppingPower(), and G4VEmModel::SetFluctuationFlag().

◆ lowLimit

G4double G4VEmModel::lowLimit

privateinherited

Definition at line 436 of file G4VEmModel.hh.

Referenced by G4VEmModel::InitialiseElementSelectors(), G4VEmModel::LowEnergyLimit(), and G4VEmModel::SetLowEnergyLimit().

◆ name

const G4String G4VEmModel::name

privateinherited

Definition at line 465 of file G4VEmModel.hh.

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

◆ nsec

G4int G4VEmModel::nsec = 5

privateinherited

Definition at line 444 of file G4VEmModel.hh.

Referenced by G4VEmModel::CrossSectionPerVolume(), and G4VEmModel::G4VEmModel().

◆ nSelectors

G4int G4VEmModel::nSelectors = 0

privateinherited

Definition at line 443 of file G4VEmModel.hh.

Referenced by G4VEmModel::InitialiseElementSelectors(), G4VEmModel::SelectRandomAtom(), G4VEmModel::SelectTargetAtom(), G4VEmModel::SetElementSelectors(), and G4VEmModel::~G4VEmModel().

◆ pBaseMaterial

const G4Material* G4VEmModel::pBaseMaterial = nullptr

protectedinherited

Definition at line 427 of file G4VEmModel.hh.

Referenced by G4VEmModel::ComputeDEDX(), G4VEmModel::CrossSection(), G4VMscModel::GetTransportMeanFreePath(), G4VEmModel::SelectRandomAtom(), G4VEmModel::SelectTargetAtom(), G4VEmModel::SetCurrentCouple(), and G4VEmModel::Value().

◆ pFactor

G4double G4VEmModel::pFactor = 1.0

protectedinherited

Definition at line 432 of file G4VEmModel.hh.

Referenced by G4VEmModel::ComputeDEDX(), G4VEmModel::CrossSection(), G4VMscModel::GetTransportMeanFreePath(), G4VEmModel::SetCurrentCouple(), and G4VEmModel::Value().

◆ polarAngleLimit

G4double G4VEmModel::polarAngleLimit

privateinherited

Definition at line 441 of file G4VEmModel.hh.

Referenced by G4VEmModel::PolarAngleLimit(), and G4VEmModel::SetPolarAngleLimit().

◆ pParticleChange

G4VParticleChange* G4VEmModel::pParticleChange = nullptr

protectedinherited

Definition at line 425 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetParticleChangeForGamma(), G4VEmModel::GetParticleChangeForLoss(), G4VMscModel::GetParticleChangeForMSC(), G4EmMultiModel::Initialise(), and G4VEmModel::SetParticleChange().

◆ secondaryThreshold

G4double G4VEmModel::secondaryThreshold = DBL_MAX

privateinherited

Definition at line 440 of file G4VEmModel.hh.

Referenced by G4VEmModel::SecondaryThreshold(), and G4VEmModel::SetSecondaryThreshold().

◆ theDensityFactor

const std::vector<G4double>* G4VEmModel::theDensityFactor = nullptr

protectedinherited

Definition at line 428 of file G4VEmModel.hh.

Referenced by G4VEmModel::G4VEmModel().

◆ theDensityIdx

const std::vector<G4int>* G4VEmModel::theDensityIdx = nullptr

protectedinherited

Definition at line 429 of file G4VEmModel.hh.

Referenced by G4VEmModel::G4VEmModel().

◆ theLPMflag

G4bool G4VEmModel::theLPMflag = false

privateinherited

Definition at line 454 of file G4VEmModel.hh.

Referenced by G4VEmModel::LPMFlag(), and G4VEmModel::SetLPMFlag().

◆ useAngularGenerator

G4bool G4VEmModel::useAngularGenerator = false

privateinherited

Definition at line 461 of file G4VEmModel.hh.

Referenced by G4VEmModel::SetAngularGeneratorFlag(), and G4VEmModel::UseAngularGeneratorFlag().

◆ useBaseMaterials

G4bool G4VEmModel::useBaseMaterials = false

privateinherited

Definition at line 462 of file G4VEmModel.hh.

Referenced by G4VEmModel::SetCurrentCouple(), G4VEmModel::SetUseBaseMaterials(), and G4VEmModel::UseBaseMaterials().

◆ xsec

std::vector<G4double> G4VEmModel::xsec

privateinherited

Definition at line 466 of file G4VEmModel.hh.

Referenced by G4VEmModel::CrossSectionPerVolume(), G4VEmModel::G4VEmModel(), and G4VEmModel::SelectRandomAtom().

◆ xSectionTable

G4PhysicsTable* G4VEmModel::xSectionTable = nullptr

protectedinherited

Definition at line 426 of file G4VEmModel.hh.

Referenced by G4VEmModel::GetCrossSectionTable(), G4VMscModel::GetParticleChangeForMSC(), G4VMscModel::GetTransportMeanFreePath(), G4VEmModel::SetCrossSectionTable(), and G4VEmModel::~G4VEmModel().

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

source/processes/electromagnetic/lowenergy/include/G4PenelopeIonisationModel.hh
source/processes/electromagnetic/lowenergy/src/G4PenelopeIonisationModel.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4PenelopeIonisationModel() [1/2]

◆ ~G4PenelopeIonisationModel()

◆ G4PenelopeIonisationModel() [2/2]

Member Function Documentation

◆ ChargeSquareRatio()

◆ ComputeCrossSectionPerAtom() [1/2]

◆ ComputeCrossSectionPerAtom() [2/2]

◆ ComputeCrossSectionPerShell()

◆ ComputeDEDX()

◆ ComputeDEDXPerVolume()

◆ ComputeMeanFreePath()

◆ CorrectionsAlongStep()

◆ CrossSection()

◆ CrossSectionPerVolume()

◆ CurrentCouple()

◆ DeexcitationFlag()

◆ DefineForRegion()

◆ FillNumberOfSecondaries()

◆ ForceBuildTableFlag()

◆ GetAngularDistribution()

◆ GetChargeSquareRatio()

◆ GetCrossSectionTable()

◆ GetCurrentElement()

◆ GetCurrentIsotope()

◆ GetElementData()

◆ GetElementSelectors()

◆ GetModelOfFluctuations()

◆ GetName()

◆ GetPartialCrossSection()

◆ GetParticleChangeForGamma()

◆ GetParticleChangeForLoss()

◆ GetParticleCharge()

◆ GetTripletModel()

◆ GetVerbosityLevel()

◆ HighEnergyActivationLimit()

◆ HighEnergyLimit()

◆ Initialise()

◆ InitialiseElementSelectors()

◆ InitialiseForElement()

◆ InitialiseForMaterial()

◆ InitialiseLocal()

◆ IsActive()

◆ IsLocked()

◆ IsMaster()

◆ LowEnergyActivationLimit()

◆ LowEnergyLimit()

◆ LPMFlag()

◆ MaxSecondaryEnergy()

◆ MaxSecondaryKinEnergy()

◆ MinEnergyCut()

◆ MinPrimaryEnergy()

◆ ModelDescription()

◆ operator=()

◆ PolarAngleLimit()

◆ SampleFinalStateElectron()

◆ SampleFinalStatePositron()

◆ SampleSecondaries()

◆ SecondaryThreshold()

◆ SelectIsotopeNumber()

◆ SelectRandomAtom() [1/2]

◆ SelectRandomAtom() [2/2]

◆ SelectRandomAtomNumber()

◆ SelectTargetAtom()

◆ SetActivationHighEnergyLimit()

◆ SetActivationLowEnergyLimit()

◆ SetAngularDistribution()

◆ SetAngularGeneratorFlag()

◆ SetCrossSectionTable()

◆ SetCurrentCouple()

◆ SetCurrentElement()

◆ SetDeexcitationFlag()

◆ SetElementSelectors()

◆ SetFluctuationFlag()