#include <G4PenelopeBremsstrahlungModel.hh>

Inheritance diagram for G4PenelopeBremsstrahlungModel:

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 *theParticle, G4double kinEnergy, G4double Z, G4double A=0, G4double cut=0, G4double emax=DBL_MAX) 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

	G4PenelopeBremsstrahlungModel (const G4ParticleDefinition *p=nullptr, const G4String &processName="PenBrem")

	G4PenelopeBremsstrahlungModel (const G4PenelopeBremsstrahlungModel &)=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

G4PenelopeBremsstrahlungModel &	operator= (const G4PenelopeBremsstrahlungModel &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	~G4PenelopeBremsstrahlungModel ()

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	BuildXSTable (const G4Material *material, G4double cut)

void	ClearTables ()

G4PenelopeCrossSection *	GetCrossSectionTableForCouple (const G4ParticleDefinition , const G4Material , G4double cut)

G4double	GetPositronXSCorrection (const G4Material *, G4double energy)

void	SetParticle (const G4ParticleDefinition *)

Private Attributes
G4VEmAngularDistribution *	anglModel = nullptr

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

G4double	eMaxActive = DBL_MAX

G4double	eMinActive = 0.0

const G4MaterialCutsCouple *	fCurrentCouple = nullptr

const G4Element *	fCurrentElement = nullptr

const G4Isotope *	fCurrentIsotope = nullptr

G4LossTableManager *	fEmManager

G4PhysicsLogVector *	fEnergyGrid

G4double	fIntrinsicHighEnergyLimit

G4double	fIntrinsicLowEnergyLimit

G4bool	fIsInitialised

G4bool	flagDeexcitation = false

G4bool	flagForceBuildTable = false

G4bool	fLocalTable

G4VEmFluctuationModel *	flucModel = nullptr

G4PenelopeOscillatorManager *	fOscManager

G4PenelopeBremsstrahlungAngular *	fPenelopeAngular

G4PenelopeBremsstrahlungFS *	fPenelopeFSHelper

G4VEmModel *	fTripletModel = nullptr

G4int	fVerboseLevel

std::map< std::pair< const G4Material , G4double >, G4PenelopeCrossSection > *	fXSTableElectron

std::map< std::pair< const G4Material , G4double >, G4PenelopeCrossSection > *	fXSTablePositron

G4double	highLimit

G4bool	isLocked = false

G4bool	isMaster = true

G4bool	localElmSelectors = true

G4bool	localTable = true

G4double	lowLimit

const G4String	name

size_t	nBins

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 62 of file G4PenelopeBremsstrahlungModel.hh.

Constructor & Destructor Documentation

◆ G4PenelopeBremsstrahlungModel() [1/2]

G4PenelopeBremsstrahlungModel::G4PenelopeBremsstrahlungModel	(	const G4ParticleDefinition *	p = `nullptr`,
		const G4String &	processName = `"PenBrem"`
	)

explicit

Definition at line 63 of file G4PenelopeBremsstrahlungModel.cc.

  :G4VEmModel(nam),fParticleChange(nullptr),fParticle(nullptr),
   fPenelopeFSHelper(nullptr),fPenelopeAngular(nullptr),fEnergyGrid(nullptr),
   fXSTableElectron(nullptr),fXSTablePositron(nullptr),
   fIsInitialised(false),fLocalTable(false)
 
{
  fIntrinsicLowEnergyLimit = 100.0*eV;
  fIntrinsicHighEnergyLimit = 100.0*GeV;
  nBins = 200;
 
  if (part)
    SetParticle(part);
 
  SetHighEnergyLimit(fIntrinsicHighEnergyLimit);
  //
  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, fOscManager, fVerboseLevel, G4PenelopeOscillatorManager::GetOscillatorManager(), GeV, nBins, G4VEmModel::SetDeexcitationFlag(), G4VEmModel::SetHighEnergyLimit(), and SetParticle().

◆ ~G4PenelopeBremsstrahlungModel()

G4PenelopeBremsstrahlungModel::~G4PenelopeBremsstrahlungModel ( )

virtual

Definition at line 97 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (IsMaster() || fLocalTable)
    {
      ClearTables();
      if (fPenelopeFSHelper)
    delete fPenelopeFSHelper;
    }
  // This is thread-local at the moment
  if (fPenelopeAngular)
    delete fPenelopeAngular;
}

References ClearTables(), fLocalTable, fPenelopeAngular, fPenelopeFSHelper, and G4VEmModel::IsMaster().

◆ G4PenelopeBremsstrahlungModel() [2/2]

G4PenelopeBremsstrahlungModel::G4PenelopeBremsstrahlungModel ( const G4PenelopeBremsstrahlungModel & )

delete

Member Function Documentation

◆ BuildXSTable()

void G4PenelopeBremsstrahlungModel::BuildXSTable	(	const G4Material *	material,
		G4double	cut
	)

private

Definition at line 469 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (!IsMaster() && !fLocalTable)
    //Should not be here!
    G4Exception("G4PenelopeBremsstrahlungModel::BuildXSTable()",
        "em0100",FatalException,"Worker thread in this method");
 
  //The key of the map
  std::pair<const G4Material*,G4double> theKey = std::make_pair(mat,cut);
 
  //The table already exists
  if (fXSTableElectron->count(theKey) && fXSTablePositron->count(theKey))
    return;
 
  //
  //This method fills the G4PenelopeCrossSection containers for electrons or positrons
  //and for the given material/cut couple.
  //Equivalent of subroutines EBRaT and PINaT of Penelope
  //
  if (fVerboseLevel > 2)
    {
      G4cout << "G4PenelopeBremsstrahlungModel: going to build cross section table " << G4endl;
      G4cout << "for e+/e- in " << mat->GetName() << " for Ecut(gamma)= " <<
    cut/keV << " keV " << G4endl;
    }
 
  //Tables have been already created (checked by GetCrossSectionTableForCouple)
  if (fEnergyGrid->GetVectorLength() != nBins)
    {
      G4ExceptionDescription ed;
      ed << "Energy Grid looks not initialized" << G4endl;
      ed << nBins << " " << fEnergyGrid->GetVectorLength() << G4endl;
      G4Exception("G4PenelopeBremsstrahlungModel::BuildXSTable()",
          "em2016",FatalException,ed);
    }
 
  G4PenelopeCrossSection* XSEntry = new G4PenelopeCrossSection(nBins);
  G4PenelopeCrossSection* XSEntry2 = new G4PenelopeCrossSection(nBins);
  const G4PhysicsTable* table = fPenelopeFSHelper->GetScaledXSTable(mat,cut);
 
  //loop on the energy grid
  for (size_t bin=0;bin<nBins;bin++)
    {
       G4double energy = fEnergyGrid->GetLowEdgeEnergy(bin);
       G4double XH0=0, XH1=0, XH2=0;
       G4double XS0=0, XS1=0, XS2=0;
 
       //Global xs factor
       G4double fact = fPenelopeFSHelper->GetEffectiveZSquared(mat)*
     ((energy+electron_mass_c2)*(energy+electron_mass_c2)/
      (energy*(energy+2.0*electron_mass_c2)));
 
       G4double restrictedCut = cut/energy;
 
       //Now I need the dSigma/dX profile - interpolated on energy - for
       //the 32-point x grid. Interpolation is log-log
       size_t nBinsX = fPenelopeFSHelper->GetNBinsX();
       G4double* tempData = new G4double[nBinsX];
       G4double logene = G4Log(energy);
       for (size_t ix=0;ix<nBinsX;ix++)
     {
       //find dSigma/dx for the given E. X belongs to the 32-point grid.
       G4double val = (*table)[ix]->Value(logene);
       tempData[ix] = G4Exp(val); //back to the real value!
     }
 
       G4double XH0A = 0.;
       if (restrictedCut <= 1) //calculate only if we are above threshold!
     XH0A = fPenelopeFSHelper->GetMomentumIntegral(tempData,1.0,-1) -
       fPenelopeFSHelper->GetMomentumIntegral(tempData,restrictedCut,-1);
       G4double XS1A = fPenelopeFSHelper->GetMomentumIntegral(tempData,
                                  restrictedCut,0);
       G4double XS2A = fPenelopeFSHelper->GetMomentumIntegral(tempData,
                                  restrictedCut,1);
       G4double XH1A=0, XH2A=0;
       if (restrictedCut <=1)
     {
       XH1A = fPenelopeFSHelper->GetMomentumIntegral(tempData,1.0,0) -
         XS1A;
       XH2A = fPenelopeFSHelper->GetMomentumIntegral(tempData,1.0,1) -
         XS2A;
     }
       delete[] tempData;
 
       XH0 = XH0A*fact;
       XS1 = XS1A*fact*energy;
       XH1 = XH1A*fact*energy;
       XS2 = XS2A*fact*energy*energy;
       XH2 = XH2A*fact*energy*energy;
 
       XSEntry->AddCrossSectionPoint(bin,energy,XH0,XH1,XH2,XS0,XS1,XS2);
 
       //take care of positrons
       G4double posCorrection = GetPositronXSCorrection(mat,energy);
       XSEntry2->AddCrossSectionPoint(bin,energy,XH0*posCorrection,
                      XH1*posCorrection,
                      XH2*posCorrection,
                      XS0,
                      XS1*posCorrection,
                      XS2*posCorrection);
    }
 
  //Insert in the appropriate table
  fXSTableElectron->insert(std::make_pair(theKey,XSEntry));
  fXSTablePositron->insert(std::make_pair(theKey,XSEntry2));
 
  return;
}

References G4PenelopeCrossSection::AddCrossSectionPoint(), source.hepunit::electron_mass_c2, G4INCL::KinematicsUtils::energy(), FatalException, fEnergyGrid, fLocalTable, fPenelopeFSHelper, fVerboseLevel, fXSTableElectron, fXSTablePositron, G4cout, G4endl, G4Exception(), G4Exp(), G4Log(), G4PenelopeBremsstrahlungFS::GetEffectiveZSquared(), G4PhysicsVector::GetLowEdgeEnergy(), G4PenelopeBremsstrahlungFS::GetMomentumIntegral(), G4Material::GetName(), G4PenelopeBremsstrahlungFS::GetNBinsX(), GetPositronXSCorrection(), G4PenelopeBremsstrahlungFS::GetScaledXSTable(), G4PhysicsVector::GetVectorLength(), G4VEmModel::IsMaster(), keV, and nBins.

Referenced by GetCrossSectionTableForCouple(), and Initialise().

◆ 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().

◆ ClearTables()

void G4PenelopeBremsstrahlungModel::ClearTables ( )

private

Definition at line 425 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (!IsMaster() && !fLocalTable)
    //Should not be here!
    G4Exception("G4PenelopeBremsstrahlungModel::ClearTables()",
        "em0100",FatalException,"Worker thread in this method");
 
  if (fXSTableElectron)
    {
      for (auto& item : (*fXSTableElectron))        
    delete item.second;        
      delete fXSTableElectron;
      fXSTableElectron = nullptr;
    }
  if (fXSTablePositron)
    {
      for (auto& item : (*fXSTablePositron))                
    delete item.second;    
      delete fXSTablePositron;
      fXSTablePositron = nullptr;
    }
  /*
  if (fEnergyGrid)
    delete fEnergyGrid;
  */
  if (fPenelopeFSHelper)
    fPenelopeFSHelper->ClearTables(IsMaster());
 
  if (fVerboseLevel > 2)
    G4cout << "G4PenelopeBremsstrahlungModel: cleared tables" << G4endl;
 
 return;
}

References G4PenelopeBremsstrahlungFS::ClearTables(), FatalException, fLocalTable, fPenelopeFSHelper, fVerboseLevel, fXSTableElectron, fXSTablePositron, G4cout, G4endl, G4Exception(), and G4VEmModel::IsMaster().

Referenced by Initialise(), and ~G4PenelopeBremsstrahlungModel().

◆ 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 G4PenelopeBremsstrahlungModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	theParticle,
		G4double	kinEnergy,
		G4double	Z,
		G4double	A = `0`,
		G4double	cut = `0`,
		G4double	emax = `DBL_MAX`
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 271 of file G4PenelopeBremsstrahlungModel.cc.

{
  G4cout << "*** G4PenelopeBremsstrahlungModel -- WARNING ***" << G4endl;
  G4cout << "Penelope Bremsstrahlung 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 G4PenelopeBremsstrahlungModel::ComputeDEDXPerVolume	(	const G4Material *	material,
		const G4ParticleDefinition *	theParticle,
		G4double	kineticEnergy,
		G4double	cutEnergy
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 287 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (fVerboseLevel > 3)
    G4cout << "Calling ComputeDEDX() of G4PenelopeBremsstrahlungModel" << G4endl;
 
  G4PenelopeCrossSection* theXS = 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 << "G4PenelopeBremsstrahlungModel " << G4endl;
      G4cout << "Stopping power < " << cutEnergy/keV << " keV at " <<
        kineticEnergy/keV << " keV = " <<
        sPowerPerVolume/(keV/mm) << " keV/mm" << G4endl;
    }
  return sPowerPerVolume;
}

References fOscManager, fVerboseLevel, G4cout, G4endl, G4PenelopeOscillatorManager::GetAtomsPerMolecule(), GetCrossSectionTableForCouple(), G4PenelopeCrossSection::GetSoftStoppingPower(), keV, eplot::material, and mm.

◆ 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 G4PenelopeBremsstrahlungModel::CrossSectionPerVolume	(	const G4Material *	material,
		const G4ParticleDefinition *	theParticle,
		G4double	kineticEnergy,
		G4double	cutEnergy,
		G4double	maxEnergy = `DBL_MAX`
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 225 of file G4PenelopeBremsstrahlungModel.cc.

{
  //
  if (fVerboseLevel > 3)
    G4cout << "Calling CrossSectionPerVolume() of G4PenelopeBremsstrahlungModel" << G4endl;
 
  SetupForMaterial(theParticle, material, energy);
  G4double crossPerMolecule = 0.;
 
  G4PenelopeCrossSection* theXS = 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 << "G4PenelopeBremsstrahlungModel " << G4endl;
    G4cout << "Mean free path for gamma emission > " << cutEnergy/keV << " keV at " <<
      energy/keV << " keV = " << (1./crossPerVolume)/mm << " mm" << G4endl;
  }
 
  return crossPerVolume;
}

References G4INCL::KinematicsUtils::energy(), fOscManager, fVerboseLevel, G4cout, G4endl, G4PenelopeOscillatorManager::GetAtomsPerMolecule(), GetCrossSectionTableForCouple(), G4PenelopeCrossSection::GetHardCrossSection(), keV, eplot::material, mm, and G4VEmModel::SetupForMaterial().

◆ 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().

◆ GetCrossSectionTableForCouple()

G4PenelopeCrossSection * G4PenelopeBremsstrahlungModel::GetCrossSectionTableForCouple	(	const G4ParticleDefinition *	part,
		const G4Material *	mat,
		G4double	cut
	)

private

Definition at line 582 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (part != G4Electron::Electron() && part != G4Positron::Positron())
    {
      G4ExceptionDescription ed;
      ed << "Invalid particle: " << part->GetParticleName() << G4endl;
      G4Exception("G4PenelopeBremsstrahlungModel::GetCrossSectionTableForCouple()",
          "em0001",FatalException,ed);
      return nullptr;
    }
 
  if (part == G4Electron::Electron())
    {
      //Either Initialize() was not called, or we are in a slave and InitializeLocal() was
      //not invoked
      if (!fXSTableElectron)
        {
      //create a **thread-local** version of the table. Used only for G4EmCalculator and
      //Unit Tests
          G4String excep = "The Cross Section Table for e- was not initialized correctly!";
          G4Exception("G4PenelopeBremsstrahlungModel::GetCrossSectionTableForCouple()",
              "em2013",JustWarning,excep);
      fLocalTable = true;
      fXSTableElectron = new
        std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*>;
      if (!fEnergyGrid)
        fEnergyGrid = new G4PhysicsLogVector(LowEnergyLimit(),
                        HighEnergyLimit(),
                        nBins-1); //one hidden bin is added
      if (!fPenelopeFSHelper)
        fPenelopeFSHelper = new G4PenelopeBremsstrahlungFS(fVerboseLevel);
        }
      std::pair<const G4Material*,G4double> theKey = std::make_pair(mat,cut);
      if (fXSTableElectron->count(theKey)) //table already built
        return fXSTableElectron->find(theKey)->second;
      else
    {
      //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)
        {
          //G4Exception (warning) is issued only in verbose mode
          G4ExceptionDescription ed;
          ed << "Unable to find e- table for " << mat->GetName() << " at Ecut(gamma)= "
         << cut/keV << " keV " << G4endl;
          ed << "This can happen only in Unit Tests or via G4EmCalculator" << G4endl;
          G4Exception("G4PenelopeBremsstrahlungModel::GetCrossSectionTableForCouple()",
              "em2009",JustWarning,ed);
        }
      //protect file reading via autolock
      G4AutoLock lock(&PenelopeBremsstrahlungModelMutex);
          fPenelopeFSHelper->BuildScaledXSTable(mat,cut,true); //pretend to be a master
      BuildXSTable(mat,cut);
      lock.unlock();
      //now it should be ok
      return fXSTableElectron->find(theKey)->second;
    }
    }
  if (part == G4Positron::Positron())
    {
      //Either Initialize() was not called, or we are in a slave and InitializeLocal() was
      //not invoked
      if (!fXSTablePositron)
        {
      G4String excep = "The Cross Section Table for e+ was not initialized correctly!";
          G4Exception("G4PenelopeBremsstrahlungModel::GetCrossSectionTableForCouple()",
              "em2013",JustWarning,excep);
      fLocalTable = true;
      fXSTablePositron = new
        std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*>;
      if (!fEnergyGrid)
        fEnergyGrid = new G4PhysicsLogVector(LowEnergyLimit(),
                        HighEnergyLimit(),
                        nBins-1); //one hidden bin is added
      if (!fPenelopeFSHelper)
            fPenelopeFSHelper = new G4PenelopeBremsstrahlungFS(fVerboseLevel);
        }
      std::pair<const G4Material*,G4double> theKey = std::make_pair(mat,cut);
      if (fXSTablePositron->count(theKey)) //table already built
        return fXSTablePositron->find(theKey)->second;
      else
        {
      //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 find e+ table for " << mat->GetName() << " at Ecut(gamma)= "
         << cut/keV << " keV " << G4endl;
          ed << "This can happen only in Unit Tests or via G4EmCalculator" << G4endl;
          G4Exception("G4PenelopeBremsstrahlungModel::GetCrossSectionTableForCouple()",
              "em2009",JustWarning,ed);
        }
      //protect file reading via autolock
      G4AutoLock lock(&PenelopeBremsstrahlungModelMutex);
          fPenelopeFSHelper->BuildScaledXSTable(mat,cut,true); //pretend to be a master
      BuildXSTable(mat,cut);
      lock.unlock();
      //now it should be ok
      return fXSTablePositron->find(theKey)->second;
        }
    }
  return nullptr;
}

References G4PenelopeBremsstrahlungFS::BuildScaledXSTable(), BuildXSTable(), G4Electron::Electron(), FatalException, fEnergyGrid, fLocalTable, fPenelopeFSHelper, fVerboseLevel, fXSTableElectron, fXSTablePositron, G4endl, G4Exception(), G4Material::GetName(), G4ParticleDefinition::GetParticleName(), G4VEmModel::HighEnergyLimit(), JustWarning, keV, G4VEmModel::LowEnergyLimit(), nBins, anonymous_namespace{G4PenelopeBremsstrahlungModel.cc}::PenelopeBremsstrahlungModelMutex, G4Positron::Positron(), and G4TemplateAutoLock< _Mutex_t >::unlock().

Referenced by ComputeDEDXPerVolume(), and CrossSectionPerVolume().

◆ 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().

◆ GetPositronXSCorrection()

G4double G4PenelopeBremsstrahlungModel::GetPositronXSCorrection	(	const G4Material *	mat,
		G4double	energy
	)

private

Definition at line 692 of file G4PenelopeBremsstrahlungModel.cc.

{
  //The electron-to-positron correction factor is set equal to the ratio of the
  //radiative stopping powers for positrons and electrons, which has been calculated
  //by Kim et al. (1986) (cf. Berger and Seltzer, 1982). Here, it is used an
  //analytical approximation which reproduces the tabulated values with 0.5%
  //accuracy
  G4double t=G4Log(1.0+1e6*energy/
              (electron_mass_c2*fPenelopeFSHelper->GetEffectiveZSquared(mat)));
  G4double corr = 1.0-G4Exp(-t*(1.2359e-1-t*(6.1274e-2-t*
                       (3.1516e-2-t*(7.7446e-3-t*(1.0595e-3-t*
                                      (7.0568e-5-t*
                                       1.8080e-6)))))));
  return corr;
}

References e6, source.hepunit::electron_mass_c2, G4INCL::KinematicsUtils::energy(), fPenelopeFSHelper, G4Exp(), G4Log(), and G4PenelopeBremsstrahlungFS::GetEffectiveZSquared().

Referenced by BuildXSTable().

◆ 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 G4PenelopeBremsstrahlungModel::GetVerbosityLevel ( )

inline

Definition at line 104 of file G4PenelopeBremsstrahlungModel.hh.

104{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 G4PenelopeBremsstrahlungModel::Initialise	(	const G4ParticleDefinition *	part,
		const G4DataVector &	theCuts
	)

overridevirtual

Implements G4VEmModel.

Definition at line 112 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (fVerboseLevel > 3)
    G4cout << "Calling G4PenelopeBremsstrahlungModel::Initialise()" << G4endl;
 
  SetParticle(part);
 
  if (IsMaster() && part == fParticle)
    {
      if (!fPenelopeFSHelper)
    fPenelopeFSHelper = new G4PenelopeBremsstrahlungFS(fVerboseLevel);
      if (!fPenelopeAngular)
    fPenelopeAngular = new G4PenelopeBremsstrahlungAngular();
      //Clear and re-build the tables
      ClearTables();
 
      //forces the cleaning of tables, in this specific case
      if (fPenelopeAngular)
    fPenelopeAngular->Initialize();
 
      //Set the number of bins for the tables. 20 points per decade
      nBins = (size_t) (20*std::log10(HighEnergyLimit()/LowEnergyLimit()));
      nBins = std::max(nBins,(size_t)100);
      fEnergyGrid = new G4PhysicsLogVector(LowEnergyLimit(),
                                      HighEnergyLimit(),
                                      nBins-1); //one hidden bin is added
 
      fXSTableElectron = new
    std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*>;
      fXSTablePositron = new
    std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*>;
 
      G4ProductionCutsTable* theCoupleTable =
    G4ProductionCutsTable::GetProductionCutsTable();
 
      //Build tables for all materials
      for (size_t i=0;i<theCoupleTable->GetTableSize();i++)
    {
      const G4Material* theMat =
        theCoupleTable->GetMaterialCutsCouple(i)->GetMaterial();
      //Forces the building of the helper tables
      fPenelopeFSHelper->BuildScaledXSTable(theMat,theCuts.at(i),IsMaster());
      fPenelopeAngular->PrepareTables(theMat,IsMaster());
      BuildXSTable(theMat,theCuts.at(i));
 
    }
 
      if (fVerboseLevel > 2) {
    G4cout << "Penelope Bremsstrahlung model v2008 is initialized " << G4endl
           << "Energy range: "
           << LowEnergyLimit() / keV << " keV - "
           << HighEnergyLimit() / GeV << " GeV."
           << G4endl;
      }
    }
 
  if(fIsInitialised) return;
  fParticleChange = GetParticleChangeForLoss();
  fIsInitialised = true;
}

References G4PenelopeBremsstrahlungFS::BuildScaledXSTable(), BuildXSTable(), ClearTables(), fEnergyGrid, fIsInitialised, fParticle, fParticleChange, fPenelopeAngular, fPenelopeFSHelper, fVerboseLevel, fXSTableElectron, fXSTablePositron, G4cout, G4endl, G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4VEmModel::GetParticleChangeForLoss(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), GeV, G4VEmModel::HighEnergyLimit(), G4PenelopeBremsstrahlungAngular::Initialize(), G4VEmModel::IsMaster(), keV, G4VEmModel::LowEnergyLimit(), G4INCL::Math::max(), nBins, G4PenelopeBremsstrahlungAngular::PrepareTables(), and SetParticle().

◆ 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 G4PenelopeBremsstrahlungModel::InitialiseLocal	(	const G4ParticleDefinition *	part,
		G4VEmModel *	masterModel
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 176 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (fVerboseLevel > 3)
    G4cout << "Calling  G4PenelopeBremsstrahlungModel::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 G4PenelopeBremsstrahlungModel* theModel =
    static_cast<G4PenelopeBremsstrahlungModel*> (masterModel);
 
      //Copy pointers to the data tables
      fEnergyGrid = theModel->fEnergyGrid;
      fXSTableElectron = theModel->fXSTableElectron;
      fXSTablePositron = theModel->fXSTablePositron;
      fPenelopeFSHelper = theModel->fPenelopeFSHelper;
 
      //created in each thread and initialized.
      if (!fPenelopeAngular)
    fPenelopeAngular = new G4PenelopeBremsstrahlungAngular();
      //forces the cleaning of tables, in this specific case
      if (fPenelopeAngular)
    fPenelopeAngular->Initialize();
 
      G4ProductionCutsTable* theCoupleTable =
    G4ProductionCutsTable::GetProductionCutsTable();
      //Build tables for all materials
      for (size_t i=0;i<theCoupleTable->GetTableSize();i++)
    {
      const G4Material* theMat =
        theCoupleTable->GetMaterialCutsCouple(i)->GetMaterial();
      fPenelopeAngular->PrepareTables(theMat,IsMaster());
    }
 
      //copy the data
      nBins = theModel->nBins;
 
      //Same verbosity for all workers, as the master
      fVerboseLevel = theModel->fVerboseLevel;
    }
  return;
}

References fEnergyGrid, fParticle, fPenelopeAngular, fPenelopeFSHelper, fVerboseLevel, fXSTableElectron, fXSTablePositron, G4cout, G4endl, G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), G4PenelopeBremsstrahlungAngular::Initialize(), G4VEmModel::IsMaster(), nBins, and G4PenelopeBremsstrahlungAngular::PrepareTables().

◆ 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 G4PenelopeBremsstrahlungModel::MinEnergyCut	(	const G4ParticleDefinition *	,
		const G4MaterialCutsCouple *
	)

overridevirtual

Reimplemented from G4VEmModel.

Definition at line 461 of file G4PenelopeBremsstrahlungModel.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=()

G4PenelopeBremsstrahlungModel & G4PenelopeBremsstrahlungModel::operator= ( const G4PenelopeBremsstrahlungModel & 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().

◆ SampleSecondaries()

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

overridevirtual

Implements G4VEmModel.

Definition at line 322 of file G4PenelopeBremsstrahlungModel.cc.

{
  if (fVerboseLevel > 3)
    G4cout << "Calling SampleSecondaries() of G4PenelopeBremsstrahlungModel" << G4endl;
 
  G4double kineticEnergy = aDynamicParticle->GetKineticEnergy();
  const G4Material* material = couple->GetMaterial();
 
  if (kineticEnergy <= fIntrinsicLowEnergyLimit)
   {
     fParticleChange->SetProposedKineticEnergy(0.);
     fParticleChange->ProposeLocalEnergyDeposit(kineticEnergy);
     return ;
   }
 
  G4ParticleMomentum particleDirection0 = aDynamicParticle->GetMomentumDirection();
  //This is the momentum
  G4ThreeVector initialMomentum =  aDynamicParticle->GetMomentum();
 
  //Not enough energy to produce a secondary! Return with nothing happened
  if (kineticEnergy < cutG)
    return;
 
  if (fVerboseLevel > 3)
    G4cout << "Going to sample gamma energy for: " <<material->GetName() << " " <<
      "energy = " << kineticEnergy/keV << ", cut = " << cutG/keV << G4endl;
 
   //Sample gamma's energy according to the spectrum
  G4double gammaEnergy =
    fPenelopeFSHelper->SampleGammaEnergy(kineticEnergy,material,cutG);
 
  if (fVerboseLevel > 3)
    G4cout << "Sampled gamma energy: " << gammaEnergy/keV << " keV" << G4endl;
 
  //Now sample the direction for the Gamma. Notice that the rotation is already done
  //Z is unused here, I plug 0. The information is in the material pointer
   G4ThreeVector gammaDirection1 =
     fPenelopeAngular->SampleDirection(aDynamicParticle,gammaEnergy,0,material);
 
  if (fVerboseLevel > 3)
    G4cout << "Sampled cosTheta for e-: " << gammaDirection1.cosTheta() << G4endl;
 
  G4double residualPrimaryEnergy = kineticEnergy-gammaEnergy;
  if (residualPrimaryEnergy < 0)
    {
      //Ok we have a problem, all energy goes with the gamma
      gammaEnergy += residualPrimaryEnergy;
      residualPrimaryEnergy = 0.0;
    }
 
  //Produce final state according to momentum conservation
  G4ThreeVector particleDirection1 = initialMomentum - gammaEnergy*gammaDirection1;
  particleDirection1 = particleDirection1.unit(); //normalize
 
  //Update the primary particle
  if (residualPrimaryEnergy > 0.)
    {
      fParticleChange->ProposeMomentumDirection(particleDirection1);
      fParticleChange->SetProposedKineticEnergy(residualPrimaryEnergy);
    }
  else
    {
      fParticleChange->SetProposedKineticEnergy(0.);
    }
 
  //Now produce the photon
  G4DynamicParticle* theGamma = new G4DynamicParticle(G4Gamma::Gamma(),
                                                      gammaDirection1,
                                                      gammaEnergy);
  fvect->push_back(theGamma);
 
  if (fVerboseLevel > 1)
    {
      G4cout << "-----------------------------------------------------------" << G4endl;
      G4cout << "Energy balance from G4PenelopeBremsstrahlung" << G4endl;
      G4cout << "Incoming primary energy: " << kineticEnergy/keV << " keV" << G4endl;
      G4cout << "-----------------------------------------------------------" << G4endl;
      G4cout << "Outgoing primary energy: " << residualPrimaryEnergy/keV << " keV" << G4endl;
      G4cout << "Bremsstrahlung photon " << gammaEnergy/keV << " keV" << G4endl;
      G4cout << "Total final state: " << (residualPrimaryEnergy+gammaEnergy)/keV
             << " keV" << G4endl;
      G4cout << "-----------------------------------------------------------" << G4endl;
    }
 
  if (fVerboseLevel > 0)
    {
      G4double energyDiff = std::fabs(residualPrimaryEnergy+gammaEnergy-kineticEnergy);
      if (energyDiff > 0.05*keV)
        G4cout << "Warning from G4PenelopeBremsstrahlung: problem with energy conservation: "
           <<
          (residualPrimaryEnergy+gammaEnergy)/keV <<
          " keV (final) vs. " <<
          kineticEnergy/keV << " keV (initial)" << G4endl;
    }
  return;
}

References CLHEP::Hep3Vector::cosTheta(), fIntrinsicLowEnergyLimit, fParticleChange, fPenelopeAngular, fPenelopeFSHelper, fVerboseLevel, G4cout, G4endl, G4Gamma::Gamma(), G4DynamicParticle::GetKineticEnergy(), G4MaterialCutsCouple::GetMaterial(), G4DynamicParticle::GetMomentum(), G4DynamicParticle::GetMomentumDirection(), keV, eplot::material, G4VParticleChange::ProposeLocalEnergyDeposit(), G4ParticleChangeForLoss::ProposeMomentumDirection(), G4PenelopeBremsstrahlungAngular::SampleDirection(), G4PenelopeBremsstrahlungFS::SampleGammaEnergy(), G4ParticleChangeForLoss::SetProposedKineticEnergy(), and CLHEP::Hep3Vector::unit().

◆ 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(), G4PenelopeComptonModel::G4PenelopeComptonModel(), G4PenelopeIonisationModel::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 G4PenelopeBremsstrahlungModel::SetParticle ( const G4ParticleDefinition * p )

private

Definition at line 711 of file G4PenelopeBremsstrahlungModel.cc.

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

References fParticle.

Referenced by G4PenelopeBremsstrahlungModel(), 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(), CrossSectionPerVolume(), G4PenelopeIonisationModel::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 G4PenelopeBremsstrahlungModel::SetVerbosityLevel ( G4int lev )

inline

Definition at line 103 of file G4PenelopeBremsstrahlungModel.hh.

103{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().

◆ 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().

◆ fEnergyGrid

G4PhysicsLogVector* G4PenelopeBremsstrahlungModel::fEnergyGrid

private

Definition at line 130 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), GetCrossSectionTableForCouple(), Initialise(), and InitialiseLocal().

◆ fIntrinsicHighEnergyLimit

G4double G4PenelopeBremsstrahlungModel::fIntrinsicHighEnergyLimit

private

Definition at line 138 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by G4PenelopeBremsstrahlungModel().

◆ fIntrinsicLowEnergyLimit

G4double G4PenelopeBremsstrahlungModel::fIntrinsicLowEnergyLimit

private

Definition at line 137 of file G4PenelopeBremsstrahlungModel.hh.

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

◆ fIsInitialised

G4bool G4PenelopeBremsstrahlungModel::fIsInitialised

private

Definition at line 141 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by Initialise().

◆ 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 G4PenelopeBremsstrahlungModel::fLocalTable

private

Definition at line 143 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), ClearTables(), GetCrossSectionTableForCouple(), and ~G4PenelopeBremsstrahlungModel().

◆ flucModel

G4VEmFluctuationModel* G4VEmModel::flucModel = nullptr

privateinherited

Definition at line 413 of file G4VEmModel.hh.

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

◆ fOscManager

G4PenelopeOscillatorManager* G4PenelopeBremsstrahlungModel::fOscManager

private

Definition at line 125 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by ComputeDEDXPerVolume(), CrossSectionPerVolume(), and G4PenelopeBremsstrahlungModel().

◆ fParticle

const G4ParticleDefinition* G4PenelopeBremsstrahlungModel::fParticle

protected

Definition at line 112 of file G4PenelopeBremsstrahlungModel.hh.

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

◆ fParticleChange

G4ParticleChangeForLoss* G4PenelopeBremsstrahlungModel::fParticleChange

protected

Definition at line 111 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by Initialise(), and SampleSecondaries().

◆ fPenelopeAngular

G4PenelopeBremsstrahlungAngular* G4PenelopeBremsstrahlungModel::fPenelopeAngular

private

Definition at line 127 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by Initialise(), InitialiseLocal(), SampleSecondaries(), and ~G4PenelopeBremsstrahlungModel().

◆ fPenelopeFSHelper

G4PenelopeBremsstrahlungFS* G4PenelopeBremsstrahlungModel::fPenelopeFSHelper

private

Definition at line 126 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), ClearTables(), GetCrossSectionTableForCouple(), GetPositronXSCorrection(), Initialise(), InitialiseLocal(), SampleSecondaries(), and ~G4PenelopeBremsstrahlungModel().

◆ 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 G4PenelopeBremsstrahlungModel::fVerboseLevel

private

Definition at line 140 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), ClearTables(), ComputeDEDXPerVolume(), CrossSectionPerVolume(), G4PenelopeBremsstrahlungModel(), GetCrossSectionTableForCouple(), GetVerbosityLevel(), Initialise(), InitialiseLocal(), SampleSecondaries(), and SetVerbosityLevel().

◆ fXSTableElectron

std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*>* G4PenelopeBremsstrahlungModel::fXSTableElectron

private

Definition at line 133 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), ClearTables(), GetCrossSectionTableForCouple(), Initialise(), and InitialiseLocal().

◆ fXSTablePositron

std::map< std::pair<const G4Material*,G4double>, G4PenelopeCrossSection*>* G4PenelopeBremsstrahlungModel::fXSTablePositron

private

Definition at line 134 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), ClearTables(), GetCrossSectionTableForCouple(), Initialise(), and InitialiseLocal().

◆ 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().

◆ nBins

size_t G4PenelopeBremsstrahlungModel::nBins

private

Definition at line 131 of file G4PenelopeBremsstrahlungModel.hh.

Referenced by BuildXSTable(), G4PenelopeBremsstrahlungModel(), GetCrossSectionTableForCouple(), Initialise(), and InitialiseLocal().

◆ 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/G4PenelopeBremsstrahlungModel.hh
source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungModel.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4PenelopeBremsstrahlungModel() [1/2]

◆ ~G4PenelopeBremsstrahlungModel()

◆ G4PenelopeBremsstrahlungModel() [2/2]

Member Function Documentation

◆ BuildXSTable()

◆ ChargeSquareRatio()

◆ ClearTables()

◆ ComputeCrossSectionPerAtom() [1/2]

◆ ComputeCrossSectionPerAtom() [2/2]

◆ ComputeCrossSectionPerShell()

◆ ComputeDEDX()

◆ ComputeDEDXPerVolume()

◆ ComputeMeanFreePath()

◆ CorrectionsAlongStep()

◆ CrossSection()

◆ CrossSectionPerVolume()

◆ CurrentCouple()

◆ DeexcitationFlag()

◆ DefineForRegion()

◆ FillNumberOfSecondaries()

◆ ForceBuildTableFlag()

◆ GetAngularDistribution()

◆ GetChargeSquareRatio()

◆ GetCrossSectionTable()

◆ GetCrossSectionTableForCouple()

◆ GetCurrentElement()

◆ GetCurrentIsotope()

◆ GetElementData()

◆ GetElementSelectors()

◆ GetModelOfFluctuations()

◆ GetName()

◆ GetPartialCrossSection()

◆ GetParticleChangeForGamma()

◆ GetParticleChangeForLoss()

◆ GetParticleCharge()

◆ GetPositronXSCorrection()

◆ GetTripletModel()

◆ GetVerbosityLevel()

◆ HighEnergyActivationLimit()

◆ HighEnergyLimit()

◆ Initialise()

◆ InitialiseElementSelectors()

◆ InitialiseForElement()

◆ InitialiseForMaterial()

◆ InitialiseLocal()

◆ IsActive()

◆ IsLocked()

◆ IsMaster()

◆ LowEnergyActivationLimit()

◆ LowEnergyLimit()

◆ LPMFlag()

◆ MaxSecondaryEnergy()

◆ MaxSecondaryKinEnergy()

◆ MinEnergyCut()

◆ MinPrimaryEnergy()

◆ ModelDescription()

◆ operator=()

◆ PolarAngleLimit()

◆ SampleSecondaries()

◆ SecondaryThreshold()

◆ SelectIsotopeNumber()

◆ SelectRandomAtom() [1/2]

◆ SelectRandomAtom() [2/2]

◆ SelectRandomAtomNumber()

◆ SelectTargetAtom()

◆ SetActivationHighEnergyLimit()

◆ SetActivationLowEnergyLimit()

◆ SetAngularDistribution()

◆ SetAngularGeneratorFlag()

◆ SetCrossSectionTable()

◆ SetCurrentCouple()

◆ SetCurrentElement()

◆ SetDeexcitationFlag()