#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 (const G4ParticleDefinition &part)

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

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

G4double	GetCurrentInteractionLength () const

void	SetPILfactor (G4double value)

G4double	GetPILfactor () const

G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)

G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)

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

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

const G4String &	GetProcessName () const

G4ProcessType	GetProcessType () const

void	SetProcessType (G4ProcessType)

G4int	GetProcessSubType () const

void	SetProcessSubType (G4int)

virtual void	EndTracking ()

virtual void	SetProcessManager (const G4ProcessManager *)

virtual const G4ProcessManager *	GetProcessManager ()

virtual void	ResetNumberOfInteractionLengthLeft ()

G4double	GetNumberOfInteractionLengthLeft () const

G4double	GetTotalNumberOfInteractionLengthTraversed () const

G4bool	isAtRestDoItIsEnabled () const

G4bool	isAlongStepDoItIsEnabled () const

G4bool	isPostStepDoItIsEnabled () const

virtual void	DumpInfo () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

virtual void	SetMasterProcess (G4VProcess *masterP)

const G4VProcess *	GetMasterProcess () const

virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)

virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Protected Member Functions
virtual void	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

void	SelectModel (G4double kinEnergy)

void	SetParticle (const G4ParticleDefinition *p)

void	SetSecondaryParticle (const G4ParticleDefinition *p)

Protected Member Functions inherited from G4VContinuousDiscreteProcess
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 122 of file G4VEnergyLossProcess.hh.

Constructor & Destructor Documentation

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

Definition at line 165 of file G4VEnergyLossProcess.cc.

References CandidateForSelection, DBL_MAX, G4Electron::Electron(), python.hepunit::eV, fParticleChange, G4Gamma::Gamma(), G4TransportationManager::GetSafetyHelper(), G4TransportationManager::GetTransportationManager(), python.hepunit::GeV, G4LossTableManager::Instance(), python.hepunit::keV, python.hepunit::mm, G4Positron::Positron(), G4VProcess::pParticleChange, G4LossTableManager::Register(), G4VParticleChange::SetSecondaryWeightByProcess(), SetStepFunction(), G4VProcess::SetVerboseLevel(), and python.hepunit::TeV.

                                           : 
   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;
   preStepKinEnergy = 0.0;
   preStepRangeEnergy = 0.0;
   computedRange = DBL_MAX;
 
   // 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
   lManager = G4LossTableManager::Instance();
   lManager->Register(this);
   fluctModel = 0;
   atomDeexcitation = 0;
 
   biasManager  = 0;
   biasFlag     = false; 
   weightFlag   = false; 
   isMaster     = true;
   lastIdx      = 0;
 
   idxDEDX = idxDEDXSub = idxDEDXunRestricted = idxIonisation =
     idxIonisationSub = idxRange = idxCSDA = idxSecRange =
     idxInverseRange = idxLambda = idxSubLambda = 0;
 
   scTracks.reserve(5);
   secParticles.reserve(5);
 
   secID = biasID = subsecID = -1;
 }

G4VEnergyLossProcess::~G4VEnergyLossProcess ( )

virtual

Definition at line 265 of file G4VEnergyLossProcess.cc.

References G4LossTableManager::DeRegister().

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

Member Function Documentation

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

Definition at line 2182 of file G4VEnergyLossProcess.cc.

References G4EmBiasingManager::ActivateForcedInteraction(), G4cout, G4endl, G4VProcess::GetProcessName(), python.hepunit::mm, and G4VProcess::verboseLevel.

Referenced by G4EmProcessOptions::ActivateForcedInteraction().

 {
   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);
 }

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

Definition at line 2202 of file G4VEnergyLossProcess.cc.

References G4EmBiasingManager::ActivateSecondaryBiasing(), G4Electron::Electron(), G4cout, G4endl, G4VProcess::GetProcessName(), python.hepunit::MeV, and G4VProcess::verboseLevel.

Referenced by G4EmProcessOptions::ActivateSecondaryBiasing().

 {
   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; 
     }
   }
 }

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

Definition at line 953 of file G4VEnergyLossProcess.cc.

References G4RegionStore::GetInstance(), and G4RegionStore::GetRegion().

Referenced by G4EmManager::Register(), and G4LossTableManager::Register().

 {
   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;
   }
 }

void G4VEnergyLossProcess::AddCollaborativeProcess ( G4VEnergyLossProcess * p )

Definition at line 1935 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

 {
   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; 
     }
   }
 }

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

Definition at line 344 of file G4VEnergyLossProcess.cc.

References G4EmModelManager::AddEmModel(), G4VProcess::pParticleChange, and G4VEmModel::SetParticleChange().

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

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1158 of file G4VEnergyLossProcess.cc.

References G4VAtomDeexcitation::AlongStepDeexcitation(), G4SafetyHelper::ComputeSafety(), G4VEmModel::CorrectionsAlongStep(), fGeomBoundary, fParticleChange, G4Track::GetDynamicParticle(), G4VEmModel::GetModelOfFluctuations(), G4VParticleChange::GetParentWeight(), G4VEmModel::GetParticleCharge(), G4Track::GetParticleDefinition(), G4StepPoint::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4ProductionCuts::GetProductionCut(), G4MaterialCutsCouple::GetProductionCuts(), G4StepPoint::GetSafety(), G4Step::GetStepLength(), G4StepPoint::GetStepStatus(), G4ParticleChangeForLoss::InitializeForAlongStep(), G4VEmModel::IsActive(), G4VEmModel::MaxSecondaryKinEnergy(), G4INCL::Math::min(), G4VParticleChange::ProposeLocalEnergyDeposit(), G4VParticleChange::ProposeWeight(), G4VEmFluctuationModel::SampleFluctuations(), SampleSubCutSecondaries(), G4ParticleChangeForLoss::SetProposedCharge(), G4ParticleChangeForLoss::SetProposedKineticEnergy(), and test::x.

 {
   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) {
     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;
   }
   //G4cout << theDEDXTable << "  idx= " << basedCoupleIndex 
   // << "  " << GetProcessName() << "  "<< currentMaterial->GetName()<<G4endl;
   //if(particle->GetParticleName() == "e-")G4cout << (*theDEDXTable) <<G4endl;
   // Short step
   eloss = GetDEDXForScaledEnergy(preStepScaledEnergy)*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;
 
   // 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);
       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(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;
 }

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1021 of file G4VEnergyLossProcess.cc.

References DBL_MAX, G4ProductionCuts::GetProductionCut(), G4MaterialCutsCouple::GetProductionCuts(), G4INCL::Math::min(), and test::x.

Referenced by ContinuousStepLimit().

 {
   G4double x = DBL_MAX;
   *selection = aGPILSelection;
   if(isIonisation) {
     fRange = GetScaledRangeForScaledEnergy(preStepScaledEnergy)*reduceFactor;
     x = fRange;
     G4double finR = finalRange;
     if(rndmStepFlag) { 
       finR = std::min(finR,
               currentCouple->GetProductionCuts()->GetProductionCut(1));
     }
     if(fRange > finR) { 
       x = fRange*dRoverRange + finR*(1.0 - dRoverRange)*(2.0 - finR/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;
 }

const G4ParticleDefinition * G4VEnergyLossProcess::BaseParticle ( ) const

inline

Definition at line 916 of file G4VEnergyLossProcess.hh.

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

 {
   return baseParticle;
 }

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

Definition at line 706 of file G4VEnergyLossProcess.cc.

References plottest35::bin, G4EmModelManager::FillDEDXVector(), G4PhysicsVector::FillSecondDerivatives(), fRestricted, fSubRestricted, fTotal, G4cout, G4endl, G4LossTableBuilder::GetFlag(), G4PhysicsTable::GetFlag(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4ProductionCutsTable::GetProductionCutsTable(), G4LossTableManager::GetTableBuilder(), G4ProductionCutsTable::GetTableSize(), G4PhysicsTableHelper::SetPhysicsVector(), G4PhysicsVector::SetSpline(), G4LossTableManager::SplineFlag(), and G4VProcess::verboseLevel.

 {
   G4bool verb = false;
   if(1 < verboseLevel || verb) {
     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;
   } 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 || verb) {
     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 = lManager->SplineFlag();
   G4PhysicsLogVector* aVector = 0;
   G4PhysicsLogVector* bVector = 0;
 
   for(size_t i=0; i<numOfCouples; ++i) {
 
     if(1 < verboseLevel || verb) {
       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);
       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 || verb) {
     G4cout << "G4VEnergyLossProcess::BuildDEDXTable(): table is built for "
            << particle->GetParticleName()
            << " and process " << GetProcessName()
            << G4endl;
     //if(2 < verboseLevel) G4cout << (*table) << G4endl;
   }
 
   return table;
 }

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

Definition at line 795 of file G4VEnergyLossProcess.cc.

References plottest35::bin, G4EmModelManager::FillLambdaVector(), G4PhysicsVector::FillSecondDerivatives(), fRestricted, fSubRestricted, G4cout, G4endl, G4Log(), G4lrint(), G4LossTableBuilder::GetCoupleIndexes(), G4LossTableBuilder::GetDensityFactors(), G4LossTableBuilder::GetFlag(), G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4ProductionCutsTable::GetProductionCutsTable(), G4LossTableManager::GetTableBuilder(), G4ProductionCutsTable::GetTableSize(), MinPrimaryEnergy(), G4PhysicsTableHelper::SetPhysicsVector(), G4PhysicsVector::SetSpline(), G4LossTableManager::SplineFlag(), and G4VProcess::verboseLevel.

 {
   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 = lManager->GetTableBuilder();
   theDensityFactor = bld->GetDensityFactors();
   theDensityIdx = bld->GetCoupleIndexes();
 
   G4bool splineFlag = lManager->SplineFlag();
   G4PhysicsLogVector* aVector = 0;
   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);
       if(bin < 3) { 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;
 }

void G4VEnergyLossProcess::BuildPhysicsTable ( const G4ParticleDefinition & part )

virtual

Reimplemented from G4VProcess.

Definition at line 603 of file G4VEnergyLossProcess.cc.

References G4LossTableManager::AtomDeexcitation(), G4LossTableManager::BuildPhysicsTable(), CSDARangeTable(), DEDXTable(), DEDXTableForSubsec(), DEDXunRestrictedTable(), fIsIonisation, fIsSubIonisation, fParticleChange, fRestricted, fSubRestricted, fTotal, G4cout, G4endl, G4LossTableBuilder::GetCoupleIndexes(), G4LossTableBuilder::GetDensityFactors(), G4VProcess::GetMasterProcess(), GetModelByIndex(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4LossTableManager::GetTableBuilder(), G4LossTableBuilder::InitialiseBaseMaterials(), G4SafetyHelper::InitialiseHelper(), G4VEmModel::InitialiseLocal(), InverseRangeTable(), IonisationTable(), IonisationTableForSubsec(), IsIonisationProcess(), G4VAtomDeexcitation::IsPIXEActive(), LambdaTable(), G4LossTableManager::LocalPhysicsTables(), G4EmModelManager::NumberOfModels(), PrintInfoDefinition(), RangeTableForLoss(), SecondaryRangeTable(), SetCSDARangeTable(), SetDEDXTable(), SetInverseRangeTable(), SetLambdaTable(), G4ParticleChangeForLoss::SetLowEnergyLimit(), SetRangeTableForLoss(), SetSecondaryRangeTable(), SetSubLambdaTable(), SubLambdaTable(), and G4VProcess::verboseLevel.

Referenced by G4ePolarizedIonisation::BuildPhysicsTable().

 {
   G4bool verb =  false;
   if(1 < verboseLevel || verb) {
   
     //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;
   }
 
   G4bool master = true;
   const G4VEnergyLossProcess* masterProcess = 
     static_cast<const G4VEnergyLossProcess*>(GetMasterProcess());
   if(masterProcess != this) { master = false; }
 
   if(&part == particle) {
 
     G4LossTableBuilder* bld = lManager->GetTableBuilder();
     if(master) {
       theDensityFactor = bld->GetDensityFactors();
       theDensityIdx = bld->GetCoupleIndexes();
       lManager->BuildPhysicsTable(particle, this);
 
     } else {
 
       // define density factors for worker thread
       bld->InitialiseBaseMaterials(masterProcess->DEDXTable()); 
       theDensityFactor = bld->GetDensityFactors();
       theDensityIdx = bld->GetCoupleIndexes();
 
       // 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")))
     { 
       PrintInfoDefinition(part); 
     }
 
   // Added tracking cut to avoid tracking artifacts
   // identify deexcitation flag
   if(isIonisation) { 
     fParticleChange.SetLowEnergyLimit(lowestKinEnergy); 
     atomDeexcitation = lManager->AtomDeexcitation();
     if(atomDeexcitation) { 
       if(atomDeexcitation->IsPIXEActive()) { useDeexcitation = true; } 
     }
   }
 
   //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;
   }
 }

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

Definition at line 1892 of file G4VEnergyLossProcess.cc.

References AlongStepGetPhysicalInteractionLength().

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

G4double G4VEnergyLossProcess::CrossSectionBiasingFactor ( ) const

inline

Definition at line 1075 of file G4VEnergyLossProcess.hh.

 {
   return biasFactor;
 }

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

Definition at line 1860 of file G4VEnergyLossProcess.cc.

References G4VEmModel::CrossSectionPerVolume(), and SelectModel().

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

G4PhysicsTable * G4VEnergyLossProcess::CSDARangeTable ( ) const

inline

Definition at line 1124 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theCSDARangeTable;
 }

size_t G4VEnergyLossProcess::CurrentMaterialCutsCoupleIndex ( ) const

inlineprotected

Definition at line 605 of file G4VEnergyLossProcess.hh.

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

 {
   return currentCoupleIndex;
 }

G4PhysicsTable * G4VEnergyLossProcess::DEDXTable ( ) const

inline

Definition at line 1089 of file G4VEnergyLossProcess.hh.

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

 {
   return theDEDXTable;
 }

G4PhysicsTable * G4VEnergyLossProcess::DEDXTableForSubsec ( ) const

inline

Definition at line 1096 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theDEDXSubTable;
 }

G4PhysicsTable * G4VEnergyLossProcess::DEDXunRestrictedTable ( ) const

inline

Definition at line 1103 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theDEDXunRestrictedTable;
 }

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

Definition at line 371 of file G4VEnergyLossProcess.cc.

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

G4VEmFluctuationModel * G4VEnergyLossProcess::FluctModel ( )

inline

Definition at line 881 of file G4VEnergyLossProcess.hh.

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

 {
   return fluctModel;
 }

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

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 1914 of file G4VEnergyLossProcess.cc.

References DBL_MAX.

 {
   return DBL_MAX;
 }

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

inline

Definition at line 801 of file G4VEnergyLossProcess.hh.

References DBL_MAX, and test::x.

Referenced by G4LossTableManager::GetCSDARange().

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

const G4Element * G4VEnergyLossProcess::GetCurrentElement ( ) const

Definition at line 2155 of file G4VEnergyLossProcess.cc.

References G4VEmModel::GetCurrentElement().

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

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

inline

Definition at line 764 of file G4VEnergyLossProcess.hh.

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

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

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

Definition at line 1842 of file G4VEnergyLossProcess.cc.

References G4VEmFluctuationModel::Dispersion(), fm, G4DynamicParticle::GetKineticEnergy(), G4VEmModel::GetModelOfFluctuations(), G4VEmModel::MaxSecondaryKinEnergy(), G4INCL::Math::min(), and SelectModel().

Referenced by G4LossTableManager::GetDEDXDispersion().

 {
   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;
 }

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

inline

Definition at line 774 of file G4VEnergyLossProcess.hh.

Referenced by G4LossTableManager::GetSubDEDX().

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

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

inline

Definition at line 830 of file G4VEnergyLossProcess.hh.

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

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

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

inline

Definition at line 840 of file G4VEnergyLossProcess.hh.

References test::x.

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

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

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 1902 of file G4VEnergyLossProcess.cc.

References MeanFreePath(), and NotForced.

Referenced by G4ePolarizedIonisation::GetMeanFreePath().

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

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

Definition at line 380 of file G4VEnergyLossProcess.cc.

References G4EmModelManager::GetModel().

Referenced by BuildPhysicsTable().

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

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

inline

Definition at line 784 of file G4VEnergyLossProcess.hh.

References test::x.

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

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

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

inline

Definition at line 815 of file G4VEnergyLossProcess.hh.

References test::x.

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

 {
   //  G4cout << "GetRangeForLoss: Range from " << GetProcessName() << G4endl;
   DefineMaterial(couple);
   G4double x = 
     GetScaledRangeForScaledEnergy(kineticEnergy*massRatio)*reduceFactor;
   //G4cout << "GetRangeForLoss: Range from " << GetProcessName() 
   //         << "  e= " << kineticEnergy << " r= " << x << G4endl;
   return x;
 }

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

protectedpure virtual

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

Referenced by PreparePhysicsTable().

G4PhysicsTable * G4VEnergyLossProcess::InverseRangeTable ( ) const

inline

Definition at line 1145 of file G4VEnergyLossProcess.hh.

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

 {
   return theInverseRangeTable;
 }

G4PhysicsTable * G4VEnergyLossProcess::IonisationTable ( ) const

inline

Definition at line 1110 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theIonisationTable;
 }

G4PhysicsTable * G4VEnergyLossProcess::IonisationTableForSubsec ( ) const

inline

Definition at line 1117 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theIonisationSubTable;
 }

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

pure virtual

Reimplemented from G4VProcess.

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

G4bool G4VEnergyLossProcess::IsIntegral ( ) const

inline

Definition at line 951 of file G4VEnergyLossProcess.hh.

 {
   return integral;
 }

G4bool G4VEnergyLossProcess::IsIonisationProcess ( ) const

inline

Definition at line 967 of file G4VEnergyLossProcess.hh.

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

 {
   return isIonisation;
 }

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

protected

Definition at line 1924 of file G4VEnergyLossProcess.cc.

References G4PhysicsVector::SetSpline(), G4LossTableManager::SplineFlag(), and test::v.

Referenced by G4ePolarizedIonisation::BuildPhysicsTable().

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

G4PhysicsTable * G4VEnergyLossProcess::LambdaTable ( ) const

inline

Definition at line 1152 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theLambdaTable;
 }

G4double G4VEnergyLossProcess::MaxKinEnergy ( ) const

inline

Definition at line 1061 of file G4VEnergyLossProcess.hh.

 {
   return maxKinEnergy;
 }

G4double G4VEnergyLossProcess::MeanFreePath ( const G4Track & track )

Definition at line 1881 of file G4VEnergyLossProcess.cc.

References DBL_MAX, G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), and test::x.

Referenced by GetMeanFreePath().

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

G4double G4VEnergyLossProcess::MinKinEnergy ( ) const

inline

Definition at line 1046 of file G4VEnergyLossProcess.hh.

 {
   return minKinEnergy;
 }

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

protectedvirtual

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

Definition at line 335 of file G4VEnergyLossProcess.cc.

Referenced by BuildLambdaTable().

 {
   return cut;
 }

G4int G4VEnergyLossProcess::NumberOfModels ( ) const

Definition at line 387 of file G4VEnergyLossProcess.cc.

References G4EmModelManager::NumberOfModels().

 {
   return modelManager->NumberOfModels();
 }

G4int G4VEnergyLossProcess::NumberOfSubCutoffRegions ( ) const

inline

Definition at line 1011 of file G4VEnergyLossProcess.hh.

 {
   return nSCoffRegions;
 }

const G4ParticleDefinition * G4VEnergyLossProcess::Particle ( ) const

inline

Definition at line 909 of file G4VEnergyLossProcess.hh.

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

 {
   return particle;
 }

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1523 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<< "Eprimary: "<<dynParticle->GetKineticEnergy()/MeV<<G4endl;
   currentModel->SampleSecondaries(&secParticles, currentCouple, 
                   dynParticle, 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());
     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;
 }

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 1056 of file G4VEnergyLossProcess.cc.

References G4VEmModel::ChargeSquareRatio(), python.hepunit::cm, G4VProcess::currentInteractionLength, DBL_MAX, G4EmBiasingManager::ForcedInteractionRegion(), G4cout, G4endl, G4Log(), G4UniformRand, G4Track::GetDefinition(), G4Track::GetKineticEnergy(), G4Track::GetMaterial(), G4Track::GetMaterialCutsCouple(), G4Material::GetName(), G4Track::GetParentID(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4EmBiasingManager::GetStepLimit(), G4VEmModel::IsActive(), python.hepunit::MeV, NotForced, SelectModel(), G4VProcess::theInitialNumberOfInteractionLength, G4VProcess::theNumberOfInteractionLengthLeft, G4VProcess::verboseLevel, and test::x.

Referenced by G4ePolarizedIonisation::PostStepGetPhysicalInteractionLength().

 {
   // 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();
   preStepScaledEnergy = preStepKinEnergy*massRatio;
   SelectModel(preStepScaledEnergy);
 
   if(!currentModel->IsActive(preStepScaledEnergy)) { 
     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; 
   // 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();
       theNumberOfInteractionLengthLeft =  -G4Log( G4UniformRand() );
       theInitialNumberOfInteractionLength = theNumberOfInteractionLengthLeft; 
 
     } 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 
        << "  " << track.GetMaterial()->GetName()
        <<G4endl;
     G4cout << "MeanFreePath = " << currentInteractionLength/cm << "[cm]" 
        << "InteractionLength= " << x/cm <<"[cm] " <<G4endl;
   }
 #endif
   return x;
 }

void G4VEnergyLossProcess::PreparePhysicsTable ( const G4ParticleDefinition & part )

virtual

Reimplemented from G4VProcess.

Definition at line 395 of file G4VEnergyLossProcess.cc.

References G4LossTableManager::BuildCSDARange(), G4PhysicsTable::clearAndDestroy(), DBL_MAX, G4ParticleTable::FindParticle(), G4cout, G4endl, G4ProcessManager::GetAlongStepProcessVector(), G4VProcess::GetMasterProcess(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4EmModelManager::GetModel(), G4Region::GetName(), G4ParticleDefinition::GetParticleName(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), G4ParticleDefinition::GetProcessManager(), G4VProcess::GetProcessName(), G4MaterialCutsCouple::GetProductionCuts(), G4ProductionCutsTable::GetProductionCutsTable(), G4LossTableManager::GetTableBuilder(), G4ProductionCutsTable::GetTableSize(), G4VEmModel::HighEnergyLimit(), G4EmBiasingManager::Initialise(), G4EmModelManager::Initialise(), G4LossTableBuilder::InitialiseBaseMaterials(), InitialiseEnergyLossProcess(), n, G4EmModelManager::NumberOfModels(), eplot::pname, G4PhysicsTableHelper::PreparePhysicsTable(), G4LossTableManager::PreparePhysicsTable(), G4PhysicsModelCatalog::Register(), G4LossTableManager::RegisterExtraParticle(), G4VEmModel::SetHighEnergyLimit(), G4VEmModel::SetMasterThread(), G4ProcessVector::size(), G4EmModelManager::SubCutoff(), test::v, and G4VProcess::verboseLevel.

 {
   if(1 < verboseLevel) {
     G4cout << "G4VEnergyLossProcess::PreparePhysicsTable for "
            << GetProcessName() << " for " << part.GetParticleName() 
        << "  " << this << G4endl;
   }
 
   if(GetMasterProcess() != this) { isMaster = false; }
 
   currentCouple = 0;
   preStepLambda = 0.0;
   mfpKinEnergy  = DBL_MAX;
   fRange        = DBL_MAX;
   preStepKinEnergy = 0.0;
   preStepRangeEnergy = 0.0;
   chargeSqRatio = 1.0;
   massRatio = 1.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);
 
   // Tables preparation
   if (isMaster && !baseParticle) {
 
     if(theDEDXTable && isIonisation) {
       if(theIonisationTable && theDEDXTable != theIonisationTable) {
     theDEDXTable->clearAndDestroy();
     theDEDXTable = theIonisationTable;
       }   
       if(theDEDXSubTable && theIonisationSubTable && 
      theDEDXSubTable != theIonisationSubTable) {
     theDEDXSubTable->clearAndDestroy();
     theDEDXSubTable = theIonisationSubTable;
       }   
     }
     
     theDEDXTable = G4PhysicsTableHelper::PreparePhysicsTable(theDEDXTable);
     bld->InitialiseBaseMaterials(theDEDXTable);
 
     if(theDEDXSubTable) {
       theDEDXSubTable = 
     G4PhysicsTableHelper::PreparePhysicsTable(theDEDXSubTable);
     }
 
     if (lManager->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) {
       theDEDXSubTable = 
     G4PhysicsTableHelper::PreparePhysicsTable(theDEDXSubTable);
       theSubLambdaTable = 
     G4PhysicsTableHelper::PreparePhysicsTable(theSubLambdaTable);
     }
   }
 
   // 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); }
   }
 
   // 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);
     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::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;
       }
     }
   }
 }

virtual void G4VEnergyLossProcess::PrintInfo ( )

pure virtual

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

Referenced by PrintInfoDefinition().

void G4VEnergyLossProcess::PrintInfoDefinition ( const G4ParticleDefinition & part )

Definition at line 876 of file G4VEnergyLossProcess.cc.

References G4EmModelManager::DumpModelList(), G4BestUnit, G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4VProcess::GetProcessSubType(), python.hepunit::mm, PrintInfo(), G4LossTableManager::SplineFlag(), and G4VProcess::verboseLevel.

Referenced by BuildPhysicsTable().

 {
   if(0 < verboseLevel) {
     G4cout << std::setprecision(6);
     G4cout << G4endl << GetProcessName() << ":   for  "
            << part.GetParticleName()
        << "    SubType= " << GetProcessSubType() 
            << G4endl;
     G4cout << "      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: " 
        << lManager->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;
       }
     }
   }
 }

G4PhysicsTable * G4VEnergyLossProcess::RangeTableForLoss ( ) const

inline

Definition at line 1138 of file G4VEnergyLossProcess.hh.

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

 {
   return theRangeTableForLoss;
 }

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

virtual

Reimplemented from G4VProcess.

Definition at line 1720 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

 {
   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;
 }

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

Definition at line 1445 of file G4VEnergyLossProcess.cc.

References python.hepunit::electron_mass_c2, G4UniformRand, G4Track::GetDynamicParticle(), G4StepPoint::GetGlobalTime(), G4DynamicParticle::GetKineticEnergy(), G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), G4Track::GetParticleDefinition(), G4StepPoint::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4Step::GetStepLength(), G4Track::GetTouchableHandle(), G4Step::GetTrack(), python.hepunit::perMillion, G4VEmModel::SampleSecondaries(), G4Track::SetCreatorModelIndex(), and G4Track::SetTouchableHandle().

Referenced by AlongStepDoIt().

 {
   // 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 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) {
 
       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;
     */
     }
   } while (fragment <= 1.0);
   return esec;
 } 

const G4ParticleDefinition * G4VEnergyLossProcess::SecondaryParticle ( ) const

inline

Definition at line 923 of file G4VEnergyLossProcess.hh.

 {
   return secondaryParticle;
 }

G4PhysicsTable * G4VEnergyLossProcess::SecondaryRangeTable ( ) const

inline

Definition at line 1131 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theSecondaryRangeTable;
 }

void G4VEnergyLossProcess::SelectModel ( G4double kinEnergy )

inlineprotected

Definition at line 612 of file G4VEnergyLossProcess.hh.

References G4EmModelManager::SelectModel(), and G4VEmModel::SetCurrentCouple().

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

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

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

inline

Definition at line 620 of file G4VEnergyLossProcess.hh.

References G4EmModelManager::SelectModel().

Referenced by G4ContinuousGainOfEnergy::GetContinuousStepLimit().

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

void G4VEnergyLossProcess::SetBaseParticle ( const G4ParticleDefinition * p )

inline

Definition at line 902 of file G4VEnergyLossProcess.hh.

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

 {
   baseParticle = p;
 }

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

Definition at line 2164 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

Referenced by G4EmProcessOptions::SetProcessBiasingFactor().

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

void G4VEnergyLossProcess::SetCSDARangeTable ( G4PhysicsTable * pRange )

Definition at line 2020 of file G4VEnergyLossProcess.cc.

References G4PhysicsTable::length(), n, and G4PhysicsVector::Value().

Referenced by BuildPhysicsTable().

 {
   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;
     }
   }
 }

void G4VEnergyLossProcess::SetDEDXBinning ( G4int nbins )

inline

Definition at line 1018 of file G4VEnergyLossProcess.hh.

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

 {
   nBins = nbins;
 }

void G4VEnergyLossProcess::SetDEDXBinningForCSDARange ( G4int nbins )

inline

Definition at line 1032 of file G4VEnergyLossProcess.hh.

 {
   nBinsCSDA = nbins;
 }

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

Definition at line 1962 of file G4VEnergyLossProcess.cc.

References fIsIonisation, fIsSubIonisation, fRestricted, fSubRestricted, fTotal, G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4PhysicsTable::length(), n, and G4PhysicsVector::Value().

Referenced by BuildPhysicsTable().

 {
   G4bool verb = false;
   if(fTotal == tType) {
     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, 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) {
     if(verb) {
       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) {
     if(verb) {
       G4cout<< "G4VEnergyLossProcess::SetIonisationTable "
         << particle->GetParticleName()
             << " oldTable " << theDEDXTable << " newTable " << p 
             << " ion " << theIonisationTable 
         << " IsMaster " << isMaster 
         << " " << GetProcessName() << G4endl;
     }
     theIonisationTable = p;
   } else if(fIsSubIonisation == tType) {
     theIonisationSubTable = p;
   }
 }

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

inline

Definition at line 645 of file G4VEnergyLossProcess.hh.

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

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

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

Definition at line 362 of file G4VEnergyLossProcess.cc.

References n.

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

void G4VEnergyLossProcess::SetFluctModel ( G4VEmFluctuationModel * p )

inline

Definition at line 874 of file G4VEnergyLossProcess.hh.

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

 {
   fluctModel = p;
 }

void G4VEnergyLossProcess::SetIntegral ( G4bool val )

inline

Definition at line 944 of file G4VEnergyLossProcess.hh.

Referenced by G4EmManager::Register(), and G4LossTableManager::Register().

 {
   integral = val;
 }

void G4VEnergyLossProcess::SetInverseRangeTable ( G4PhysicsTable * p )

Definition at line 2070 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

Referenced by BuildPhysicsTable().

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

void G4VEnergyLossProcess::SetIonisation ( G4bool val )

inline

Definition at line 958 of file G4VEnergyLossProcess.hh.

References CandidateForSelection, and NotCandidateForSelection.

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

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

void G4VEnergyLossProcess::SetLambdaBinning ( G4int nbins )

inline

Definition at line 1025 of file G4VEnergyLossProcess.hh.

 {
   nBins = nbins;
 }

void G4VEnergyLossProcess::SetLambdaFactor ( G4double val )

inline

Definition at line 988 of file G4VEnergyLossProcess.hh.

Referenced by G4EmProcessOptions::SetLambdaFactor().

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

void G4VEnergyLossProcess::SetLambdaTable ( G4PhysicsTable * p )

Definition at line 2082 of file G4VEnergyLossProcess.cc.

References DBL_MAX, G4PhysicsVector::Energy(), G4cout, G4endl, G4LossTableBuilder::GetCoupleIndexes(), G4LossTableBuilder::GetDensityFactors(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4LossTableManager::GetTableBuilder(), G4PhysicsVector::GetVectorLength(), G4PhysicsTable::length(), python.hepunit::MeV, n, smax, and G4VProcess::verboseLevel.

Referenced by BuildPhysicsTable().

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

void G4VEnergyLossProcess::SetLinearLossLimit ( G4double val )

inline

Definition at line 974 of file G4VEnergyLossProcess.hh.

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

 {
   linLossLimit = val;
 }

void G4VEnergyLossProcess::SetLossFluctuations ( G4bool val )

inline

Definition at line 930 of file G4VEnergyLossProcess.hh.

Referenced by G4AdjointPhysicsList::ConstructEM(), G4EmManager::Register(), G4LossTableManager::Register(), and G4AdjointPhysicsList::SetLossFluctuationFlag().

 {
   lossFluctuationFlag = val;
 }

void G4VEnergyLossProcess::SetLowestEnergyLimit ( G4double val )

inline

Definition at line 1004 of file G4VEnergyLossProcess.hh.

Referenced by G4EmUserPhysics::ConstructProcess().

 {
   lowestKinEnergy = val;
 }

void G4VEnergyLossProcess::SetMaxKinEnergy ( G4double e )

inline

Definition at line 1053 of file G4VEnergyLossProcess.hh.

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

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

void G4VEnergyLossProcess::SetMaxKinEnergyForCSDARange ( G4double e )

inline

Definition at line 1068 of file G4VEnergyLossProcess.hh.

 {
   maxKinEnergyCSDA = e;
 }

void G4VEnergyLossProcess::SetMinKinEnergy ( G4double e )

inline

Definition at line 1039 of file G4VEnergyLossProcess.hh.

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

 {
   minKinEnergy = e;
 }

void G4VEnergyLossProcess::SetMinSubRange ( G4double val )

inline

Definition at line 981 of file G4VEnergyLossProcess.hh.

 {
   minSubRange = val;
 }

void G4VEnergyLossProcess::SetParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 888 of file G4VEnergyLossProcess.hh.

 {
   particle = p;
 }

void G4VEnergyLossProcess::SetRandomStep ( G4bool val )

inline

Definition at line 937 of file G4VEnergyLossProcess.hh.

Referenced by G4EmManager::Register(), and G4LossTableManager::Register().

 {
   rndmStepFlag = val;
 }

void G4VEnergyLossProcess::SetRangeTableForLoss ( G4PhysicsTable * p )

Definition at line 2046 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

Referenced by BuildPhysicsTable().

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

void G4VEnergyLossProcess::SetSecondaryParticle ( const G4ParticleDefinition * p )

inlineprotected

Definition at line 895 of file G4VEnergyLossProcess.hh.

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

 {
   secondaryParticle = p;
 }

void G4VEnergyLossProcess::SetSecondaryRangeTable ( G4PhysicsTable * p )

Definition at line 2058 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

Referenced by BuildPhysicsTable().

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

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

inline

Definition at line 995 of file G4VEnergyLossProcess.hh.

 {
   dRoverRange = v1;
   finalRange = v2;
   if (dRoverRange > 0.999) { dRoverRange = 1.0; }
 }

void G4VEnergyLossProcess::SetSubLambdaTable ( G4PhysicsTable * p )

Definition at line 2143 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

Referenced by BuildPhysicsTable().

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

void G4VEnergyLossProcess::StartTracking ( G4Track * track )

virtual

Reimplemented from G4VProcess.

Definition at line 982 of file G4VEnergyLossProcess.cc.

References DBL_MAX, G4Track::GetDefinition(), G4Track::GetParentID(), G4ParticleDefinition::GetPDGMass(), python.hepunit::proton_mass_c2, G4EmBiasingManager::ResetForcedInteraction(), and G4VProcess::theNumberOfInteractionLengthLeft.

 {
   /*    
     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;
   mfpKinEnergy = DBL_MAX; 
   preStepRangeEnergy = 0.0;
 
   // reset ion
   if(isIon) {
     chargeSqRatio = 0.5;
 
     G4double newmass = track->GetDefinition()->GetPDGMass();
     if(baseParticle) {
       massRatio = baseParticle->GetPDGMass()/newmass;
     } else if(theGenericIon) {
       massRatio = proton_mass_c2/newmass;
     } else {
       massRatio = 1.0;
     }
   }  
   // forced biasing only for primary particles
   if(biasManager) {
     if(0 == track->GetParentID()) {
       // primary particle
       biasFlag = true; 
       biasManager->ResetForcedInteraction(); 
     }
   }
 }

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

virtual

Reimplemented from G4VProcess.

Definition at line 1659 of file G4VEnergyLossProcess.cc.

References G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

 {
   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;
 }

G4PhysicsTable * G4VEnergyLossProcess::SubLambdaTable ( ) const

inline

Definition at line 1159 of file G4VEnergyLossProcess.hh.

Referenced by BuildPhysicsTable().

 {
   return theSubLambdaTable;
 }

G4bool G4VEnergyLossProcess::TablesAreBuilt ( ) const

inline

Definition at line 1082 of file G4VEnergyLossProcess.hh.

 {
   return  tablesAreBuilt;
 }

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

Definition at line 354 of file G4VEnergyLossProcess.cc.

References G4EmModelManager::UpdateEmModel().

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

Field Documentation

G4ParticleChangeForLoss G4VEnergyLossProcess::fParticleChange

protected

Definition at line 563 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

Member Function Documentation

Field Documentation