G4ionIonisation Class Reference

#include <G4ionIonisation.hh>

Inheritance diagram for G4ionIonisation:

Public Member Functions
	G4ionIonisation (const G4String &name="ionIoni")
virtual	~G4ionIonisation ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &p)
virtual void	PrintInfo ()
void	AddStoppingData (G4int Z, G4int A, const G4String &materialName, G4PhysicsVector *dVector)
void	ActivateStoppingData (G4bool)
Public Member Functions inherited from G4VEnergyLossProcess
	G4VEnergyLossProcess (const G4String &name="EnergyLoss", G4ProcessType type=fElectromagnetic)
virtual	~G4VEnergyLossProcess ()
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 *)
virtual G4double	MinPrimaryEnergy (const G4ParticleDefinition *p, const G4Material *, G4double cut)
G4double	BetheBlochEnergyThreshold ()
Protected Member Functions inherited from G4VEnergyLossProcess
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)
Protected Attributes inherited from G4VEnergyLossProcess
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

Detailed Description

Definition at line 78 of file G4ionIonisation.hh.

Constructor & Destructor Documentation

G4ionIonisation::G4ionIonisation

(

const G4String &

name = "ionIoni"

)

Definition at line 83 of file G4ionIonisation.cc.

References G4Electron::Electron(), G4LossTableManager::EmCorrections(), fIonisation, G4LossTableManager::Instance(), python.hepunit::MeV, python.hepunit::mm, G4VEnergyLossProcess::SetLinearLossLimit(), G4VProcess::SetProcessSubType(), G4VEnergyLossProcess::SetSecondaryParticle(), and G4VEnergyLossProcess::SetStepFunction().

  : G4VEnergyLossProcess(name),
    theParticle(0),
    isInitialised(false),
    stopDataActive(false)
{
  SetLinearLossLimit(0.02);
  SetStepFunction(0.1, 0.01*mm);
  SetProcessSubType(fIonisation);
  SetSecondaryParticle(G4Electron::Electron());
  corr = G4LossTableManager::Instance()->EmCorrections();
  eth = 2*MeV;
}

G4ionIonisation::~G4ionIonisation ( )

virtual

Definition at line 99 of file G4ionIonisation.cc.

{}

Member Function Documentation

void G4ionIonisation::ActivateStoppingData ( G4bool  val )

inline

Definition at line 125 of file G4ionIonisation.hh.

{
  stopDataActive = val;
}

void G4ionIonisation::AddStoppingData	(	G4int	Z,
		G4int	A,
		const G4String &	materialName,
		G4PhysicsVector *	dVector
	)

Definition at line 184 of file G4ionIonisation.cc.

References G4EmCorrections::AddStoppingData().

{
  corr->AddStoppingData(Z, A, mname, dVector);
}

G4double G4ionIonisation::BetheBlochEnergyThreshold ( )

inlineprotected

Definition at line 132 of file G4ionIonisation.hh.

{
  return eth;
}

void G4ionIonisation::InitialiseEnergyLossProcess ( const G4ParticleDefinition *  part, const G4ParticleDefinition *  bpart  )

protectedvirtual

Implements G4VEnergyLossProcess.

Definition at line 122 of file G4ionIonisation.cc.

References G4VEnergyLossProcess::AddEmModel(), G4VEnergyLossProcess::EmModel(), G4VEnergyLossProcess::FluctModel(), G4GenericIon::GenericIon(), G4ParticleDefinition::GetPDGMass(), G4VEmModel::HighEnergyLimit(), G4EmCorrections::InitialiseForNewRun(), G4VEnergyLossProcess::MaxKinEnergy(), G4VEnergyLossProcess::MinKinEnergy(), python.hepunit::proton_mass_c2, G4VEnergyLossProcess::SetBaseParticle(), G4VEnergyLossProcess::SetEmModel(), G4VEnergyLossProcess::SetFluctModel(), G4VEmModel::SetHighEnergyLimit(), G4EmCorrections::SetIonisationModels(), G4VEmModel::SetLowEnergyLimit(), and python.hepunit::TeV.

{
  const G4ParticleDefinition* ion = G4GenericIon::GenericIon();

  if(!isInitialised) {

    theParticle = part;

    // define base particle
    const G4ParticleDefinition* theBaseParticle = 0;

    if(part == ion)     { theBaseParticle = 0; }
    else if(bpart == 0) { theBaseParticle = ion; }
    else                { theBaseParticle = bpart; }

    SetBaseParticle(theBaseParticle);

    if (!EmModel(1)) { SetEmModel(new G4BraggIonModel(), 1); }
    EmModel(1)->SetLowEnergyLimit(MinKinEnergy());

    // model limit defined for protons
    eth = (EmModel(1)->HighEnergyLimit())*part->GetPDGMass()/proton_mass_c2;
    EmModel(1)->SetHighEnergyLimit(eth);

    if (!FluctModel()) { SetFluctModel(new G4IonFluctuations()); }
    AddEmModel(1, EmModel(1), FluctModel());

    if(eth < 10*TeV) {
      if (!EmModel(2)) { SetEmModel(new G4BetheBlochModel(),2); }  
      EmModel(2)->SetLowEnergyLimit(eth);
      EmModel(2)->SetHighEnergyLimit(MaxKinEnergy());
      AddEmModel(2, EmModel(2), FluctModel());    

      // Add ion stoping tables for Generic Ion
      if(part == ion) {
    stopDataActive = true;
    G4WaterStopping  ws(corr);
    corr->SetIonisationModels(EmModel(1),EmModel(2));
      }
    }
    isInitialised = true;
  }
  // reinitialisation of corrections for the new run
  if(part == ion) { corr->InitialiseForNewRun(); }
}

G4bool G4ionIonisation::IsApplicable ( const G4ParticleDefinition &  p )

virtual

Implements G4VEnergyLossProcess.

Definition at line 104 of file G4ionIonisation.cc.

References G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetPDGCharge(), and G4ParticleDefinition::IsShortLived().

{
  return (p.GetPDGCharge() != 0.0 && !p.IsShortLived() &&
          p.GetParticleType() == "nucleus");
}

G4double G4ionIonisation::MinPrimaryEnergy ( const G4ParticleDefinition *  p, const G4Material *  , G4double  cut  )

protectedvirtual

Reimplemented from G4VEnergyLossProcess.

Definition at line 112 of file G4ionIonisation.cc.

References G4ParticleDefinition::GetPDGMass().

{
  return 
    p->GetPDGMass()*(std::sqrt(1. + 0.5*cut/CLHEP::electron_mass_c2) - 1.0);
}

void G4ionIonisation::PrintInfo ( )

virtual

Implements G4VEnergyLossProcess.

Definition at line 172 of file G4ionIonisation.cc.

References G4cout, G4endl, G4GenericIon::GenericIon(), and G4EmCorrections::GetNumberOfStoppingVectors().

{
  if (stopDataActive && G4GenericIon::GenericIon() == theParticle) {
    G4cout << "      Stopping Power data for " 
           << corr->GetNumberOfStoppingVectors()
       << " ion/material pairs "
           << G4endl;
  }
}

---

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

• G4ionIonisation.hh
• G4ionIonisation.cc