#include <G4hRDEnergyLoss.hh>

Inheritance 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.

References RecorderOfpbarProcess, and RecorderOfpProcess.

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

References G4PhysicsTable::clearAndDestroy(), and theLossTable.

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

Member Function Documentation

void G4hRDEnergyLoss::BuildDEDXTable ( const G4ParticleDefinition & aParticleType )

staticprotected

Definition at line 251 of file G4hRDEnergyLoss.cc.

References Charge, G4PhysicsTable::clearAndDestroy(), CounterOfpbarProcess, CounterOfpProcess, python.hepunit::eplus, G4PhysicsVector::GetLowEdgeEnergy(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), HighestKineticEnergy, G4PhysicsTable::insert(), LowestKineticEnergy, ParticleMass, python.hepunit::proton_mass_c2, G4PhysicsVector::PutValue(), RecorderOfpbarProcess, RecorderOfpProcess, G4EnergyLossTables::Register(), theDEDXpbarTable, theDEDXpTable, theInverseRangepbarTable, theInverseRangepTable, theLabTimepbarTable, theLabTimepTable, theProperTimepbarTable, theProperTimepTable, theRangepbarTable, theRangepTable, and TotBin.

Referenced by G4hImpactIonisation::BuildPhysicsTable().

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

G4bool G4hRDEnergyLoss::CutsWhereModified ( )

protected

Definition at line 1155 of file G4hRDEnergyLoss.cc.

References G4ProductionCutsTable::GetMaterialCutsCouple(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), and G4MaterialCutsCouple::IsRecalcNeeded().

Referenced by G4hImpactIonisation::BuildPhysicsTable().

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

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.

References dRoverRange.

 {
   dRoverRange = value;
 }

void G4hRDEnergyLoss::SetEnlossFluc ( G4bool value )

static

Definition at line 233 of file G4hRDEnergyLoss.cc.

References EnlossFlucFlag.

 {
   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.

References rndmStepFlag.

 {
   rndmStepFlag = value;
 }

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

static

Definition at line 240 of file G4hRDEnergyLoss.cc.

References plottest35::c1, c1lim, c2lim, c3lim, dRoverRange, and finalRange.

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