#include <G4VEnergyLossProcess.hh>

Inheritance diagram for G4VEnergyLossProcess:

Public Member Functions
	G4VEnergyLossProcess (const G4String &name="EnergyLoss", G4ProcessType type=fElectromagnetic)

virtual	~G4VEnergyLossProcess ()

virtual G4bool	IsApplicable (const G4ParticleDefinition &p) override=0

virtual void	PrintInfo ()

virtual void	ProcessDescription (std::ostream &outFile) const override

virtual void	PreparePhysicsTable (const G4ParticleDefinition &) override

virtual void	BuildPhysicsTable (const G4ParticleDefinition &) override

G4PhysicsTable *	BuildDEDXTable (G4EmTableType tType=fRestricted)

G4PhysicsTable *	BuildLambdaTable (G4EmTableType tType=fRestricted)

virtual void	StartTracking (G4Track *) override

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety, G4GPILSelection *selection) override

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &) override

G4double	SampleSubCutSecondaries (std::vector< G4Track * > &, const G4Step &, G4VEmModel *model, G4int matIdx)

virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &) override

virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii=false) override

virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii) override

G4double	GetDEDXDispersion (const G4MaterialCutsCouple couple, const G4DynamicParticle dp, G4double length)

G4double	CrossSectionPerVolume (G4double kineticEnergy, const G4MaterialCutsCouple *couple)

G4double	CrossSectionPerVolume (G4double kineticEnergy, const G4MaterialCutsCouple *couple, G4double logKineticEnergy)

G4double	MeanFreePath (const G4Track &track)

G4double	ContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)

G4VEmModel *	SelectModelForMaterial (G4double kinEnergy, size_t &idx) const

void	AddEmModel (G4int, G4VEmModel , G4VEmFluctuationModel fluc=0, const G4Region *region=nullptr)

void	UpdateEmModel (const G4String &, G4double, G4double)

void	SetEmModel (G4VEmModel *, G4int index=0)

G4VEmModel *	EmModel (size_t index=0) const

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

G4int	NumberOfModels () const

void	SetFluctModel (G4VEmFluctuationModel *)

G4VEmFluctuationModel *	FluctModel ()

void	SetBaseParticle (const G4ParticleDefinition *p)

const G4ParticleDefinition *	Particle () const

const G4ParticleDefinition *	BaseParticle () const

const G4ParticleDefinition *	SecondaryParticle () const

void	ActivateSubCutoff (G4bool val, const G4Region *region=nullptr)

void	SetCrossSectionBiasingFactor (G4double f, G4bool flag=true)

void	ActivateForcedInteraction (G4double length, const G4String &region, G4bool flag=true)

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

void	AddCollaborativeProcess (G4VEnergyLossProcess *)

void	SetLossFluctuations (G4bool val)

void	SetIntegral (G4bool val)

G4bool	IsIntegral () const

void	SetIonisation (G4bool val)

G4bool	IsIonisationProcess () const

void	SetLinearLossLimit (G4double val)

void	SetStepFunction (G4double v1, G4double v2)

void	SetLowestEnergyLimit (G4double)

G4int	NumberOfSubCutoffRegions () const

void	SetDEDXTable (G4PhysicsTable *p, G4EmTableType tType)

void	SetCSDARangeTable (G4PhysicsTable *pRange)

void	SetRangeTableForLoss (G4PhysicsTable *p)

void	SetSecondaryRangeTable (G4PhysicsTable *p)

void	SetInverseRangeTable (G4PhysicsTable *p)

void	SetLambdaTable (G4PhysicsTable *p)

void	SetSubLambdaTable (G4PhysicsTable *p)

void	SetDEDXBinning (G4int nbins)

void	SetMinKinEnergy (G4double e)

G4double	MinKinEnergy () const

void	SetMaxKinEnergy (G4double e)

G4double	MaxKinEnergy () const

G4double	CrossSectionBiasingFactor () const

G4double	GetDEDX (G4double kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetDEDX (G4double kineticEnergy, const G4MaterialCutsCouple *, G4double logKineticEnergy)

G4double	GetDEDXForSubsec (G4double kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetRange (G4double kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetRange (G4double kineticEnergy, const G4MaterialCutsCouple *, G4double logKineticEnergy)

G4double	GetCSDARange (G4double kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetRangeForLoss (G4double kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetRangeForLoss (G4double kineticEnergy, const G4MaterialCutsCouple *, G4double logKineticEnergy)

G4double	GetKineticEnergy (G4double range, const G4MaterialCutsCouple *)

G4double	GetLambda (G4double kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetLambda (G4double kineticEnergy, const G4MaterialCutsCouple *, G4double logKineticEnergy)

G4bool	TablesAreBuilt () const

G4PhysicsTable *	DEDXTable () const

G4PhysicsTable *	DEDXTableForSubsec () const

G4PhysicsTable *	DEDXunRestrictedTable () const

G4PhysicsTable *	IonisationTable () const

G4PhysicsTable *	IonisationTableForSubsec () const

G4PhysicsTable *	CSDARangeTable () const

G4PhysicsTable *	SecondaryRangeTable () const

G4PhysicsTable *	RangeTableForLoss () const

G4PhysicsTable *	InverseRangeTable () const

G4PhysicsTable *	LambdaTable () const

G4PhysicsTable *	SubLambdaTable () const

const G4Element *	GetCurrentElement () const

void	SetDynamicMassCharge (G4double massratio, G4double charge2ratio)

Public Member Functions inherited from G4VContinuousDiscreteProcess
	G4VContinuousDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)

	G4VContinuousDiscreteProcess (G4VContinuousDiscreteProcess &)

virtual	~G4VContinuousDiscreteProcess ()

G4VContinuousDiscreteProcess &	operator= (const G4VContinuousDiscreteProcess &)=delete

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

virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety, G4GPILSelection *selection)

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)

virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)

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

virtual G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4VParticleChange *	AtRestDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0

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)

virtual G4bool	IsApplicable (const G4ParticleDefinition &)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	PreparePhysicsTable (const G4ParticleDefinition &)

virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

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 void	StartTracking (G4Track *)

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

virtual void	ProcessDescription (std::ostream &outfile) 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	InitialiseEnergyLossProcess (const G4ParticleDefinition , const G4ParticleDefinition )=0

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

virtual G4double	GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override

virtual G4double	GetContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety) override

G4PhysicsVector *	LambdaPhysicsVector (const G4MaterialCutsCouple *, G4double cut)

size_t	CurrentMaterialCutsCoupleIndex () const

void	SelectModel (G4double kinEnergy)

void	SetParticle (const G4ParticleDefinition *p)

void	SetSecondaryParticle (const G4ParticleDefinition *p)

Protected Member Functions inherited from G4VContinuousDiscreteProcess
virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0

virtual G4double	GetContinuousStepLimit (const G4Track &aTrack, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)=0

void	SetGPILSelection (G4GPILSelection selection)

G4GPILSelection	GetGPILSelection () const

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)

void	ClearNumberOfInteractionLengthLeft ()

Protected Attributes
G4ParticleChangeForLoss	fParticleChange

const G4Material *	currentMaterial

const G4MaterialCutsCouple *	currentCouple

size_t	currentCoupleIndex

G4double	preStepLambda

G4double	fRange

G4double	computedRange

G4double	preStepKinEnergy

G4double	preStepLogKinEnergy

G4double	preStepScaledEnergy

G4double	preStepLogScaledEnergy

G4double	preStepRangeEnergy

G4double	mfpKinEnergy

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 81 of file G4VEnergyLossProcess.hh.

Constructor & Destructor Documentation

◆ G4VEnergyLossProcess()

G4VEnergyLossProcess::G4VEnergyLossProcess	(	const G4String &	name = `"EnergyLoss"`,
		G4ProcessType	type = `fElectromagnetic`
	)

Definition at line 86 of file G4VEnergyLossProcess.cc.

                                                              : 
  G4VContinuousDiscreteProcess(name, type),
  secondaryParticle(nullptr),
  nSCoffRegions(0),
  idxSCoffRegions(nullptr),
  nProcesses(0),
  theDEDXTable(nullptr),
  theDEDXSubTable(nullptr),
  theDEDXunRestrictedTable(nullptr),
  theIonisationTable(nullptr),
  theIonisationSubTable(nullptr),
  theRangeTableForLoss(nullptr),
  theCSDARangeTable(nullptr),
  theSecondaryRangeTable(nullptr),
  theInverseRangeTable(nullptr),
  theLambdaTable(nullptr),
  theSubLambdaTable(nullptr),
  baseParticle(nullptr),
  lossFluctuationFlag(true),
  rndmStepFlag(false),
  tablesAreBuilt(false),
  integral(true),
  isIon(false),
  isIonisation(true),
  useSubCutoff(false),
  useDeexcitation(false),
  currentCouple(nullptr),
  mfpKinEnergy(0.0),
  particle(nullptr)
{
  theParameters = G4EmParameters::Instance();
  SetVerboseLevel(1);
 
  // low energy limit
  lowestKinEnergy     = theParameters->LowestElectronEnergy();
  preStepKinEnergy    = 0.0;
  preStepLogKinEnergy = LOG_EKIN_MIN;
  preStepRangeEnergy  = 0.0;
  computedRange       = DBL_MAX;
 
  // Size of tables assuming spline
  minKinEnergy     = 0.1*keV;
  maxKinEnergy     = 100.0*TeV;
  nBins            = 84;
  maxKinEnergyCSDA = 1.0*GeV;
  nBinsCSDA        = 35;
  actMinKinEnergy = actMaxKinEnergy = actBinning = actLinLossLimit 
    = actLossFluc = actIntegral = false;
 
  // default linear loss limit for spline
  linLossLimit = 0.01;
  dRoverRange  = 0.2;
  finalRange   = CLHEP::mm;
 
  // default lambda factor
  lambdaFactor    = 0.8;
  logLambdafactor = G4Log(lambdaFactor);
 
  // cross section biasing
  biasFactor = 1.0;
 
  // particle types
  theElectron   = G4Electron::Electron();
  thePositron   = G4Positron::Positron();
  theGamma      = G4Gamma::Gamma();
  theGenericIon = nullptr;
 
  // run time objects
  pParticleChange = &fParticleChange;
  fParticleChange.SetSecondaryWeightByProcess(true);
  modelManager = new G4EmModelManager();
  safetyHelper = G4TransportationManager::GetTransportationManager()
    ->GetSafetyHelper();
  aGPILSelection = CandidateForSelection;
 
  // initialise model
  lManager = G4LossTableManager::Instance();
  lManager->Register(this);
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  theDensityFactor = bld->GetDensityFactors();
  theDensityIdx = bld->GetCoupleIndexes();
 
  fluctModel = nullptr;
  currentModel = nullptr;
  atomDeexcitation = nullptr;
  subcutProducer = nullptr;
 
  biasManager  = nullptr;
  biasFlag     = false; 
  weightFlag   = false; 
  isMaster     = true;
  lastIdx      = 0;
 
  scTracks.reserve(5);
  secParticles.reserve(5);
 
  theCuts = theSubCuts = nullptr;
  currentMaterial = nullptr;
  currentCoupleIndex  = basedCoupleIndex = 0;
  massRatio = fFactor = reduceFactor = chargeSqRatio = 1.0;
  preStepLambda = preStepScaledEnergy = fRange = logMassRatio = 0.0;
  preStepLogScaledEnergy = LOG_EKIN_MIN;
 
  secID = biasID = subsecID = -1;
}

◆ ~G4VEnergyLossProcess()

G4VEnergyLossProcess::~G4VEnergyLossProcess ( )

virtual

Definition at line 195 of file G4VEnergyLossProcess.cc.

{
  /*
  G4cout << "** G4VEnergyLossProcess::~G4VEnergyLossProcess() for " 
         << GetProcessName() << " isMaster: " << isMaster
         << "  basePart: " << baseParticle 
         << G4endl;
  */
  Clean();
 
  // G4cout << " isIonisation " << isIonisation << "  " 
  //   << theDEDXTable << "  " <<  theIonisationTable << G4endl;
 
  if (isMaster && !baseParticle) {
    if(theDEDXTable) {
 
      //G4cout << " theIonisationTable " << theIonisationTable << G4endl;
      if(theIonisationTable == theDEDXTable) { theIonisationTable = nullptr; }
      //G4cout << " delete theDEDXTable " << theDEDXTable << G4endl;
      theDEDXTable->clearAndDestroy();
      delete theDEDXTable;
      theDEDXTable = nullptr;
      if(theDEDXSubTable) {
        if(theIonisationSubTable == theDEDXSubTable) 
          { theIonisationSubTable = nullptr; }
        theDEDXSubTable->clearAndDestroy();
        delete theDEDXSubTable;
        theDEDXSubTable = nullptr;
      }
    }
    //G4cout << " theIonisationTable " << theIonisationTable << G4endl;
    if(theIonisationTable) {
      //G4cout << " delete theIonisationTable " << theIonisationTable << G4endl;
      theIonisationTable->clearAndDestroy();
      delete theIonisationTable;
      theIonisationTable = nullptr;
    }
    if(theIonisationSubTable) {
      theIonisationSubTable->clearAndDestroy();
      delete theIonisationSubTable;
      theIonisationSubTable = nullptr;
    }
    if(theDEDXunRestrictedTable && isIonisation) {
      theDEDXunRestrictedTable->clearAndDestroy();
      delete theDEDXunRestrictedTable;
      theDEDXunRestrictedTable = nullptr;
    }
    if(theCSDARangeTable && isIonisation) {
      theCSDARangeTable->clearAndDestroy();
      delete theCSDARangeTable;
      theCSDARangeTable = nullptr;
    }
    //G4cout << "delete RangeTable: " << theRangeTableForLoss << G4endl;
    if(theRangeTableForLoss && isIonisation) {
      theRangeTableForLoss->clearAndDestroy();
      delete theRangeTableForLoss;
      theRangeTableForLoss = nullptr;
    }
    //G4cout << "delete InvRangeTable: " << theInverseRangeTable << G4endl;
    if(theInverseRangeTable && isIonisation /*&& !isIon*/) {
      theInverseRangeTable->clearAndDestroy();
      delete theInverseRangeTable;
      theInverseRangeTable = nullptr;
    }
    //G4cout << "delete LambdaTable: " << theLambdaTable << G4endl;
    if(theLambdaTable) {
      theLambdaTable->clearAndDestroy();
      delete theLambdaTable;
      theLambdaTable = nullptr;
    }
    if(theSubLambdaTable) {
      theSubLambdaTable->clearAndDestroy();
      delete theSubLambdaTable;
      theSubLambdaTable = nullptr;
    }
  }
 
  delete modelManager;
  delete biasManager;
  lManager->DeRegister(this);
  //G4cout << "** all removed" << G4endl;
}

Member Function Documentation

◆ ActivateForcedInteraction()

void G4VEnergyLossProcess::ActivateForcedInteraction	(	G4double	length,
		const G4String &	region,
		G4bool	flag = `true`
	)

Definition at line 2236 of file G4VEnergyLossProcess.cc.

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

Referenced by G4EmExtraParameters::DefineRegParamForLoss().

◆ ActivateSecondaryBiasing()

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

Definition at line 2256 of file G4VEnergyLossProcess.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::DefineRegParamForLoss().

◆ ActivateSubCutoff()

void G4VEnergyLossProcess::ActivateSubCutoff	(	G4bool	val,
		const G4Region *	region = `nullptr`
	)

Definition at line 956 of file G4VEnergyLossProcess.cc.

{
  G4RegionStore* regionStore = G4RegionStore::GetInstance();
  const G4Region* reg = r;
  if (!reg) {
    reg = regionStore->GetRegion("DefaultRegionForTheWorld", false);
  }
 
  // the region is in the list
  if (nSCoffRegions > 0) {
    for (G4int i=0; i<nSCoffRegions; ++i) {
      if (reg == scoffRegions[i]) {
        return;
      }
    }
  }
  // new region 
  if(val) {
    scoffRegions.push_back(reg);
    ++nSCoffRegions;
  }
}

Referenced by G4EmExtraParameters::DefineRegParamForLoss().

◆ AddCollaborativeProcess()

void G4VEnergyLossProcess::AddCollaborativeProcess ( G4VEnergyLossProcess * p )

Definition at line 1985 of file G4VEnergyLossProcess.cc.

{
  G4bool add = true;
  if(p->GetProcessName() != "eBrem") { add = false; }
  if(add && nProcesses > 0) {
    for(G4int i=0; i<nProcesses; ++i) {
      if(p == scProcesses[i]) {
        add = false;
        break;
      }
    }
  }
  if(add) {
    scProcesses.push_back(p);
    ++nProcesses;
    if (1 < verboseLevel) { 
      G4cout << "### The process " << p->GetProcessName() 
             << " is added to the list of collaborative processes of "
             << GetProcessName() << G4endl; 
    }
  }
}

◆ AddEmModel()

void G4VEnergyLossProcess::AddEmModel	(	G4int	order,
		G4VEmModel *	p,
		G4VEmFluctuationModel *	fluc = `0`,
		const G4Region *	region = `nullptr`
	)

Definition at line 312 of file G4VEnergyLossProcess.cc.

{
  modelManager->AddEmModel(order, p, fluc, region);
  if(p) { p->SetParticleChange(pParticleChange, fluc); }
}

◆ AlongStepDoIt()

G4VParticleChange * G4VEnergyLossProcess::AlongStepDoIt	(	const G4Track &	track,
		const G4Step &	step
	)

overridevirtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1186 of file G4VEnergyLossProcess.cc.

{
  fParticleChange.InitializeForAlongStep(track);
  // The process has range table - calculate energy loss
  if(!isIonisation || !currentModel->IsActive(preStepScaledEnergy)) {
    return &fParticleChange;
  }
 
  // Get the actual (true) Step length
  G4double length = step.GetStepLength();
  if(length <= 0.0) { return &fParticleChange; }
  G4double eloss  = 0.0;
 
  /*      
  if(-1 < verboseLevel) {
    const G4ParticleDefinition* d = track.GetParticleDefinition();
    G4cout << "AlongStepDoIt for "
           << GetProcessName() << " and particle "
           << d->GetParticleName()
           << "  eScaled(MeV)= " << preStepScaledEnergy/MeV
           << "  range(mm)= " << fRange/mm
           << "  s(mm)= " << length/mm
           << "  rf= " << reduceFactor
           << "  q^2= " << chargeSqRatio
           << " md= " << d->GetPDGMass()
           << "  status= " << track.GetTrackStatus()
           << "  " << track.GetMaterial()->GetName()
           << G4endl;
  }
  */
 
  const G4DynamicParticle* dynParticle = track.GetDynamicParticle();
 
  // define new weight for primary and secondaries
  G4double weight = fParticleChange.GetParentWeight();
  if(weightFlag) {
    weight /= biasFactor;
    fParticleChange.ProposeWeight(weight);
  }
 
  // stopping
  if (length >= fRange || preStepKinEnergy <= lowestKinEnergy) {
    eloss = preStepKinEnergy;
    if (useDeexcitation) {
      atomDeexcitation->AlongStepDeexcitation(scTracks, step, 
                                              eloss, currentCoupleIndex);
      if(scTracks.size() > 0) { FillSecondariesAlongStep(eloss, weight); }
      eloss = std::max(eloss, 0.0);
    }
    fParticleChange.SetProposedKineticEnergy(0.0);
    fParticleChange.ProposeLocalEnergyDeposit(eloss);
    return &fParticleChange;
  }
  //G4cout << theDEDXTable << "  idx= " << basedCoupleIndex 
  // << "  " << GetProcessName() << "  "<< currentMaterial->GetName()<<G4endl;
  //if(particle->GetParticleName() == "e-")G4cout << (*theDEDXTable) <<G4endl;
  // Short step
  eloss  = GetDEDXForScaledEnergy(preStepScaledEnergy, preStepLogScaledEnergy);
  eloss *= length;
 
  //G4cout << "eloss= " << eloss << G4endl;
 
  // Long step
  if(eloss > preStepKinEnergy*linLossLimit) {
 
    G4double x = (fRange - length)/reduceFactor;
    //G4cout << "x= " << x << "  " << theInverseRangeTable << G4endl;
    eloss = preStepKinEnergy - ScaledKinEnergyForLoss(x)/massRatio;
   
    /*
    if(-1 < verboseLevel) 
      G4cout << "Long STEP: rPre(mm)= " 
             << GetScaledRangeForScaledEnergy(preStepScaledEnergy)/mm
             << " rPost(mm)= " << x/mm
             << " ePre(MeV)= " << preStepScaledEnergy/MeV
             << " eloss(MeV)= " << eloss/MeV
             << " eloss0(MeV)= "
             << GetDEDXForScaledEnergy(preStepScaledEnergy)*length/MeV
             << " lim(MeV)= " << preStepKinEnergy*linLossLimit/MeV
             << G4endl;
    */
  }
 
  /*   
  G4double eloss0 = eloss;
  if(-1 < verboseLevel ) {
    G4cout << "Before fluct: eloss(MeV)= " << eloss/MeV
           << " e-eloss= " << preStepKinEnergy-eloss
           << " step(mm)= " << length/mm
           << " range(mm)= " << fRange/mm
           << " fluct= " << lossFluctuationFlag
           << G4endl;
  }
  */
 
  G4double cut  = (*theCuts)[currentCoupleIndex];
  G4double esec = 0.0;
 
  //G4cout << "cut= " << cut << " useSubCut= " << useSubCutoff << G4endl;
 
  // SubCutOff 
  if(useSubCutoff && !subcutProducer) {
    if(idxSCoffRegions[currentCoupleIndex]) {
 
      G4bool yes = false;
      const G4StepPoint* prePoint = step.GetPreStepPoint();
 
      // Check boundary
      if(prePoint->GetStepStatus() == fGeomBoundary) { yes = true; }
 
      // Check PrePoint
      else {
        G4double preSafety  = prePoint->GetSafety();
        G4double rcut = 
          currentCouple->GetProductionCuts()->GetProductionCut(1);
 
        // recompute presafety
        if(preSafety < rcut) {
          preSafety = safetyHelper->ComputeSafety(prePoint->GetPosition(),
                                                  rcut);
        }
 
        if(preSafety < rcut) { yes = true; }
 
        // Check PostPoint
        else {
          G4double postSafety = preSafety - length; 
          if(postSafety < rcut) {
            postSafety = safetyHelper->ComputeSafety(
              step.GetPostStepPoint()->GetPosition(), rcut);
            if(postSafety < rcut) { yes = true; }
          }
        }
      }
  
      // Decided to start subcut sampling
      if(yes) {
 
        cut = (*theSubCuts)[currentCoupleIndex];
         eloss -= GetSubDEDXForScaledEnergy(preStepScaledEnergy)*length;
        esec = SampleSubCutSecondaries(scTracks, step, 
                                       currentModel,currentCoupleIndex);
        // add bremsstrahlung sampling
        /*
        if(nProcesses > 0) {
          for(G4int i=0; i<nProcesses; ++i) {
            (scProcesses[i])->SampleSubCutSecondaries(
                scTracks, step, (scProcesses[i])->
                SelectModelForMaterial(preStepKinEnergy, currentCoupleIndex),
                currentCoupleIndex);
          }
        } 
        */
      }   
    }
  }
 
  // Corrections, which cannot be tabulated
  if(isIon) {
    G4double eadd = 0.0;
    G4double eloss_before = eloss;
    currentModel->CorrectionsAlongStep(currentCouple, dynParticle, 
                                       eloss, eadd, length);
    if(eloss < 0.0) { eloss = 0.5*eloss_before; }
  }
 
  // Sample fluctuations
  if (lossFluctuationFlag) {
    G4VEmFluctuationModel* fluc = currentModel->GetModelOfFluctuations();
    if(eloss + esec < preStepKinEnergy) {
 
      G4double tmax = 
        std::min(currentModel->MaxSecondaryKinEnergy(dynParticle),cut);
      eloss = fluc->SampleFluctuations(currentCouple,dynParticle,
                                       tmax,length,eloss);
      /*                                  
      if(-1 < verboseLevel) 
      G4cout << "After fluct: eloss(MeV)= " << eloss/MeV
             << " fluc= " << (eloss-eloss0)/MeV
             << " ChargeSqRatio= " << chargeSqRatio
             << " massRatio= " << massRatio
             << " tmax= " << tmax
             << G4endl;
      */
    }
  }
 
  // deexcitation
  if (useDeexcitation) {
    G4double esecfluo = preStepKinEnergy - esec;
    G4double de = esecfluo;
    //G4double eloss0 = eloss;
    /*
    G4cout << "### 1: E(keV)= " << preStepKinEnergy/keV
           << " Efluomax(keV)= " << de/keV
           << " Eloss(keV)= " << eloss/keV << G4endl; 
    */
    atomDeexcitation->AlongStepDeexcitation(scTracks, step, 
                                            de, currentCoupleIndex);
 
    // sum of de-excitation energies
    esecfluo -= de;
 
    // subtracted from energy loss
    if(eloss >= esecfluo) {
      esec  += esecfluo;
      eloss -= esecfluo;
    } else {
      esec += esecfluo;
      eloss = 0.0; 
    } 
    /*    
    if(esecfluo > 0.0) {
      G4cout << "### 2: E(keV)= " << preStepKinEnergy/keV
             << " Esec(keV)= " << esec/keV
             << " Esecf(kV)= " << esecfluo/keV
             << " Eloss0(kV)= " << eloss0/keV
             << " Eloss(keV)= " << eloss/keV 
             << G4endl; 
    } 
    */   
  }
  if(subcutProducer && idxSCoffRegions[currentCoupleIndex]) {
    subcutProducer->SampleSecondaries(step, scTracks, eloss, cut);
  }
  if(scTracks.size() > 0) { FillSecondariesAlongStep(eloss, weight); }
 
  // Energy balance
  G4double finalT = preStepKinEnergy - eloss - esec;
  if (finalT <= lowestKinEnergy) {
    eloss += finalT;
    finalT = 0.0;
  } else if(isIon) {
    fParticleChange.SetProposedCharge(
      currentModel->GetParticleCharge(track.GetParticleDefinition(),
                                      currentMaterial,finalT));
  }
 
  eloss = std::max(eloss, 0.0);
 
  fParticleChange.SetProposedKineticEnergy(finalT);
  fParticleChange.ProposeLocalEnergyDeposit(eloss);
  /*
  if(-1 < verboseLevel) {
    G4double del = finalT + eloss + esec - preStepKinEnergy;
    G4cout << "Final value eloss(MeV)= " << eloss/MeV
           << " preStepKinEnergy= " << preStepKinEnergy
           << " postStepKinEnergy= " << finalT
           << " de(keV)= " << del/keV
           << " lossFlag= " << lossFluctuationFlag
           << "  status= " << track.GetTrackStatus()
           << G4endl;
  }
  */
  return &fParticleChange;
}

◆ AlongStepGetPhysicalInteractionLength()

G4double G4VEnergyLossProcess::AlongStepGetPhysicalInteractionLength	(	const G4Track &	,
		G4double	previousStepSize,
		G4double	currentMinimumStep,
		G4double &	currentSafety,
		G4GPILSelection *	selection
	)

overridevirtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1023 of file G4VEnergyLossProcess.cc.

{
  G4double x = DBL_MAX;
  *selection = aGPILSelection;
  if(isIonisation && currentModel->IsActive(preStepScaledEnergy)) {
    fRange = reduceFactor*GetScaledRangeForScaledEnergy(preStepScaledEnergy,
                                                        preStepLogScaledEnergy);
    G4double finR = (rndmStepFlag) ? std::min(finalRange,
      currentCouple->GetProductionCuts()->GetProductionCut(1)) : finalRange;
    x = (fRange > finR) ? 
      fRange*dRoverRange + finR*(1.0-dRoverRange)*(2.0-finR/fRange) : fRange; 
   // if(particle->GetPDGMass() > 0.9*GeV)
    /*
    G4cout<<GetProcessName()<<": e= "<<preStepKinEnergy
          <<" range= "<<fRange << " idx= " << basedCoupleIndex
              << " finR= " << finR
          << " limit= " << x <<G4endl;
    G4cout << "massRatio= " << massRatio << " Q^2= " << chargeSqRatio 
           << " finR= " << finR << " dRoverRange= " << dRoverRange 
           << " finalRange= " << finalRange << G4endl;
    */
  }
  //G4cout<<GetProcessName()<<": e= "<<preStepKinEnergy 
  //<<" stepLimit= "<<x<<G4endl;
  return x;
}

Referenced by ContinuousStepLimit().

◆ BaseParticle()

const G4ParticleDefinition * G4VEnergyLossProcess::BaseParticle ( ) const

inline

Definition at line 980 of file G4VEnergyLossProcess.hh.

{
  return baseParticle;
}

Referenced by G4LossTableManager::BuildPhysicsTable(), G4LossTableManager::LocalPhysicsTables(), and G4LossTableManager::RegisterExtraParticle().

◆ BuildDEDXTable()

G4PhysicsTable * G4VEnergyLossProcess::BuildDEDXTable ( G4EmTableType tType = fRestricted )

Definition at line 713 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel ) {
    G4cout << "G4VEnergyLossProcess::BuildDEDXTable() of type " << tType
           << " for " << GetProcessName()
           << " and particle " << particle->GetParticleName()
           << G4endl;
  }
  G4PhysicsTable* table = nullptr;
  G4double emax = maxKinEnergy;
  G4int bin = nBins;
 
  if(fTotal == tType) {
    emax  = maxKinEnergyCSDA;
    bin   = nBinsCSDA;
    table = theDEDXunRestrictedTable;
  } else if(fRestricted == tType) {
    table = theDEDXTable;
  } else if(fSubRestricted == tType) {    
    table = theDEDXSubTable;
  } else {
    G4cout << "G4VEnergyLossProcess::BuildDEDXTable WARNING: wrong type "
           << tType << G4endl;
  }
 
  // Access to materials
  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();
 
  if(1 < verboseLevel) {
    G4cout << numOfCouples << " materials"
           << " minKinEnergy= " << minKinEnergy
           << " maxKinEnergy= " << emax
           << " nbin= " << bin
           << " EmTableType= " << tType
           << " table= " << table << "  " << this 
           << G4endl;
  }
  if(!table) { return table; }
 
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  G4bool splineFlag = theParameters->Spline();
  G4PhysicsLogVector* aVector = nullptr;
  G4PhysicsLogVector* bVector = nullptr;
 
  for(size_t i=0; i<numOfCouples; ++i) {
 
    if(1 < verboseLevel) {
      G4cout << "G4VEnergyLossProcess::BuildDEDXVector Idx= " << i 
             << "  flagTable=  " << table->GetFlag(i) 
             << " Flag= " << bld->GetFlag(i) << G4endl;
    }
    if(bld->GetFlag(i)) {
 
      // create physics vector and fill it
      const G4MaterialCutsCouple* couple = 
        theCoupleTable->GetMaterialCutsCouple(i);
      if((*table)[i]) { delete (*table)[i]; }
      if(bVector) {
        aVector = new G4PhysicsLogVector(*bVector);
      } else {
        bVector = new G4PhysicsLogVector(minKinEnergy, emax, bin);
        aVector = bVector;
      }
      aVector->SetSpline(splineFlag);
 
      modelManager->FillDEDXVector(aVector, couple, tType);
      if(splineFlag) { aVector->FillSecondDerivatives(); }
 
      // Insert vector for this material into the table
      G4PhysicsTableHelper::SetPhysicsVector(table, i, aVector);
    }
  }
 
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::BuildDEDXTable(): table is built for "
           << particle->GetParticleName()
           << " and process " << GetProcessName()
           << G4endl;
    if(2 < verboseLevel) G4cout << (*table) << G4endl;
  }
 
  return table;
}

◆ BuildLambdaTable()

G4PhysicsTable * G4VEnergyLossProcess::BuildLambdaTable ( G4EmTableType tType = fRestricted )

Definition at line 801 of file G4VEnergyLossProcess.cc.

{
  G4PhysicsTable* table = nullptr;
 
  if(fRestricted == tType) {
    table = theLambdaTable;
  } else if(fSubRestricted == tType) {    
    table = theSubLambdaTable;
  } else {
    G4cout << "G4VEnergyLossProcess::BuildLambdaTable WARNING: wrong type "
           << tType << G4endl;
  }
 
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::BuildLambdaTable() of type "
           << tType << " for process "
           << GetProcessName() << " and particle "
           << particle->GetParticleName()
           << " EmTableType= " << tType
           << " table= " << table
           << G4endl;
  }
  if(!table) {return table;}
 
  // Access to materials
  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();
 
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  theDensityFactor = bld->GetDensityFactors();
  theDensityIdx = bld->GetCoupleIndexes();
 
  G4bool splineFlag = theParameters->Spline();
  G4PhysicsLogVector* aVector = nullptr;
  G4double scale = G4Log(maxKinEnergy/minKinEnergy);
 
  for(size_t i=0; i<numOfCouples; ++i) {
 
    if (bld->GetFlag(i)) {
 
      // create physics vector and fill it
      const G4MaterialCutsCouple* couple = 
        theCoupleTable->GetMaterialCutsCouple(i);
      delete (*table)[i];
 
      G4bool startNull = true;
      G4double emin = 
        MinPrimaryEnergy(particle,couple->GetMaterial(),(*theCuts)[i]);
      if(minKinEnergy > emin) { 
        emin = minKinEnergy; 
        startNull = false;
      }
 
      G4double emax = maxKinEnergy;
      if(emax <= emin) { emax = 2*emin; }
      G4int bin = G4lrint(nBins*G4Log(emax/emin)/scale);
      bin = std::max(bin, 3);
      aVector = new G4PhysicsLogVector(emin, emax, bin);
      aVector->SetSpline(splineFlag);
 
      modelManager->FillLambdaVector(aVector, couple, startNull, tType);
      if(splineFlag) { aVector->FillSecondDerivatives(); }
 
      // Insert vector for this material into the table
      G4PhysicsTableHelper::SetPhysicsVector(table, i, aVector);
    }
  }
 
  if(1 < verboseLevel) {
    G4cout << "Lambda table is built for "
           << particle->GetParticleName()
           << G4endl;
  }
 
  return table;
}

◆ BuildPhysicsTable()

void G4VEnergyLossProcess::BuildPhysicsTable ( const G4ParticleDefinition & part )

overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4ePolarizedIonisation.

Definition at line 613 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {  
    G4cout << "### G4VEnergyLossProcess::BuildPhysicsTable() for "
           << GetProcessName()
           << " and particle " << part.GetParticleName()
           << "; local: " << particle->GetParticleName();
    if(baseParticle) { 
      G4cout << "; base: " << baseParticle->GetParticleName(); 
    }
    G4cout << " TablesAreBuilt= " << tablesAreBuilt
           << " isIon= " << isIon << "  " << this << G4endl;
  }
 
  if(&part == particle) {
 
    if(isMaster) {
      lManager->BuildPhysicsTable(particle, this);
 
    } else {
 
      const G4VEnergyLossProcess* masterProcess = 
        static_cast<const G4VEnergyLossProcess*>(GetMasterProcess());
 
      // copy table pointers from master thread
      SetDEDXTable(masterProcess->DEDXTable(),fRestricted);
      SetDEDXTable(masterProcess->DEDXTableForSubsec(),fSubRestricted);
      SetDEDXTable(masterProcess->DEDXunRestrictedTable(),fTotal);
      SetDEDXTable(masterProcess->IonisationTable(),fIsIonisation);
      SetDEDXTable(masterProcess->IonisationTableForSubsec(),fIsSubIonisation);
      SetRangeTableForLoss(masterProcess->RangeTableForLoss());
      SetCSDARangeTable(masterProcess->CSDARangeTable());
      SetSecondaryRangeTable(masterProcess->SecondaryRangeTable());
      SetInverseRangeTable(masterProcess->InverseRangeTable());
      SetLambdaTable(masterProcess->LambdaTable());
      SetSubLambdaTable(masterProcess->SubLambdaTable());
      isIonisation = masterProcess->IsIonisationProcess();
 
      tablesAreBuilt = true;  
      // local initialisation of models
      G4bool printing = true;
      G4int numberOfModels = modelManager->NumberOfModels();
      for(G4int i=0; i<numberOfModels; ++i) {
        G4VEmModel* mod = GetModelByIndex(i, printing);
        G4VEmModel* mod0= masterProcess->GetModelByIndex(i,printing);
        mod->InitialiseLocal(particle, mod0);
      }
 
      lManager->LocalPhysicsTables(particle, this);
    }
   
    // needs to be done only once
    safetyHelper->InitialiseHelper();
  }
  // explicitly defined printout by particle name
  G4String num = part.GetParticleName();
  if(1 < verboseLevel || 
     (0 < verboseLevel && (num == "e-" || 
                           num == "e+"    || num == "mu+" || 
                           num == "mu-"   || num == "proton"|| 
                           num == "pi+"   || num == "pi-" || 
                           num == "kaon+" || num == "kaon-" || 
                           num == "alpha" || num == "anti_proton" || 
                           num == "GenericIon"|| num == "alpha++" ||
                           num == "alpha+" )))
    { 
      StreamInfo(G4cout, part); 
    }
 
  // Added tracking cut to avoid tracking artifacts
  // identify deexcitation flag
  if(isIonisation) { 
    atomDeexcitation = lManager->AtomDeexcitation();
    if(nSCoffRegions > 0) { subcutProducer = lManager->SubCutProducer(); }
    if(atomDeexcitation) { 
      if(atomDeexcitation->IsPIXEActive()) { useDeexcitation = true; } 
    }
  }
  /*  
  G4cout << "** G4VEnergyLossProcess::BuildPhysicsTable() for " 
         << GetProcessName() << " and " << particle->GetParticleName()
         << " isMaster: " << isMaster << " isIonisation: " 
         << isIonisation << G4endl;
  G4cout << " theDEDX: " << theDEDXTable 
         << " theRange: " << theRangeTableForLoss
         << " theInverse: " << theInverseRangeTable
         << " theLambda: " << theLambdaTable << G4endl;
  */
  //if(1 < verboseLevel || verb) {
  if(1 < verboseLevel) {
    G4cout << "### G4VEnergyLossProcess::BuildPhysicsTable() done for "
           << GetProcessName()
           << " and particle " << part.GetParticleName();
    if(isIonisation) { G4cout << "  isIonisation  flag = 1"; }
    G4cout << G4endl;
  }
}

Referenced by G4ePolarizedIonisation::BuildPhysicsTable().

◆ ContinuousStepLimit()

G4double G4VEnergyLossProcess::ContinuousStepLimit	(	const G4Track &	track,
		G4double	previousStepSize,
		G4double	currentMinimumStep,
		G4double &	currentSafety
	)

Definition at line 1942 of file G4VEnergyLossProcess.cc.

{
  G4GPILSelection sel;
  return AlongStepGetPhysicalInteractionLength(track, x, y, z, &sel);
}

◆ CrossSectionBiasingFactor()

G4double G4VEnergyLossProcess::CrossSectionBiasingFactor ( ) const

inline

Definition at line 1046 of file G4VEnergyLossProcess.hh.

{
  return biasFactor;
}

◆ CrossSectionPerVolume() [1/2]

G4double G4VEnergyLossProcess::CrossSectionPerVolume	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

◆ CrossSectionPerVolume() [2/2]

G4double G4VEnergyLossProcess::CrossSectionPerVolume	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple,
		G4double	logKineticEnergy
	)

Definition at line 1908 of file G4VEnergyLossProcess.cc.

{
  // Cross section per volume is calculated
  DefineMaterial(couple);
  G4double cross = 0.0;
  if (theLambdaTable) {
    cross = GetLambdaForScaledEnergy(kineticEnergy * massRatio,
                                     logKineticEnergy + logMassRatio);
  } else {
    SelectModel(kineticEnergy*massRatio);
    cross = biasFactor*(*theDensityFactor)[currentCoupleIndex]
      *(currentModel->CrossSectionPerVolume(currentMaterial,
                                            particle, kineticEnergy,
                                            (*theCuts)[currentCoupleIndex]));
  }
  return std::max(cross, 0.0);
}

◆ CSDARangeTable()

G4PhysicsTable * G4VEnergyLossProcess::CSDARangeTable ( ) const

inline

Definition at line 1095 of file G4VEnergyLossProcess.hh.

{
  return theCSDARangeTable;
}

Referenced by BuildPhysicsTable().

◆ CurrentMaterialCutsCoupleIndex()

size_t G4VEnergyLossProcess::CurrentMaterialCutsCoupleIndex ( ) const

inlineprotected

Definition at line 595 of file G4VEnergyLossProcess.hh.

{
  return currentCoupleIndex;
}

◆ DEDXTable()

G4PhysicsTable * G4VEnergyLossProcess::DEDXTable ( ) const

inline

Definition at line 1060 of file G4VEnergyLossProcess.hh.

{
  return theDEDXTable;
}

Referenced by BuildPhysicsTable(), G4LossTableManager::LocalPhysicsTables(), and G4EmCalculator::PrintDEDXTable().

◆ DEDXTableForSubsec()

G4PhysicsTable * G4VEnergyLossProcess::DEDXTableForSubsec ( ) const

inline

Definition at line 1067 of file G4VEnergyLossProcess.hh.

{
  return theDEDXSubTable;
}

Referenced by BuildPhysicsTable().

◆ DEDXunRestrictedTable()

G4PhysicsTable * G4VEnergyLossProcess::DEDXunRestrictedTable ( ) const

inline

Definition at line 1074 of file G4VEnergyLossProcess.hh.

{
  return theDEDXunRestrictedTable;
}

Referenced by BuildPhysicsTable().

◆ EmModel()

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

Definition at line 338 of file G4VEnergyLossProcess.cc.

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

◆ FluctModel()

G4VEmFluctuationModel * G4VEnergyLossProcess::FluctModel ( )

inline

Definition at line 943 of file G4VEnergyLossProcess.hh.

{
  return fluctModel;
}

Referenced by G4MuIonisation::InitialiseEnergyLossProcess(), G4ePolarizedIonisation::InitialiseEnergyLossProcess(), G4alphaIonisation::InitialiseEnergyLossProcess(), G4eIonisation::InitialiseEnergyLossProcess(), G4hIonisation::InitialiseEnergyLossProcess(), and G4ionIonisation::InitialiseEnergyLossProcess().

◆ GetContinuousStepLimit()

G4double G4VEnergyLossProcess::GetContinuousStepLimit	(	const G4Track &	track,
		G4double	previousStepSize,
		G4double	currentMinimumStep,
		G4double &	currentSafety
	)

overrideprotectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 1964 of file G4VEnergyLossProcess.cc.

{
  return DBL_MAX;
}

◆ GetCSDARange()

G4double G4VEnergyLossProcess::GetCSDARange	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 871 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return (theCSDARangeTable) ? 
    GetLimitScaledRangeForScaledEnergy(kineticEnergy*massRatio)*reduceFactor
    : DBL_MAX;
}

Referenced by G4LossTableManager::GetCSDARange().

◆ GetCurrentElement()

const G4Element * G4VEnergyLossProcess::GetCurrentElement ( ) const

Definition at line 2209 of file G4VEnergyLossProcess.cc.

{
  const G4Element* elm = nullptr;
  if(currentModel) { elm = currentModel->GetCurrentElement(); }
  return elm;
}

◆ GetDEDX() [1/2]

G4double G4VEnergyLossProcess::GetDEDX	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 809 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return GetDEDXForScaledEnergy(kinEnergy*massRatio);
}

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4LossTableManager::GetDEDX(), and G4VMscModel::GetDEDX().

◆ GetDEDX() [2/2]

G4double G4VEnergyLossProcess::GetDEDX	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple,
		G4double	logKineticEnergy
	)

inline

Definition at line 817 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return GetDEDXForScaledEnergy(kinEnergy*massRatio, logKinEnergy+logMassRatio);
}

◆ GetDEDXDispersion()

G4double G4VEnergyLossProcess::GetDEDXDispersion	(	const G4MaterialCutsCouple *	couple,
		const G4DynamicParticle *	dp,
		G4double	length
	)

Definition at line 1889 of file G4VEnergyLossProcess.cc.

{
  DefineMaterial(couple);
  G4double ekin = dp->GetKineticEnergy();
  SelectModel(ekin*massRatio);
  G4double tmax = currentModel->MaxSecondaryKinEnergy(dp);
  tmax = std::min(tmax,(*theCuts)[currentCoupleIndex]);
  G4double d = 0.0;
  G4VEmFluctuationModel* fm = currentModel->GetModelOfFluctuations();
  if(fm) { d = fm->Dispersion(currentMaterial,dp,tmax,length); }
  return d;
}

Referenced by G4LossTableManager::GetDEDXDispersion().

◆ GetDEDXForSubsec()

G4double G4VEnergyLossProcess::GetDEDXForSubsec	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 828 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return GetSubDEDXForScaledEnergy(kineticEnergy*massRatio);
}

Referenced by G4LossTableManager::GetSubDEDX().

◆ GetKineticEnergy()

G4double G4VEnergyLossProcess::GetKineticEnergy	(	G4double	range,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 906 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return ScaledKinEnergyForLoss(range/reduceFactor)/massRatio;
}

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4LossTableManager::GetEnergy(), and G4VMscModel::GetEnergy().

◆ GetLambda() [1/2]

G4double G4VEnergyLossProcess::GetLambda	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 916 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return theLambdaTable ? GetLambdaForScaledEnergy(kinEnergy*massRatio) : 0.0;
}

◆ GetLambda() [2/2]

G4double G4VEnergyLossProcess::GetLambda	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple,
		G4double	logKineticEnergy
	)

inline

Definition at line 924 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  return theLambdaTable
    ?  GetLambdaForScaledEnergy(kinEnergy*massRatio, logKinEnergy+logMassRatio)
    :  0.0;
}

◆ GetMeanFreePath()

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

overrideprotectedvirtual

Implements G4VContinuousDiscreteProcess.

Reimplemented in G4ePolarizedIonisation.

Definition at line 1952 of file G4VEnergyLossProcess.cc.

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

Referenced by G4ePolarizedIonisation::GetMeanFreePath().

◆ GetModelByIndex()

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

Definition at line 345 of file G4VEnergyLossProcess.cc.

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

Referenced by BuildPhysicsTable().

◆ GetRange() [1/2]

G4double G4VEnergyLossProcess::GetRange	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 838 of file G4VEnergyLossProcess.hh.

{
  G4double x = fRange;
  DefineMaterial(couple);
  if(theCSDARangeTable) {
    x  = reduceFactor * GetLimitScaledRangeForScaledEnergy(kinEnergy*massRatio);
  } else if(theRangeTableForLoss) {
    x  = reduceFactor * GetScaledRangeForScaledEnergy(kinEnergy*massRatio);
  }
  return x;
}

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), G4ContinuousGainOfEnergy::GetContinuousStepLimit(), and G4LossTableManager::GetRange().

◆ GetRange() [2/2]

G4double G4VEnergyLossProcess::GetRange	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple,
		G4double	logKineticEnergy
	)

inline

Definition at line 852 of file G4VEnergyLossProcess.hh.

{
  G4double x = fRange;
  DefineMaterial(couple);
  if(theCSDARangeTable) {
    x  = reduceFactor * GetLimitScaledRangeForScaledEnergy(kinEnergy*massRatio,
                                                     logKinEnergy+logMassRatio);
  } else if(theRangeTableForLoss) {
    x  = reduceFactor * GetScaledRangeForScaledEnergy(kinEnergy*massRatio,
                                                     logKinEnergy+logMassRatio);
  }
  return x;
}

◆ GetRangeForLoss() [1/2]

G4double G4VEnergyLossProcess::GetRangeForLoss	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple
	)

inline

Definition at line 883 of file G4VEnergyLossProcess.hh.

{
  //  G4cout << "GetRangeForLoss: Range from " << GetProcessName() << G4endl;
  DefineMaterial(couple);
  return reduceFactor * GetScaledRangeForScaledEnergy(kinEnergy*massRatio);
}

Referenced by G4VMscModel::GetRange(), and G4LossTableManager::GetRangeFromRestricteDEDX().

◆ GetRangeForLoss() [2/2]

G4double G4VEnergyLossProcess::GetRangeForLoss	(	G4double	kineticEnergy,
		const G4MaterialCutsCouple *	couple,
		G4double	logKineticEnergy
	)

inline

Definition at line 893 of file G4VEnergyLossProcess.hh.

{
  //  G4cout << "GetRangeForLoss: Range from " << GetProcessName() << G4endl;
  DefineMaterial(couple);
  return reduceFactor * GetScaledRangeForScaledEnergy(kinEnergy*massRatio,
                                                     logKinEnergy+logMassRatio);
}

◆ InitialiseEnergyLossProcess()

virtual void G4VEnergyLossProcess::InitialiseEnergyLossProcess	(	const G4ParticleDefinition *	,
		const G4ParticleDefinition *
	)

protectedpure virtual

Implemented in G4hBremsstrahlung, G4hhIonisation, G4hPairProduction, G4mplIonisation, G4ePairProduction, G4MuBremsstrahlung, G4MuIonisation, G4MuPairProduction, G4ePolarizedBremsstrahlung, G4ePolarizedIonisation, G4alphaIonisation, G4eBremsstrahlung, G4eIonisation, G4hIonisation, and G4ionIonisation.

Referenced by PreparePhysicsTable().

◆ InverseRangeTable()

G4PhysicsTable * G4VEnergyLossProcess::InverseRangeTable ( ) const

inline

Definition at line 1116 of file G4VEnergyLossProcess.hh.

{
  return theInverseRangeTable;
}

Referenced by BuildPhysicsTable(), G4LossTableManager::LocalPhysicsTables(), and G4EmCalculator::PrintInverseRangeTable().

◆ IonisationTable()

G4PhysicsTable * G4VEnergyLossProcess::IonisationTable ( ) const

inline

Definition at line 1081 of file G4VEnergyLossProcess.hh.

{
  return theIonisationTable;
}

Referenced by BuildPhysicsTable().

◆ IonisationTableForSubsec()

G4PhysicsTable * G4VEnergyLossProcess::IonisationTableForSubsec ( ) const

inline

Definition at line 1088 of file G4VEnergyLossProcess.hh.

{
  return theIonisationSubTable;
}

Referenced by BuildPhysicsTable().

◆ IsApplicable()

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

overridepure virtual

Reimplemented from G4VProcess.

Implemented in G4alphaIonisation, G4eBremsstrahlung, G4eIonisation, G4ionIonisation, G4hBremsstrahlung, G4hhIonisation, G4hPairProduction, G4mplIonisation, G4ePairProduction, G4MuBremsstrahlung, G4MuIonisation, G4MuPairProduction, G4ePolarizedIonisation, and G4hIonisation.

◆ IsIntegral()

G4bool G4VEnergyLossProcess::IsIntegral ( ) const

inline

Definition at line 1011 of file G4VEnergyLossProcess.hh.

{
  return integral;
}

◆ IsIonisationProcess()

G4bool G4VEnergyLossProcess::IsIonisationProcess ( ) const

inline

Definition at line 1018 of file G4VEnergyLossProcess.hh.

{
  return isIonisation;
}

Referenced by BuildPhysicsTable(), G4LossTableManager::BuildPhysicsTable(), and G4LossTableManager::LocalPhysicsTables().

◆ LambdaPhysicsVector()

G4PhysicsVector * G4VEnergyLossProcess::LambdaPhysicsVector	(	const G4MaterialCutsCouple *	,
		G4double	cut
	)

protected

Definition at line 1974 of file G4VEnergyLossProcess.cc.

{
  G4PhysicsVector* v = 
    new G4PhysicsLogVector(minKinEnergy, maxKinEnergy, nBins);
  v->SetSpline(theParameters->Spline());
  return v;
}

◆ LambdaTable()

G4PhysicsTable * G4VEnergyLossProcess::LambdaTable ( ) const

inline

Definition at line 1123 of file G4VEnergyLossProcess.hh.

{
  return theLambdaTable;
}

Referenced by BuildPhysicsTable().

◆ MaxKinEnergy()

G4double G4VEnergyLossProcess::MaxKinEnergy ( ) const

inline

Definition at line 1039 of file G4VEnergyLossProcess.hh.

{
  return maxKinEnergy;
}

◆ MeanFreePath()

G4double G4VEnergyLossProcess::MeanFreePath ( const G4Track & track )

Definition at line 1930 of file G4VEnergyLossProcess.cc.

{
  DefineMaterial(track.GetMaterialCutsCouple());
  const G4double kinEnergy    = track.GetKineticEnergy();
  const G4double logKinEnergy = track.GetDynamicParticle()->GetLogKineticEnergy();
  const G4double cs = GetLambdaForScaledEnergy(kinEnergy * massRatio, 
                                               logKinEnergy + logMassRatio);
  return (0.0 < cs) ? 1.0/cs : DBL_MAX;
}

Referenced by GetMeanFreePath().

◆ MinKinEnergy()

G4double G4VEnergyLossProcess::MinKinEnergy ( ) const

inline

Definition at line 1032 of file G4VEnergyLossProcess.hh.

{
  return minKinEnergy;
}

◆ MinPrimaryEnergy()

G4double G4VEnergyLossProcess::MinPrimaryEnergy	(	const G4ParticleDefinition *	,
		const G4Material *	,
		G4double	cut
	)

protectedvirtual

Reimplemented in G4eIonisation, G4ePolarizedIonisation, G4mplIonisation, G4alphaIonisation, G4hIonisation, G4ionIonisation, G4hhIonisation, G4ePairProduction, G4MuBremsstrahlung, G4MuIonisation, and G4MuPairProduction.

Definition at line 303 of file G4VEnergyLossProcess.cc.

{
  return cut;
}

Referenced by BuildLambdaTable().

◆ NumberOfModels()

G4int G4VEnergyLossProcess::NumberOfModels ( ) const

Definition at line 352 of file G4VEnergyLossProcess.cc.

{
  return modelManager->NumberOfModels();
}

◆ NumberOfSubCutoffRegions()

G4int G4VEnergyLossProcess::NumberOfSubCutoffRegions ( ) const

inline

Definition at line 1025 of file G4VEnergyLossProcess.hh.

{
  return nSCoffRegions;
}

◆ Particle()

const G4ParticleDefinition * G4VEnergyLossProcess::Particle ( ) const

inline

Definition at line 973 of file G4VEnergyLossProcess.hh.

{
  return particle;
}

Referenced by G4LossTableManager::BuildPhysicsTable(), and G4LossTableManager::LocalPhysicsTables().

◆ PostStepDoIt()

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

overridevirtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1557 of file G4VEnergyLossProcess.cc.

{
  // In all cases clear number of interaction lengths
  theNumberOfInteractionLengthLeft = -1.0;
  mfpKinEnergy = currentInteractionLength = DBL_MAX; 
 
  fParticleChange.InitializeForPostStep(track);
  G4double finalT = track.GetKineticEnergy();
  if(finalT <= lowestKinEnergy) { return &fParticleChange; }
 
  G4double postStepScaledEnergy = finalT*massRatio;
  SelectModel(postStepScaledEnergy);
 
  if(!currentModel->IsActive(postStepScaledEnergy)) { 
    return &fParticleChange; 
  }
  /*
  if(-1 < verboseLevel) {
    G4cout << GetProcessName()
           << "::PostStepDoIt: E(MeV)= " << finalT/MeV
           << G4endl;
  }
  */
 
  // forced process - should happen only once per track
  if(biasFlag) {
    if(biasManager->ForcedInteractionRegion(currentCoupleIndex)) {
      biasFlag = false;
    }
  }
 
  const G4DynamicParticle* dp = track.GetDynamicParticle();
  const G4double logFinalT    = dp->GetLogKineticEnergy();
  // postStepLogScaledEnergy = logFinalT + logMassRatio;
 
  // Integral approach
  if (integral) {
    const G4double lx = GetLambdaForScaledEnergy(postStepScaledEnergy,
                                                 logFinalT + logMassRatio);
    /*
    if(preStepLambda<lx && 1 < verboseLevel) {
      G4cout << "WARNING: for " << particle->GetParticleName()
             << " and " << GetProcessName()
             << " E(MeV)= " << finalT/MeV
             << " preLambda= " << preStepLambda 
             << " < " << lx << " (postLambda) "
             << G4endl;
    }
    */
    if(lx <= 0.0 || preStepLambda*G4UniformRand() > lx) {
      return &fParticleChange;
    }
  }
 
  SelectModel(postStepScaledEnergy);
 
  // define new weight for primary and secondaries
  G4double weight = fParticleChange.GetParentWeight();
  if(weightFlag) {
    weight /= biasFactor;
    fParticleChange.ProposeWeight(weight);
  }
 
  G4double tcut = (*theCuts)[currentCoupleIndex];
 
  // sample secondaries
  secParticles.clear();
  //G4cout<< "@@@ Eprimary= "<<dynParticle->GetKineticEnergy()/MeV
  //        << " cut= " << tcut/MeV << G4endl;
  currentModel->SampleSecondaries(&secParticles, currentCouple, dp, tcut);
 
  G4int num0 = secParticles.size();
 
  // bremsstrahlung splitting or Russian roulette  
  if(biasManager) {
    if(biasManager->SecondaryBiasingRegion(currentCoupleIndex)) {
      G4double eloss = 0.0;
      weight *= biasManager->ApplySecondaryBiasing(
                                      secParticles,
                                      track, currentModel, 
                                      &fParticleChange, eloss,
                                      currentCoupleIndex, tcut, 
                                      step.GetPostStepPoint()->GetSafety());
      if(eloss > 0.0) {
        eloss += fParticleChange.GetLocalEnergyDeposit();
        fParticleChange.ProposeLocalEnergyDeposit(eloss);
      }
    }
  }
 
  // save secondaries
  G4int num = secParticles.size();
  if(num > 0) {
 
    fParticleChange.SetNumberOfSecondaries(num);
    G4double time = track.GetGlobalTime();
 
    for (G4int i=0; i<num; ++i) {
      if(secParticles[i]) {
        G4Track* t = new G4Track(secParticles[i], time, track.GetPosition());
        t->SetTouchableHandle(track.GetTouchableHandle());
        if (biasManager) {
          t->SetWeight(weight * biasManager->GetWeight(i));
        } else {
          t->SetWeight(weight);
        }
        if(i < num0) { t->SetCreatorModelIndex(secID); }
        else         { t->SetCreatorModelIndex(biasID); }
 
        //G4cout << "Secondary(post step) has weight " << t->GetWeight() 
        //       << ", kenergy " << t->GetKineticEnergy()/MeV << " MeV" 
        //       << " time= " << time/ns << " ns " << G4endl;
        pParticleChange->AddSecondary(t);
      }
    }
  }
 
  if(0.0 == fParticleChange.GetProposedKineticEnergy() &&
     fAlive == fParticleChange.GetTrackStatus()) {
    if(particle->GetProcessManager()->GetAtRestProcessVector()->size() > 0)
         { fParticleChange.ProposeTrackStatus(fStopButAlive); }
    else { fParticleChange.ProposeTrackStatus(fStopAndKill); }
  }
 
  /*
  if(-1 < verboseLevel) {
    G4cout << "::PostStepDoIt: Sample secondary; Efin= " 
    << fParticleChange.GetProposedKineticEnergy()/MeV
           << " MeV; model= (" << currentModel->LowEnergyLimit()
           << ", " <<  currentModel->HighEnergyLimit() << ")"
           << "  preStepLambda= " << preStepLambda
           << "  dir= " << track.GetMomentumDirection()
           << "  status= " << track.GetTrackStatus()
           << G4endl;
  }
  */
  return &fParticleChange;
}

◆ PostStepGetPhysicalInteractionLength()

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

overridevirtual

Reimplemented from G4VContinuousDiscreteProcess.

Reimplemented in G4ePolarizedIonisation.

Definition at line 1054 of file G4VEnergyLossProcess.cc.

{
  // condition is set to "Not Forced"
  *condition = NotForced;
  G4double x = DBL_MAX;
 
  // initialisation of material, mass, charge, model 
  // at the beginning of the step
  DefineMaterial(track.GetMaterialCutsCouple());
  preStepKinEnergy       = track.GetKineticEnergy();
  preStepLogKinEnergy    = track.GetDynamicParticle()->GetLogKineticEnergy();
  preStepScaledEnergy    = preStepKinEnergy*massRatio;
  preStepLogScaledEnergy = preStepLogKinEnergy + logMassRatio;
  SelectModel(preStepScaledEnergy);
 
  if(!currentModel->IsActive(preStepScaledEnergy)) { 
    theNumberOfInteractionLengthLeft = -1.0;
    currentInteractionLength = DBL_MAX;
    return x; 
  }
 
  // change effective charge of an ion on fly
  if(isIon) {
    G4double q2 = currentModel->ChargeSquareRatio(track);
    if(q2 != chargeSqRatio && q2 > 0.0) {
      chargeSqRatio = q2;
      fFactor = q2*biasFactor*(*theDensityFactor)[currentCoupleIndex];
      reduceFactor = 1.0/(fFactor*massRatio);
    }
  }
  //  if(particle->GetPDGMass() > 0.9*GeV)
  //G4cout << "q2= "<<chargeSqRatio << " massRatio= " << massRatio << G4endl; 
 
  // forced biasing only for primary particles
  if(biasManager) {
    if(0 == track.GetParentID() && biasFlag && 
       biasManager->ForcedInteractionRegion(currentCoupleIndex)) {
      return biasManager->GetStepLimit(currentCoupleIndex, previousStepSize);
    }
  }
 
  // compute mean free path
  if(preStepScaledEnergy < mfpKinEnergy) {
    if (integral) {
      ComputeLambdaForScaledEnergy(preStepScaledEnergy, preStepLogScaledEnergy);
    } else {
      preStepLambda = 
          GetLambdaForScaledEnergy(preStepScaledEnergy, preStepLogScaledEnergy);
    }
 
    // zero cross section
    if(preStepLambda <= 0.0) { 
      theNumberOfInteractionLengthLeft = -1.0;
      currentInteractionLength = DBL_MAX;
    }
  }
 
  // non-zero cross section
  if(preStepLambda > 0.0) { 
    if (theNumberOfInteractionLengthLeft < 0.0) {
 
      // beggining of tracking (or just after DoIt of this process)
      theNumberOfInteractionLengthLeft =  -G4Log( G4UniformRand() );
      theInitialNumberOfInteractionLength = theNumberOfInteractionLengthLeft; 
 
    } else if(currentInteractionLength < DBL_MAX) {
 
      // subtract NumberOfInteractionLengthLeft using previous step
      theNumberOfInteractionLengthLeft -= 
        previousStepSize/currentInteractionLength;
 
      theNumberOfInteractionLengthLeft = 
        std::max(theNumberOfInteractionLengthLeft, 0.0);
    }
 
    // new mean free path and step limit
    currentInteractionLength = 1.0/preStepLambda;
    x = theNumberOfInteractionLengthLeft * currentInteractionLength;
  }
#ifdef G4VERBOSE
  if (verboseLevel>2){
    //  if(particle->GetPDGMass() > 0.9*GeV){
    G4cout << "G4VEnergyLossProcess::PostStepGetPhysicalInteractionLength ";
    G4cout << "[ " << GetProcessName() << "]" << G4endl; 
    G4cout << " for " << track.GetDefinition()->GetParticleName() 
           << " in Material  " <<  currentMaterial->GetName()
           << " Ekin(MeV)= " << preStepKinEnergy/MeV 
           << "  " << track.GetMaterial()->GetName()
           <<G4endl;
    G4cout << "MeanFreePath = " << currentInteractionLength/cm << "[cm]" 
           << "InteractionLength= " << x/cm <<"[cm] " <<G4endl;
  }
#endif
  return x;
}

Referenced by G4ePolarizedIonisation::PostStepGetPhysicalInteractionLength().

◆ PreparePhysicsTable()

void G4VEnergyLossProcess::PreparePhysicsTable ( const G4ParticleDefinition & part )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 360 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::PreparePhysicsTable for "
           << GetProcessName() << " for " << part.GetParticleName() 
           << "  " << this << G4endl;
  }
  isMaster = lManager->IsMaster();
 
  currentCouple = nullptr;
  preStepLambda = 0.0;
  mfpKinEnergy  = DBL_MAX;
  fRange        = DBL_MAX;
  preStepKinEnergy    = 0.0;
  preStepLogKinEnergy = LOG_EKIN_MIN;
  preStepRangeEnergy  = 0.0;
  chargeSqRatio = 1.0;
  massRatio     = 1.0;
  logMassRatio  = 0.;
  reduceFactor  = 1.0;
  fFactor = 1.0;
  lastIdx = 0;
 
  // Are particle defined?
  if( !particle ) { particle = &part; }
 
  if(part.GetParticleType() == "nucleus") {
 
    G4String pname = part.GetParticleName();
    if(pname != "deuteron" && pname != "triton" &&
       pname != "alpha+"   && pname != "helium" &&
       pname != "hydrogen") {
 
      if(!theGenericIon) {
        theGenericIon = 
          G4ParticleTable::GetParticleTable()->FindParticle("GenericIon");
      }
      isIon = true; 
      if(theGenericIon && particle != theGenericIon) {
        G4ProcessManager* pm =  theGenericIon->GetProcessManager();
        G4ProcessVector* v = pm->GetAlongStepProcessVector();
        size_t n = v->size();
        for(size_t j=0; j<n; ++j) {
          if((*v)[j] == this) {
            particle = theGenericIon;
            break;
          } 
        }
      }
    }
  }
 
  if( particle != &part ) {
    if(!isIon) {
      lManager->RegisterExtraParticle(&part, this);
    }
    if(1 < verboseLevel) {
      G4cout << "### G4VEnergyLossProcess::PreparePhysicsTable()"
             << " interrupted for "
             << part.GetParticleName() << "  isIon= " << isIon 
             << "  particle " << particle << "  GenericIon " << theGenericIon 
             << G4endl;
    }
    return;
  }
 
  Clean();
  lManager->PreparePhysicsTable(&part, this, isMaster);
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
 
  // Base particle and set of models can be defined here
  InitialiseEnergyLossProcess(particle, baseParticle);
 
  const G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  size_t n = theCoupleTable->GetTableSize();
 
  theDEDXAtMaxEnergy.resize(n, 0.0);
  theRangeAtMaxEnergy.resize(n, 0.0);
  theEnergyOfCrossSectionMax.resize(n, 0.0);
  theCrossSectionMax.resize(n, DBL_MAX);
 
  // parameters of the process
  if(!actIntegral) { integral = theParameters->Integral(); }
  if(!actLossFluc) { lossFluctuationFlag = theParameters->LossFluctuation(); }
  rndmStepFlag = theParameters->UseCutAsFinalRange();
  if(!actMinKinEnergy) { minKinEnergy = theParameters->MinKinEnergy(); }
  if(!actMaxKinEnergy) { maxKinEnergy = theParameters->MaxKinEnergy(); }
  if(!actBinning) { 
    nBins = theParameters->NumberOfBinsPerDecade()
      *G4lrint(std::log10(maxKinEnergy/minKinEnergy));
  }
  maxKinEnergyCSDA = theParameters->MaxEnergyForCSDARange();
  nBinsCSDA = theParameters->NumberOfBinsPerDecade()
    *G4lrint(std::log10(maxKinEnergyCSDA/minKinEnergy));
  if(!actLinLossLimit) { linLossLimit = theParameters->LinearLossLimit(); }
  lambdaFactor    = theParameters->LambdaFactor();
  logLambdafactor = G4Log(lambdaFactor);
  if(isMaster) { SetVerboseLevel(theParameters->Verbose()); }
  else {  SetVerboseLevel(theParameters->WorkerVerbose()); }
 
  theParameters->DefineRegParamForLoss(this);
 
  G4double initialCharge = particle->GetPDGCharge();
  G4double initialMass   = particle->GetPDGMass();
 
  theParameters->FillStepFunction(particle, this);
 
  if (baseParticle) {
    massRatio    = (baseParticle->GetPDGMass())/initialMass;
    logMassRatio = G4Log(massRatio);
    G4double q = initialCharge/baseParticle->GetPDGCharge();
    chargeSqRatio = q*q;
    if(chargeSqRatio > 0.0) { reduceFactor = 1.0/(chargeSqRatio*massRatio); }
  }
  lowestKinEnergy = (initialMass < MeV) ? theParameters->LowestElectronEnergy()
    : theParameters->LowestMuHadEnergy();
 
  // Tables preparation
  if (isMaster && !baseParticle) {
 
    if(theDEDXTable && isIonisation) {
      if(theIonisationTable && theDEDXTable != theIonisationTable) {
        theDEDXTable->clearAndDestroy();
        delete theDEDXTable;
        theDEDXTable = theIonisationTable;
      }   
      if(theDEDXSubTable && theIonisationSubTable && 
         theDEDXSubTable != theIonisationSubTable) {
        theDEDXSubTable->clearAndDestroy();
        delete theDEDXSubTable;
        theDEDXSubTable = theIonisationSubTable;
      }   
    }
    
    theDEDXTable = G4PhysicsTableHelper::PreparePhysicsTable(theDEDXTable);
    bld->InitialiseBaseMaterials(theDEDXTable);
 
    if(theDEDXSubTable) {
      theDEDXSubTable = 
        G4PhysicsTableHelper::PreparePhysicsTable(theDEDXSubTable);
    }
 
    if (theParameters->BuildCSDARange()) {
      theDEDXunRestrictedTable = 
        G4PhysicsTableHelper::PreparePhysicsTable(theDEDXunRestrictedTable);
      theCSDARangeTable = 
        G4PhysicsTableHelper::PreparePhysicsTable(theCSDARangeTable);
    }
 
    theLambdaTable = G4PhysicsTableHelper::PreparePhysicsTable(theLambdaTable);
 
    if(isIonisation) {
      theRangeTableForLoss = 
        G4PhysicsTableHelper::PreparePhysicsTable(theRangeTableForLoss);
      theInverseRangeTable = 
        G4PhysicsTableHelper::PreparePhysicsTable(theInverseRangeTable);  
    }
 
    if (nSCoffRegions && !lManager->SubCutProducer()) {
      theDEDXSubTable = 
        G4PhysicsTableHelper::PreparePhysicsTable(theDEDXSubTable);
      theSubLambdaTable = 
        G4PhysicsTableHelper::PreparePhysicsTable(theSubLambdaTable);
    }
  }
  /*
  G4cout << "** G4VEnergyLossProcess::PreparePhysicsTable() for " 
         << GetProcessName() << " and " << particle->GetParticleName()
         << " isMaster: " << isMaster << " isIonisation: " 
         << isIonisation << G4endl;
  G4cout << " theDEDX: " << theDEDXTable 
         << " theRange: " << theRangeTableForLoss
         << " theInverse: " << theInverseRangeTable
         << " theLambda: " << theLambdaTable << G4endl;
  */
  // forced biasing
  if(biasManager) { 
    biasManager->Initialise(part,GetProcessName(),verboseLevel); 
    biasFlag = false; 
  }
 
  // defined ID of secondary particles
  if(isMaster) {
    G4String nam1 = GetProcessName();
    G4String nam4 = nam1 + "_split";
    G4String nam5 = nam1 + "_subcut";
    secID   = G4PhysicsModelCatalog::Register(nam1); 
    biasID  = G4PhysicsModelCatalog::Register(nam4); 
    subsecID= G4PhysicsModelCatalog::Register(nam5);
  } 
 
  // initialisation of models
  G4int nmod = modelManager->NumberOfModels();
  for(G4int i=0; i<nmod; ++i) {
    G4VEmModel* mod = modelManager->GetModel(i);
    mod->SetMasterThread(isMaster);
    mod->SetAngularGeneratorFlag(
      theParameters->UseAngularGeneratorForIonisation());
    if(mod->HighEnergyLimit() > maxKinEnergy) {
      mod->SetHighEnergyLimit(maxKinEnergy);
    }
  }
  theCuts = modelManager->Initialise(particle, secondaryParticle, 
                                     theParameters->MinSubRange(), 
                                     verboseLevel);
 
  // Sub Cutoff 
  if(nSCoffRegions > 0) {
    if(theParameters->MinSubRange() < 1.0) { useSubCutoff = true; }
 
    theSubCuts = modelManager->SubCutoff();
 
    idxSCoffRegions = new G4bool[n]; 
    for (size_t j=0; j<n; ++j) {
 
      const G4MaterialCutsCouple* couple = 
        theCoupleTable->GetMaterialCutsCouple(j);
      const G4ProductionCuts* pcuts = couple->GetProductionCuts();
      
      G4bool reg = false;
      for(G4int i=0; i<nSCoffRegions; ++i) {
        if( pcuts == scoffRegions[i]->GetProductionCuts()) { 
          reg = true;
          break; 
        }
      }
      idxSCoffRegions[j] = reg;
    }
  }
 
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::PrepearPhysicsTable() is done "
           << " for local " << particle->GetParticleName()
           << " isIon= " << isIon;
    if(baseParticle) { 
      G4cout << "; base: " << baseParticle->GetParticleName(); 
    }
    G4cout << " chargeSqRatio= " << chargeSqRatio
           << " massRatio= " << massRatio
           << " reduceFactor= " << reduceFactor << G4endl;
    if (nSCoffRegions) {
      G4cout << " SubCutoff Regime is ON for regions: " << G4endl;
      for (G4int i=0; i<nSCoffRegions; ++i) {
        const G4Region* r = scoffRegions[i];
        G4cout << "           " << r->GetName() << G4endl;
      }
    }
  }
}

◆ PrintInfo()

virtual void G4VEnergyLossProcess::PrintInfo ( )

inlinevirtual

Reimplemented in G4hIonisation, G4hhIonisation, G4mplIonisation, G4MuBremsstrahlung, G4MuIonisation, G4ePolarizedIonisation, G4alphaIonisation, and G4eIonisation.

Definition at line 102 of file G4VEnergyLossProcess.hh.

102{};

◆ ProcessDescription()

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

overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4hBremsstrahlung, G4hhIonisation, G4hPairProduction, G4mplIonisation, G4ePairProduction, G4MuBremsstrahlung, G4MuIonisation, G4MuPairProduction, G4alphaIonisation, G4eBremsstrahlung, G4eIonisation, G4hIonisation, and G4ionIonisation.

Definition at line 2368 of file G4VEnergyLossProcess.cc.

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

◆ RangeTableForLoss()

G4PhysicsTable * G4VEnergyLossProcess::RangeTableForLoss ( ) const

inline

Definition at line 1109 of file G4VEnergyLossProcess.hh.

{
  return theRangeTableForLoss;
}

Referenced by BuildPhysicsTable(), G4LossTableManager::LocalPhysicsTables(), and G4EmCalculator::PrintRangeTable().

◆ RetrievePhysicsTable()

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

overridevirtual

Reimplemented from G4VProcess.

Definition at line 1762 of file G4VEnergyLossProcess.cc.

{
  G4bool res = true;
  if (!isMaster) return res;
  const G4String particleName = part->GetParticleName();
 
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::RetrievePhysicsTable() for "
           << particleName << " and process " << GetProcessName()
           << "; tables_are_built= " << tablesAreBuilt
           << G4endl;
  }
  if(particle == part) {
 
    if ( !baseParticle ) {
 
      G4bool fpi = true;
      if(!RetrieveTable(part,theDEDXTable,ascii,directory,"DEDX",fpi)) 
        {fpi = false;}
 
      // ionisation table keeps individual dEdx and not sum of sub-processes
      if(!RetrieveTable(part,theDEDXTable,ascii,directory,"Ionisation",false)) 
        {fpi = false;}
 
      if(!RetrieveTable(part,theRangeTableForLoss,ascii,directory,"Range",fpi)) 
        {res = false;}
 
      if(!RetrieveTable(part,theDEDXunRestrictedTable,ascii,directory,
                        "DEDXnr",false)) 
        {res = false;}
 
      if(!RetrieveTable(part,theCSDARangeTable,ascii,directory,
                        "CSDARange",false)) 
        {res = false;}
 
      if(!RetrieveTable(part,theInverseRangeTable,ascii,directory,
                        "InverseRange",fpi)) 
        {res = false;}
 
      if(!RetrieveTable(part,theLambdaTable,ascii,directory,"Lambda",true)) 
        {res = false;}
 
      G4bool yes = false;
      if(nSCoffRegions > 0) {yes = true;}
 
      if(!RetrieveTable(part,theDEDXSubTable,ascii,directory,"SubDEDX",yes)) 
        {res = false;}
 
      if(!RetrieveTable(part,theSubLambdaTable,ascii,directory,
                        "SubLambda",yes)) 
        {res = false;}
 
      if(!fpi) yes = false;
      if(!RetrieveTable(part,theIonisationSubTable,ascii,directory,
                        "SubIonisation",yes))
        {res = false;}
    }
  }
 
  return res;
}

◆ SampleSubCutSecondaries()

G4double G4VEnergyLossProcess::SampleSubCutSecondaries	(	std::vector< G4Track * > &	tracks,
		const G4Step &	step,
		G4VEmModel *	model,
		G4int	matIdx
	)

Definition at line 1479 of file G4VEnergyLossProcess.cc.

{
  // Fast check weather subcutoff can work
  G4double esec = 0.0;
  G4double subcut = (*theSubCuts)[idx];
  G4double cut = (*theCuts)[idx];
  if(cut <= subcut) { return esec; }
 
  const G4Track* track = step.GetTrack();
  const G4DynamicParticle* dp = track->GetDynamicParticle();
  G4double e = dp->GetKineticEnergy()*massRatio;
  G4double cross = (*theDensityFactor)[idx]*chargeSqRatio
    *(((*theSubLambdaTable)[(*theDensityIdx)[idx]])->Value(e, idxSubLambda));
  G4double length = step.GetStepLength();
 
  // negligible probability to get any interaction
  if(length*cross < perMillion) { return esec; }
  /*      
  if(-1 < verboseLevel) 
    G4cout << "<<< Subcutoff for " << GetProcessName()
           << " cross(1/mm)= " << cross*mm << ">>>"
           << " e(MeV)= " << preStepScaledEnergy
           << " matIdx= " << currentCoupleIndex
           << G4endl;
  */
 
  // Sample subcutoff secondaries
  G4StepPoint* preStepPoint = step.GetPreStepPoint();
  G4StepPoint* postStepPoint = step.GetPostStepPoint();
  G4ThreeVector prepoint = preStepPoint->GetPosition();
  G4ThreeVector dr = postStepPoint->GetPosition() - prepoint;
  G4double pretime = preStepPoint->GetGlobalTime();
  G4double dt = postStepPoint->GetGlobalTime() - pretime;
  G4double fragment = 0.0;
 
  do {
    G4double del = -G4Log(G4UniformRand())/cross;
    fragment += del/length;
    if (fragment > 1.0) { break; }
 
    // sample secondaries
    secParticles.clear();
    model->SampleSecondaries(&secParticles,track->GetMaterialCutsCouple(),
                             dp,subcut,cut);
 
    // position of subcutoff particles
    G4ThreeVector r = prepoint + fragment*dr;
    std::vector<G4DynamicParticle*>::iterator it;
    for(it=secParticles.begin(); it!=secParticles.end(); ++it) {
 
      G4Track* t = new G4Track((*it), pretime + fragment*dt, r);
      t->SetTouchableHandle(track->GetTouchableHandle());
      t->SetCreatorModelIndex(subsecID);
      tracks.push_back(t);
      esec += t->GetKineticEnergy();
      if (t->GetParticleDefinition() == thePositron) { 
        esec += 2.0*electron_mass_c2; 
      }
 
        /*        
        if(-1 < verboseLevel) 
          G4cout << "New track " 
                 << t->GetParticleDefinition()->GetParticleName()
                 << " e(keV)= " << t->GetKineticEnergy()/keV
                 << " fragment= " << fragment
                 << G4endl;
        */
    }
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  } while (fragment <= 1.0);
  return esec;
} 

Referenced by AlongStepDoIt().

◆ SecondaryParticle()

const G4ParticleDefinition * G4VEnergyLossProcess::SecondaryParticle ( ) const

inline

Definition at line 988 of file G4VEnergyLossProcess.hh.

{
  return secondaryParticle;
}

◆ SecondaryRangeTable()

G4PhysicsTable * G4VEnergyLossProcess::SecondaryRangeTable ( ) const

inline

Definition at line 1102 of file G4VEnergyLossProcess.hh.

{
  return theSecondaryRangeTable;
}

Referenced by BuildPhysicsTable().

◆ SelectModel()

void G4VEnergyLossProcess::SelectModel ( G4double kinEnergy )

inlineprotected

Definition at line 602 of file G4VEnergyLossProcess.hh.

{
  currentModel = modelManager->SelectModel(kinEnergy, currentCoupleIndex);
  currentModel->SetCurrentCouple(currentCouple);
}

Referenced by CrossSectionPerVolume(), GetDEDXDispersion(), PostStepDoIt(), and PostStepGetPhysicalInteractionLength().

◆ SelectModelForMaterial()

G4VEmModel * G4VEnergyLossProcess::SelectModelForMaterial	(	G4double	kinEnergy,
		size_t &	idx
	)		const

inline

Definition at line 610 of file G4VEnergyLossProcess.hh.

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

Referenced by G4ContinuousGainOfEnergy::GetContinuousStepLimit().

◆ SetBaseParticle()

void G4VEnergyLossProcess::SetBaseParticle ( const G4ParticleDefinition * p )

inline

Definition at line 966 of file G4VEnergyLossProcess.hh.

{
  baseParticle = p;
}

Referenced by G4hhIonisation::InitialiseEnergyLossProcess(), G4mplIonisation::InitialiseEnergyLossProcess(), G4alphaIonisation::InitialiseEnergyLossProcess(), G4hIonisation::InitialiseEnergyLossProcess(), and G4ionIonisation::InitialiseEnergyLossProcess().

◆ SetCrossSectionBiasingFactor()

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

Definition at line 2218 of file G4VEnergyLossProcess.cc.

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

Referenced by G4EmExtraParameters::DefineRegParamForLoss().

◆ SetCSDARangeTable()

void G4VEnergyLossProcess::SetCSDARangeTable ( G4PhysicsTable * pRange )

Definition at line 2073 of file G4VEnergyLossProcess.cc.

{
  theCSDARangeTable = p; 
 
  if(p) {
    size_t n = p->length();
    G4PhysicsVector* pv;
    G4double emax = maxKinEnergyCSDA;
 
    for (size_t i=0; i<n; ++i) {
      pv = (*p)[i];
      G4double rmax = 0.0;
      if(pv) { rmax = pv->Value(emax, idxCSDA); }
      else {
        pv = (*p)[(*theDensityIdx)[i]];
        if(pv) { rmax = pv->Value(emax, idxCSDA)/(*theDensityFactor)[i]; }
      }
      theRangeAtMaxEnergy[i] = rmax;
      //G4cout << "i= " << i << " Emax(MeV)= " << emax/MeV << " Rmax= " 
      //<< rmax<< G4endl;
    }
  }
}

Referenced by BuildPhysicsTable().

◆ SetDEDXBinning()

void G4VEnergyLossProcess::SetDEDXBinning ( G4int nbins )

Definition at line 2320 of file G4VEnergyLossProcess.cc.

{
  if(2 < n && n < 1000000000) { 
    nBins = n; 
    actBinning = true;
  } else {
    G4double e = (G4double)n;
    PrintWarning("SetDEDXBinning", e); 
  } 
}

Referenced by G4hhIonisation::InitialiseEnergyLossProcess(), G4mplIonisation::InitialiseEnergyLossProcess(), and G4hIonisation::InitialiseEnergyLossProcess().

◆ SetDEDXTable()

void G4VEnergyLossProcess::SetDEDXTable	(	G4PhysicsTable *	p,
		G4EmTableType	tType
	)

Definition at line 2012 of file G4VEnergyLossProcess.cc.

{
  if(fTotal == tType) {
    theDEDXunRestrictedTable = p;
    if(p) {
      size_t n = p->length();
      G4PhysicsVector* pv = (*p)[0];
      G4double emax = maxKinEnergyCSDA;
 
      G4LossTableBuilder* bld = lManager->GetTableBuilder();
      theDensityFactor = bld->GetDensityFactors();
      theDensityIdx = bld->GetCoupleIndexes();
 
      for (size_t i=0; i<n; ++i) {
        G4double dedx = 0.0; 
        pv = (*p)[i];
        if(pv) { 
          dedx = pv->Value(emax, idxDEDXunRestricted); 
        } else {
          pv = (*p)[(*theDensityIdx)[i]];
          if(pv) { 
            dedx = 
              pv->Value(emax, idxDEDXunRestricted)*(*theDensityFactor)[i]; 
          }
        }
        theDEDXAtMaxEnergy[i] = dedx;
        //G4cout << "i= " << i << " emax(MeV)= " << emax/MeV<< " dedx= " 
        //     << dedx << G4endl;
      }
    }
 
  } else if(fRestricted == tType) {
    /*
      G4cout<< "G4VEnergyLossProcess::SetDEDXTable "
            << particle->GetParticleName()
            << " oldTable " << theDEDXTable << " newTable " << p 
            << " ion " << theIonisationTable 
            << " IsMaster " << isMaster 
            << " " << GetProcessName() << G4endl;
      G4cout << (*p) << G4endl;
    */
    theDEDXTable = p;
  } else if(fSubRestricted == tType) {
      theDEDXSubTable = p;
  } else if(fIsIonisation == tType) {
    /*
      G4cout<< "G4VEnergyLossProcess::SetIonisationTable "
            << particle->GetParticleName()
            << " oldTable " << theDEDXTable << " newTable " << p 
            << " ion " << theIonisationTable 
            << " IsMaster " << isMaster 
            << " " << GetProcessName() << G4endl;
    */
    theIonisationTable = p;
  } else if(fIsSubIonisation == tType) {
    theIonisationSubTable = p;
  }
}

Referenced by BuildPhysicsTable().

◆ SetDynamicMassCharge()

void G4VEnergyLossProcess::SetDynamicMassCharge	(	G4double	massratio,
		G4double	charge2ratio
	)

inline

Definition at line 635 of file G4VEnergyLossProcess.hh.

{
  massRatio     = massratio;
  logMassRatio  = G4Log(massRatio);
  fFactor = charge2ratio*biasFactor*(*theDensityFactor)[currentCoupleIndex];
  chargeSqRatio = charge2ratio;
  reduceFactor  = 1.0/(fFactor*massRatio);
}

Referenced by G4ContinuousGainOfEnergy::AlongStepDoIt(), and G4ContinuousGainOfEnergy::GetContinuousStepLimit().

◆ SetEmModel()

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

Definition at line 330 of file G4VEnergyLossProcess.cc.

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

◆ SetFluctModel()

void G4VEnergyLossProcess::SetFluctModel ( G4VEmFluctuationModel * p )

inline

Definition at line 936 of file G4VEnergyLossProcess.hh.

{
  fluctModel = p;
}

Referenced by G4EmStandardPhysicsWVI::ConstructProcess(), G4MuIonisation::InitialiseEnergyLossProcess(), G4ePolarizedIonisation::InitialiseEnergyLossProcess(), G4alphaIonisation::InitialiseEnergyLossProcess(), G4eIonisation::InitialiseEnergyLossProcess(), G4hIonisation::InitialiseEnergyLossProcess(), G4ionIonisation::InitialiseEnergyLossProcess(), and G4EmConfigurator::PrepareModels().

◆ SetIntegral()

void G4VEnergyLossProcess::SetIntegral ( G4bool val )

inline

Definition at line 1003 of file G4VEnergyLossProcess.hh.

{
  integral = val;
  actIntegral = true;
}

◆ SetInverseRangeTable()

void G4VEnergyLossProcess::SetInverseRangeTable ( G4PhysicsTable * p )

Definition at line 2123 of file G4VEnergyLossProcess.cc.

{
  theInverseRangeTable = p;
  if(1 < verboseLevel) {
    G4cout << "### Set InverseRange table " << p 
           << " for " << particle->GetParticleName()
           << " and process " << GetProcessName() << G4endl;
  }
}

Referenced by BuildPhysicsTable().

◆ SetIonisation()

void G4VEnergyLossProcess::SetIonisation ( G4bool val )

Definition at line 2281 of file G4VEnergyLossProcess.cc.

{
  isIonisation = val;
  aGPILSelection = (val) ? CandidateForSelection : NotCandidateForSelection;
}

Referenced by G4LossTableManager::BuildPhysicsTable(), G4eBremsstrahlung::G4eBremsstrahlung(), G4ePairProduction::G4ePairProduction(), G4MuBremsstrahlung::G4MuBremsstrahlung(), G4MuPairProduction::G4MuPairProduction(), and G4ePolarizedBremsstrahlung::InitialiseEnergyLossProcess().

◆ SetLambdaTable()

void G4VEnergyLossProcess::SetLambdaTable ( G4PhysicsTable * p )

Definition at line 2135 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {
    G4cout << "### Set Lambda table " << p 
           << " for " << particle->GetParticleName()
           << " and process " << GetProcessName() << G4endl;
    //G4cout << *p << G4endl;
  }
  theLambdaTable = p; 
  tablesAreBuilt = true;
 
  G4LossTableBuilder* bld = lManager->GetTableBuilder();
  theDensityFactor = bld->GetDensityFactors();
  theDensityIdx = bld->GetCoupleIndexes();
 
  if(theLambdaTable) {
    size_t n = theLambdaTable->length();
    G4PhysicsVector* pv = (*theLambdaTable)[0];
    G4double e, ss, smax, emax;
 
    size_t i;
 
    // first loop on existing vectors
    for (i=0; i<n; ++i) {
      pv = (*theLambdaTable)[i];
      if(pv) {
        size_t nb = pv->GetVectorLength();
        emax = DBL_MAX;
        smax = 0.0;
        if(nb > 0) {
          for (size_t j=0; j<nb; ++j) {
            e = pv->Energy(j);
            ss = (*pv)(j);
            if(ss > smax) {
              smax = ss;
              emax = e;
            }
          }
        }
        theEnergyOfCrossSectionMax[i] = emax;
        theCrossSectionMax[i] = smax;
        if(1 < verboseLevel) {
          G4cout << "For " << particle->GetParticleName() 
                 << " Max CS at i= " << i << " emax(MeV)= " << emax/MeV
                 << " lambda= " << smax << G4endl;
        }
      }
    }
    // second loop using base materials
    for (i=0; i<n; ++i) {
      pv = (*theLambdaTable)[i];
      if(!pv){
        G4int j = (*theDensityIdx)[i];
        theEnergyOfCrossSectionMax[i] = theEnergyOfCrossSectionMax[j];
        theCrossSectionMax[i] = (*theDensityFactor)[i]*theCrossSectionMax[j];
      }
    }
  }
}

Referenced by BuildPhysicsTable().

◆ SetLinearLossLimit()

void G4VEnergyLossProcess::SetLinearLossLimit ( G4double val )

Definition at line 2289 of file G4VEnergyLossProcess.cc.

{
  if(0.0 < val && val < 1.0) { 
    linLossLimit = val;
    actLinLossLimit = true; 
  } else { PrintWarning("SetLinearLossLimit", val); }
}

Referenced by G4alphaIonisation::G4alphaIonisation(), and G4ionIonisation::G4ionIonisation().

◆ SetLossFluctuations()

void G4VEnergyLossProcess::SetLossFluctuations ( G4bool val )

inline

Definition at line 995 of file G4VEnergyLossProcess.hh.

{
  lossFluctuationFlag = val;
  actLossFluc = true;
}

◆ SetLowestEnergyLimit()

void G4VEnergyLossProcess::SetLowestEnergyLimit ( G4double val )

Definition at line 2312 of file G4VEnergyLossProcess.cc.

{
  if(1.e-18 < val && val < 1.e+50) { lowestKinEnergy = val; }
  else { PrintWarning("SetLowestEnergyLimit", val); }
}

◆ SetMaxKinEnergy()

void G4VEnergyLossProcess::SetMaxKinEnergy ( G4double e )

Definition at line 2343 of file G4VEnergyLossProcess.cc.

{
  if(minKinEnergy < e && e < 1.e+50) { 
    maxKinEnergy = e;
    actMaxKinEnergy = true;
    if(e < maxKinEnergyCSDA) { maxKinEnergyCSDA = e; }
  } else { PrintWarning("SetMaxKinEnergy", e); } 
}

Referenced by G4hhIonisation::InitialiseEnergyLossProcess(), G4mplIonisation::InitialiseEnergyLossProcess(), and G4hIonisation::InitialiseEnergyLossProcess().

◆ SetMinKinEnergy()

void G4VEnergyLossProcess::SetMinKinEnergy ( G4double e )

Definition at line 2333 of file G4VEnergyLossProcess.cc.

{
  if(1.e-18 < e && e < maxKinEnergy) { 
    minKinEnergy = e; 
    actMinKinEnergy = true;
  } else { PrintWarning("SetMinKinEnergy", e); } 
}

Referenced by G4hhIonisation::InitialiseEnergyLossProcess(), G4mplIonisation::InitialiseEnergyLossProcess(), and G4hIonisation::InitialiseEnergyLossProcess().

◆ SetParticle()

void G4VEnergyLossProcess::SetParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 950 of file G4VEnergyLossProcess.hh.

{
  particle = p;
}

◆ SetRangeTableForLoss()

void G4VEnergyLossProcess::SetRangeTableForLoss ( G4PhysicsTable * p )

Definition at line 2099 of file G4VEnergyLossProcess.cc.

{
  theRangeTableForLoss = p;
  if(1 < verboseLevel) {
    G4cout << "### Set Range table " << p 
           << " for " << particle->GetParticleName()
           << " and process " << GetProcessName() << G4endl;
  }
}

Referenced by BuildPhysicsTable().

◆ SetSecondaryParticle()

void G4VEnergyLossProcess::SetSecondaryParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 958 of file G4VEnergyLossProcess.hh.

{
  secondaryParticle = p;
}

Referenced by G4eBremsstrahlung::G4eBremsstrahlung(), G4eIonisation::G4eIonisation(), G4ePairProduction::G4ePairProduction(), G4ePolarizedIonisation::G4ePolarizedIonisation(), G4hhIonisation::G4hhIonisation(), G4hIonisation::G4hIonisation(), G4ionIonisation::G4ionIonisation(), G4mplIonisation::G4mplIonisation(), G4MuBremsstrahlung::G4MuBremsstrahlung(), G4MuIonisation::G4MuIonisation(), G4MuPairProduction::G4MuPairProduction(), G4ePolarizedBremsstrahlung::InitialiseEnergyLossProcess(), and G4alphaIonisation::InitialiseEnergyLossProcess().

◆ SetSecondaryRangeTable()

void G4VEnergyLossProcess::SetSecondaryRangeTable ( G4PhysicsTable * p )

Definition at line 2111 of file G4VEnergyLossProcess.cc.

{
  theSecondaryRangeTable = p;
  if(1 < verboseLevel) {
    G4cout << "### Set SecondaryRange table " << p 
           << " for " << particle->GetParticleName()
           << " and process " << GetProcessName() << G4endl;
  }
}

Referenced by BuildPhysicsTable().

◆ SetStepFunction()

void G4VEnergyLossProcess::SetStepFunction	(	G4double	v1,
		G4double	v2
	)

Definition at line 2299 of file G4VEnergyLossProcess.cc.

{
  if(0.0 < v1 && 0.0 < v2) { 
    dRoverRange = std::min(1.0, v1);
    finalRange = std::min(v2, 1.e+50);
  } else {
    PrintWarning("SetStepFunctionV1", v1); 
    PrintWarning("SetStepFunctionV2", v2); 
  }
}

◆ SetSubLambdaTable()

void G4VEnergyLossProcess::SetSubLambdaTable ( G4PhysicsTable * p )

Definition at line 2197 of file G4VEnergyLossProcess.cc.

{
  theSubLambdaTable = p;
  if(1 < verboseLevel) {
    G4cout << "### Set SebLambda table " << p 
           << " for " << particle->GetParticleName()
           << " and process " << GetProcessName() << G4endl;
  }
}

Referenced by BuildPhysicsTable().

◆ StartTracking()

void G4VEnergyLossProcess::StartTracking ( G4Track * track )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 981 of file G4VEnergyLossProcess.cc.

{
  /*      
    G4cout << track->GetDefinition()->GetParticleName() 
           << " e(MeV)= " << track->GetKineticEnergy()
           << "  baseParticle " << baseParticle << " proc " << this;
    if(particle) G4cout << "  " << particle->GetParticleName();
    G4cout << " isIon= " << isIon << " dedx " << theDEDXTable <<G4endl;
  */
  // reset parameters for the new track
  theNumberOfInteractionLengthLeft = -1.0;
  currentInteractionLength = mfpKinEnergy = DBL_MAX; 
  preStepRangeEnergy = 0.0;
 
  // reset ion
  if(isIon) {
    chargeSqRatio = 0.5;
 
    G4double newmass = track->GetDefinition()->GetPDGMass();
    if(baseParticle) {
      massRatio    = baseParticle->GetPDGMass()/newmass;
      logMassRatio = G4Log(massRatio);
    } else if(theGenericIon) {
      massRatio    = proton_mass_c2/newmass;
      logMassRatio = G4Log(massRatio);
    } else {
      massRatio    = 1.0;
      logMassRatio = 0.0;
    }
  }  
  // forced biasing only for primary particles
  if(biasManager) {
    if(0 == track->GetParentID()) {
      // primary particle
      biasFlag = true; 
      biasManager->ResetForcedInteraction(); 
    }
  }
}

◆ StorePhysicsTable()

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

overridevirtual

Reimplemented from G4VProcess.

Definition at line 1699 of file G4VEnergyLossProcess.cc.

{
  G4bool res = true;
  //G4cout << "G4VEnergyLossProcess::StorePhysicsTable: " << part->GetParticleName()
  //         << "  " << directory << "  " << ascii << G4endl;
  if (!isMaster || baseParticle || part != particle ) return res;
 
  if(!StoreTable(part,theDEDXTable,ascii,directory,"DEDX")) 
    {res = false;}
 
  if(!StoreTable(part,theDEDXunRestrictedTable,ascii,directory,"DEDXnr")) 
    {res = false;}
 
  if(!StoreTable(part,theDEDXSubTable,ascii,directory,"SubDEDX")) 
    {res = false;}
 
  if(!StoreTable(part,theIonisationTable,ascii,directory,"Ionisation")) 
    {res = false;}
 
  if(!StoreTable(part,theIonisationSubTable,ascii,directory,"SubIonisation")) 
    {res = false;}
 
  if(isIonisation &&
     !StoreTable(part,theCSDARangeTable,ascii,directory,"CSDARange")) 
    {res = false;}
 
  if(isIonisation &&
     !StoreTable(part,theRangeTableForLoss,ascii,directory,"Range")) 
    {res = false;}
  
  if(isIonisation &&
     !StoreTable(part,theInverseRangeTable,ascii,directory,"InverseRange")) 
    {res = false;}
  
  if(!StoreTable(part,theLambdaTable,ascii,directory,"Lambda")) 
    {res = false;}
 
  if(!StoreTable(part,theSubLambdaTable,ascii,directory,"SubLambda")) 
    {res = false;}
 
  if ( !res ) {
    if(1 < verboseLevel) {
      G4cout << "Physics tables are stored for " 
             << particle->GetParticleName()
             << " and process " << GetProcessName()
             << " in the directory <" << directory
             << "> " << G4endl;
    }
  } else {
    G4cout << "Fail to store Physics Tables for " 
           << particle->GetParticleName()
           << " and process " << GetProcessName()
           << " in the directory <" << directory
           << "> " << G4endl;
  }
  return res;
}

◆ StreamProcessInfo()

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

inlineprotectedvirtual

Reimplemented in G4ePairProduction, G4MuPairProduction, G4eBremsstrahlung, and G4ionIonisation.

Definition at line 108 of file G4VEnergyLossProcess.hh.

108{};

◆ SubLambdaTable()

G4PhysicsTable * G4VEnergyLossProcess::SubLambdaTable ( ) const

inline

Definition at line 1130 of file G4VEnergyLossProcess.hh.

{
  return theSubLambdaTable;
}

Referenced by BuildPhysicsTable().

◆ TablesAreBuilt()

G4bool G4VEnergyLossProcess::TablesAreBuilt ( ) const

inline

Definition at line 1053 of file G4VEnergyLossProcess.hh.

{
  return tablesAreBuilt;
}

◆ UpdateEmModel()

void G4VEnergyLossProcess::UpdateEmModel	(	const G4String &	nam,
		G4double	emin,
		G4double	emax
	)

Definition at line 322 of file G4VEnergyLossProcess.cc.

{
  modelManager->UpdateEmModel(nam, emin, emax);
}

Member Data Documentation

◆ computedRange

G4double G4VEnergyLossProcess::computedRange

protected

Definition at line 559 of file G4VEnergyLossProcess.hh.

Referenced by G4VEnergyLossProcess().

◆ currentCouple

const G4MaterialCutsCouple* G4VEnergyLossProcess::currentCouple

protected

Definition at line 554 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), AlongStepGetPhysicalInteractionLength(), PostStepDoIt(), PreparePhysicsTable(), and SelectModel().

◆ currentCoupleIndex

size_t G4VEnergyLossProcess::currentCoupleIndex

protected

Definition at line 555 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), CrossSectionPerVolume(), CurrentMaterialCutsCoupleIndex(), G4VEnergyLossProcess(), GetDEDXDispersion(), PostStepDoIt(), PostStepGetPhysicalInteractionLength(), SelectModel(), and SetDynamicMassCharge().

◆ currentMaterial

const G4Material* G4VEnergyLossProcess::currentMaterial

protected

Definition at line 553 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), CrossSectionPerVolume(), G4VEnergyLossProcess(), GetDEDXDispersion(), and PostStepGetPhysicalInteractionLength().

◆ fParticleChange

G4ParticleChangeForLoss G4VEnergyLossProcess::fParticleChange

protected

Definition at line 552 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), G4VEnergyLossProcess(), and PostStepDoIt().

◆ fRange

G4double G4VEnergyLossProcess::fRange

protected

Definition at line 558 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), AlongStepGetPhysicalInteractionLength(), G4VEnergyLossProcess(), GetRange(), and PreparePhysicsTable().

◆ mfpKinEnergy

G4double G4VEnergyLossProcess::mfpKinEnergy

protected

Definition at line 565 of file G4VEnergyLossProcess.hh.

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

◆ preStepKinEnergy

G4double G4VEnergyLossProcess::preStepKinEnergy

protected

Definition at line 560 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), G4VEnergyLossProcess(), PostStepGetPhysicalInteractionLength(), and PreparePhysicsTable().

◆ preStepLambda

G4double G4VEnergyLossProcess::preStepLambda

protected

Definition at line 557 of file G4VEnergyLossProcess.hh.

Referenced by G4VEnergyLossProcess(), PostStepDoIt(), PostStepGetPhysicalInteractionLength(), and PreparePhysicsTable().

◆ preStepLogKinEnergy

G4double G4VEnergyLossProcess::preStepLogKinEnergy

protected

Definition at line 561 of file G4VEnergyLossProcess.hh.

Referenced by G4VEnergyLossProcess(), PostStepGetPhysicalInteractionLength(), and PreparePhysicsTable().

◆ preStepLogScaledEnergy

G4double G4VEnergyLossProcess::preStepLogScaledEnergy

protected

Definition at line 563 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), AlongStepGetPhysicalInteractionLength(), G4VEnergyLossProcess(), and PostStepGetPhysicalInteractionLength().

◆ preStepRangeEnergy

G4double G4VEnergyLossProcess::preStepRangeEnergy

protected

Definition at line 564 of file G4VEnergyLossProcess.hh.

Referenced by G4VEnergyLossProcess(), PreparePhysicsTable(), and StartTracking().

◆ preStepScaledEnergy

G4double G4VEnergyLossProcess::preStepScaledEnergy

protected

Definition at line 562 of file G4VEnergyLossProcess.hh.

Referenced by AlongStepDoIt(), AlongStepGetPhysicalInteractionLength(), G4VEnergyLossProcess(), and PostStepGetPhysicalInteractionLength().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4VEnergyLossProcess()

◆ ~G4VEnergyLossProcess()

Member Function Documentation

◆ ActivateForcedInteraction()

◆ ActivateSecondaryBiasing()

◆ ActivateSubCutoff()

◆ AddCollaborativeProcess()

◆ AddEmModel()

◆ AlongStepDoIt()

◆ AlongStepGetPhysicalInteractionLength()

◆ BaseParticle()

◆ BuildDEDXTable()

◆ BuildLambdaTable()

◆ BuildPhysicsTable()

◆ ContinuousStepLimit()

◆ CrossSectionBiasingFactor()

◆ CrossSectionPerVolume() [1/2]

◆ CrossSectionPerVolume() [2/2]

◆ CSDARangeTable()

◆ CurrentMaterialCutsCoupleIndex()

◆ DEDXTable()

◆ DEDXTableForSubsec()

◆ DEDXunRestrictedTable()

◆ EmModel()

◆ FluctModel()

◆ GetContinuousStepLimit()

◆ GetCSDARange()

◆ GetCurrentElement()

◆ GetDEDX() [1/2]

◆ GetDEDX() [2/2]

◆ GetDEDXDispersion()

◆ GetDEDXForSubsec()

◆ GetKineticEnergy()

◆ GetLambda() [1/2]

◆ GetLambda() [2/2]

◆ GetMeanFreePath()

◆ GetModelByIndex()

◆ GetRange() [1/2]

◆ GetRange() [2/2]

◆ GetRangeForLoss() [1/2]

◆ GetRangeForLoss() [2/2]

◆ InitialiseEnergyLossProcess()

◆ InverseRangeTable()

◆ IonisationTable()

◆ IonisationTableForSubsec()

◆ IsApplicable()

◆ IsIntegral()

◆ IsIonisationProcess()

◆ LambdaPhysicsVector()

◆ LambdaTable()

◆ MaxKinEnergy()

◆ MeanFreePath()

◆ MinKinEnergy()

◆ MinPrimaryEnergy()

◆ NumberOfModels()

◆ NumberOfSubCutoffRegions()

◆ Particle()

◆ PostStepDoIt()

◆ PostStepGetPhysicalInteractionLength()

◆ PreparePhysicsTable()

◆ PrintInfo()

◆ ProcessDescription()

◆ RangeTableForLoss()

◆ RetrievePhysicsTable()

◆ SampleSubCutSecondaries()

◆ SecondaryParticle()

◆ SecondaryRangeTable()

◆ SelectModel()

◆ SelectModelForMaterial()

◆ SetBaseParticle()

◆ SetCrossSectionBiasingFactor()

◆ SetCSDARangeTable()

◆ SetDEDXBinning()

◆ SetDEDXTable()