#include <G4EnergyLossForExtrapolator.hh>

Collaboration diagram for G4EnergyLossForExtrapolator:

Public Member Functions
	G4EnergyLossForExtrapolator (G4int verb=1)

	~G4EnergyLossForExtrapolator ()

G4double	ComputeDEDX (G4double kinEnergy, const G4ParticleDefinition *)

G4double	ComputeRange (G4double kinEnergy, const G4ParticleDefinition *)

G4double	ComputeEnergy (G4double range, const G4ParticleDefinition *)

G4double	EnergyAfterStep (G4double kinEnergy, G4double step, const G4Material , const G4ParticleDefinition )

G4double	EnergyBeforeStep (G4double kinEnergy, G4double step, const G4Material , const G4ParticleDefinition )

G4double	TrueStepLength (G4double kinEnergy, G4double step, const G4Material , const G4ParticleDefinition part)

G4double	EnergyAfterStep (G4double kinEnergy, G4double step, const G4Material *, const G4String &particleName)

G4double	EnergyBeforeStep (G4double kinEnergy, G4double step, const G4Material *, const G4String &particleName)

G4double	AverageScatteringAngle (G4double kinEnergy, G4double step, const G4Material , const G4ParticleDefinition part)

G4double	AverageScatteringAngle (G4double kinEnergy, G4double step, const G4Material *, const G4String &particleName)

G4double	ComputeTrueStep (const G4Material , const G4ParticleDefinition part, G4double kinEnergy, G4double stepLength)

G4double	EnergyDispersion (G4double kinEnergy, G4double step, const G4Material , const G4ParticleDefinition )

G4double	EnergyDispersion (G4double kinEnergy, G4double step, const G4Material *, const G4String &particleName)

void	SetVerbose (G4int val)

void	SetMinKinEnergy (G4double)

void	SetMaxKinEnergy (G4double)

void	SetMaxEnergyTransfer (G4double)

Private Member Functions
void	Initialisation ()

void	BuildTables ()

G4bool	SetupKinematics (const G4ParticleDefinition , const G4Material , G4double kinEnergy)

const G4ParticleDefinition *	FindParticle (const G4String &name)

G4double	ComputeValue (G4double x, const G4PhysicsTable *table, size_t idx)

const G4PhysicsTable *	GetPhysicsTable (ExtTableType type) const

G4EnergyLossForExtrapolator &	operator= (const G4EnergyLossForExtrapolator &right)

	G4EnergyLossForExtrapolator (const G4EnergyLossForExtrapolator &)

Private Attributes
const G4ParticleDefinition *	currentParticle

const G4ParticleDefinition *	electron

const G4ParticleDefinition *	positron

const G4ParticleDefinition *	muonPlus

const G4ParticleDefinition *	muonMinus

const G4ParticleDefinition *	proton

G4String	currentParticleName

size_t	idxDedxElectron

size_t	idxDedxPositron

size_t	idxDedxMuon

size_t	idxDedxProton

size_t	idxRangeElectron

size_t	idxRangePositron

size_t	idxRangeMuon

size_t	idxRangeProton

size_t	idxInvRangeElectron

size_t	idxInvRangePositron

size_t	idxInvRangeMuon

size_t	idxInvRangeProton

size_t	idxMscElectron

const G4Material *	currentMaterial

G4int	index

G4double	electronDensity

G4double	radLength

G4double	mass

G4double	charge2

G4double	kineticEnergy

G4double	gam

G4double	bg2

G4double	beta2

G4double	tmax

G4double	linLossLimit

G4double	emin

G4double	emax

G4double	maxEnergyTransfer

G4int	nbins

G4int	nmat

G4int	verbose

Static Private Attributes
static G4TablesForExtrapolator *	tables = 0

Detailed Description

Definition at line 67 of file G4EnergyLossForExtrapolator.hh.

Constructor & Destructor Documentation

◆ G4EnergyLossForExtrapolator() [1/2]

G4EnergyLossForExtrapolator::G4EnergyLossForExtrapolator ( G4int verb = 1 )

Definition at line 67 of file G4EnergyLossForExtrapolator.cc.

   : maxEnergyTransfer(DBL_MAX), verbose(verb)
 {
   currentParticle = 0;
   currentMaterial = 0;
 
   linLossLimit = 0.01;
   emin         = 1.*MeV;
   emax         = 10.*TeV;
   nbins        = 70;
 
   nmat = index = 0;
 
   mass = charge2 = electronDensity = radLength = bg2 = beta2 
     = kineticEnergy = tmax = 0.0;
   gam = 1.0;
 
   idxDedxElectron = idxDedxPositron = idxDedxMuon = idxDedxProton 
     = idxRangeElectron = idxRangePositron = idxRangeMuon = idxRangeProton
     = idxInvRangeElectron = idxInvRangePositron = idxInvRangeMuon
     = idxInvRangeProton = idxMscElectron = 0;
  
   electron = positron = proton = muonPlus = muonMinus = 0;
 }

◆ ~G4EnergyLossForExtrapolator()

G4EnergyLossForExtrapolator::~G4EnergyLossForExtrapolator ( )

Definition at line 346 of file G4EnergyLossForExtrapolator.cc.

 {
   G4AutoLock l(&G4EnergyLossForExtrapolatorMutex);
   delete tables;
   tables = 0;
 }

◆ G4EnergyLossForExtrapolator() [2/2]

G4EnergyLossForExtrapolator::G4EnergyLossForExtrapolator ( const G4EnergyLossForExtrapolator & )

private

Member Function Documentation

◆ AverageScatteringAngle() [1/2]

G4double G4EnergyLossForExtrapolator::AverageScatteringAngle	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4ParticleDefinition *	part
	)

inline

Definition at line 248 of file G4EnergyLossForExtrapolator.hh.

 {
   G4double theta = 0.0;
   if(SetupKinematics(part, mat, kinEnergy)) {
     G4double t = stepLength/radLength;
     G4double y = std::max(0.001, t); 
     theta = 19.23*CLHEP::MeV*std::sqrt(charge2*t)*(1.0 + 0.038*G4Log(y))
       /(beta2*gam*mass);
   }
   return theta;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ AverageScatteringAngle() [2/2]

G4double G4EnergyLossForExtrapolator::AverageScatteringAngle	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4String &	particleName
	)

inline

Definition at line 227 of file G4EnergyLossForExtrapolator.hh.

 {
   return AverageScatteringAngle(kinEnergy,step,mat,FindParticle(name));
 }

Here is the call graph for this function:

◆ BuildTables()

void G4EnergyLossForExtrapolator::BuildTables ( )

private

Definition at line 330 of file G4EnergyLossForExtrapolator.cc.

 {
   G4AutoLock l(&G4EnergyLossForExtrapolatorMutex);
 
   if(!tables) {
     if(verbose > 0) {
       G4cout << "### G4EnergyLossForExtrapolator::BuildTables for "
          << nmat << " materials Nbins= " << nbins 
          << " Emin(MeV)= " << emin << "  Emax(MeV)= " << emax << G4endl;
     }
     tables = new G4TablesForExtrapolator(verbose, nbins, emin, emax);
   }
 }

Here is the caller graph for this function:

◆ ComputeDEDX()

G4double G4EnergyLossForExtrapolator::ComputeDEDX	(	G4double	kinEnergy,
		const G4ParticleDefinition *	part
	)

Definition at line 240 of file G4EnergyLossForExtrapolator.cc.

 {
   G4double x = 0.0;
   if(part == electron)  { 
     x = ComputeValue(ekin, GetPhysicsTable(fDedxElectron), idxDedxElectron);
   } else if(part == positron) {
     x = ComputeValue(ekin, GetPhysicsTable(fDedxPositron), idxDedxPositron);
   } else if(part == muonPlus || part == muonMinus) {
     x = ComputeValue(ekin, GetPhysicsTable(fDedxMuon), idxDedxMuon);
   } else {
     G4double e = ekin*proton_mass_c2/mass;
     x = ComputeValue(e, GetPhysicsTable(fDedxProton), idxDedxProton)*charge2;
   }
   return x;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeEnergy()

G4double G4EnergyLossForExtrapolator::ComputeEnergy	(	G4double	range,
		const G4ParticleDefinition *	part
	)

Definition at line 282 of file G4EnergyLossForExtrapolator.cc.

 {
   G4double x = 0.0;
   if(part == electron) {
     x = ComputeValue(range, GetPhysicsTable(fInvRangeElectron), 
              idxInvRangeElectron);
   } else if(part == positron) {
     x = ComputeValue(range, GetPhysicsTable(fInvRangePositron), 
              idxInvRangePositron);
   } else if(part == muonPlus || part == muonMinus) {
     x = ComputeValue(range, GetPhysicsTable(fInvRangeMuon), idxInvRangeMuon);
   } else {
     G4double massratio = proton_mass_c2/mass;
     G4double r = range*massratio*charge2;
     x = ComputeValue(r, GetPhysicsTable(fInvRangeProton), 
              idxInvRangeProton)/massratio;
   }
   return x;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeRange()

G4double G4EnergyLossForExtrapolator::ComputeRange	(	G4double	kinEnergy,
		const G4ParticleDefinition *	part
	)

Definition at line 260 of file G4EnergyLossForExtrapolator.cc.

 {
   G4double x = 0.0;
   if(part == electron) { 
     x = ComputeValue(ekin, GetPhysicsTable(fRangeElectron), idxRangeElectron);
   } else if(part == positron) {
     x = ComputeValue(ekin, GetPhysicsTable(fRangePositron), idxRangePositron);
   } else if(part == muonPlus || part == muonMinus) { 
     x = ComputeValue(ekin, GetPhysicsTable(fRangeMuon), idxRangeMuon);
   } else {
     G4double massratio = proton_mass_c2/mass;
     G4double e = ekin*massratio;
     x = ComputeValue(e, GetPhysicsTable(fRangeProton), idxRangeProton)
       /(charge2*massratio);
   }
   return x;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeTrueStep()

G4double G4EnergyLossForExtrapolator::ComputeTrueStep	(	const G4Material *	mat,
		const G4ParticleDefinition *	part,
		G4double	kinEnergy,
		G4double	stepLength
	)

inline

Definition at line 267 of file G4EnergyLossForExtrapolator.hh.

 {
   G4double theta = AverageScatteringAngle(kinEnergy,stepLength,mat,part);
   return stepLength*std::sqrt(1.0 + 0.625*theta*theta);
 } 

Here is the call graph for this function:

Here is the caller graph for this function:

◆ ComputeValue()

G4double G4EnergyLossForExtrapolator::ComputeValue	(	G4double	x,
		const G4PhysicsTable *	table,
		size_t	idx
	)

inlineprivate

Definition at line 296 of file G4EnergyLossForExtrapolator.hh.

 {
   G4double res = 0.0;
   if(table) { res = ((*table)[index])->Value(x, idx); }
   return res;
 }

Here is the caller graph for this function:

◆ EnergyAfterStep() [1/2]

G4double G4EnergyLossForExtrapolator::EnergyAfterStep	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4ParticleDefinition *	part
	)

Definition at line 95 of file G4EnergyLossForExtrapolator.cc.

 {
   if(0 == nmat) { Initialisation(); }
   G4double kinEnergyFinal = kinEnergy;
   if(SetupKinematics(part, mat, kinEnergy)) {
     G4double step = TrueStepLength(kinEnergy,stepLength,mat,part);
     G4double r  = ComputeRange(kinEnergy,part);
     if(r <= step) {
       kinEnergyFinal = 0.0;
     } else if(step < linLossLimit*r) {
       kinEnergyFinal -= step*ComputeDEDX(kinEnergy,part);
     } else {  
       G4double r1 = r - step;
       kinEnergyFinal = ComputeEnergy(r1,part);
     }
   }
   return kinEnergyFinal;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ EnergyAfterStep() [2/2]

G4double G4EnergyLossForExtrapolator::EnergyAfterStep	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4String &	particleName
	)

inline

Definition at line 205 of file G4EnergyLossForExtrapolator.hh.

 {
   return EnergyAfterStep(kinEnergy,step,mat,FindParticle(name));
 }

Here is the call graph for this function:

◆ EnergyBeforeStep() [1/2]

G4double G4EnergyLossForExtrapolator::EnergyBeforeStep	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4ParticleDefinition *	part
	)

Definition at line 120 of file G4EnergyLossForExtrapolator.cc.

 {
   //G4cout << "G4EnergyLossForExtrapolator::EnergyBeforeStep" << G4endl;
   if(0 == nmat) { Initialisation(); }
   G4double kinEnergyFinal = kinEnergy;
 
   if(SetupKinematics(part, mat, kinEnergy)) {
     G4double step = TrueStepLength(kinEnergy,stepLength,mat,part);
     G4double r  = ComputeRange(kinEnergy,part);
 
     if(step < linLossLimit*r) {
       kinEnergyFinal += step*ComputeDEDX(kinEnergy,part);
     } else {  
       G4double r1 = r + step;
       kinEnergyFinal = ComputeEnergy(r1,part);
     }
   }
   return kinEnergyFinal;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ EnergyBeforeStep() [2/2]

G4double G4EnergyLossForExtrapolator::EnergyBeforeStep	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4String &	particleName
	)

inline

Definition at line 216 of file G4EnergyLossForExtrapolator.hh.

 {
   return EnergyBeforeStep(kinEnergy,step,mat,FindParticle(name));
 }

Here is the call graph for this function:

◆ EnergyDispersion() [1/2]

G4double G4EnergyLossForExtrapolator::EnergyDispersion	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4ParticleDefinition *	part
	)

inline

Definition at line 279 of file G4EnergyLossForExtrapolator.hh.

 {
   G4double sig2 = 0.0;
   if(SetupKinematics(part, mat, kinEnergy)) {
     G4double step = ComputeTrueStep(mat,part,kinEnergy,stepLength);
     sig2 = (1.0/beta2 - 0.5)
       *CLHEP::twopi_mc2_rcl2*tmax*step*electronDensity*charge2;
   }
   return sig2;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ EnergyDispersion() [2/2]

G4double G4EnergyLossForExtrapolator::EnergyDispersion	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4String &	particleName
	)

inline

Definition at line 238 of file G4EnergyLossForExtrapolator.hh.

 {
   return EnergyDispersion(kinEnergy,step,mat,FindParticle(name));
 }

Here is the call graph for this function:

◆ FindParticle()

const G4ParticleDefinition * G4EnergyLossForExtrapolator::FindParticle ( const G4String & name )

private

Definition at line 221 of file G4EnergyLossForExtrapolator.cc.

 {
   const G4ParticleDefinition* p = 0;
   if(name != currentParticleName) {
     p = G4ParticleTable::GetParticleTable()->FindParticle(name);
     if(!p) {
       G4cout << "### G4EnergyLossForExtrapolator WARNING: "
          << "FindParticle() fails to find " 
          << name << G4endl;
     }
   } else {
     p = currentParticle;
   }
   return p;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetPhysicsTable()

const G4PhysicsTable * G4EnergyLossForExtrapolator::GetPhysicsTable ( ExtTableType type ) const

inlineprivate

Definition at line 197 of file G4EnergyLossForExtrapolator.hh.

 {
   return tables->GetPhysicsTable(type);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ Initialisation()

void G4EnergyLossForExtrapolator::Initialisation ( )

private

Definition at line 305 of file G4EnergyLossForExtrapolator.cc.

 {
   if(verbose>1) {
     G4cout << "### G4EnergyLossForExtrapolator::Initialisation" << G4endl;
   }
   currentParticle = 0;
   currentMaterial = 0;
   kineticEnergy   = 0.0;
 
   electron = G4Electron::Electron();
   positron = G4Positron::Positron();
   proton   = G4Proton::Proton();
   muonPlus = G4MuonPlus::MuonPlus();
   muonMinus= G4MuonMinus::MuonMinus();
 
   nmat = G4Material::GetNumberOfMaterials();
   currentParticleName = "";
   BuildTables(); 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ operator=()

G4EnergyLossForExtrapolator& G4EnergyLossForExtrapolator::operator= ( const G4EnergyLossForExtrapolator & right )

private

◆ SetMaxEnergyTransfer()

void G4EnergyLossForExtrapolator::SetMaxEnergyTransfer ( G4double val )

inline

Definition at line 328 of file G4EnergyLossForExtrapolator.hh.

 {
   maxEnergyTransfer = val;
 }

◆ SetMaxKinEnergy()

void G4EnergyLossForExtrapolator::SetMaxKinEnergy ( G4double val )

inline

Definition at line 321 of file G4EnergyLossForExtrapolator.hh.

 {
   emax = val;
 }

◆ SetMinKinEnergy()

void G4EnergyLossForExtrapolator::SetMinKinEnergy ( G4double val )

inline

Definition at line 314 of file G4EnergyLossForExtrapolator.hh.

 {
   emin = val;
 }

◆ SetupKinematics()

G4bool G4EnergyLossForExtrapolator::SetupKinematics	(	const G4ParticleDefinition *	part,
		const G4Material *	mat,
		G4double	kinEnergy
	)

private

Definition at line 173 of file G4EnergyLossForExtrapolator.cc.

 {
   if(0 == nmat) { Initialisation(); }
   if(!part || !mat || kinEnergy < keV) { return false; }
   G4bool flag = false;
   if(part != currentParticle) {
     flag = true;
     currentParticle = part;
     mass = part->GetPDGMass();
     G4double q = part->GetPDGCharge()/eplus;
     charge2 = q*q;
   }
   if(mat != currentMaterial) {
     G4int i = mat->GetIndex();
     if(i >= nmat) {
       G4cout << "### G4EnergyLossForExtrapolator WARNING:index i= " 
          << i << " is out of table - NO extrapolation" << G4endl;
     } else {
       flag = true;
       currentMaterial = mat;
       electronDensity = mat->GetElectronDensity();
       radLength       = mat->GetRadlen();
       index           = i;
     }
   }
   if(flag || kinEnergy != kineticEnergy) {
     kineticEnergy = kinEnergy;
     G4double tau  = kinEnergy/mass;
 
     gam   = tau + 1.0;
     bg2   = tau * (tau + 2.0);
     beta2 = bg2/(gam*gam);
     tmax  = kinEnergy;
     if(part == electron) tmax *= 0.5;
     else if(part != positron) {
       G4double r = electron_mass_c2/mass;
       tmax = 2.0*bg2*electron_mass_c2/(1.0 + 2.0*gam*r + r*r);
     }
     if(tmax > maxEnergyTransfer) { tmax = maxEnergyTransfer; }
   }
   return true;
 } 

Here is the call graph for this function:

Here is the caller graph for this function:

◆ SetVerbose()

void G4EnergyLossForExtrapolator::SetVerbose ( G4int val )

inline

Definition at line 307 of file G4EnergyLossForExtrapolator.hh.

 {
   verbose = val;
 }

◆ TrueStepLength()

G4double G4EnergyLossForExtrapolator::TrueStepLength	(	G4double	kinEnergy,
		G4double	step,
		const G4Material *	mat,
		const G4ParticleDefinition *	part
	)

Definition at line 146 of file G4EnergyLossForExtrapolator.cc.

 {
   G4double res = stepLength;
   if(0 == nmat) { Initialisation(); }
   //G4cout << "## G4EnergyLossForExtrapolator::TrueStepLength L= " << res 
   //     <<  "  " << part->GetParticleName() << G4endl;
   if(SetupKinematics(part, mat, kinEnergy)) {
     if(part == electron || part == positron) {
       G4double x = stepLength*
     ComputeValue(kinEnergy, GetPhysicsTable(fMscElectron), idxMscElectron);
       //G4cout << " x= " << x << G4endl;
       if(x < 0.2)         { res *= (1.0 + 0.5*x + x*x/3.0); }
       else if(x < 0.9999) { res = -G4Log(1.0 - x)*stepLength/x; }
       else                { res = ComputeRange(kinEnergy,part); }
     } else {
       res = ComputeTrueStep(mat,part,kinEnergy,stepLength);
     }
   }
   return res;
 }