#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:

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation