#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)=0

virtual void	PrintInfo ()=0

void	PreparePhysicsTable (const G4ParticleDefinition &)

void	BuildPhysicsTable (const G4ParticleDefinition &)

G4PhysicsTable *	BuildDEDXTable (G4EmTableType tType=fRestricted)

G4PhysicsTable *	BuildLambdaTable (G4EmTableType tType=fRestricted)

void	PrintInfoDefinition ()

void	StartTracking (G4Track *)

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

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

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

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

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

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

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

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

G4double	CrossSectionPerVolume (G4double kineticEnergy, const G4MaterialCutsCouple *couple)

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=0)

void	UpdateEmModel (const G4String &, G4double, G4double)

void	SetEmModel (G4VEmModel *, G4int index=1)

G4VEmModel *	EmModel (G4int index=1)

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

G4int	NumberOfModels ()

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=0)

void	SetCrossSectionBiasingFactor (G4double f, G4bool flag=true)

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

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

void	AddCollaborativeProcess (G4VEnergyLossProcess *)

void	SetLossFluctuations (G4bool val)

void	SetRandomStep (G4bool val)

void	SetIntegral (G4bool val)

G4bool	IsIntegral () const

void	SetIonisation (G4bool val)

G4bool	IsIonisationProcess () const

void	SetLinearLossLimit (G4double val)

void	SetMinSubRange (G4double val)

void	SetLambdaFactor (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	SetLambdaBinning (G4int nbins)

void	SetDEDXBinningForCSDARange (G4int nbins)

void	SetMinKinEnergy (G4double e)

G4double	MinKinEnergy () const

void	SetMaxKinEnergy (G4double e)

G4double	MaxKinEnergy () const

void	SetMaxKinEnergyForCSDARange (G4double e)

G4double	CrossSectionBiasingFactor () const

G4double	GetDEDX (G4double &kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetDEDXForSubsec (G4double &kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetRange (G4double &kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetCSDARange (G4double &kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetRangeForLoss (G4double &kineticEnergy, const G4MaterialCutsCouple *)

G4double	GetKineticEnergy (G4double &range, const G4MaterialCutsCouple *)

G4double	GetLambda (G4double &kineticEnergy, const G4MaterialCutsCouple *)

G4bool	TablesAreBuilt () const

G4PhysicsTable *	DEDXTable () const

G4PhysicsTable *	DEDXTableForSubsec () const

G4PhysicsTable *	DEDXunRestrictedTable () const

G4PhysicsTable *	IonisationTable () const

G4PhysicsTable *	IonisationTableForSubsec () const

G4PhysicsTable *	CSDARangeTable () const

G4PhysicsTable *	RangeTableForLoss () const

G4PhysicsTable *	InverseRangeTable () const

G4PhysicsTable *	LambdaTable ()

G4PhysicsTable *	SubLambdaTable ()

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

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

G4int	operator== (const G4VProcess &right) const

G4int	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

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

Protected Member Functions
virtual void	InitialiseEnergyLossProcess (const G4ParticleDefinition , const G4ParticleDefinition )=0

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

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

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

G4PhysicsVector *	LambdaPhysicsVector (const G4MaterialCutsCouple *, G4double cut)

size_t	CurrentMaterialCutsCoupleIndex () const

G4double	GetCurrentRange () 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 previousStepSize)

void	ClearNumberOfInteractionLengthLeft ()

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

Detailed Description

Definition at line 121 of file G4VEnergyLossProcess.hh.

Constructor & Destructor Documentation

◆ G4VEnergyLossProcess()

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

Definition at line 157 of file G4VEnergyLossProcess.cc.

                                          : 
  G4VContinuousDiscreteProcess(name, type),
  secondaryParticle(0),
  nSCoffRegions(0),
  idxSCoffRegions(0),
  nProcesses(0),
  theDEDXTable(0),
  theDEDXSubTable(0),
  theDEDXunRestrictedTable(0),
  theIonisationTable(0),
  theIonisationSubTable(0),
  theRangeTableForLoss(0),
  theCSDARangeTable(0),
  theSecondaryRangeTable(0),
  theInverseRangeTable(0),
  theLambdaTable(0),
  theSubLambdaTable(0),
  theDensityFactor(0),
  theDensityIdx(0),
  baseParticle(0),
  minSubRange(0.1),
  lossFluctuationFlag(true),
  rndmStepFlag(false),
  tablesAreBuilt(false),
  integral(true),
  isIon(false),
  isIonisation(true),
  useSubCutoff(false),
  useDeexcitation(false),
  particle(0),
  currentCouple(0),
  nWarnings(0),
  mfpKinEnergy(0.0)
{
  SetVerboseLevel(1);
 
  // low energy limit
  lowestKinEnergy  = 1.*eV;
 
  // Size of tables assuming spline
  minKinEnergy     = 0.1*keV;
  maxKinEnergy     = 10.0*TeV;
  nBins            = 77;
  maxKinEnergyCSDA = 1.0*GeV;
  nBinsCSDA        = 35;
 
  // default linear loss limit for spline
  linLossLimit  = 0.01;
 
  // default dRoverRange and finalRange
  SetStepFunction(0.2, 1.0*mm);
 
  // default lambda factor
  lambdaFactor  = 0.8;
 
  // cross section biasing
  biasFactor = 1.0;
 
  // particle types
  theElectron   = G4Electron::Electron();
  thePositron   = G4Positron::Positron();
  theGamma      = G4Gamma::Gamma();
  theGenericIon = 0;
 
  // run time objects
  pParticleChange = &fParticleChange;
  fParticleChange.SetSecondaryWeightByProcess(true);
  modelManager = new G4EmModelManager();
  safetyHelper = G4TransportationManager::GetTransportationManager()
    ->GetSafetyHelper();
  aGPILSelection = CandidateForSelection;
 
  // initialise model
  (G4LossTableManager::Instance())->Register(this);
  fluctModel = 0;
  atomDeexcitation = 0;
 
  biasManager  = 0;
  biasFlag     = false; 
  weightFlag   = false; 
 
  scTracks.reserve(5);
  secParticles.reserve(5);
}

◆ ~G4VEnergyLossProcess()

G4VEnergyLossProcess::~G4VEnergyLossProcess ( )

virtual

Definition at line 245 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess destruct " << GetProcessName()
       << "  " << this << "  " << baseParticle << G4endl;
  }
  Clean();
 
  if ( !baseParticle ) {
    if(theDEDXTable) {
      if(theIonisationTable == theDEDXTable) { theIonisationTable = 0; }
      theDEDXTable->clearAndDestroy();
      delete theDEDXTable;
      if(theDEDXSubTable) {
    if(theIonisationSubTable == theDEDXSubTable) 
      theIonisationSubTable = 0;
    theDEDXSubTable->clearAndDestroy();
        delete theDEDXSubTable;
      }
    }
    if(theIonisationTable) {
      theIonisationTable->clearAndDestroy(); 
      delete theIonisationTable;
    }
    if(theIonisationSubTable) {
      theIonisationSubTable->clearAndDestroy(); 
      delete theIonisationSubTable;
    }
    if(theDEDXunRestrictedTable && theCSDARangeTable) {
       theDEDXunRestrictedTable->clearAndDestroy();
       delete theDEDXunRestrictedTable;
    }
    if(theCSDARangeTable) {
      theCSDARangeTable->clearAndDestroy();
      delete theCSDARangeTable;
    }
    if(theRangeTableForLoss) {
      theRangeTableForLoss->clearAndDestroy();
      delete theRangeTableForLoss;
    }
    if(theInverseRangeTable) {
      theInverseRangeTable->clearAndDestroy();
      delete theInverseRangeTable;
    }
    if(theLambdaTable) {
      theLambdaTable->clearAndDestroy();
      delete theLambdaTable;
    }
    if(theSubLambdaTable) {
      theSubLambdaTable->clearAndDestroy();
      delete theSubLambdaTable;
    }
  }
 
  delete modelManager;
  delete biasManager;
  (G4LossTableManager::Instance())->DeRegister(this);
}

Member Function Documentation

◆ ActivateForcedInteraction()

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

Definition at line 2058 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 G4EmProcessOptions::ActivateForcedInteraction().

◆ ActivateSecondaryBiasing()

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

Definition at line 2078 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 G4EmProcessOptions::ActivateSecondaryBiasing().

◆ ActivateSubCutoff()

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

Definition at line 847 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) {
    for (G4int i=0; i<nSCoffRegions; ++i) {
      if (reg == scoffRegions[i]) {
        return;
      }
    }
  }
 
  // new region 
  if(val) {
    useSubCutoff = true;
    scoffRegions.push_back(reg);
    ++nSCoffRegions;
  } else {
    useSubCutoff = false;
  }
}

Referenced by G4LossTableManager::Register().

◆ AddCollaborativeProcess()

void G4VEnergyLossProcess::AddCollaborativeProcess ( G4VEnergyLossProcess * p )

Definition at line 1805 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 = `0`
	)

Definition at line 331 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
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1038 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()
           << 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) {
    eloss = preStepKinEnergy;
    if (useDeexcitation) {
      atomDeexcitation->AlongStepDeexcitation(scTracks, step, 
                          eloss, currentCoupleIndex);
      if(scTracks.size() > 0) { FillSecondariesAlongStep(eloss, weight); }
      if(eloss < 0.0) { eloss = 0.0; }
    }
    fParticleChange.SetProposedKineticEnergy(0.0);
    fParticleChange.ProposeLocalEnergyDeposit(eloss);
    return &fParticleChange;
  }
 
  // Short step
  eloss = GetDEDXForScaledEnergy(preStepScaledEnergy)*length;
 
  // Long step
  if(eloss > preStepKinEnergy*linLossLimit) {
 
    //G4double x = GetScaledRangeForScaledEnergy(preStepScaledEnergy) 
    //  - length/reduceFactor;
    G4double x = (fRange - length)/reduceFactor;
    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;
 
  // SubCutOff 
  if(useSubCutoff) {
    if(idxSCoffRegions[currentCoupleIndex]) {
 
      G4bool yes = false;
      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());
    }
 
        if(preSafety < rcut) { yes = true; }
 
    // Check PostPoint
    else {
      G4double postSafety = preSafety - length; 
      if(postSafety < rcut) {
        postSafety = 
          safetyHelper->ComputeSafety(step.GetPostStepPoint()->GetPosition());
        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(fluc && 
      (eloss + esec + lowestKinEnergy) < preStepKinEnergy) {
 
      G4double tmax = 
    std::min(currentModel->MaxSecondaryKinEnergy(dynParticle),cut);
      G4double emean = eloss;
      eloss = fluc->SampleFluctuations(currentMaterial,dynParticle,
                       tmax,length,emean);
      /*                            
      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(scTracks.size() > 0) { FillSecondariesAlongStep(eloss, weight); }
 
  // Energy balanse
  G4double finalT = preStepKinEnergy - eloss - esec;
  if (finalT <= lowestKinEnergy) {
    eloss += finalT;
    finalT = 0.0;
  } else if(isIon) {
    fParticleChange.SetProposedCharge(
      currentModel->GetParticleCharge(track.GetParticleDefinition(),
                      currentMaterial,finalT));
  }
 
  if(eloss < 0.0) { 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
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 903 of file G4VEnergyLossProcess.cc.

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

Referenced by ContinuousStepLimit().

◆ BaseParticle()

const G4ParticleDefinition * G4VEnergyLossProcess::BaseParticle ( ) const

inline

Definition at line 887 of file G4VEnergyLossProcess.hh.

{
  return baseParticle;
}

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

◆ BuildDEDXTable()

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

Definition at line 612 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::BuildDEDXTable() of type " << tType
       << " for " << GetProcessName()
           << " and particle " << particle->GetParticleName()
           << G4endl;
  }
  G4PhysicsTable* table = 0;
  G4double emax = maxKinEnergy;
  G4int bin = nBins;
 
  if(fTotal == tType) {
    emax  = maxKinEnergyCSDA;
    bin   = nBinsCSDA;
    table = theDEDXunRestrictedTable;
  } else if(fRestricted == tType) {
    table = theDEDXTable;
    if(theIonisationTable) 
      table = G4PhysicsTableHelper::PreparePhysicsTable(theIonisationTable); 
  } else if(fSubRestricted == tType) {    
    table = theDEDXSubTable;
    if(theIonisationSubTable) 
      table = G4PhysicsTableHelper::PreparePhysicsTable(theIonisationSubTable); 
  } 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 = (G4LossTableManager::Instance())->GetTableBuilder();
  G4bool splineFlag = (G4LossTableManager::Instance())->SplineFlag();
  G4PhysicsLogVector* aVector = 0;
  G4PhysicsLogVector* bVector = 0;
 
  for(size_t i=0; i<numOfCouples; ++i) {
 
    if(1 < verboseLevel) {
      G4cout << "G4VEnergyLossProcess::BuildDEDXVector 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);
      delete (*table)[i];
      if(!bVector) {
    aVector = new G4PhysicsLogVector(minKinEnergy, emax, bin);
        bVector = aVector;
      } else {
        aVector = new G4PhysicsLogVector(*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 703 of file G4VEnergyLossProcess.cc.

{
  G4PhysicsTable* table = 0;
 
  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 = (G4LossTableManager::Instance())->GetTableBuilder();
  G4bool splineFlag = (G4LossTableManager::Instance())->SplineFlag();
  G4PhysicsLogVector* aVector = 0;
  G4double scale = std::log(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];
      G4double emin = MinPrimaryEnergy(particle,couple->GetMaterial(),(*theCuts)[i]);
      if(0.0 >= emin) { emin = eV; }
      else if(maxKinEnergy <= emin) { emin = 0.5*maxKinEnergy; }
      G4int bin = G4int(nBins*std::log(maxKinEnergy/emin)/scale + 0.5);
      if(bin < 3) { bin = 3; }
      aVector = new G4PhysicsLogVector(emin, maxKinEnergy, bin);
      aVector->SetSpline(splineFlag);
 
      modelManager->FillLambdaVector(aVector, couple, true, 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 )

virtual

Reimplemented from G4VProcess.

Definition at line 554 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(!tablesAreBuilt) {
      G4LossTableManager::Instance()->BuildPhysicsTable(particle, this);
    }
    if(!baseParticle) {
      // 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")))
    { 
      particle = &part;
      PrintInfoDefinition(); 
    }
 
  // Added tracking cut to avoid tracking artifacts
  // identify deexcitation flag
  if(isIonisation) { 
    fParticleChange.SetLowEnergyLimit(lowestKinEnergy); 
    atomDeexcitation = G4LossTableManager::Instance()->AtomDeexcitation();
    if(atomDeexcitation) { 
      if(atomDeexcitation->IsPIXEActive()) { useDeexcitation = true; } 
    }
  }
 
  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 1763 of file G4VEnergyLossProcess.cc.

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

◆ CrossSectionBiasingFactor()

G4double G4VEnergyLossProcess::CrossSectionBiasingFactor ( ) const

inline

Definition at line 1048 of file G4VEnergyLossProcess.hh.

{
  return biasFactor;
}

◆ CrossSectionPerVolume()

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

Definition at line 1731 of file G4VEnergyLossProcess.cc.

{
  // Cross section per volume is calculated
  DefineMaterial(couple);
  G4double cross = 0.0;
  if(theLambdaTable) {
    cross = (*theDensityFactor)[currentCoupleIndex]*
      ((*theLambdaTable)[basedCoupleIndex])->Value(kineticEnergy);
  } else {
    SelectModel(kineticEnergy);
    cross = currentModel->CrossSectionPerVolume(currentMaterial,
                        particle, kineticEnergy,
                        (*theCuts)[currentCoupleIndex]);
  }
  if(cross < 0.0) { cross = 0.0; }
  return cross;
}

◆ CSDARangeTable()

G4PhysicsTable * G4VEnergyLossProcess::CSDARangeTable ( ) const

inline

Definition at line 1101 of file G4VEnergyLossProcess.hh.

{
  return theCSDARangeTable;
}

◆ CurrentMaterialCutsCoupleIndex()

size_t G4VEnergyLossProcess::CurrentMaterialCutsCoupleIndex ( ) const

inlineprotected

Definition at line 584 of file G4VEnergyLossProcess.hh.

{
  return currentCoupleIndex;
}

Referenced by G4ePolarizedIonisation::GetMeanFreePath(), and G4ePolarizedIonisation::PostStepGetPhysicalInteractionLength().

◆ DEDXTable()

G4PhysicsTable * G4VEnergyLossProcess::DEDXTable ( ) const

inline

Definition at line 1062 of file G4VEnergyLossProcess.hh.

{
  return theDEDXTable;
}

Referenced by G4EmCalculator::PrintDEDXTable().

◆ DEDXTableForSubsec()

G4PhysicsTable * G4VEnergyLossProcess::DEDXTableForSubsec ( ) const

inline

Definition at line 1069 of file G4VEnergyLossProcess.hh.

{
  return theDEDXSubTable;
}

◆ DEDXunRestrictedTable()

G4PhysicsTable * G4VEnergyLossProcess::DEDXunRestrictedTable ( ) const

inline

Definition at line 1076 of file G4VEnergyLossProcess.hh.

{
  return theDEDXunRestrictedTable;
}

◆ EmModel()

G4VEmModel * G4VEnergyLossProcess::EmModel ( G4int index = 1 )

Definition at line 358 of file G4VEnergyLossProcess.cc.

{
  G4VEmModel* p = 0;
  if(index >= 0 && index <  G4int(emModels.size())) { p = emModels[index]; }
  return p;
}

◆ FluctModel()

G4VEmFluctuationModel * G4VEnergyLossProcess::FluctModel ( )

inline

Definition at line 852 of file G4VEnergyLossProcess.hh.

{
  return fluctModel;
}

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

◆ GetContinuousStepLimit()

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

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 1785 of file G4VEnergyLossProcess.cc.

{
  return DBL_MAX;
}

◆ GetCSDARange()

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

inline

Definition at line 771 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  G4double x = DBL_MAX;
  if(theCSDARangeTable) {
    x = GetLimitScaledRangeForScaledEnergy(kineticEnergy*massRatio)*reduceFactor;
  }
  return x;
}

Referenced by G4LossTableManager::GetCSDARange().

◆ GetCurrentElement()

const G4Element * G4VEnergyLossProcess::GetCurrentElement ( ) const

Definition at line 2031 of file G4VEnergyLossProcess.cc.

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

◆ GetCurrentRange()

G4double G4VEnergyLossProcess::GetCurrentRange ( ) const

inlineprotected

Definition at line 591 of file G4VEnergyLossProcess.hh.

{
  return fRange;
}

◆ GetDEDX()

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

inline

Definition at line 732 of file G4VEnergyLossProcess.hh.

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

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

◆ GetDEDXDispersion()

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

Definition at line 1713 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 742 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 801 of file G4VEnergyLossProcess.hh.

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

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

◆ GetLambda()

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

inline

Definition at line 811 of file G4VEnergyLossProcess.hh.

{
  DefineMaterial(couple);
  G4double x = 0.0;
  if(theLambdaTable) { x = GetLambdaForScaledEnergy(kineticEnergy*massRatio); }
  return x;
}

◆ GetMeanFreePath()

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

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 1773 of file G4VEnergyLossProcess.cc.

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

Referenced by G4ePolarizedIonisation::GetMeanFreePath().

◆ GetModelByIndex()

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

Definition at line 367 of file G4VEnergyLossProcess.cc.

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

◆ GetRange()

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

inline

Definition at line 752 of file G4VEnergyLossProcess.hh.

{
  G4double x = fRange;
  if(couple != currentCouple || kineticEnergy != preStepKinEnergy) { 
    DefineMaterial(couple);
    if(theCSDARangeTable) {
      x = GetLimitScaledRangeForScaledEnergy(kineticEnergy*massRatio)
    * reduceFactor;
    } else if(theRangeTableForLoss) {
      x = GetScaledRangeForScaledEnergy(kineticEnergy*massRatio)*reduceFactor;
    }
  }
  return x;
}

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

◆ GetRangeForLoss()

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

inline

Definition at line 785 of file G4VEnergyLossProcess.hh.

{
  G4double x = fRange;
  if(couple != currentCouple || kineticEnergy != preStepKinEnergy) {
    DefineMaterial(couple);
    x = GetScaledRangeForScaledEnergy(kineticEnergy*massRatio)*reduceFactor;
  }
  //  G4cout << "Range from " << GetProcessName() 
  //         << "  e= " << kineticEnergy << " r= " << x << G4endl;
  return x;
}

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

◆ InitialiseEnergyLossProcess()

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

protectedpure virtual

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

Referenced by PreparePhysicsTable().

◆ InverseRangeTable()

G4PhysicsTable * G4VEnergyLossProcess::InverseRangeTable ( ) const

inline

Definition at line 1115 of file G4VEnergyLossProcess.hh.

{
  return theInverseRangeTable;
}

Referenced by G4EmCalculator::PrintInverseRangeTable().

◆ IonisationTable()

G4PhysicsTable * G4VEnergyLossProcess::IonisationTable ( ) const

inline

Definition at line 1083 of file G4VEnergyLossProcess.hh.

{
  G4PhysicsTable* t = theDEDXTable;
  if(theIonisationTable) { t = theIonisationTable; } 
  return t;
}

◆ IonisationTableForSubsec()

G4PhysicsTable * G4VEnergyLossProcess::IonisationTableForSubsec ( ) const

inline

Definition at line 1092 of file G4VEnergyLossProcess.hh.

{
  G4PhysicsTable* t = theDEDXSubTable;
  if(theIonisationSubTable) { t = theIonisationSubTable; } 
  return t;
}

◆ IsApplicable()

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

pure virtual

Reimplemented from G4VProcess.

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

◆ IsIntegral()

G4bool G4VEnergyLossProcess::IsIntegral ( ) const

inline

Definition at line 922 of file G4VEnergyLossProcess.hh.

{
  return integral;
}

◆ IsIonisationProcess()

G4bool G4VEnergyLossProcess::IsIonisationProcess ( ) const

inline

Definition at line 938 of file G4VEnergyLossProcess.hh.

{
  return isIonisation;
}

Referenced by G4LossTableManager::BuildPhysicsTable().

◆ LambdaPhysicsVector()

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

protected

Definition at line 1795 of file G4VEnergyLossProcess.cc.

{
  G4PhysicsVector* v = new G4PhysicsLogVector(minKinEnergy, maxKinEnergy, nBins);
  v->SetSpline((G4LossTableManager::Instance())->SplineFlag());
  return v;
}

Referenced by G4ePolarizedIonisation::BuildPhysicsTable().

◆ LambdaTable()

G4PhysicsTable * G4VEnergyLossProcess::LambdaTable ( )

inline

Definition at line 1122 of file G4VEnergyLossProcess.hh.

{
  return theLambdaTable;
}

◆ MaxKinEnergy()

G4double G4VEnergyLossProcess::MaxKinEnergy ( ) const

inline

Definition at line 1034 of file G4VEnergyLossProcess.hh.

{
  return maxKinEnergy;
}

◆ MeanFreePath()

G4double G4VEnergyLossProcess::MeanFreePath ( const G4Track & track )

Definition at line 1752 of file G4VEnergyLossProcess.cc.

{
  DefineMaterial(track.GetMaterialCutsCouple());
  preStepLambda = GetLambdaForScaledEnergy(track.GetKineticEnergy()*massRatio);
  G4double x = DBL_MAX;
  if(0.0 < preStepLambda) { x = 1.0/preStepLambda; }
  return x;
}

Referenced by GetMeanFreePath().

◆ MinKinEnergy()

G4double G4VEnergyLossProcess::MinKinEnergy ( ) const

inline

Definition at line 1019 of file G4VEnergyLossProcess.hh.

{
  return minKinEnergy;
}

◆ MinPrimaryEnergy()

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

protectedvirtual

Reimplemented in G4ePolarizedIonisation, G4eIonisation, G4hBremsstrahlung, G4hhIonisation, G4hPairProduction, G4MuBremsstrahlung, G4MuIonisation, G4MuPairProduction, G4alphaIonisation, G4hIonisation, and G4ionIonisation.

Definition at line 322 of file G4VEnergyLossProcess.cc.

{
  return cut;
}

Referenced by BuildLambdaTable().

◆ NumberOfModels()

G4int G4VEnergyLossProcess::NumberOfModels ( )

Definition at line 374 of file G4VEnergyLossProcess.cc.

{
  return modelManager->NumberOfModels();
}

◆ NumberOfSubCutoffRegions()

G4int G4VEnergyLossProcess::NumberOfSubCutoffRegions ( ) const

inline

Definition at line 984 of file G4VEnergyLossProcess.hh.

{
  return nSCoffRegions;
}

◆ Particle()

const G4ParticleDefinition * G4VEnergyLossProcess::Particle ( ) const

inline

Definition at line 880 of file G4VEnergyLossProcess.hh.

{
  return particle;
}

Referenced by G4LossTableManager::BuildPhysicsTable().

◆ PostStepDoIt()

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1410 of file G4VEnergyLossProcess.cc.

{
  // In all cases clear number of interaction lengths
  theNumberOfInteractionLengthLeft = -1.0;
  mfpKinEnergy = DBL_MAX; 
 
  fParticleChange.InitializeForPostStep(track);
  G4double finalT = track.GetKineticEnergy();
  if(finalT <= lowestKinEnergy) { return &fParticleChange; }
 
  G4double postStepScaledEnergy = finalT*massRatio;
 
  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;
    }
  }
 
  // Integral approach
  if (integral) {
    G4double lx = GetLambdaForScaledEnergy(postStepScaledEnergy);
    /*
    if(preStepLambda<lx && 1 < verboseLevel && nWarnings<200) {
      G4cout << "WARNING: for " << particle->GetParticleName()
             << " and " << GetProcessName()
             << " E(MeV)= " << finalT/MeV
             << " preLambda= " << preStepLambda 
         << " < " << lx << " (postLambda) "
         << G4endl;
      ++nWarnings;
    }
    */
    if(lx <= 0.0) {
      return &fParticleChange;
    } else if(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);
  }
 
  const G4DynamicParticle* dynParticle = track.GetDynamicParticle();
  G4double tcut = (*theCuts)[currentCoupleIndex];
 
  // sample secondaries
  secParticles.clear();
  //G4cout << "Energy of primary: " << dynParticle->GetKineticEnergy()/MeV<<G4endl;
  currentModel->SampleSecondaries(&secParticles, currentCouple, dynParticle, tcut);
 
  // 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);
 
    for (G4int i=0; i<num; ++i) {
      if(secParticles[i]) {
    G4Track* t = new G4Track(secParticles[i], track.GetGlobalTime(), 
                 track.GetPosition());
    t->SetTouchableHandle(track.GetTouchableHandle());
    t->SetWeight(weight); 
    //G4cout << "Secondary(post step) has weight " << t->GetWeight() 
    //<< ", kenergy " << t->GetKineticEnergy()/MeV << " MeV" <<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
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 933 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
  /*
  if(!theDensityFactor || !theDensityIdx) {
    G4cout << "G4VEnergyLossProcess::PostStepGetPhysicalInteractionLength 1: "
       <<  theDensityFactor << "  " << theDensityIdx
       << G4endl;
    G4cout << track.GetDefinition()->GetParticleName() 
       << " e(MeV)= " << track.GetKineticEnergy()
       << " mat " << track.GetMaterialCutsCouple()->GetMaterial()->GetName()
       << G4endl;
  }
  */
  DefineMaterial(track.GetMaterialCutsCouple());
  preStepKinEnergy    = track.GetKineticEnergy();
  preStepScaledEnergy = preStepKinEnergy*massRatio;
  SelectModel(preStepScaledEnergy);
 
  if(!currentModel->IsActive(preStepScaledEnergy)) { return x; }
 
  // change effective charge of an ion on fly
  if(isIon) {
    G4double q2 = currentModel->ChargeSquareRatio(track);
    if(q2 != chargeSqRatio) {
      chargeSqRatio = q2;
      fFactor = q2*biasFactor*(*theDensityFactor)[currentCoupleIndex];
      reduceFactor = 1.0/(fFactor*massRatio);
    }
  }
  //  if(particle->GetPDGMass() > 0.9*GeV)
  //G4cout << "q2= " << chargeSqRatio << " massRatio= " << massRatio << G4endl; 
  // initialisation for sampling of the interaction length 
  //if(previousStepSize <= 0.0) { theNumberOfInteractionLengthLeft = -1.0; }
  //if(theNumberOfInteractionLengthLeft < 0.0) { mfpKinEnergy = DBL_MAX; }
 
  // forced biasing only for primary particles
  if(biasManager) {
    if(0 == track.GetParentID()) {
      if(biasFlag && biasManager->ForcedInteractionRegion(currentCoupleIndex)) {
        return biasManager->GetStepLimit(currentCoupleIndex, previousStepSize);
      }
    }
  }
 
  // compute mean free path
  if(preStepScaledEnergy < mfpKinEnergy) {
    if (integral) { ComputeLambdaForScaledEnergy(preStepScaledEnergy); }
    else  { preStepLambda = GetLambdaForScaledEnergy(preStepScaledEnergy); }
 
    // 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)
      ResetNumberOfInteractionLengthLeft();
 
    } else if(currentInteractionLength < DBL_MAX) {
 
      // subtract NumberOfInteractionLengthLeft using previous step
      theNumberOfInteractionLengthLeft -= previousStepSize/currentInteractionLength;
      //    SubtractNumberOfInteractionLengthLeft(previousStepSize);
      if(theNumberOfInteractionLengthLeft < 0.) {
    theNumberOfInteractionLengthLeft = 0.0;
    //theNumberOfInteractionLengthLeft = perMillion;
      }
    }
 
    // 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 
         <<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 )

virtual

Reimplemented from G4VProcess.

Definition at line 382 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::PreparePhysicsTable for "
           << GetProcessName() << " for " << part.GetParticleName() 
       << "  " << this << G4endl;
  }
 
  currentCouple = 0;
  preStepLambda = 0.0;
  mfpKinEnergy  = DBL_MAX;
  fRange        = DBL_MAX;
  preStepKinEnergy = 0.0;
  chargeSqRatio = 1.0;
  massRatio = 1.0;
  reduceFactor = 1.0;
  fFactor = 1.0;
 
  G4LossTableManager* lManager = G4LossTableManager::Instance();
 
  // 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") {
 
      theGenericIon = G4GenericIon::GenericIon();
      isIon = true; 
      // process is shared between all ions inheriting G4GenericIon
      // for all excluding He3 and alpha
      if(pname != "He3" && pname != "alpha") { particle = theGenericIon; }
    }
  }
 
  if( particle != &part ) {
    if(isIon) {
      lManager->RegisterIon(&part, this);
    } else { 
      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);
  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);
 
  // Tables preparation
  if (!baseParticle) {
    
    theDEDXTable = G4PhysicsTableHelper::PreparePhysicsTable(theDEDXTable);
    bld->InitialiseBaseMaterials(theDEDXTable);
 
    if (lManager->BuildCSDARange()) {
      theDEDXunRestrictedTable = 
    G4PhysicsTableHelper::PreparePhysicsTable(theDEDXunRestrictedTable);
      theCSDARangeTable = 
    G4PhysicsTableHelper::PreparePhysicsTable(theCSDARangeTable);
    }
 
    theLambdaTable = G4PhysicsTableHelper::PreparePhysicsTable(theLambdaTable);
    bld->InitialiseBaseMaterials(theLambdaTable);  
 
    if(isIonisation) {
      theRangeTableForLoss = 
    G4PhysicsTableHelper::PreparePhysicsTable(theRangeTableForLoss);
      theInverseRangeTable = 
    G4PhysicsTableHelper::PreparePhysicsTable(theInverseRangeTable);  
    }
 
    if (nSCoffRegions) {
      theDEDXSubTable = 
    G4PhysicsTableHelper::PreparePhysicsTable(theDEDXSubTable);
      theSubLambdaTable = 
    G4PhysicsTableHelper::PreparePhysicsTable(theSubLambdaTable);
    }
  }
 
  theDensityFactor = bld->GetDensityFactors();
  theDensityIdx = bld->GetCoupleIndexes();
 
  // forced biasing
  if(biasManager) { 
    biasManager->Initialise(part,GetProcessName(),verboseLevel); 
    biasFlag = false; 
  }
 
  G4double initialCharge = particle->GetPDGCharge();
  G4double initialMass   = particle->GetPDGMass();
 
  if (baseParticle) {
    massRatio = (baseParticle->GetPDGMass())/initialMass;
    G4double q = initialCharge/baseParticle->GetPDGCharge();
    chargeSqRatio = q*q;
    if(chargeSqRatio > 0.0) { reduceFactor = 1.0/(chargeSqRatio*massRatio); }
  }
 
  // initialisation of models
  G4int nmod = modelManager->NumberOfModels();
  for(G4int i=0; i<nmod; ++i) {
    G4VEmModel* mod = modelManager->GetModel(i);
    if(mod->HighEnergyLimit() > maxKinEnergy) {
      mod->SetHighEnergyLimit(maxKinEnergy);
    }
  }
 
  theCuts = modelManager->Initialise(particle, secondaryParticle, 
                     minSubRange, verboseLevel);
 
  // Sub Cutoff 
  if (nSCoffRegions>0) {
    theSubCuts = modelManager->SubCutoff();
 
    if(nSCoffRegions>0) { idxSCoffRegions = new G4bool[n]; }
    for (size_t j=0; j<n; ++j) {
 
      const G4MaterialCutsCouple* couple = 
    theCoupleTable->GetMaterialCutsCouple(j);
      const G4ProductionCuts* pcuts = couple->GetProductionCuts();
      
      if(nSCoffRegions>0) {
    G4bool reg = false;
    for(G4int i=0; i<nSCoffRegions; ++i) {
      if( pcuts == scoffRegions[i]->GetProductionCuts()) { reg = true; }
    }
    idxSCoffRegions[j] = reg;
      }
    }
  }
 
  if(1 < verboseLevel) {
    G4cout << "G4VEnergyLossProcess::Initialise() 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 ( )

pure virtual

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

Referenced by PrintInfoDefinition().

◆ PrintInfoDefinition()

void G4VEnergyLossProcess::PrintInfoDefinition ( )

Definition at line 772 of file G4VEnergyLossProcess.cc.

{
  if(0 < verboseLevel) {
    G4cout << G4endl << GetProcessName() << ":   for  "
           << particle->GetParticleName()
       << "    SubType= " << GetProcessSubType() 
           << G4endl
           << "      dE/dx and range tables from "
       << G4BestUnit(minKinEnergy,"Energy")
           << " to " << G4BestUnit(maxKinEnergy,"Energy")
           << " in " << nBins << " bins" << G4endl
           << "      Lambda tables from threshold to "
           << G4BestUnit(maxKinEnergy,"Energy")
           << " in " << nBins << " bins, spline: " 
       << (G4LossTableManager::Instance())->SplineFlag()
           << G4endl;
    if(theRangeTableForLoss && isIonisation) {
      G4cout << "      finalRange(mm)= " << finalRange/mm
             << ", dRoverRange= " << dRoverRange
             << ", integral: " << integral
             << ", fluct: " << lossFluctuationFlag
         << ", linLossLimit= " << linLossLimit
             << G4endl;
    }
    PrintInfo();
    modelManager->DumpModelList(verboseLevel);
    if(theCSDARangeTable && isIonisation) {
      G4cout << "      CSDA range table up"
             << " to " << G4BestUnit(maxKinEnergyCSDA,"Energy")
             << " in " << nBinsCSDA << " bins" << G4endl;
    }
    if(nSCoffRegions>0 && isIonisation) {
      G4cout << "      Subcutoff sampling in " << nSCoffRegions 
         << " regions" << G4endl;
    }
    if(2 < verboseLevel) {
      G4cout << "      DEDXTable address= " << theDEDXTable << G4endl;
      if(theDEDXTable && isIonisation) G4cout << (*theDEDXTable) << G4endl;
      G4cout << "non restricted DEDXTable address= " 
         << theDEDXunRestrictedTable << G4endl;
      if(theDEDXunRestrictedTable && isIonisation) {
           G4cout << (*theDEDXunRestrictedTable) << G4endl;
      }
      if(theDEDXSubTable && isIonisation) {
    G4cout << (*theDEDXSubTable) << G4endl;
      }
      G4cout << "      CSDARangeTable address= " << theCSDARangeTable 
         << G4endl;
      if(theCSDARangeTable && isIonisation) {
    G4cout << (*theCSDARangeTable) << G4endl;
      }
      G4cout << "      RangeTableForLoss address= " << theRangeTableForLoss 
         << G4endl;
      if(theRangeTableForLoss && isIonisation) {
             G4cout << (*theRangeTableForLoss) << G4endl;
      }
      G4cout << "      InverseRangeTable address= " << theInverseRangeTable 
         << G4endl;
      if(theInverseRangeTable && isIonisation) {
             G4cout << (*theInverseRangeTable) << G4endl;
      }
      G4cout << "      LambdaTable address= " << theLambdaTable << G4endl;
      if(theLambdaTable && isIonisation) {
    G4cout << (*theLambdaTable) << G4endl;
      }
      G4cout << "      SubLambdaTable address= " << theSubLambdaTable << G4endl;
      if(theSubLambdaTable && isIonisation) {
    G4cout << (*theSubLambdaTable) << G4endl;
      }
    }
  }
}

Referenced by BuildPhysicsTable().

◆ RangeTableForLoss()

G4PhysicsTable * G4VEnergyLossProcess::RangeTableForLoss ( ) const

inline

Definition at line 1108 of file G4VEnergyLossProcess.hh.

{
  return theRangeTableForLoss;
}

Referenced by G4EmCalculator::PrintRangeTable().

◆ RetrievePhysicsTable()

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

virtual

Reimplemented from G4VProcess.

Definition at line 1596 of file G4VEnergyLossProcess.cc.

{
  G4bool res = true;
  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 1318 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));
  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 dt = length/preStepPoint->GetVelocity();
  G4double fragment = 0.0;
 
  do {
    G4double del = -std::log(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) {
 
      G4bool addSec = true;
      /*
      // do not track very low-energy delta-electrons
      if(theSecondaryRangeTable && (*it)->GetParticleDefinition() == theElectron) {
    G4double ekin = (*it)->GetKineticEnergy();
    G4double rg = ((*theSecondaryRangeTable)[idx]->Value(ekin));
    //          if(rg < currentMinSafety) {
    if(rg < safetyHelper->ComputeSafety(r)) {
      extraEdep += ekin;
      delete (*it);
      addSec = false;
    }
      }
      */
      if(addSec) {
    G4Track* t = new G4Track((*it), pretime + fragment*dt, r);
    t->SetTouchableHandle(track->GetTouchableHandle());
    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;
    */
      }
    }
  } while (fragment <= 1.0);
  return esec;
} 

Referenced by AlongStepDoIt().

◆ SecondaryParticle()

const G4ParticleDefinition * G4VEnergyLossProcess::SecondaryParticle ( ) const

inline

Definition at line 894 of file G4VEnergyLossProcess.hh.

{
  return secondaryParticle;
}

◆ SelectModel()

void G4VEnergyLossProcess::SelectModel ( G4double kinEnergy )

inlineprotected

Definition at line 598 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 606 of file G4VEnergyLossProcess.hh.

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

Referenced by G4ContinuousGainOfEnergy::GetContinuousStepLimit().

◆ SetBaseParticle()

void G4VEnergyLossProcess::SetBaseParticle ( const G4ParticleDefinition * p )

inline

Definition at line 873 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 2040 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 G4EmProcessOptions::SetProcessBiasingFactor().

◆ SetCSDARangeTable()

void G4VEnergyLossProcess::SetCSDARangeTable ( G4PhysicsTable * pRange )

Definition at line 1892 of file G4VEnergyLossProcess.cc.

{
  if(theCSDARangeTable != p) { 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); }
      else {
    pv = (*p)[(*theDensityIdx)[i]];
        if(pv) { rmax = pv->Value(emax)/(*theDensityFactor)[i]; }
      }
      theRangeAtMaxEnergy[i] = rmax;
      //G4cout << "i= " << i << " Emax(MeV)= " << emax/MeV << " Rmax= " 
      //<< rmax<< G4endl;
    }
  }
}

◆ SetDEDXBinning()

void G4VEnergyLossProcess::SetDEDXBinning ( G4int nbins )

inline

Definition at line 991 of file G4VEnergyLossProcess.hh.

{
  nBins = nbins;
}

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

◆ SetDEDXBinningForCSDARange()

void G4VEnergyLossProcess::SetDEDXBinningForCSDARange ( G4int nbins )

inline

Definition at line 1005 of file G4VEnergyLossProcess.hh.

{
  nBinsCSDA = nbins;
}

◆ SetDEDXTable()

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

Definition at line 1831 of file G4VEnergyLossProcess.cc.

{
  if(fTotal == tType && theDEDXunRestrictedTable != p && !baseParticle) {
    if(theDEDXunRestrictedTable) {
      theDEDXunRestrictedTable->clearAndDestroy();
      delete theDEDXunRestrictedTable;
    } 
    theDEDXunRestrictedTable = p;
    if(p) {
      size_t n = p->length();
      G4PhysicsVector* pv = (*p)[0];
      G4double emax = maxKinEnergyCSDA;
 
      for (size_t i=0; i<n; ++i) {
    G4double dedx = 0.0; 
    pv = (*p)[i];
    if(pv) { dedx = pv->Value(emax); }
    else {
      pv = (*p)[(*theDensityIdx)[i]];
      if(pv) { dedx = pv->Value(emax)*(*theDensityFactor)[i]; }
    }
    theDEDXAtMaxEnergy[i] = dedx;
    //G4cout << "i= " << i << " emax(MeV)= " << emax/MeV<< " dedx= " 
    //<< dedx << G4endl;
      }
    }
 
  } else if(fRestricted == tType && theDEDXTable != p) {
    //G4cout << "G4VEnergyLossProcess::SetDEDXTable " << particle->GetParticleName()
    //     << " old table " << theDEDXTable << " new table " << p 
    //     << " ion " << theIonisationTable << " bp " << baseParticle << G4endl;
    if(theDEDXTable && !baseParticle) {
      if(theDEDXTable == theIonisationTable) { theIonisationTable = 0; }
      theDEDXTable->clearAndDestroy();
      delete theDEDXTable;
    }
    theDEDXTable = p;
  } else if(fSubRestricted == tType && theDEDXSubTable != p) {    
    if(theDEDXSubTable && !baseParticle) {
      if(theDEDXSubTable == theIonisationSubTable) { theIonisationSubTable = 0; }
      theDEDXSubTable->clearAndDestroy();
      delete theDEDXSubTable;
    }
    theDEDXSubTable = p;
  } else if(fIsIonisation == tType && theIonisationTable != p) {    
    if(theIonisationTable && theIonisationTable != theDEDXTable && !baseParticle) {
      theIonisationTable->clearAndDestroy();
      delete theIonisationTable;
    }
    theIonisationTable = p;
  } else if(fIsSubIonisation == tType && theIonisationSubTable != p) {    
    if(theIonisationSubTable && theIonisationSubTable != theDEDXSubTable && !baseParticle) {
      theIonisationSubTable->clearAndDestroy();
      delete theIonisationSubTable;
    }
    theIonisationSubTable = p;
  }
}

◆ SetDynamicMassCharge()

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

inline

Definition at line 630 of file G4VEnergyLossProcess.hh.

{
  massRatio     = massratio;
  fFactor      *= charge2ratio/chargeSqRatio;
  chargeSqRatio = charge2ratio;
  reduceFactor  = 1.0/(fFactor*massRatio);
}

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

◆ SetEmModel()

void G4VEnergyLossProcess::SetEmModel	(	G4VEmModel *	p,
		G4int	index = `1`
	)

Definition at line 349 of file G4VEnergyLossProcess.cc.

{
  G4int n = emModels.size();
  if(index >= n) { for(G4int i=n; i<=index; ++i) {emModels.push_back(0);} }
  emModels[index] = p;
}

◆ SetFluctModel()

void G4VEnergyLossProcess::SetFluctModel ( G4VEmFluctuationModel * p )

inline

Definition at line 845 of file G4VEnergyLossProcess.hh.

{
  fluctModel = p;
}

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

◆ SetIntegral()

void G4VEnergyLossProcess::SetIntegral ( G4bool val )

inline

Definition at line 915 of file G4VEnergyLossProcess.hh.

{
  integral = val;
}

Referenced by G4LossTableManager::Register().

◆ SetInverseRangeTable()

void G4VEnergyLossProcess::SetInverseRangeTable ( G4PhysicsTable * p )

Definition at line 1946 of file G4VEnergyLossProcess.cc.

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

◆ SetIonisation()

void G4VEnergyLossProcess::SetIonisation ( G4bool val )

inline

Definition at line 929 of file G4VEnergyLossProcess.hh.

{
  isIonisation = val;
  if(val) { aGPILSelection = CandidateForSelection; }
  else    { aGPILSelection = NotCandidateForSelection; }
}

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

◆ SetLambdaBinning()

void G4VEnergyLossProcess::SetLambdaBinning ( G4int nbins )

inline

Definition at line 998 of file G4VEnergyLossProcess.hh.

{
  nBins = nbins;
}

◆ SetLambdaFactor()

void G4VEnergyLossProcess::SetLambdaFactor ( G4double val )

inline

Definition at line 959 of file G4VEnergyLossProcess.hh.

{
  if(val > 0.0 && val <= 1.0) { lambdaFactor = val; }
}

Referenced by G4EmProcessOptions::SetLambdaFactor().

◆ SetLambdaTable()

void G4VEnergyLossProcess::SetLambdaTable ( G4PhysicsTable * p )

Definition at line 1960 of file G4VEnergyLossProcess.cc.

{
  if(1 < verboseLevel) {
    G4cout << "### Set Lambda table " << p 
       << " for " << particle->GetParticleName()
           << " and process " << GetProcessName() << G4endl;
  }
  if(theLambdaTable != p) { theLambdaTable = p; }
  tablesAreBuilt = true;
 
  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];
      }
    }
  }
}

◆ SetLinearLossLimit()

void G4VEnergyLossProcess::SetLinearLossLimit ( G4double val )

inline

Definition at line 945 of file G4VEnergyLossProcess.hh.

{
  linLossLimit = val;
}

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

◆ SetLossFluctuations()

void G4VEnergyLossProcess::SetLossFluctuations ( G4bool val )

inline

Definition at line 901 of file G4VEnergyLossProcess.hh.

{
  lossFluctuationFlag = val;
}

Referenced by G4LossTableManager::Register().

◆ SetLowestEnergyLimit()

void G4VEnergyLossProcess::SetLowestEnergyLimit ( G4double val )

inline

Definition at line 977 of file G4VEnergyLossProcess.hh.

{
  lowestKinEnergy = val;
}

◆ SetMaxKinEnergy()

void G4VEnergyLossProcess::SetMaxKinEnergy ( G4double e )

inline

Definition at line 1026 of file G4VEnergyLossProcess.hh.

{
  maxKinEnergy = e;
  if(e < maxKinEnergyCSDA) { maxKinEnergyCSDA = e; }
}

Referenced by G4hhIonisation::InitialiseEnergyLossProcess(), G4mplIonisation::InitialiseEnergyLossProcess(), and G4LossTableManager::Register().

◆ SetMaxKinEnergyForCSDARange()

void G4VEnergyLossProcess::SetMaxKinEnergyForCSDARange ( G4double e )

inline

Definition at line 1041 of file G4VEnergyLossProcess.hh.

{
  maxKinEnergyCSDA = e;
}

◆ SetMinKinEnergy()

void G4VEnergyLossProcess::SetMinKinEnergy ( G4double e )

inline

Definition at line 1012 of file G4VEnergyLossProcess.hh.

{
  minKinEnergy = e;
}

Referenced by G4hhIonisation::InitialiseEnergyLossProcess(), G4mplIonisation::InitialiseEnergyLossProcess(), and G4LossTableManager::Register().

◆ SetMinSubRange()

void G4VEnergyLossProcess::SetMinSubRange ( G4double val )

inline

Definition at line 952 of file G4VEnergyLossProcess.hh.

{
  minSubRange = val;
}

◆ SetParticle()

void G4VEnergyLossProcess::SetParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 859 of file G4VEnergyLossProcess.hh.

{
  particle = p;
}

◆ SetRandomStep()

void G4VEnergyLossProcess::SetRandomStep ( G4bool val )

inline

Definition at line 908 of file G4VEnergyLossProcess.hh.

{
  rndmStepFlag = val;
}

Referenced by G4LossTableManager::Register().

◆ SetRangeTableForLoss()

void G4VEnergyLossProcess::SetRangeTableForLoss ( G4PhysicsTable * p )

Definition at line 1918 of file G4VEnergyLossProcess.cc.

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

◆ SetSecondaryParticle()

void G4VEnergyLossProcess::SetSecondaryParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 866 of file G4VEnergyLossProcess.hh.

{
  secondaryParticle = p;
}

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

◆ SetSecondaryRangeTable()

void G4VEnergyLossProcess::SetSecondaryRangeTable ( G4PhysicsTable * p )

Definition at line 1932 of file G4VEnergyLossProcess.cc.

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

◆ SetStepFunction()

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

inline

Definition at line 966 of file G4VEnergyLossProcess.hh.

{
  dRoverRange = v1;
  finalRange = v2;
  if (dRoverRange > 0.999) { dRoverRange = 1.0; }
  currentCouple = 0;
  mfpKinEnergy  = DBL_MAX;
}

◆ SetSubLambdaTable()

void G4VEnergyLossProcess::SetSubLambdaTable ( G4PhysicsTable * p )

Definition at line 2017 of file G4VEnergyLossProcess.cc.

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

◆ StartTracking()

void G4VEnergyLossProcess::StartTracking ( G4Track * track )

virtual

Reimplemented from G4VProcess.

Definition at line 874 of file G4VEnergyLossProcess.cc.

{
  // reset parameters for the new track
  theNumberOfInteractionLengthLeft = -1.0;
  mfpKinEnergy = DBL_MAX; 
 
  // reset ion
  if(isIon) {
    chargeSqRatio = 0.5;
 
    G4double newmass = track->GetDefinition()->GetPDGMass();
    if(baseParticle) {
      massRatio = baseParticle->GetPDGMass()/newmass;
    } else {
      massRatio = proton_mass_c2/newmass;
    }
  }  
  // 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`
	)

virtual

Reimplemented from G4VProcess.

Definition at line 1536 of file G4VEnergyLossProcess.cc.

{
  G4bool res = true;
  if ( 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(0 < 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;
}

◆ SubLambdaTable()

G4PhysicsTable * G4VEnergyLossProcess::SubLambdaTable ( )

inline

Definition at line 1129 of file G4VEnergyLossProcess.hh.

{
  return theSubLambdaTable;
}

◆ TablesAreBuilt()

G4bool G4VEnergyLossProcess::TablesAreBuilt ( ) const

inline

Definition at line 1055 of file G4VEnergyLossProcess.hh.

{
  return  tablesAreBuilt;
}

◆ UpdateEmModel()

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

Definition at line 341 of file G4VEnergyLossProcess.cc.

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

Member Data Documentation

◆ fParticleChange

G4ParticleChangeForLoss G4VEnergyLossProcess::fParticleChange

protected

Definition at line 548 of file G4VEnergyLossProcess.hh.

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

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

Protected Attributes
G4ParticleChangeForLoss	fParticleChange

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager

G4VParticleChange *	pParticleChange

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft

G4double	currentInteractionLength

G4double	theInitialNumberOfInteractionLength

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType

G4int	theProcessSubType

G4double	thePILfactor

G4bool	enableAtRestDoIt

G4bool	enableAlongStepDoIt

G4bool	enablePostStepDoIt

G4int	verboseLevel

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

◆ CSDARangeTable()

◆ CurrentMaterialCutsCoupleIndex()

◆ DEDXTable()

◆ DEDXTableForSubsec()

◆ DEDXunRestrictedTable()

◆ EmModel()

◆ FluctModel()

◆ GetContinuousStepLimit()

◆ GetCSDARange()

◆ GetCurrentElement()

◆ GetCurrentRange()

◆ GetDEDX()

◆ GetDEDXDispersion()

◆ GetDEDXForSubsec()

◆ GetKineticEnergy()

◆ GetLambda()

◆ GetMeanFreePath()

◆ GetModelByIndex()

◆ GetRange()

◆ GetRangeForLoss()

◆ InitialiseEnergyLossProcess()

◆ InverseRangeTable()

◆ IonisationTable()

◆ IonisationTableForSubsec()

◆ IsApplicable()

◆ IsIntegral()

◆ IsIonisationProcess()

◆ LambdaPhysicsVector()

◆ LambdaTable()

◆ MaxKinEnergy()

◆ MeanFreePath()

◆ MinKinEnergy()

◆ MinPrimaryEnergy()

◆ NumberOfModels()

◆ NumberOfSubCutoffRegions()

◆ Particle()

◆ PostStepDoIt()

◆ PostStepGetPhysicalInteractionLength()

◆ PreparePhysicsTable()

◆ PrintInfo()

◆ PrintInfoDefinition()

◆ RangeTableForLoss()

◆ RetrievePhysicsTable()

◆ SampleSubCutSecondaries()

◆ SecondaryParticle()

◆ SelectModel()

◆ SelectModelForMaterial()

◆ SetBaseParticle()

◆ SetCrossSectionBiasingFactor()

◆ SetCSDARangeTable()

◆ SetDEDXBinning()

◆ SetDEDXBinningForCSDARange()

◆ SetDEDXTable()

◆ SetDynamicMassCharge()

◆ SetEmModel()

◆ SetFluctModel()

◆ SetIntegral()