G4BraggModel Class Reference

#include <G4BraggModel.hh>

Inheritance diagram for G4BraggModel:

Public Member Functions
	G4BraggModel (const G4ParticleDefinition *p=0, const G4String &nam="Bragg")
virtual	~G4BraggModel ()
virtual void	Initialise (const G4ParticleDefinition *, const G4DataVector &)
virtual G4double	ComputeCrossSectionPerElectron (const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual G4double	ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy, G4double maxEnergy)
virtual G4double	CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual G4double	ComputeDEDXPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy)
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
virtual G4double	GetChargeSquareRatio (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
virtual G4double	GetParticleCharge (const G4ParticleDefinition *p, const G4Material *mat, G4double kineticEnergy)
Protected Member Functions
virtual G4double	MaxSecondaryEnergy (const G4ParticleDefinition *, G4double kinEnergy)
G4double	GetChargeSquareRatio () const
void	SetChargeSquareRatio (G4double val)

---

Detailed Description

Definition at line 71 of file G4BraggModel.hh.

---

Constructor & Destructor Documentation

G4BraggModel::G4BraggModel	(	const G4ParticleDefinition *	p = 0,
		const G4String &	nam = "Bragg"
	)

Definition at line 79 of file G4BraggModel.cc.

References G4Electron::Electron(), G4LossTableManager::EmCorrections(), G4LossTableManager::Instance(), and G4VEmModel::SetHighEnergyLimit().

  : G4VEmModel(nam),
    particle(0),
    currentMaterial(0),
    protonMassAMU(1.007276),
    iMolecula(-1),
    iPSTAR(-1),
    isIon(false),
    isInitialised(false)
{
  fParticleChange = 0;
  SetHighEnergyLimit(2.0*MeV);

  lowestKinEnergy  = 1.0*keV;
  theZieglerFactor = eV*cm2*1.0e-15;
  theElectron = G4Electron::Electron();
  expStopPower125 = 0.0;

  corr = G4LossTableManager::Instance()->EmCorrections();
  if(p) { SetParticle(p); }
  else  { SetParticle(theElectron); }
}

G4BraggModel::~G4BraggModel

(

)

[virtual]

Definition at line 104 of file G4BraggModel.cc.

{}

---

Member Function Documentation

G4double G4BraggModel::ComputeCrossSectionPerAtom	(	const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	Z,
		G4double	A,
		G4double	cutEnergy,
		G4double	maxEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 180 of file G4BraggModel.cc.

References ComputeCrossSectionPerElectron().

{
  G4double cross = Z*ComputeCrossSectionPerElectron
                                         (p,kineticEnergy,cutEnergy,maxEnergy);
  return cross;
}

G4double G4BraggModel::ComputeCrossSectionPerElectron	(	const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	cutEnergy,
		G4double	maxEnergy
	)			[virtual]

Definition at line 154 of file G4BraggModel.cc.

References MaxSecondaryEnergy().

Referenced by ComputeCrossSectionPerAtom(), and CrossSectionPerVolume().

{
  G4double cross     = 0.0;
  G4double tmax      = MaxSecondaryEnergy(p, kineticEnergy);
  G4double maxEnergy = std::min(tmax,maxKinEnergy);
  if(cutEnergy < maxEnergy) {

    G4double energy  = kineticEnergy + mass;
    G4double energy2 = energy*energy;
    G4double beta2   = kineticEnergy*(kineticEnergy + 2.0*mass)/energy2;
    cross = 1.0/cutEnergy - 1.0/maxEnergy - beta2*log(maxEnergy/cutEnergy)/tmax;

    cross *= twopi_mc2_rcl2*chargeSquare/beta2;
  }
 //   G4cout << "BR: e= " << kineticEnergy << " tmin= " << cutEnergy 
 //          << " tmax= " << tmax << " cross= " << cross << G4endl;
 
  return cross;
}

G4double G4BraggModel::ComputeDEDXPerVolume	(	const G4Material *	,
		const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	cutEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 209 of file G4BraggModel.cc.

References G4InuclParticleNames::gam, G4Material::GetElectronDensity(), and MaxSecondaryEnergy().

{
  G4double tmax  = MaxSecondaryEnergy(p, kineticEnergy);
  G4double tkin  = kineticEnergy/massRate;
  G4double dedx  = 0.0;

  if(tkin < lowestKinEnergy) {
    dedx = DEDX(material, lowestKinEnergy)*sqrt(tkin/lowestKinEnergy);
  } else {
    dedx = DEDX(material, tkin); 
  }

  if (cutEnergy < tmax) {

    G4double tau   = kineticEnergy/mass;
    G4double gam   = tau + 1.0;
    G4double bg2   = tau * (tau+2.0);
    G4double beta2 = bg2/(gam*gam);
    G4double x     = cutEnergy/tmax;

    dedx += (log(x) + (1.0 - x)*beta2) * twopi_mc2_rcl2
          * (material->GetElectronDensity())/beta2;
  }

  // now compute the total ionization loss

  if (dedx < 0.0) { dedx = 0.0; }

  dedx *= chargeSquare;

  //G4cout << "E(MeV)= " << tkin/MeV << " dedx= " << dedx 
  //     << "  " << material->GetName() << G4endl;

  return dedx;
}

G4double G4BraggModel::CrossSectionPerVolume	(	const G4Material *	,
		const G4ParticleDefinition *	,
		G4double	kineticEnergy,
		G4double	cutEnergy,
		G4double	maxEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 194 of file G4BraggModel.cc.

References ComputeCrossSectionPerElectron(), and G4Material::GetElectronDensity().

{
  G4double eDensity = material->GetElectronDensity();
  G4double cross = eDensity*ComputeCrossSectionPerElectron
                                         (p,kineticEnergy,cutEnergy,maxEnergy);
  return cross;
}

G4double G4BraggModel::GetChargeSquareRatio

(

)

const [inline, protected]

Definition at line 194 of file G4BraggModel.hh.

{
  return chargeSquare;
}

G4double G4BraggModel::GetChargeSquareRatio	(	const G4ParticleDefinition *	,
		const G4Material *	,
		G4double	kineticEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Definition at line 132 of file G4BraggModel.cc.

References G4EmCorrections::EffectiveChargeCorrection(), G4EmCorrections::EffectiveChargeSquareRatio(), G4VEmModel::GetModelOfFluctuations(), and G4VEmFluctuationModel::SetParticleAndCharge().

{
  // this method is called only for ions
  G4double q2 = corr->EffectiveChargeSquareRatio(p,mat,kineticEnergy);
  GetModelOfFluctuations()->SetParticleAndCharge(p, q2);
  return q2*corr->EffectiveChargeCorrection(p,mat,kineticEnergy);
}

G4double G4BraggModel::GetParticleCharge	(	const G4ParticleDefinition *	p,
		const G4Material *	mat,
		G4double	kineticEnergy
	)			[virtual]

Reimplemented from G4VEmModel.

Reimplemented in G4BraggIonGasModel.

Definition at line 144 of file G4BraggModel.cc.

References G4EmCorrections::GetParticleCharge().

{
  // this method is called only for ions, so no check if it is an ion 
  return corr->GetParticleCharge(p,mat,kineticEnergy);
}

void G4BraggModel::Initialise	(	const G4ParticleDefinition *	,
		const G4DataVector &
	)			[virtual]

Implements G4VEmModel.

Definition at line 109 of file G4BraggModel.cc.

References G4VEmModel::GetParticleChangeForLoss(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetParticleType(), and G4VEmModel::SetDeexcitationFlag().

{
  if(p != particle) { SetParticle(p); }

  // always false before the run
  SetDeexcitationFlag(false);

  if(!isInitialised) {
    isInitialised = true;

    G4String pname = particle->GetParticleName();
    if(particle->GetParticleType() == "nucleus" && 
       pname != "deuteron" && pname != "triton" &&
       pname != "alpha+"   && pname != "helium" &&
       pname != "hydrogen") { isIon = true; }

    fParticleChange = GetParticleChangeForLoss();
  }
}

G4double G4BraggModel::MaxSecondaryEnergy	(	const G4ParticleDefinition *	,
		G4double	kinEnergy
	)			[protected, virtual]

Reimplemented from G4VEmModel.

Definition at line 315 of file G4BraggModel.cc.

Referenced by ComputeCrossSectionPerElectron(), and ComputeDEDXPerVolume().

{
  if(pd != particle) { SetParticle(pd); }
  G4double tau  = kinEnergy/mass;
  G4double tmax = 2.0*electron_mass_c2*tau*(tau + 2.) /
                  (1. + 2.0*(tau + 1.)*ratio + ratio*ratio);
  return tmax;
}

virtual void G4BraggModel::SampleSecondaries	(	std::vector< G4DynamicParticle * > *	,
		const G4MaterialCutsCouple *	,
		const G4DynamicParticle *	,
		G4double	tmin,
		G4double	maxEnergy
	)			[virtual]

Implements G4VEmModel.

void G4BraggModel::SetChargeSquareRatio

(

G4double

val

)

[inline, protected]

Definition at line 199 of file G4BraggModel.hh.

Referenced by G4BraggIonGasModel::ChargeSquareRatio().

{
  chargeSquare = val;
}

---

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

• G4BraggModel.hh
• G4BraggModel.cc

---