G4hRDEnergyLoss Class Reference

#include <G4hRDEnergyLoss.hh>

Inheritance diagram for G4hRDEnergyLoss:

Collaboration diagram for G4hRDEnergyLoss:

Public Member Functions
	G4hRDEnergyLoss (const G4String &)
	~G4hRDEnergyLoss ()
virtual G4double	GetMeanFreePath (const G4Track &track, G4double previousStepSize, enum G4ForceCondition *condition)=0
virtual G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &Step)=0
Public Member Functions inherited from G4VContinuousDiscreteProcess
	G4VContinuousDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)
	G4VContinuousDiscreteProcess (G4VContinuousDiscreteProcess &)
virtual	~G4VContinuousDiscreteProcess ()
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety, G4GPILSelection *selection)
virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)
Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
	G4VProcess (const G4VProcess &right)
virtual	~G4VProcess ()
G4int	operator== (const G4VProcess &right) const
G4int	operator!= (const G4VProcess &right) const
G4double	GetCurrentInteractionLength () const
void	SetPILfactor (G4double value)
G4double	GetPILfactor () const
G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	PreparePhysicsTable (const G4ParticleDefinition &)
virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
const G4String &	GetProcessName () const
G4ProcessType	GetProcessType () const
void	SetProcessType (G4ProcessType)
G4int	GetProcessSubType () const
void	SetProcessSubType (G4int)
virtual void	StartTracking (G4Track *)
virtual void	EndTracking ()
virtual void	SetProcessManager (const G4ProcessManager *)
virtual const G4ProcessManager *	GetProcessManager ()
virtual void	ResetNumberOfInteractionLengthLeft ()
G4double	GetNumberOfInteractionLengthLeft () const
G4double	GetTotalNumberOfInteractionLengthTraversed () const
G4bool	isAtRestDoItIsEnabled () const
G4bool	isAlongStepDoItIsEnabled () const
G4bool	isPostStepDoItIsEnabled () const
virtual void	DumpInfo () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
virtual void	SetMasterProcess (G4VProcess *masterP)
const G4VProcess *	GetMasterProcess () const
virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Static Public Member Functions
static G4int	GetNumberOfProcesses ()
static void	SetNumberOfProcesses (G4int number)
static void	PlusNumberOfProcesses ()
static void	MinusNumberOfProcesses ()
static void	SetdRoverRange (G4double value)
static void	SetRndmStep (G4bool value)
static void	SetEnlossFluc (G4bool value)
static void	SetStepFunction (G4double c1, G4double c2)
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Protected Member Functions
G4bool	CutsWhereModified ()
Protected Member Functions inherited from G4VContinuousDiscreteProcess
virtual G4double	GetContinuousStepLimit (const G4Track &aTrack, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)=0
void	SetGPILSelection (G4GPILSelection selection)
G4GPILSelection	GetGPILSelection () const
Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
void	ClearNumberOfInteractionLengthLeft ()

Static Protected Member Functions
static void	BuildDEDXTable (const G4ParticleDefinition &aParticleType)

Protected Attributes
const G4double	MaxExcitationNumber
const G4double	probLimFluct
const long	nmaxDirectFluct
const long	nmaxCont1
const long	nmaxCont2
G4PhysicsTable *	theLossTable
G4double	linLossLimit
G4double	MinKineticEnergy
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

Static Protected Attributes
static G4ThreadLocal G4PhysicsTable *	theDEDXpTable = 0
static G4ThreadLocal G4PhysicsTable *	theDEDXpbarTable = 0
static G4ThreadLocal G4PhysicsTable *	theRangepTable = 0
static G4ThreadLocal G4PhysicsTable *	theRangepbarTable = 0
static G4ThreadLocal G4PhysicsTable *	theInverseRangepTable = 0
static G4ThreadLocal G4PhysicsTable *	theInverseRangepbarTable = 0
static G4ThreadLocal G4PhysicsTable *	theLabTimepTable = 0
static G4ThreadLocal G4PhysicsTable *	theLabTimepbarTable = 0
static G4ThreadLocal G4PhysicsTable *	theProperTimepTable = 0
static G4ThreadLocal G4PhysicsTable *	theProperTimepbarTable = 0
static G4ThreadLocal G4PhysicsTable **	RecorderOfpProcess = 0
static G4ThreadLocal G4PhysicsTable **	RecorderOfpbarProcess = 0
static G4ThreadLocal G4int	CounterOfpProcess = 0
static G4ThreadLocal G4int	CounterOfpbarProcess = 0
static G4ThreadLocal G4double	ParticleMass
static G4ThreadLocal G4double	ptableElectronCutInRange = 0.0
static G4ThreadLocal G4double	pbartableElectronCutInRange = 0.0
static G4ThreadLocal G4double	Charge
static G4ThreadLocal G4double	LowestKineticEnergy = 1e-05
static G4ThreadLocal G4double	HighestKineticEnergy = 1.e5
static G4ThreadLocal G4int	TotBin = 360
static G4ThreadLocal G4double	RTable
static G4ThreadLocal G4double	LOGRTable
static G4ThreadLocal G4double	dRoverRange = 0.20
static G4ThreadLocal G4double	finalRange = 0.2
static G4ThreadLocal G4double	c1lim = 0.20
static G4ThreadLocal G4double	c2lim = 0.32
static G4ThreadLocal G4double	c3lim = -0.032
static G4ThreadLocal G4bool	rndmStepFlag = false
static G4ThreadLocal G4bool	EnlossFlucFlag = true

Detailed Description

Definition at line 93 of file G4hRDEnergyLoss.hh.

Constructor & Destructor Documentation

G4hRDEnergyLoss::G4hRDEnergyLoss

(

const G4String &

processName

)

Definition at line 162 of file G4hRDEnergyLoss.cc.

  : G4VContinuousDiscreteProcess (processName),
    MaxExcitationNumber (1.e6),
    probLimFluct (0.01),
    nmaxDirectFluct (100),
    nmaxCont1(4),
    nmaxCont2(16),
    theLossTable(0),
    linLossLimit(0.05),
    MinKineticEnergy(0.0) 
{
  if (!RecorderOfpbarProcess) RecorderOfpbarProcess = new G4PhysicsTable*[100];
  if (!RecorderOfpProcess) RecorderOfpProcess = new G4PhysicsTable*[100];
  if (!RecorderOfProcess) RecorderOfProcess = new G4PhysicsTable*[100];
}

G4hRDEnergyLoss::~G4hRDEnergyLoss

(

)

Definition at line 180 of file G4hRDEnergyLoss.cc.

{
  if(theLossTable)
  {
    theLossTable->clearAndDestroy();
    delete theLossTable;
  }
}

Here is the call graph for this function:

Member Function Documentation

void G4hRDEnergyLoss::BuildDEDXTable ( const G4ParticleDefinition &  aParticleType )

staticprotected

Definition at line 251 of file G4hRDEnergyLoss.cc.

{
  //  calculate data members TotBin,LOGRTable,RTable first

  if (!RecorderOfpbarProcess) RecorderOfpbarProcess = new G4PhysicsTable*[100];
  if (!RecorderOfpProcess) RecorderOfpProcess = new G4PhysicsTable*[100];
  if (!RecorderOfProcess) RecorderOfProcess = new G4PhysicsTable*[100];

  const G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();

  // create/fill proton or antiproton tables depending on the charge
  Charge = aParticleType.GetPDGCharge()/eplus;
  ParticleMass = aParticleType.GetPDGMass() ;

  if (Charge>0.) {theDEDXTable= theDEDXpTable;}
  else           {theDEDXTable= theDEDXpbarTable;}

  if( ((Charge>0.) && (theDEDXTable==0)) ||
      ((Charge<0.) && (theDEDXTable==0)) 
      )
    {

      // Build energy loss table as a sum of the energy loss due to the
      //              different processes.
      if( Charge >0.)
    {
      RecorderOfProcess=RecorderOfpProcess;
      CounterOfProcess=CounterOfpProcess;

      if(CounterOfProcess == NumberOfProcesses)
        {
          if(theDEDXpTable)
        { theDEDXpTable->clearAndDestroy();
          delete theDEDXpTable; }
          theDEDXpTable = new G4PhysicsTable(numOfCouples);
          theDEDXTable = theDEDXpTable;
        }
    }
      else
    {
      RecorderOfProcess=RecorderOfpbarProcess;
      CounterOfProcess=CounterOfpbarProcess;

      if(CounterOfProcess == NumberOfProcesses)
        {
          if(theDEDXpbarTable)
        { theDEDXpbarTable->clearAndDestroy();
          delete theDEDXpbarTable; }
          theDEDXpbarTable = new G4PhysicsTable(numOfCouples);
          theDEDXTable = theDEDXpbarTable;
        }
    }

      if(CounterOfProcess == NumberOfProcesses)
    {
      //  loop for materials
      G4double LowEdgeEnergy , Value ;
      G4bool isOutRange ;
      G4PhysicsTable* pointer ;

      for (size_t J=0; J<numOfCouples; J++)
        {
          // create physics vector and fill it
          G4PhysicsLogVector* aVector =
                new G4PhysicsLogVector(LowestKineticEnergy,
                                       HighestKineticEnergy, TotBin);

          // loop for the kinetic energy
          for (G4int i=0; i<TotBin; i++)
        {
          LowEdgeEnergy = aVector->GetLowEdgeEnergy(i) ;
          Value = 0. ;

          // loop for the contributing processes
          for (G4int process=0; process < NumberOfProcesses; process++)
            {
              pointer= RecorderOfProcess[process];
              Value += (*pointer)[J]->
            GetValue(LowEdgeEnergy,isOutRange) ;
            }

          aVector->PutValue(i,Value) ;
        }

          theDEDXTable->insert(aVector) ;
        }

      //  reset counter to zero ..................
      if( Charge >0.)
        CounterOfpProcess=0 ;
      else
        CounterOfpbarProcess=0 ;

      // Build range table
      BuildRangeTable( aParticleType);

      // Build lab/proper time tables
      BuildTimeTables( aParticleType) ;

      // Build coeff tables for the energy loss calculation
      BuildRangeCoeffATable( aParticleType);
      BuildRangeCoeffBTable( aParticleType);
      BuildRangeCoeffCTable( aParticleType);

      // invert the range table

      BuildInverseRangeTable(aParticleType);
    }
    }
  // make the energy loss and the range table available

  G4EnergyLossTables::Register(&aParticleType,
                   (Charge>0)?
                   theDEDXpTable: theDEDXpbarTable,
                   (Charge>0)?
                   theRangepTable: theRangepbarTable,
                   (Charge>0)?
                   theInverseRangepTable: theInverseRangepbarTable,
                   (Charge>0)?
                   theLabTimepTable: theLabTimepbarTable,
                   (Charge>0)?
                   theProperTimepTable: theProperTimepbarTable,
                   LowestKineticEnergy, HighestKineticEnergy,
                   proton_mass_c2/aParticleType.GetPDGMass(),
                   TotBin);

}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4hRDEnergyLoss::CutsWhereModified ( )

protected

Definition at line 1155 of file G4hRDEnergyLoss.cc.

{
  G4bool wasModified = false;
  const G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();

  for (size_t j=0; j<numOfCouples; j++){
    if (theCoupleTable->GetMaterialCutsCouple(j)->IsRecalcNeeded()) {
      wasModified = true;
      break;
    }
  }
  return wasModified;
}

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4double G4hRDEnergyLoss::GetMeanFreePath ( const G4Track &  track, G4double  previousStepSize, enum G4ForceCondition *  condition  )

pure virtual

Implements G4VContinuousDiscreteProcess.

Implemented in G4hImpactIonisation.

G4int G4hRDEnergyLoss::GetNumberOfProcesses ( )

static

Definition at line 191 of file G4hRDEnergyLoss.cc.

{
  return NumberOfProcesses; 
} 

void G4hRDEnergyLoss::MinusNumberOfProcesses ( )

static

Definition at line 212 of file G4hRDEnergyLoss.cc.

{
  NumberOfProcesses--; 
} 

void G4hRDEnergyLoss::PlusNumberOfProcesses ( )

static

Definition at line 205 of file G4hRDEnergyLoss.cc.

{
  NumberOfProcesses++; 
} 

virtual G4VParticleChange* G4hRDEnergyLoss::PostStepDoIt ( const G4Track &  track, const G4Step &  Step  )

pure virtual

Reimplemented from G4VContinuousDiscreteProcess.

Implemented in G4hImpactIonisation.

void G4hRDEnergyLoss::SetdRoverRange ( G4double  value )

static

Definition at line 219 of file G4hRDEnergyLoss.cc.

{
  dRoverRange = value;
}

void G4hRDEnergyLoss::SetEnlossFluc ( G4bool  value )

static

Definition at line 233 of file G4hRDEnergyLoss.cc.

{
  EnlossFlucFlag = value;
}

void G4hRDEnergyLoss::SetNumberOfProcesses ( G4int  number )

static

Definition at line 198 of file G4hRDEnergyLoss.cc.

{
  NumberOfProcesses=number; 
} 

void G4hRDEnergyLoss::SetRndmStep ( G4bool  value )

static

Definition at line 226 of file G4hRDEnergyLoss.cc.

{
  rndmStepFlag = value;
}

void G4hRDEnergyLoss::SetStepFunction ( G4double  c1, G4double  c2  )

static

Definition at line 240 of file G4hRDEnergyLoss.cc.

{
  dRoverRange = c1; 
  finalRange = c2;
  c1lim=dRoverRange;
  c2lim=2.*(1-dRoverRange)*finalRange;
  c3lim=-(1.-dRoverRange)*finalRange*finalRange;
}

Member Data Documentation

G4ThreadLocal G4double G4hRDEnergyLoss::c1lim = 0.20

staticprotected

Definition at line 194 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::c2lim = 0.32

staticprotected

Definition at line 194 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::c3lim = -0.032

staticprotected

Definition at line 194 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::Charge

staticprotected

Definition at line 174 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4int G4hRDEnergyLoss::CounterOfpbarProcess = 0

staticprotected

Definition at line 165 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4int G4hRDEnergyLoss::CounterOfpProcess = 0

staticprotected

Definition at line 164 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::dRoverRange = 0.20

staticprotected

Definition at line 191 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4bool G4hRDEnergyLoss::EnlossFlucFlag = true

staticprotected

Definition at line 197 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::finalRange = 0.2

staticprotected

Definition at line 193 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::HighestKineticEnergy = 1.e5

staticprotected

Definition at line 177 of file G4hRDEnergyLoss.hh.

G4double G4hRDEnergyLoss::linLossLimit

protected

Definition at line 187 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::LOGRTable

staticprotected

Definition at line 181 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::LowestKineticEnergy = 1e-05

staticprotected

Definition at line 176 of file G4hRDEnergyLoss.hh.

const G4double G4hRDEnergyLoss::MaxExcitationNumber

protected

Definition at line 136 of file G4hRDEnergyLoss.hh.

G4double G4hRDEnergyLoss::MinKineticEnergy

protected

Definition at line 189 of file G4hRDEnergyLoss.hh.

const long G4hRDEnergyLoss::nmaxCont1

protected

Definition at line 138 of file G4hRDEnergyLoss.hh.

const long G4hRDEnergyLoss::nmaxCont2

protected

Definition at line 138 of file G4hRDEnergyLoss.hh.

const long G4hRDEnergyLoss::nmaxDirectFluct

protected

Definition at line 138 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::ParticleMass

staticprotected

Definition at line 168 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::pbartableElectronCutInRange = 0.0

staticprotected

Definition at line 172 of file G4hRDEnergyLoss.hh.

const G4double G4hRDEnergyLoss::probLimFluct

protected

Definition at line 137 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::ptableElectronCutInRange = 0.0

staticprotected

Definition at line 171 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable ** G4hRDEnergyLoss::RecorderOfpbarProcess = 0

staticprotected

Definition at line 163 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable ** G4hRDEnergyLoss::RecorderOfpProcess = 0

staticprotected

Definition at line 162 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4bool G4hRDEnergyLoss::rndmStepFlag = false

staticprotected

Definition at line 196 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4double G4hRDEnergyLoss::RTable

staticprotected

Definition at line 181 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theDEDXpbarTable = 0

staticprotected

Definition at line 145 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theDEDXpTable = 0

staticprotected

Definition at line 144 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theInverseRangepbarTable = 0

staticprotected

Definition at line 151 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theInverseRangepTable = 0

staticprotected

Definition at line 150 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theLabTimepbarTable = 0

staticprotected

Definition at line 155 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theLabTimepTable = 0

staticprotected

Definition at line 154 of file G4hRDEnergyLoss.hh.

G4PhysicsTable* G4hRDEnergyLoss::theLossTable

protected

Definition at line 185 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theProperTimepbarTable = 0

staticprotected

Definition at line 158 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theProperTimepTable = 0

staticprotected

Definition at line 157 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theRangepbarTable = 0

staticprotected

Definition at line 147 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4PhysicsTable * G4hRDEnergyLoss::theRangepTable = 0

staticprotected

Definition at line 146 of file G4hRDEnergyLoss.hh.

G4ThreadLocal G4int G4hRDEnergyLoss::TotBin = 360

staticprotected

Definition at line 178 of file G4hRDEnergyLoss.hh.

---

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

• source/geant4.10.03.p03/source/processes/electromagnetic/pii/include/G4hRDEnergyLoss.hh
• source/geant4.10.03.p03/source/processes/electromagnetic/pii/src/G4hRDEnergyLoss.cc