G4VEmProcess Class Reference

#include <G4VEmProcess.hh>

Inheritance diagram for G4VEmProcess:

Public Member Functions
	G4VEmProcess (const G4String &name, G4ProcessType type=fElectromagnetic)
	~G4VEmProcess () override
virtual G4bool	IsApplicable (const G4ParticleDefinition &p) override=0
void	ProcessDescription (std::ostream &outFile) const override
void	PreparePhysicsTable (const G4ParticleDefinition &) override
void	BuildPhysicsTable (const G4ParticleDefinition &) override
void	StartTracking (G4Track *) override
G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &) override
G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii=false) override
G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii) override
virtual G4VEmProcess *	GetEmProcess (const G4String &name)
G4double	GetLambda (G4double kinEnergy, const G4MaterialCutsCouple *couple, G4double logKinEnergy)
G4double	GetCrossSection (const G4double kinEnergy, const G4MaterialCutsCouple *couple) override
G4double	ComputeCrossSectionPerAtom (G4double kineticEnergy, G4double Z, G4double A=0., G4double cut=0.0)
G4double	MeanFreePath (const G4Track &track)
void	SetLambdaBinning (G4int nbins)
void	SetMinKinEnergy (G4double e)
void	SetMinKinEnergyPrim (G4double e)
void	SetMaxKinEnergy (G4double e)
G4PhysicsTable *	LambdaTable () const
G4PhysicsTable *	LambdaTablePrim () const
void	SetLambdaTable (G4PhysicsTable *)
void	SetLambdaTablePrim (G4PhysicsTable *)
std::vector< G4double > *	EnergyOfCrossSectionMax () const
void	SetEnergyOfCrossSectionMax (std::vector< G4double > *)
G4CrossSectionType	CrossSectionType () const
void	SetCrossSectionType (G4CrossSectionType val)
const G4ParticleDefinition *	Particle () const
const G4ParticleDefinition *	SecondaryParticle () const
G4VEmModel *	SelectModelForMaterial (G4double kinEnergy, size_t idxCouple) const
void	AddEmModel (G4int, G4VEmModel *, const G4Region *region=nullptr)
void	SetEmModel (G4VEmModel *, G4int index=0)
G4int	NumberOfModels () const
G4VEmModel *	EmModel (size_t index=0) const
const G4VEmModel *	GetCurrentModel () const
G4VEmModel *	GetModelByIndex (G4int idx=0, G4bool ver=false) const
const G4Element *	GetCurrentElement () const
void	SetCrossSectionBiasingFactor (G4double f, G4bool flag=true)
G4double	CrossSectionBiasingFactor () const
void	ActivateForcedInteraction (G4double length=0.0, const G4String &r="", G4bool flag=true)
void	ActivateSecondaryBiasing (const G4String &region, G4double factor, G4double energyLimit)
void	SetEmMasterProcess (const G4VEmProcess *)
void	SetBuildTableFlag (G4bool val)
void	CurrentSetup (const G4MaterialCutsCouple *, G4double energy)
G4bool	UseBaseMaterial () const
void	BuildLambdaTable ()
void	StreamInfo (std::ostream &outFile, const G4ParticleDefinition &, G4bool rst=false) const
	G4VEmProcess (G4VEmProcess &)=delete
G4VEmProcess &	operator= (const G4VEmProcess &right)=delete
Public Member Functions inherited from G4VDiscreteProcess
	G4VDiscreteProcess (const G4String &aName, G4ProcessType aType=fNotDefined)
	G4VDiscreteProcess (G4VDiscreteProcess &)
virtual	~G4VDiscreteProcess ()
G4VDiscreteProcess &	operator= (const G4VDiscreteProcess &)=delete
virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)
virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)
virtual G4double	MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *)
Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
	G4VProcess (const G4VProcess &right)
virtual	~G4VProcess ()
G4VProcess &	operator= (const G4VProcess &)=delete
G4bool	operator== (const G4VProcess &right) const
G4bool	operator!= (const G4VProcess &right) const
G4double	GetCurrentInteractionLength () const
void	SetPILfactor (G4double value)
G4double	GetPILfactor () const
G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
const G4String &	GetProcessName () const
G4ProcessType	GetProcessType () const
void	SetProcessType (G4ProcessType)
G4int	GetProcessSubType () const
void	SetProcessSubType (G4int)
virtual const G4VProcess *	GetCreatorProcess () const
virtual void	EndTracking ()
virtual void	SetProcessManager (const G4ProcessManager *)
virtual const G4ProcessManager *	GetProcessManager ()
virtual void	ResetNumberOfInteractionLengthLeft ()
G4double	GetNumberOfInteractionLengthLeft () const
G4double	GetTotalNumberOfInteractionLengthTraversed () const
G4bool	isAtRestDoItIsEnabled () const
G4bool	isAlongStepDoItIsEnabled () const
G4bool	isPostStepDoItIsEnabled () const
virtual void	DumpInfo () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
virtual void	SetMasterProcess (G4VProcess *masterP)
const G4VProcess *	GetMasterProcess () const
virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Protected Member Functions
virtual void	StreamProcessInfo (std::ostream &) const
virtual void	InitialiseProcess (const G4ParticleDefinition *)=0
G4VEmModel *	SelectModel (G4double kinEnergy, size_t)
G4double	GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
G4PhysicsVector *	LambdaPhysicsVector (const G4MaterialCutsCouple *)
void	DefineMaterial (const G4MaterialCutsCouple *couple)
G4int	LambdaBinning () const
G4double	MinKinEnergy () const
G4double	MaxKinEnergy () const
G4double	PolarAngleLimit () const
G4ParticleChangeForGamma *	GetParticleChange ()
void	SetParticle (const G4ParticleDefinition *p)
void	SetSecondaryParticle (const G4ParticleDefinition *p)
size_t	CurrentMaterialCutsCoupleIndex () const
const G4MaterialCutsCouple *	MaterialCutsCouple () const
G4bool	ApplyCuts () const
G4double	GetGammaEnergyCut ()
G4double	GetElectronEnergyCut ()
void	SetStartFromNullFlag (G4bool val)
void	SetSplineFlag (G4bool val)
const G4Element *	GetTargetElement () const
const G4Isotope *	GetTargetIsotope () const
G4int	DensityIndex (G4int idx) const
G4double	DensityFactor (G4int idx) const
Protected Member Functions inherited from G4VDiscreteProcess
Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)
void	ClearNumberOfInteractionLengthLeft ()

Protected Attributes
const G4MaterialCutsCouple *	currentCouple = nullptr
const G4Material *	currentMaterial = nullptr
G4EmBiasingManager *	biasManager = nullptr
std::vector< G4double > *	theEnergyOfCrossSectionMax = nullptr
G4double	mfpKinEnergy = DBL_MAX
G4double	preStepKinEnergy = 0.0
G4double	preStepLambda = 0.0
G4int	mainSecondaries = 1
G4int	secID = _EM
G4int	fluoID = _Fluorescence
G4int	augerID = _AugerElectron
G4int	biasID = _EM
G4int	tripletID = _TripletElectron
size_t	currentCoupleIndex = 0
size_t	basedCoupleIndex = 0
size_t	coupleIdxLambda = 0
size_t	idxLambda = 0
G4bool	isTheMaster = true
G4bool	baseMat = false
std::vector< G4DynamicParticle * >	secParticles
G4ParticleChangeForGamma	fParticleChange
Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager = nullptr
G4VParticleChange *	pParticleChange = nullptr
G4ParticleChange	aParticleChange
G4double	theNumberOfInteractionLengthLeft = -1.0
G4double	currentInteractionLength = -1.0
G4double	theInitialNumberOfInteractionLength = -1.0
G4String	theProcessName
G4String	thePhysicsTableFileName
G4ProcessType	theProcessType = fNotDefined
G4int	theProcessSubType = -1
G4double	thePILfactor = 1.0
G4int	verboseLevel = 0
G4bool	enableAtRestDoIt = true
G4bool	enableAlongStepDoIt = true
G4bool	enablePostStepDoIt = true

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Detailed Description

Definition at line 76 of file G4VEmProcess.hh.

Constructor & Destructor Documentation

G4VEmProcess() [1/2]

G4VEmProcess::G4VEmProcess	(	const G4String &	name,
		G4ProcessType	type = fElectromagnetic )

Definition at line 79 of file G4VEmProcess.cc.

                                                                  :
  G4VDiscreteProcess(name, type)
{
  theParameters = G4EmParameters::Instance();
  SetVerboseLevel(1);
 
  // Size of tables
  minKinEnergy = 0.1*CLHEP::keV;
  maxKinEnergy = 100.0*CLHEP::TeV;
 
  // default lambda factor
  invLambdaFactor = 1.0/lambdaFactor;
 
  // particle types
  theGamma = G4Gamma::Gamma();
  theElectron = G4Electron::Electron();
  thePositron = G4Positron::Positron();
 
  pParticleChange = &fParticleChange;
  fParticleChange.SetSecondaryWeightByProcess(true);
  secParticles.reserve(5);
 
  modelManager = new G4EmModelManager();
  lManager = G4LossTableManager::Instance();
  lManager->Register(this);
  isTheMaster = lManager->IsMaster();
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  theDensityFactor = bld->GetDensityFactors();
  theDensityIdx = bld->GetCoupleIndexes();
}

~G4VEmProcess()

G4VEmProcess::~G4VEmProcess ( )

override

Definition at line 112 of file G4VEmProcess.cc.

{
  if(isTheMaster) {
    delete theData;
    delete theEnergyOfCrossSectionMax;
  }
  delete modelManager;
  delete biasManager;
  lManager->DeRegister(this);
}

G4VEmProcess() [2/2]

G4VEmProcess::G4VEmProcess ( G4VEmProcess & )

delete

Member Function Documentation

ActivateForcedInteraction()

void G4VEmProcess::ActivateForcedInteraction	(	G4double	length = 0.0,
		const G4String &	r = "",
		G4bool	flag = true )

Definition at line 764 of file G4VEmProcess.cc.

{
  if(nullptr == biasManager) { biasManager = new G4EmBiasingManager(); }
  if(1 < verboseLevel) {
    G4cout << "### ActivateForcedInteraction: for " 
           << particle->GetParticleName() 
           << " and process " << GetProcessName()
           << " length(mm)= " << length/mm
           << " in G4Region <" << r 
           << "> weightFlag= " << flag 
           << G4endl; 
  }
  weightFlag = flag;
  biasManager->ActivateForcedInteraction(length, r);
}

Referenced by G4EmExtraParameters::DefineRegParamForEM().

ActivateSecondaryBiasing()

void G4VEmProcess::ActivateSecondaryBiasing	(	const G4String &	region,
		G4double	factor,
		G4double	energyLimit )

Definition at line 784 of file G4VEmProcess.cc.

{
  if (0.0 <= factor) {
 
    // Range cut can be applied only for e-
    if(0.0 == factor && secondaryParticle != G4Electron::Electron())
      { return; }
 
    if(!biasManager) { biasManager = new G4EmBiasingManager(); }
    biasManager->ActivateSecondaryBiasing(region, factor, energyLimit);
    if(1 < verboseLevel) {
      G4cout << "### ActivateSecondaryBiasing: for "
       << " process " << GetProcessName()
       << " factor= " << factor
       << " in G4Region <" << region
       << "> energyLimit(MeV)= " << energyLimit/MeV
       << G4endl;
    }
  }
}

Referenced by G4EmExtraParameters::DefineRegParamForEM().

AddEmModel()

void G4VEmProcess::AddEmModel	(	G4int	order,
		G4VEmModel *	ptr,
		const G4Region *	region = nullptr )

Definition at line 125 of file G4VEmProcess.cc.

{
  if(nullptr == ptr) { return; }
  G4VEmFluctuationModel* fm = nullptr;
  modelManager->AddEmModel(order, ptr, fm, region);
  ptr->SetParticleChange(pParticleChange);
}

Referenced by G4EmDNABuilder::ConstructDNAElectronPhysics(), LBE::ConstructEM(), G4EmLivermorePhysics::ConstructProcess(), G4EmPenelopePhysics::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), G4ComptonScattering::InitialiseProcess(), G4CoulombScattering::InitialiseProcess(), G4DNAAttachment::InitialiseProcess(), G4DNAChargeDecrease::InitialiseProcess(), G4DNAChargeIncrease::InitialiseProcess(), G4DNADissociation::InitialiseProcess(), G4DNAElastic::InitialiseProcess(), G4DNAElectronSolvation::InitialiseProcess(), G4DNAExcitation::InitialiseProcess(), G4DNAIonisation::InitialiseProcess(), G4DNAPlasmonExcitation::InitialiseProcess(), G4DNAPositronium::InitialiseProcess(), G4DNARotExcitation::InitialiseProcess(), G4DNAVibExcitation::InitialiseProcess(), G4eeToHadrons::InitialiseProcess(), G4eplusAnnihilation::InitialiseProcess(), G4GammaConversion::InitialiseProcess(), G4JAEAElasticScattering::InitialiseProcess(), G4MicroElecElastic::InitialiseProcess(), G4MicroElecInelastic::InitialiseProcess(), G4MicroElecLOPhononScattering::InitialiseProcess(), G4NuclearStopping::InitialiseProcess(), G4PhotoElectricEffect::InitialiseProcess(), G4PolarizedCompton::InitialiseProcess(), G4PolarizedGammaConversion::InitialiseProcess(), G4PolarizedPhotoElectric::InitialiseProcess(), G4RayleighScattering::InitialiseProcess(), and G4EmConfigurator::PrepareModels().

ApplyCuts()

G4bool G4VEmProcess::ApplyCuts ( ) const

inlineprotected

Definition at line 615 of file G4VEmProcess.hh.

{
  return applyCuts;
}

Referenced by G4eplusAnnihilation::AtRestDoIt().

BuildLambdaTable()

void G4VEmProcess::BuildLambdaTable

(

)

Definition at line 246 of file G4VEmProcess.cc.

{
  G4double scale = theParameters->MaxKinEnergy()/theParameters->MinKinEnergy();
  G4int nbin = 
    theParameters->NumberOfBinsPerDecade()*G4lrint(std::log10(scale));
  if(actBinning) { nbin = std::max(nbin, nLambdaBins); }
  scale = nbin/G4Log(scale);
  
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  G4EmTableUtil::BuildLambdaTable(this, particle, modelManager,
                                  bld, theLambdaTable, theLambdaTablePrim,
                                  minKinEnergy, minKinEnergyPrim,
                                  maxKinEnergy, scale, verboseLevel,
                                  startFromNull, splineFlag);
}

Referenced by G4EmTableUtil::BuildEmProcess().

BuildPhysicsTable()

void G4VEmProcess::BuildPhysicsTable ( const G4ParticleDefinition & part )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 229 of file G4VEmProcess.cc.

{
  if(nullptr == masterProc) {
    if(isTheMaster) { masterProc = this; }
    else { masterProc = static_cast<const G4VEmProcess*>(GetMasterProcess());}
  }
  G4int nModels = modelManager->NumberOfModels();
  G4bool isLocked = theParameters->IsPrintLocked();
  G4bool toBuild = (buildLambdaTable || minKinEnergyPrim < maxKinEnergy);
 
  G4EmTableUtil::BuildEmProcess(this, masterProc, particle, &part,
                                nModels, verboseLevel, isTheMaster,
                                isLocked, toBuild, baseMat);
}

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), G4PolarizedAnnihilation::BuildPhysicsTable(), and G4PolarizedCompton::BuildPhysicsTable().

ComputeCrossSectionPerAtom()

G4double G4VEmProcess::ComputeCrossSectionPerAtom	(	G4double	kineticEnergy,
		G4double	Z,
		G4double	A = 0.,
		G4double	cut = 0.0 )

Definition at line 700 of file G4VEmProcess.cc.

{
  SelectModel(kinEnergy, currentCoupleIndex);
  return (currentModel) ? 
    currentModel->ComputeCrossSectionPerAtom(currentParticle, kinEnergy,
                                             Z, A, cut) : 0.0;
}

CrossSectionBiasingFactor()

G4double G4VEmProcess::CrossSectionBiasingFactor ( ) const

inline

Definition at line 643 of file G4VEmProcess.hh.

{
  return biasFactor;
}

CrossSectionType()

G4CrossSectionType G4VEmProcess::CrossSectionType ( ) const

inline

Definition at line 714 of file G4VEmProcess.hh.

{
  return fXSType;
}

Referenced by G4EmTableUtil::BuildEmProcess().

CurrentMaterialCutsCoupleIndex()

size_t G4VEmProcess::CurrentMaterialCutsCoupleIndex ( ) const

inlineprotected

Definition at line 445 of file G4VEmProcess.hh.

{
  return currentCoupleIndex;
}

Referenced by G4eplusAnnihilation::AtRestDoIt().

CurrentSetup()

void G4VEmProcess::CurrentSetup ( const G4MaterialCutsCouple * couple, G4double energy )

inline

Definition at line 586 of file G4VEmProcess.hh.

{
  DefineMaterial(couple);
  SelectModel(energy*massRatio, currentCoupleIndex);
}

Referenced by GetCrossSection(), GetLambda(), MeanFreePath(), and G4GammaGeneralProcess::SelectEmProcess().

DefineMaterial()

void G4VEmProcess::DefineMaterial ( const G4MaterialCutsCouple * couple )

inlineprotected

Definition at line 473 of file G4VEmProcess.hh.

{
  if(couple != currentCouple) {
    currentCouple = couple;
    baseMaterial = currentMaterial = couple->GetMaterial();
    basedCoupleIndex = currentCoupleIndex = couple->GetIndex();
    fFactor = biasFactor;
    mfpKinEnergy = DBL_MAX;
    if(baseMat) {
      basedCoupleIndex = (*theDensityIdx)[currentCoupleIndex];
      if(nullptr != currentMaterial->GetBaseMaterial())
        baseMaterial = currentMaterial->GetBaseMaterial();
      fFactor *= (*theDensityFactor)[currentCoupleIndex];
    }
  }
}

Referenced by G4eplusAnnihilation::AtRestDoIt(), CurrentSetup(), LambdaPhysicsVector(), and PostStepGetPhysicalInteractionLength().

DensityFactor()

G4double G4VEmProcess::DensityFactor ( G4int idx ) const

inlineprotected

Definition at line 771 of file G4VEmProcess.hh.

{
  return (*theDensityFactor)[idx];  
}

Referenced by G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength().

DensityIndex()

G4int G4VEmProcess::DensityIndex ( G4int idx ) const

inlineprotected

Definition at line 764 of file G4VEmProcess.hh.

{
  return (*theDensityIdx)[idx];  
}

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), and G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength().

EmModel()

G4VEmModel * G4VEmProcess::EmModel ( size_t index = 0 ) const

inline

Definition at line 806 of file G4VEmProcess.hh.

{
  return (index < emModels.size()) ? emModels[index] : nullptr;
}

Referenced by G4EmCalculator::ComputeNuclearDEDX(), G4EmDNAChemistry::ConstructProcess(), G4EmDNAChemistry_option1::ConstructProcess(), G4EmDNAChemistry_option2::ConstructProcess(), G4EmDNAChemistry_option3::ConstructProcess(), G4ComptonScattering::InitialiseProcess(), G4CoulombScattering::InitialiseProcess(), G4DNAAttachment::InitialiseProcess(), G4DNAChargeDecrease::InitialiseProcess(), G4DNAChargeIncrease::InitialiseProcess(), G4DNADissociation::InitialiseProcess(), G4DNAElastic::InitialiseProcess(), G4DNAElectronSolvation::InitialiseProcess(), G4DNAExcitation::InitialiseProcess(), G4DNAIonisation::InitialiseProcess(), G4DNAPlasmonExcitation::InitialiseProcess(), G4DNAPositronium::InitialiseProcess(), G4DNARotExcitation::InitialiseProcess(), G4DNAVibExcitation::InitialiseProcess(), G4eplusAnnihilation::InitialiseProcess(), G4GammaConversion::InitialiseProcess(), G4JAEAElasticScattering::InitialiseProcess(), G4MicroElecElastic::InitialiseProcess(), G4MicroElecInelastic::InitialiseProcess(), G4MicroElecLOPhononScattering::InitialiseProcess(), G4NuclearStopping::InitialiseProcess(), G4PhotoElectricEffect::InitialiseProcess(), G4PolarizedCompton::InitialiseProcess(), G4PolarizedGammaConversion::InitialiseProcess(), G4PolarizedPhotoElectric::InitialiseProcess(), G4RayleighScattering::InitialiseProcess(), G4DNAAttachment::PrintInfo(), G4DNAChargeDecrease::PrintInfo(), G4DNAChargeIncrease::PrintInfo(), G4DNADissociation::PrintInfo(), G4DNAElastic::PrintInfo(), G4DNAIonisation::PrintInfo(), G4DNAPlasmonExcitation::PrintInfo(), and G4DNAPositronium::PrintInfo().

EnergyOfCrossSectionMax()

std::vector< G4double > * G4VEmProcess::EnergyOfCrossSectionMax ( ) const

inline

Definition at line 678 of file G4VEmProcess.hh.

{
  return theEnergyOfCrossSectionMax;
}

Referenced by G4EmTableUtil::BuildEmProcess().

GetCrossSection()

G4double G4VEmProcess::GetCrossSection ( const G4double kinEnergy, const G4MaterialCutsCouple * couple )

overridevirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 680 of file G4VEmProcess.cc.

{
  CurrentSetup(couple, kinEnergy);
  return GetCurrentLambda(kinEnergy, G4Log(kinEnergy));
}

Referenced by G4EmCalculator::GetCrossSectionPerVolume(), and G4AdjointComptonModel::RapidSampleSecondaries().

GetCurrentElement()

const G4Element * G4VEmProcess::GetCurrentElement

(

)

const

Definition at line 722 of file G4VEmProcess.cc.

{
  return (nullptr != currentModel) ?
    currentModel->GetCurrentElement(currentMaterial) : nullptr; 
}

Referenced by GetTargetIsotope().

GetCurrentModel()

const G4VEmModel * G4VEmProcess::GetCurrentModel ( ) const

inline

Definition at line 785 of file G4VEmProcess.hh.

{
  return currentModel;
}

GetElectronEnergyCut()

G4double G4VEmProcess::GetElectronEnergyCut ( )

inlineprotected

Definition at line 466 of file G4VEmProcess.hh.

{
  return (*theCutsElectron)[currentCoupleIndex];
}

Referenced by G4eplusAnnihilation::AtRestDoIt().

GetEmProcess()

G4VEmProcess * G4VEmProcess::GetEmProcess ( const G4String & name )

virtual

Reimplemented in G4GammaGeneralProcess.

Definition at line 855 of file G4VEmProcess.cc.

{
  return (nam == GetProcessName()) ? this : nullptr;
}

GetGammaEnergyCut()

G4double G4VEmProcess::GetGammaEnergyCut ( )

inlineprotected

Definition at line 459 of file G4VEmProcess.hh.

{
  return (*theCutsGamma)[currentCoupleIndex];
}

Referenced by G4eplusAnnihilation::AtRestDoIt().

GetLambda()

G4double G4VEmProcess::GetLambda ( G4double kinEnergy, const G4MaterialCutsCouple * couple, G4double logKinEnergy )

inline

Definition at line 595 of file G4VEmProcess.hh.

{
  CurrentSetup(couple, kinEnergy);
  return GetCurrentLambda(kinEnergy, logKinEnergy);
}

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), and G4GammaGeneralProcess::TotalCrossSectionPerVolume().

GetMeanFreePath()

G4double G4VEmProcess::GetMeanFreePath ( const G4Track & track, G4double previousStepSize, G4ForceCondition * condition )

overrideprotectedvirtual

Implements G4VDiscreteProcess.

Definition at line 689 of file G4VEmProcess.cc.

{
  *condition = NotForced;
  return G4VEmProcess::MeanFreePath(track);
}

Referenced by G4PolarizedAnnihilation::GetMeanFreePath(), and G4PolarizedCompton::GetMeanFreePath().

GetModelByIndex()

G4VEmModel * G4VEmProcess::GetModelByIndex ( G4int idx = 0, G4bool ver = false ) const

inline

Definition at line 813 of file G4VEmProcess.hh.

{
  return modelManager->GetModel(idx, ver);
}

Referenced by G4EmTableUtil::BuildEmProcess().

GetParticleChange()

G4ParticleChangeForGamma * G4VEmProcess::GetParticleChange ( )

inlineprotected

Definition at line 728 of file G4VEmProcess.hh.

{
  return &fParticleChange;
}

GetTargetElement()

const G4Element * G4VEmProcess::GetTargetElement ( ) const

protected

Definition at line 730 of file G4VEmProcess.cc.

{
  return (nullptr != currentModel) ?
    currentModel->GetCurrentElement(currentMaterial) : nullptr;
}

GetTargetIsotope()

const G4Isotope * G4VEmProcess::GetTargetIsotope ( ) const

protected

Definition at line 738 of file G4VEmProcess.cc.

{
  return (nullptr != currentModel) ?
    currentModel->GetCurrentIsotope(GetCurrentElement()) : nullptr;
}

InitialiseProcess()

virtual void G4VEmProcess::InitialiseProcess ( const G4ParticleDefinition * )

protectedpure virtual

Implemented in G4ComptonScattering, G4CoulombScattering, G4DNAAttachment, G4DNAChargeDecrease, G4DNAChargeIncrease, G4DNADissociation, G4DNAElastic, G4DNAElectronSolvation, G4DNAExcitation, G4DNAIonisation, G4DNAPlasmonExcitation, G4DNAPositronium, G4DNARotExcitation, G4DNAVibExcitation, G4eeToHadrons, G4eplusAnnihilation, G4GammaConversion, G4GammaGeneralProcess, G4JAEAElasticScattering, G4MicroElecElastic, G4MicroElecInelastic, G4MicroElecLOPhononScattering, G4NuclearStopping, G4PhotoElectricEffect, G4PolarizedCompton, G4PolarizedGammaConversion, G4PolarizedPhotoElectric, and G4RayleighScattering.

Referenced by PreparePhysicsTable().

IsApplicable()

virtual G4bool G4VEmProcess::IsApplicable ( const G4ParticleDefinition & p )

overridepure virtual

Reimplemented from G4VProcess.

Implemented in G4ComptonScattering, G4CoulombScattering, G4DNAAttachment, G4DNAChargeDecrease, G4DNAChargeIncrease, G4DNADissociation, G4DNAElastic, G4DNAElectronSolvation, G4DNAExcitation, G4DNAIonisation, G4DNAPlasmonExcitation, G4DNAPositronium, G4DNARotExcitation, G4DNAVibExcitation, G4eeToHadrons, G4eplusAnnihilation, G4GammaConversion, G4GammaGeneralProcess, G4JAEAElasticScattering, G4MicroElecElastic, G4MicroElecInelastic, G4MicroElecLOPhononScattering, G4NuclearStopping, G4PhotoElectricEffect, G4PolarizedCompton, G4PolarizedGammaConversion, G4PolarizedPhotoElectric, and G4RayleighScattering.

LambdaBinning()

G4int G4VEmProcess::LambdaBinning ( ) const

inlineprotected

Definition at line 622 of file G4VEmProcess.hh.

{
  return nLambdaBins;
}

LambdaPhysicsVector()

G4PhysicsVector * G4VEmProcess::LambdaPhysicsVector ( const G4MaterialCutsCouple * couple )

protected

Definition at line 712 of file G4VEmProcess.cc.

{
  DefineMaterial(couple);
  G4PhysicsVector* newv = new G4PhysicsLogVector(minKinEnergy, maxKinEnergy, 
                                                 nLambdaBins, splineFlag);
  return newv;
}

LambdaTable()

G4PhysicsTable * G4VEmProcess::LambdaTable ( ) const

inline

Definition at line 650 of file G4VEmProcess.hh.

{
  return theLambdaTable;
}

Referenced by G4EmTableUtil::BuildEmProcess().

LambdaTablePrim()

G4PhysicsTable * G4VEmProcess::LambdaTablePrim ( ) const

inline

Definition at line 657 of file G4VEmProcess.hh.

{
  return theLambdaTablePrim;
}

Referenced by G4EmTableUtil::BuildEmProcess().

MaterialCutsCouple()

const G4MaterialCutsCouple * G4VEmProcess::MaterialCutsCouple ( ) const

inlineprotected

Definition at line 452 of file G4VEmProcess.hh.

{
  return currentCouple;
}

Referenced by G4eplusAnnihilation::AtRestDoIt().

MaxKinEnergy()

G4double G4VEmProcess::MaxKinEnergy ( ) const

inlineprotected

Definition at line 636 of file G4VEmProcess.hh.

{
  return maxKinEnergy;
}

Referenced by G4eplusAnnihilation::InitialiseProcess(), and SetMinKinEnergyPrim().

MeanFreePath()

G4double G4VEmProcess::MeanFreePath ( const G4Track & track )

inline

Definition at line 604 of file G4VEmProcess.hh.

{
  const G4double kinEnergy = track.GetKineticEnergy();
  CurrentSetup(track.GetMaterialCutsCouple(), kinEnergy);
  const G4double xs = GetCurrentLambda(kinEnergy,
                             track.GetDynamicParticle()->GetLogKineticEnergy());
  return (0.0 < xs) ? 1.0/xs : DBL_MAX; 
}

Referenced by G4GammaGeneralProcess::GetMeanFreePath(), and GetMeanFreePath().

MinKinEnergy()

G4double G4VEmProcess::MinKinEnergy ( ) const

inlineprotected

Definition at line 629 of file G4VEmProcess.hh.

{
  return minKinEnergy;
}

Referenced by G4eplusAnnihilation::InitialiseProcess().

NumberOfModels()

G4int G4VEmProcess::NumberOfModels ( ) const

inline

Definition at line 799 of file G4VEmProcess.hh.

{
  return numberOfModels;
}

operator=()

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

delete

Particle()

const G4ParticleDefinition * G4VEmProcess::Particle ( ) const

inline

Definition at line 693 of file G4VEmProcess.hh.

{
  return particle;
}

PolarAngleLimit()

G4double G4VEmProcess::PolarAngleLimit ( ) const

inlineprotected

Definition at line 862 of file G4VEmProcess.cc.

{
  return theParameters->MscThetaLimit();
}

PostStepDoIt()

G4VParticleChange * G4VEmProcess::PostStepDoIt ( const G4Track & track, const G4Step & step )

overridevirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 460 of file G4VEmProcess.cc.

{
  // clear number of interaction lengths in any case
  theNumberOfInteractionLengthLeft = -1.0;
  mfpKinEnergy = DBL_MAX;
 
  fParticleChange.InitializeForPostStep(track);
 
  // Do not make anything if particle is stopped, the annihilation then
  // should be performed by the AtRestDoIt!
  if (track.GetTrackStatus() == fStopButAlive) { return &fParticleChange; }
 
  const G4double finalT = track.GetKineticEnergy();
 
  // forced process - should happen only once per track
  if(biasFlag) {
    if(biasManager->ForcedInteractionRegion((G4int)currentCoupleIndex)) {
      biasFlag = false;
    }
  }
 
  // check active and select model
  const G4double scaledEnergy = finalT*massRatio;
  SelectModel(scaledEnergy, currentCoupleIndex);
  if(!currentModel->IsActive(scaledEnergy)) { return &fParticleChange; }
 
  // Integral approach
  if (fXSType != fEmNoIntegral) {
    const G4double logFinalT = 
      track.GetDynamicParticle()->GetLogKineticEnergy();
    const G4double lx = std::max(GetCurrentLambda(finalT, logFinalT), 0.0);
#ifdef G4VERBOSE
    if(preStepLambda < lx && 1 < verboseLevel) {
      G4cout << "WARNING: for " << currentParticle->GetParticleName() 
             << " and " << GetProcessName() << " E(MeV)= " << finalT/MeV
             << " preLambda= " << preStepLambda 
             << " < " << lx << " (postLambda) " << G4endl;  
    }
#endif
    // if false interaction then use new cross section value
    // if both values are zero - no interaction
    if(preStepLambda*G4UniformRand() >= lx) {
      return &fParticleChange;
    }
  }
 
  // define new weight for primary and secondaries
  G4double weight = fParticleChange.GetParentWeight();
  if(weightFlag) { 
    weight /= biasFactor; 
    fParticleChange.ProposeWeight(weight);
  }
  
#ifdef G4VERBOSE
  if(1 < verboseLevel) {
    G4cout << "G4VEmProcess::PostStepDoIt: Sample secondary; E= "
           << finalT/MeV
           << " MeV; model= (" << currentModel->LowEnergyLimit()
           << ", " <<  currentModel->HighEnergyLimit() << ")"
           << G4endl;
  }
#endif
 
  // sample secondaries
  secParticles.clear();
  currentModel->SampleSecondaries(&secParticles, 
                                  currentCouple, 
                                  track.GetDynamicParticle(),
                                  (*theCuts)[currentCoupleIndex]);
 
  G4int num0 = (G4int)secParticles.size();
 
  // splitting or Russian roulette
  if(biasManager) {
    if(biasManager->SecondaryBiasingRegion((G4int)currentCoupleIndex)) {
      G4double eloss = 0.0;
      weight *= biasManager->ApplySecondaryBiasing(
        secParticles, track, currentModel, &fParticleChange, eloss, 
        (G4int)currentCoupleIndex, (*theCuts)[currentCoupleIndex],
        step.GetPostStepPoint()->GetSafety());
      if(eloss > 0.0) {
        eloss += fParticleChange.GetLocalEnergyDeposit();
        fParticleChange.ProposeLocalEnergyDeposit(eloss);
      }
    }
  }
 
  // save secondaries
  G4int num = (G4int)secParticles.size();
  if(num > 0) {
 
    fParticleChange.SetNumberOfSecondaries(num);
    G4double edep = fParticleChange.GetLocalEnergyDeposit();
    G4double time = track.GetGlobalTime();
 
    G4int n1(0), n2(0);
    if(num > mainSecondaries) { 
      currentModel->FillNumberOfSecondaries(n1, n2);
    }
     
    for (G4int i=0; i<num; ++i) {
      G4DynamicParticle* dp = secParticles[i];
      if (nullptr != dp) {
        const G4ParticleDefinition* p = dp->GetParticleDefinition();
        G4double e = dp->GetKineticEnergy();
        G4bool good = true;
        if(applyCuts) {
          if (p == theGamma) {
            if (e < (*theCutsGamma)[currentCoupleIndex]) { good = false; }
 
          } else if (p == theElectron) {
            if (e < (*theCutsElectron)[currentCoupleIndex]) { good = false; }
 
          } else if (p == thePositron) {
            if (electron_mass_c2 < (*theCutsGamma)[currentCoupleIndex] &&
                e < (*theCutsPositron)[currentCoupleIndex]) {
              good = false;
              e += 2.0*electron_mass_c2;
            }
          }
          // added secondary if it is good
        }
        if (good) { 
          G4Track* t = new G4Track(dp, time, track.GetPosition());
          t->SetTouchableHandle(track.GetTouchableHandle());
          if (biasManager) {
            t->SetWeight(weight * biasManager->GetWeight(i));
          } else {
            t->SetWeight(weight);
          }
          pParticleChange->AddSecondary(t);
 
          // define type of secondary
          if(i < mainSecondaries) { 
            t->SetCreatorModelID(secID);
            if(GetProcessSubType() == fComptonScattering && p == theGamma) {
              t->SetCreatorModelID(_ComptonGamma);
            }
          } else if(i < mainSecondaries + n1) {
            t->SetCreatorModelID(tripletID);
          } else if(i < mainSecondaries + n1 + n2) {
            t->SetCreatorModelID(_IonRecoil);
          } else {
            if(i < num0) {
              if(p == theGamma) { 
                t->SetCreatorModelID(fluoID);
              } else {
                t->SetCreatorModelID(augerID);
              }
            } else {
              t->SetCreatorModelID(secID);
            }
          }
          /* 
          G4cout << "Secondary(post step) has weight " << t->GetWeight() 
                 << ", Ekin= " << t->GetKineticEnergy()/MeV << " MeV "
                 << GetProcessName() << " fluoID= " << fluoID
                 << " augerID= " << augerID <<G4endl;
          */
        } else {
          delete dp;
          edep += e;
        }
      } 
    }
    fParticleChange.ProposeLocalEnergyDeposit(edep);
  }
 
  if(0.0 == fParticleChange.GetProposedKineticEnergy() &&
     fAlive == fParticleChange.GetTrackStatus()) {
    if(particle->GetProcessManager()->GetAtRestProcessVector()->size() > 0)
         { fParticleChange.ProposeTrackStatus(fStopButAlive); }
    else { fParticleChange.ProposeTrackStatus(fStopAndKill); }
  }
 
  return &fParticleChange;
}

PostStepGetPhysicalInteractionLength()

G4double G4VEmProcess::PostStepGetPhysicalInteractionLength ( const G4Track & track, G4double previousStepSize, G4ForceCondition * condition )

overridevirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 354 of file G4VEmProcess.cc.

{
  *condition = NotForced;
  G4double x = DBL_MAX;
 
  DefineMaterial(track.GetMaterialCutsCouple());
  preStepKinEnergy = track.GetKineticEnergy();
  const G4double scaledEnergy = preStepKinEnergy*massRatio;
  SelectModel(scaledEnergy, currentCoupleIndex);
  /*
  G4cout << "PostStepGetPhysicalInteractionLength: idx= " << currentCoupleIndex
         << "  couple: " << currentCouple << G4endl;
  */
  if(!currentModel->IsActive(scaledEnergy)) { 
    theNumberOfInteractionLengthLeft = -1.0;
    currentInteractionLength = DBL_MAX;
    mfpKinEnergy = DBL_MAX;
    preStepLambda = 0.0;
    return x; 
  }
 
  // forced biasing only for primary particles
  if(biasManager) {
    if(0 == track.GetParentID()) {
      if(biasFlag && 
         biasManager->ForcedInteractionRegion((G4int)currentCoupleIndex)) {
        return biasManager->GetStepLimit((G4int)currentCoupleIndex, previousStepSize);
      }
    }
  }
 
  // compute mean free path
 
  ComputeIntegralLambda(preStepKinEnergy, track);
 
  // zero cross section
  if(preStepLambda <= 0.0) { 
    theNumberOfInteractionLengthLeft = -1.0;
    currentInteractionLength = DBL_MAX;
 
  } else {
 
    // non-zero cross section
    if (theNumberOfInteractionLengthLeft < 0.0) {
 
      // beggining of tracking (or just after DoIt of this process)
      theNumberOfInteractionLengthLeft = -G4Log( G4UniformRand() );
      theInitialNumberOfInteractionLength = theNumberOfInteractionLengthLeft; 
 
    } else {
 
      theNumberOfInteractionLengthLeft -= 
        previousStepSize/currentInteractionLength;
      theNumberOfInteractionLengthLeft = 
        std::max(theNumberOfInteractionLengthLeft, 0.0);
    }
 
    // new mean free path and step limit for the next step
    currentInteractionLength = 1.0/preStepLambda;
    x = theNumberOfInteractionLengthLeft * currentInteractionLength;
  }
  return x;
}

Referenced by G4PolarizedAnnihilation::PostStepGetPhysicalInteractionLength(), and G4PolarizedCompton::PostStepGetPhysicalInteractionLength().

PreparePhysicsTable()

void G4VEmProcess::PreparePhysicsTable ( const G4ParticleDefinition & part )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 147 of file G4VEmProcess.cc.

{
  if(nullptr == particle) { SetParticle(&part); }
 
  if(part.GetParticleType() == "nucleus" && 
     part.GetParticleSubType() == "generic") {
 
    G4String pname = part.GetParticleName();
    if(pname != "deuteron" && pname != "triton" &&
       pname != "alpha" && pname != "alpha+" &&
       pname != "helium" && pname != "hydrogen") {
 
      particle = G4GenericIon::GenericIon();
      isIon = true;
    }
  }
  if(particle != &part) { return; }
 
  lManager->PreparePhysicsTable(&part, this);
 
  // for new run
  currentCouple = nullptr;
  preStepLambda = 0.0;
  fLambdaEnergy = 0.0;
 
  InitialiseProcess(particle);
 
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  const G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  theCutsGamma    = theCoupleTable->GetEnergyCutsVector(idxG4GammaCut);
  theCutsElectron = theCoupleTable->GetEnergyCutsVector(idxG4ElectronCut);
  theCutsPositron = theCoupleTable->GetEnergyCutsVector(idxG4PositronCut);
 
  // initialisation of the process  
  if(!actMinKinEnergy) { minKinEnergy = theParameters->MinKinEnergy(); }
  if(!actMaxKinEnergy) { maxKinEnergy = theParameters->MaxKinEnergy(); }
 
  applyCuts       = theParameters->ApplyCuts();
  lambdaFactor    = theParameters->LambdaFactor();
  invLambdaFactor = 1.0/lambdaFactor;
  theParameters->DefineRegParamForEM(this);
 
  // integral option may be disabled
  if(!theParameters->Integral()) { fXSType = fEmNoIntegral; }
 
  // prepare tables
  if(isTheMaster) {
    if(nullptr == theData) { theData = new G4EmDataHandler(2); }
 
    if(buildLambdaTable) {
      theLambdaTable = theData->MakeTable(0);
      bld->InitialiseBaseMaterials(theLambdaTable);
    }
    // high energy table
    if(minKinEnergyPrim < maxKinEnergy) {
      theLambdaTablePrim = theData->MakeTable(1);
      bld->InitialiseBaseMaterials(theLambdaTablePrim);
    }
  }
  // models
  baseMat = bld->GetBaseMaterialFlag();
  numberOfModels = modelManager->NumberOfModels();
  currentModel = modelManager->GetModel(0);
  if(nullptr != lManager->AtomDeexcitation()) { 
    modelManager->SetFluoFlag(true); 
  }
  // forced biasing
  if(nullptr != biasManager) { 
    biasManager->Initialise(part, GetProcessName(), verboseLevel); 
    biasFlag = false;
  }
 
  theCuts =
    G4EmTableUtil::PrepareEmProcess(this, particle, secondaryParticle,
                                    modelManager, maxKinEnergy, 
                                    secID, tripletID, mainSecondaries,
                                    verboseLevel, isTheMaster);
}

Referenced by G4GammaGeneralProcess::PreparePhysicsTable().

ProcessDescription()

void G4VEmProcess::ProcessDescription ( std::ostream & outFile ) const

overridevirtual

Reimplemented from G4VProcess.

Definition at line 883 of file G4VEmProcess.cc.

{
  if(nullptr != particle) {
    StreamInfo(out, *particle, true);
  }
}

Referenced by G4ComptonScattering::ProcessDescription(), G4CoulombScattering::ProcessDescription(), G4DNAExcitation::ProcessDescription(), G4DNARotExcitation::ProcessDescription(), G4DNAVibExcitation::ProcessDescription(), G4eeToHadrons::ProcessDescription(), G4eplusAnnihilation::ProcessDescription(), G4GammaConversion::ProcessDescription(), G4GammaGeneralProcess::ProcessDescription(), G4NuclearStopping::ProcessDescription(), G4PhotoElectricEffect::ProcessDescription(), G4PolarizedCompton::ProcessDescription(), G4PolarizedGammaConversion::ProcessDescription(), G4PolarizedPhotoElectric::ProcessDescription(), and G4RayleighScattering::ProcessDescription().

RetrievePhysicsTable()

G4bool G4VEmProcess::RetrievePhysicsTable ( const G4ParticleDefinition * part, const G4String & directory, G4bool ascii )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 659 of file G4VEmProcess.cc.

{
  if(!isTheMaster || part != particle) { return true; }
  G4bool yes = true;
  if(buildLambdaTable) {
    yes = G4EmTableUtil::RetrieveTable(this, part, theLambdaTable, dir,
                                       "Lambda", verboseLevel,
                                       ascii, splineFlag);
  }
  if(yes && minKinEnergyPrim < maxKinEnergy) {
    yes = G4EmTableUtil::RetrieveTable(this, part, theLambdaTablePrim, dir,
                                       "LambdaPrim", verboseLevel,
                                       ascii, splineFlag);
  }
  return yes;
}

Referenced by G4GammaGeneralProcess::RetrievePhysicsTable().

SecondaryParticle()

const G4ParticleDefinition * G4VEmProcess::SecondaryParticle ( ) const

inline

Definition at line 700 of file G4VEmProcess.hh.

{
  return secondaryParticle;
}

SelectModel()

G4VEmModel * G4VEmProcess::SelectModel ( G4double kinEnergy, size_t  )

inlineprotected

Definition at line 493 of file G4VEmProcess.hh.

{
  if(1 < numberOfModels) {
    currentModel = modelManager->SelectModel(kinEnergy, currentCoupleIndex);
  }
  currentModel->SetCurrentCouple(currentCouple);
  return currentModel;
}

Referenced by G4NuclearStopping::AlongStepDoIt(), G4eplusAnnihilation::AtRestDoIt(), ComputeCrossSectionPerAtom(), CurrentSetup(), PostStepDoIt(), and PostStepGetPhysicalInteractionLength().

SelectModelForMaterial()

G4VEmModel * G4VEmProcess::SelectModelForMaterial ( G4double kinEnergy, size_t idxCouple ) const

inline

Definition at line 505 of file G4VEmProcess.hh.

{
  return modelManager->SelectModel(kinEnergy, idxCouple);
}

SetBuildTableFlag()

void G4VEmProcess::SetBuildTableFlag ( G4bool val )

inline

Definition at line 721 of file G4VEmProcess.hh.

{
  buildLambdaTable = val;
}

Referenced by G4ComptonScattering::G4ComptonScattering(), G4CoulombScattering::G4CoulombScattering(), G4DNARotExcitation::G4DNARotExcitation(), G4eeToHadrons::G4eeToHadrons(), G4eplusAnnihilation::G4eplusAnnihilation(), G4GammaConversion::G4GammaConversion(), G4JAEAElasticScattering::G4JAEAElasticScattering(), G4NuclearStopping::G4NuclearStopping(), G4PhotoElectricEffect::G4PhotoElectricEffect(), G4PolarizedCompton::G4PolarizedCompton(), G4PolarizedGammaConversion::G4PolarizedGammaConversion(), G4PolarizedPhotoElectric::G4PolarizedPhotoElectric(), G4RayleighScattering::G4RayleighScattering(), G4CoulombScattering::InitialiseProcess(), G4DNAAttachment::InitialiseProcess(), G4DNAChargeDecrease::InitialiseProcess(), G4DNAChargeIncrease::InitialiseProcess(), G4DNADissociation::InitialiseProcess(), G4DNAElastic::InitialiseProcess(), G4DNAElectronSolvation::InitialiseProcess(), G4DNAExcitation::InitialiseProcess(), G4DNAIonisation::InitialiseProcess(), G4DNAPlasmonExcitation::InitialiseProcess(), G4DNAPositronium::InitialiseProcess(), G4DNAVibExcitation::InitialiseProcess(), G4MicroElecElastic::InitialiseProcess(), G4MicroElecInelastic::InitialiseProcess(), and G4MicroElecLOPhononScattering::InitialiseProcess().

SetCrossSectionBiasingFactor()

void G4VEmProcess::SetCrossSectionBiasingFactor	(	G4double	f,
		G4bool	flag = true )

Definition at line 746 of file G4VEmProcess.cc.

{
  if(f > 0.0) { 
    biasFactor = f; 
    weightFlag = flag;
    if(1 < verboseLevel) {
      G4cout << "### SetCrossSectionBiasingFactor: for " 
             << particle->GetParticleName() 
             << " and process " << GetProcessName()
             << " biasFactor= " << f << " weightFlag= " << flag 
             << G4endl; 
    }
  }
}

Referenced by G4EmExtraParameters::DefineRegParamForEM().

SetCrossSectionType()

void G4VEmProcess::SetCrossSectionType ( G4CrossSectionType val )

inline

Definition at line 707 of file G4VEmProcess.hh.

{
  fXSType = val; 
}

Referenced by G4EmTableUtil::BuildEmProcess(), G4CoulombScattering::G4CoulombScattering(), G4eeToHadrons::G4eeToHadrons(), G4eplusAnnihilation::G4eplusAnnihilation(), and G4CoulombScattering::InitialiseProcess().

SetEmMasterProcess()

void G4VEmProcess::SetEmMasterProcess ( const G4VEmProcess * ptr )

inline

Definition at line 792 of file G4VEmProcess.hh.

{ 
  masterProc = ptr;
}

Referenced by G4GammaGeneralProcess::BuildPhysicsTable().

SetEmModel()

void G4VEmProcess::SetEmModel	(	G4VEmModel *	ptr,
		G4int	index = 0 )

Definition at line 136 of file G4VEmProcess.cc.

{
  if(nullptr == ptr) { return; }
  if(!emModels.empty()) {
    for(auto & em : emModels) { if(em == ptr) { return; } }
  }
  emModels.push_back(ptr);
}

Referenced by G4EmBuilder::ConstructElectronSSProcess(), G4EmLivermorePhysics::ConstructProcess(), G4EmLowEPPhysics::ConstructProcess(), G4EmPenelopePhysics::ConstructProcess(), G4EmStandardPhysics::ConstructProcess(), G4EmStandardPhysics_option1::ConstructProcess(), G4EmStandardPhysics_option2::ConstructProcess(), G4EmStandardPhysics_option3::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), G4EmStandardPhysicsGS::ConstructProcess(), G4EmStandardPhysicsSS::ConstructProcess(), G4EmStandardPhysicsWVI::ConstructProcess(), G4EmDNABuilder::FindOrBuildAttachment(), G4EmDNABuilder::FindOrBuildChargeDecrease(), G4EmDNABuilder::FindOrBuildChargeIncrease(), G4EmDNABuilder::FindOrBuildElastic(), G4EmDNABuilder::FindOrBuildElectronSolvation(), G4EmDNABuilder::FindOrBuildExcitation(), G4EmDNABuilder::FindOrBuildIonisation(), G4EmDNABuilder::FindOrBuildVibExcitation(), G4PolarizedAnnihilation::G4PolarizedAnnihilation(), G4ComptonScattering::InitialiseProcess(), G4CoulombScattering::InitialiseProcess(), G4DNAAttachment::InitialiseProcess(), G4DNAChargeDecrease::InitialiseProcess(), G4DNAChargeIncrease::InitialiseProcess(), G4DNADissociation::InitialiseProcess(), G4DNAElastic::InitialiseProcess(), G4DNAElectronSolvation::InitialiseProcess(), G4DNAExcitation::InitialiseProcess(), G4DNAIonisation::InitialiseProcess(), G4DNAPlasmonExcitation::InitialiseProcess(), G4DNAPositronium::InitialiseProcess(), G4DNARotExcitation::InitialiseProcess(), G4DNAVibExcitation::InitialiseProcess(), G4eeToHadrons::InitialiseProcess(), G4eplusAnnihilation::InitialiseProcess(), G4GammaConversion::InitialiseProcess(), G4JAEAElasticScattering::InitialiseProcess(), G4MicroElecElastic::InitialiseProcess(), G4MicroElecInelastic::InitialiseProcess(), G4MicroElecLOPhononScattering::InitialiseProcess(), G4NuclearStopping::InitialiseProcess(), G4PhotoElectricEffect::InitialiseProcess(), G4PolarizedCompton::InitialiseProcess(), G4PolarizedGammaConversion::InitialiseProcess(), G4PolarizedPhotoElectric::InitialiseProcess(), G4RayleighScattering::InitialiseProcess(), and G4EmTableUtil::PrepareEmProcess().

SetEnergyOfCrossSectionMax()

void G4VEmProcess::SetEnergyOfCrossSectionMax ( std::vector< G4double > * ptr )

inline

Definition at line 686 of file G4VEmProcess.hh.

{
  theEnergyOfCrossSectionMax = ptr;
}

Referenced by G4EmTableUtil::BuildEmProcess().

SetLambdaBinning()

void G4VEmProcess::SetLambdaBinning

(

G4int

nbins

)

Definition at line 809 of file G4VEmProcess.cc.

{
  if(5 < n && n < 10000000) {  
    nLambdaBins = n; 
    actBinning = true;
  } else { 
    G4double e = (G4double)n;
    PrintWarning("SetLambdaBinning", e); 
  } 
}

Referenced by G4GammaConversion::G4GammaConversion(), and G4PolarizedGammaConversion::G4PolarizedGammaConversion().

SetLambdaTable()

void G4VEmProcess::SetLambdaTable ( G4PhysicsTable * ptr )

inline

Definition at line 664 of file G4VEmProcess.hh.

{
  theLambdaTable = ptr;
}

Referenced by G4EmTableUtil::BuildEmProcess().

SetLambdaTablePrim()

void G4VEmProcess::SetLambdaTablePrim ( G4PhysicsTable * ptr )

inline

Definition at line 671 of file G4VEmProcess.hh.

{
  theLambdaTablePrim = ptr;
}

Referenced by G4EmTableUtil::BuildEmProcess().

SetMaxKinEnergy()

void G4VEmProcess::SetMaxKinEnergy

(

G4double

e

)

Definition at line 834 of file G4VEmProcess.cc.

{
  if(minKinEnergy < e && e < 1.e+6*TeV) { 
    nLambdaBins = G4lrint(nLambdaBins*G4Log(e/minKinEnergy)
                          /G4Log(maxKinEnergy/minKinEnergy));
    maxKinEnergy = e;
    actMaxKinEnergy = true;
  } else { PrintWarning("SetMaxKinEnergy", e); } 
}

Referenced by G4EmLivermorePhysics::ConstructProcess(), G4EmLowEPPhysics::ConstructProcess(), G4EmPenelopePhysics::ConstructProcess(), G4EmStandardPhysics::ConstructProcess(), G4EmStandardPhysics_option3::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), and G4EmStandardPhysicsWVI::ConstructProcess().

SetMinKinEnergy()

void G4VEmProcess::SetMinKinEnergy

(

G4double

e

)

Definition at line 822 of file G4VEmProcess.cc.

{
  if(1.e-3*eV < e && e < maxKinEnergy) { 
    nLambdaBins = G4lrint(nLambdaBins*G4Log(maxKinEnergy/e)
                          /G4Log(maxKinEnergy/minKinEnergy));
    minKinEnergy = e;
    actMinKinEnergy = true;
  } else { PrintWarning("SetMinKinEnergy", e); } 
}

Referenced by G4EmLivermorePhysics::ConstructProcess(), G4EmLowEPPhysics::ConstructProcess(), G4EmPenelopePhysics::ConstructProcess(), G4EmStandardPhysics::ConstructProcess(), G4EmStandardPhysics_option1::ConstructProcess(), G4EmStandardPhysics_option2::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), G4EmStandardPhysicsGS::ConstructProcess(), G4EmStandardPhysicsWVI::ConstructProcess(), G4GammaConversion::G4GammaConversion(), G4PolarizedGammaConversion::G4PolarizedGammaConversion(), and G4GammaConversion::InitialiseProcess().

SetMinKinEnergyPrim()

void G4VEmProcess::SetMinKinEnergyPrim

(

G4double

e

)

Definition at line 846 of file G4VEmProcess.cc.

{
  if(theParameters->MinKinEnergy() <= e && 
     e <= theParameters->MaxKinEnergy()) { minKinEnergyPrim = e; } 
  else { PrintWarning("SetMinKinEnergyPrim", e); } 
}

Referenced by G4ComptonScattering::G4ComptonScattering(), G4JAEAElasticScattering::G4JAEAElasticScattering(), G4PhotoElectricEffect::G4PhotoElectricEffect(), G4PolarizedCompton::G4PolarizedCompton(), and G4RayleighScattering::G4RayleighScattering().

SetParticle()

void G4VEmProcess::SetParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 735 of file G4VEmProcess.hh.

{
  particle = p;
  currentParticle = p;
}

Referenced by G4GammaGeneralProcess::G4GammaGeneralProcess(), G4eeToHadrons::InitialiseProcess(), G4GammaGeneralProcess::PreparePhysicsTable(), and PreparePhysicsTable().

SetSecondaryParticle()

void G4VEmProcess::SetSecondaryParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 743 of file G4VEmProcess.hh.

{
  secondaryParticle = p;
}

Referenced by G4ComptonScattering::G4ComptonScattering(), G4CoulombScattering::G4CoulombScattering(), G4eeToHadrons::G4eeToHadrons(), G4eplusAnnihilation::G4eplusAnnihilation(), G4GammaConversion::G4GammaConversion(), G4PhotoElectricEffect::G4PhotoElectricEffect(), G4PolarizedCompton::G4PolarizedCompton(), G4PolarizedGammaConversion::G4PolarizedGammaConversion(), and G4PolarizedPhotoElectric::G4PolarizedPhotoElectric().

SetSplineFlag()

void G4VEmProcess::SetSplineFlag ( G4bool val )

inlineprotected

Definition at line 757 of file G4VEmProcess.hh.

{
  splineFlag = val;
}

Referenced by G4ComptonScattering::G4ComptonScattering(), G4CoulombScattering::G4CoulombScattering(), G4JAEAElasticScattering::G4JAEAElasticScattering(), G4PolarizedCompton::G4PolarizedCompton(), G4RayleighScattering::G4RayleighScattering(), and G4CoulombScattering::InitialiseProcess().

SetStartFromNullFlag()

void G4VEmProcess::SetStartFromNullFlag ( G4bool val )

inlineprotected

Definition at line 750 of file G4VEmProcess.hh.

{
  startFromNull = val;
}

Referenced by G4ComptonScattering::G4ComptonScattering(), G4CoulombScattering::G4CoulombScattering(), G4eplusAnnihilation::G4eplusAnnihilation(), G4GammaConversion::G4GammaConversion(), G4JAEAElasticScattering::G4JAEAElasticScattering(), G4PolarizedCompton::G4PolarizedCompton(), G4RayleighScattering::G4RayleighScattering(), and G4CoulombScattering::InitialiseProcess().

StartTracking()

void G4VEmProcess::StartTracking ( G4Track * track )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 332 of file G4VEmProcess.cc.

{
  // reset parameters for the new track
  currentParticle = track->GetParticleDefinition();
  theNumberOfInteractionLengthLeft = -1.0;
  mfpKinEnergy = DBL_MAX;
  preStepLambda = 0.0;
 
  if(isIon) { massRatio = proton_mass_c2/currentParticle->GetPDGMass(); }
 
  // forced biasing only for primary particles
  if(biasManager) {
    if(0 == track->GetParentID()) {
      // primary particle
      biasFlag = true; 
      biasManager->ResetForcedInteraction(); 
    }
  }
}

StorePhysicsTable()

G4bool G4VEmProcess::StorePhysicsTable ( const G4ParticleDefinition * part, const G4String & directory, G4bool ascii = false )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 641 of file G4VEmProcess.cc.

{
  if(!isTheMaster || part != particle) { return true; }
  if(G4EmTableUtil::StoreTable(this, part, theLambdaTable,
                   directory, "Lambda",
                               verboseLevel, ascii) &&
     G4EmTableUtil::StoreTable(this, part, theLambdaTablePrim,
                   directory, "LambdaPrim",
                               verboseLevel, ascii)) { 
     return true;
  }
  return false;
}

Referenced by G4GammaGeneralProcess::StorePhysicsTable().

StreamInfo()

void G4VEmProcess::StreamInfo	(	std::ostream &	outFile,
		const G4ParticleDefinition &	part,
		G4bool	rst = false ) const

Definition at line 264 of file G4VEmProcess.cc.

{
  G4String indent = (rst ? "  " : "");
  out << std::setprecision(6);
  out << G4endl << indent << GetProcessName() << ": ";
  if (!rst) {
    out << " for " << part.GetParticleName();
  }
  if(fXSType != fEmNoIntegral)  { out << " XStype:" << fXSType; }
  if(applyCuts) { out << " applyCuts:1 "; }
  out << " SubType=" << GetProcessSubType();
  if(biasFactor != 1.0) { out << "  BiasingFactor= " << biasFactor; }
  out << " BuildTable=" << buildLambdaTable << G4endl;
  if(buildLambdaTable) {
    if(particle == &part) { 
      for(auto & v : *theLambdaTable) {
        if(nullptr != v) {
          out << "      Lambda table from ";
          G4double emin = v->Energy(0);
          G4double emax = v->GetMaxEnergy();
          G4int nbin = G4int(v->GetVectorLength() - 1);
          if(emin > minKinEnergy) { out << "threshold "; }
          else { out << G4BestUnit(emin,"Energy"); } 
          out << " to "
              << G4BestUnit(emax,"Energy")
              << ", " << G4lrint(nbin/std::log10(emax/emin))
              << " bins/decade, spline: " 
              << splineFlag << G4endl;
          break;
        }
      }
    } else {
      out << "      Used Lambda table of " 
      << particle->GetParticleName() << G4endl;
    }
  }
  if(minKinEnergyPrim < maxKinEnergy) {
    if(particle == &part) {
      for(auto & v : *theLambdaTablePrim) {
        if(nullptr != v) {
          out << "      LambdaPrime table from "
              << G4BestUnit(v->Energy(0),"Energy") 
              << " to "
              << G4BestUnit(v->GetMaxEnergy(),"Energy")
              << " in " << v->GetVectorLength()-1
              << " bins " << G4endl;
          break;
        }
      }
    } else {
      out << "      Used LambdaPrime table of " 
               << particle->GetParticleName() << G4endl;
    }
  }
  StreamProcessInfo(out);
  modelManager->DumpModelList(out, verboseLevel);
 
  if(verboseLevel > 2 && buildLambdaTable) {
    out << "      LambdaTable address= " << theLambdaTable << G4endl;
    if(theLambdaTable && particle == &part) { 
      out << (*theLambdaTable) << G4endl;
    }
  }
}

Referenced by G4EmTableUtil::BuildEmProcess(), and ProcessDescription().

StreamProcessInfo()

virtual void G4VEmProcess::StreamProcessInfo ( std::ostream & ) const

inlineprotectedvirtual

Reimplemented in G4CoulombScattering, G4eeToHadrons, and G4eplusAnnihilation.

Definition at line 94 of file G4VEmProcess.hh.

{};

Referenced by StreamInfo().

UseBaseMaterial()

G4bool G4VEmProcess::UseBaseMaterial ( ) const

inline

Definition at line 778 of file G4VEmProcess.hh.

{
  return baseMat;
}

Referenced by G4EmTableUtil::BuildEmProcess(), and G4GammaGeneralProcess::BuildPhysicsTable().

Member Data Documentation

augerID

G4int G4VEmProcess::augerID = _AugerElectron

protected

Definition at line 404 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt(), and PostStepDoIt().

basedCoupleIndex

size_t G4VEmProcess::basedCoupleIndex = 0

protected

Definition at line 408 of file G4VEmProcess.hh.

Referenced by G4GammaGeneralProcess::ComputeGeneralLambda(), DefineMaterial(), G4GammaGeneralProcess::GetProbability(), and G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength().

baseMat

G4bool G4VEmProcess::baseMat = false

protected

Definition at line 413 of file G4VEmProcess.hh.

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), BuildPhysicsTable(), DefineMaterial(), G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength(), G4GammaGeneralProcess::PreparePhysicsTable(), PreparePhysicsTable(), and UseBaseMaterial().

biasID

G4int G4VEmProcess::biasID = _EM

protected

Definition at line 405 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt().

biasManager

G4EmBiasingManager* G4VEmProcess::biasManager = nullptr

protected

Definition at line 366 of file G4VEmProcess.hh.

Referenced by ActivateForcedInteraction(), ActivateSecondaryBiasing(), G4eplusAnnihilation::AtRestDoIt(), PostStepDoIt(), PostStepGetPhysicalInteractionLength(), PreparePhysicsTable(), StartTracking(), and ~G4VEmProcess().

coupleIdxLambda

size_t G4VEmProcess::coupleIdxLambda = 0

protected

Definition at line 409 of file G4VEmProcess.hh.

currentCouple

const G4MaterialCutsCouple* G4VEmProcess::currentCouple = nullptr

protected

Definition at line 364 of file G4VEmProcess.hh.

Referenced by DefineMaterial(), MaterialCutsCouple(), PostStepDoIt(), G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength(), G4GammaGeneralProcess::PreparePhysicsTable(), PreparePhysicsTable(), G4GammaGeneralProcess::SelectEmProcess(), SelectModel(), and G4GammaGeneralProcess::TotalCrossSectionPerVolume().

currentCoupleIndex

size_t G4VEmProcess::currentCoupleIndex = 0

protected

Definition at line 407 of file G4VEmProcess.hh.

Referenced by ComputeCrossSectionPerAtom(), CurrentMaterialCutsCoupleIndex(), CurrentSetup(), DefineMaterial(), GetElectronEnergyCut(), GetGammaEnergyCut(), PostStepDoIt(), G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength(), PostStepGetPhysicalInteractionLength(), and SelectModel().

currentMaterial

const G4Material* G4VEmProcess::currentMaterial = nullptr

protected

Definition at line 365 of file G4VEmProcess.hh.

Referenced by DefineMaterial(), GetCurrentElement(), GetTargetElement(), G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength(), and G4GammaGeneralProcess::SelectHadProcess().

fluoID

G4int G4VEmProcess::fluoID = _Fluorescence

protected

Definition at line 403 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt(), and PostStepDoIt().

fParticleChange

G4ParticleChangeForGamma G4VEmProcess::fParticleChange

protected

Definition at line 432 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt(), G4VEmProcess(), GetParticleChange(), G4GammaGeneralProcess::PostStepDoIt(), and PostStepDoIt().

idxLambda

size_t G4VEmProcess::idxLambda = 0

protected

Definition at line 410 of file G4VEmProcess.hh.

isTheMaster

G4bool G4VEmProcess::isTheMaster = true

protected

Definition at line 412 of file G4VEmProcess.hh.

Referenced by G4GammaGeneralProcess::BuildPhysicsTable(), G4PolarizedAnnihilation::BuildPhysicsTable(), BuildPhysicsTable(), G4VEmProcess(), G4GammaGeneralProcess::InitialiseProcess(), G4GammaGeneralProcess::PreparePhysicsTable(), PreparePhysicsTable(), RetrievePhysicsTable(), G4GammaGeneralProcess::StorePhysicsTable(), StorePhysicsTable(), G4GammaGeneralProcess::~G4GammaGeneralProcess(), and ~G4VEmProcess().

mainSecondaries

G4int G4VEmProcess::mainSecondaries = 1

protected

Definition at line 401 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt(), G4eplusAnnihilation::G4eplusAnnihilation(), PostStepDoIt(), and PreparePhysicsTable().

mfpKinEnergy

G4double G4VEmProcess::mfpKinEnergy = DBL_MAX

protected

Definition at line 388 of file G4VEmProcess.hh.

Referenced by DefineMaterial(), PostStepDoIt(), PostStepGetPhysicalInteractionLength(), and StartTracking().

preStepKinEnergy

G4double G4VEmProcess::preStepKinEnergy = 0.0

protected

Definition at line 389 of file G4VEmProcess.hh.

Referenced by G4GammaGeneralProcess::ComputeGeneralLambda(), G4GammaGeneralProcess::GetProbability(), G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength(), PostStepGetPhysicalInteractionLength(), G4GammaGeneralProcess::SelectEmProcess(), and G4GammaGeneralProcess::TotalCrossSectionPerVolume().

preStepLambda

G4double G4VEmProcess::preStepLambda = 0.0

protected

Definition at line 390 of file G4VEmProcess.hh.

Referenced by G4GammaGeneralProcess::PostStepDoIt(), PostStepDoIt(), G4GammaGeneralProcess::PostStepGetPhysicalInteractionLength(), PostStepGetPhysicalInteractionLength(), G4GammaGeneralProcess::PreparePhysicsTable(), PreparePhysicsTable(), and StartTracking().

secID

G4int G4VEmProcess::secID = _EM

protected

Definition at line 402 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt(), PostStepDoIt(), and PreparePhysicsTable().

secParticles

std::vector<G4DynamicParticle*> G4VEmProcess::secParticles

protected

Definition at line 431 of file G4VEmProcess.hh.

Referenced by G4eplusAnnihilation::AtRestDoIt(), G4VEmProcess(), and PostStepDoIt().

theEnergyOfCrossSectionMax

std::vector<G4double>* G4VEmProcess::theEnergyOfCrossSectionMax = nullptr

protected

Definition at line 367 of file G4VEmProcess.hh.

Referenced by EnergyOfCrossSectionMax(), SetEnergyOfCrossSectionMax(), and ~G4VEmProcess().

tripletID

G4int G4VEmProcess::tripletID = _TripletElectron

protected

Definition at line 406 of file G4VEmProcess.hh.

Referenced by PostStepDoIt(), and PreparePhysicsTable().

---

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

• G4VEmProcess.hh
• G4VEmProcess.cc