G4EnergyLossForExtrapolator Class Reference

#include <G4EnergyLossForExtrapolator.hh>

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)

Detailed Description

Definition at line 68 of file G4EnergyLossForExtrapolator.hh.

Constructor & Destructor Documentation

G4EnergyLossForExtrapolator::G4EnergyLossForExtrapolator

(

G4int

verb = 1

)

Definition at line 75 of file G4EnergyLossForExtrapolator.cc.

References python.hepunit::MeV, and python.hepunit::TeV.

  :maxEnergyTransfer(DBL_MAX),verbose(verb),isInitialised(false)
{
  currentParticle = 0;
  currentMaterial = 0;

  linLossLimit = 0.01;
  emin         = 1.*MeV;
  emax         = 10.*TeV;
  nbins        = 70;

  nmat = index = 0;
  pcuts = 0;

  mass = charge2 = electronDensity = radLength = bg2 = beta2 
    = kineticEnergy = tmax = 0;
  gam = 1.0;

  dedxElectron = dedxPositron = dedxProton = rangeElectron 
    = rangePositron = rangeProton = invRangeElectron = invRangePositron 
    = invRangeProton = mscElectron = dedxMuon = rangeMuon = invRangeMuon = 0;
 
  electron = positron = proton = muonPlus = muonMinus = 0;
}

G4EnergyLossForExtrapolator::~G4EnergyLossForExtrapolator

(

)

Definition at line 102 of file G4EnergyLossForExtrapolator.cc.

{
  for(G4int i=0; i<nmat; i++) {delete couples[i];}
  delete dedxElectron;
  delete dedxPositron;
  delete dedxProton;
  delete dedxMuon;
  delete rangeElectron;
  delete rangePositron;
  delete rangeProton;
  delete rangeMuon;
  delete invRangeElectron;
  delete invRangePositron;
  delete invRangeProton;
  delete invRangeMuon;
  delete mscElectron;
  delete pcuts;
}

Member Function Documentation

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

inline

Definition at line 250 of file G4EnergyLossForExtrapolator.hh.

References G4Log(), and G4INCL::Math::max().

Referenced by AverageScatteringAngle(), and ComputeTrueStep().

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

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

inline

Definition at line 228 of file G4EnergyLossForExtrapolator.hh.

References AverageScatteringAngle().

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

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

Definition at line 360 of file G4EnergyLossForExtrapolator.cc.

References python.hepunit::proton_mass_c2, and test::x.

Referenced by EnergyAfterStep(), and EnergyBeforeStep().

{
  G4double x = 0.0;
  if(part == electron)      x = ComputeValue(kinEnergy, dedxElectron);
  else if(part == positron) x = ComputeValue(kinEnergy, dedxPositron);
  else if(part == muonPlus || part == muonMinus) {
    x = ComputeValue(kinEnergy, dedxMuon);
  } else {
    G4double e = kinEnergy*proton_mass_c2/mass;
    x = ComputeValue(e, dedxProton)*charge2;
  }
  return x;
}

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

Definition at line 397 of file G4EnergyLossForExtrapolator.cc.

References python.hepunit::proton_mass_c2, and test::x.

Referenced by EnergyAfterStep(), and EnergyBeforeStep().

{
  G4double x = 0.0;
  if(part == electron)      x = ComputeValue(range, invRangeElectron);
  else if(part == positron) x = ComputeValue(range, invRangePositron);
  else if(part == muonPlus || part == muonMinus) 
    x = ComputeValue(range, invRangeMuon);
  else {
    G4double massratio = proton_mass_c2/mass;
    G4double r = range*massratio*charge2;
    x = ComputeValue(r, invRangeProton)/massratio;
  }
  return x;
}

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

Definition at line 378 of file G4EnergyLossForExtrapolator.cc.

References python.hepunit::proton_mass_c2, and test::x.

Referenced by EnergyAfterStep(), EnergyBeforeStep(), and TrueStepLength().

{
  G4double x = 0.0;
  if(part == electron)      x = ComputeValue(kinEnergy, rangeElectron);
  else if(part == positron) x = ComputeValue(kinEnergy, rangePositron);
  else if(part == muonPlus || part == muonMinus) 
    x = ComputeValue(kinEnergy, rangeMuon);
  else {
    G4double massratio = proton_mass_c2/mass;
    G4double e = kinEnergy*massratio;
    x = ComputeValue(e, rangeProton)/(charge2*massratio);
  }
  return x;
}

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

inline

Definition at line 268 of file G4EnergyLossForExtrapolator.hh.

References AverageScatteringAngle().

Referenced by EnergyDispersion(), and TrueStepLength().

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

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

Definition at line 124 of file G4EnergyLossForExtrapolator.cc.

References ComputeDEDX(), ComputeEnergy(), ComputeRange(), and TrueStepLength().

Referenced by G4ErrorEnergyLoss::AlongStepDoIt(), EnergyAfterStep(), and G4ErrorEnergyLoss::GetContinuousStepLimit().

{
  if(!isInitialised) 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;
}

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

inline

Definition at line 206 of file G4EnergyLossForExtrapolator.hh.

References EnergyAfterStep().

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

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

Definition at line 149 of file G4EnergyLossForExtrapolator.cc.

References ComputeDEDX(), ComputeEnergy(), ComputeRange(), and TrueStepLength().

Referenced by G4ErrorEnergyLoss::AlongStepDoIt(), EnergyBeforeStep(), and G4ErrorEnergyLoss::GetContinuousStepLimit().

{
  if(!isInitialised) { 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;
}

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

inline

Definition at line 217 of file G4EnergyLossForExtrapolator.hh.

References EnergyBeforeStep().

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

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

inline

Definition at line 280 of file G4EnergyLossForExtrapolator.hh.

References ComputeTrueStep().

Referenced by EnergyDispersion().

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

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

inline

Definition at line 239 of file G4EnergyLossForExtrapolator.hh.

References EnergyDispersion().

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

void G4EnergyLossForExtrapolator::SetMaxEnergyTransfer ( G4double  val )

inline

Definition at line 328 of file G4EnergyLossForExtrapolator.hh.

{
  maxEnergyTransfer = val;
}

void G4EnergyLossForExtrapolator::SetMaxKinEnergy ( G4double  val )

inline

Definition at line 321 of file G4EnergyLossForExtrapolator.hh.

{
  emax = val;
}

void G4EnergyLossForExtrapolator::SetMinKinEnergy ( G4double  val )

inline

Definition at line 314 of file G4EnergyLossForExtrapolator.hh.

{
  emin = val;
}

void G4EnergyLossForExtrapolator::SetVerbose ( G4int  val )

inline

Definition at line 307 of file G4EnergyLossForExtrapolator.hh.

{
  verbose = val;
}

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

Definition at line 174 of file G4EnergyLossForExtrapolator.cc.

References ComputeRange(), ComputeTrueStep(), G4Log(), and test::x.

Referenced by EnergyAfterStep(), EnergyBeforeStep(), and G4EnergySplitter::SplitEnergyInVolumes().

{
  G4double res = stepLength;
  if(!isInitialised) Initialisation();
  if(SetupKinematics(part, mat, kinEnergy)) {
    if(part == electron || part == positron) {
      G4double x = stepLength*ComputeValue(kinEnergy, mscElectron);
      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;
}

---

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

• G4EnergyLossForExtrapolator.hh
• G4EnergyLossForExtrapolator.cc