G4eBremsstrahlungSpectrum Class Reference

#include <G4eBremsstrahlungSpectrum.hh>

Inheritance diagram for G4eBremsstrahlungSpectrum:

Public Member Functions
	G4eBremsstrahlungSpectrum (const G4DataVector &bins, const G4String &name)
	~G4eBremsstrahlungSpectrum ()
G4double	Probability (G4int Z, G4double tMin, G4double tMax, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const
G4double	AverageEnergy (G4int Z, G4double tMin, G4double tMax, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const
G4double	SampleEnergy (G4int Z, G4double tMin, G4double tMax, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const
G4double	MaxEnergyOfSecondaries (G4double kineticEnergy, G4int Z=0, const G4ParticleDefinition *pd=0) const
G4double	Excitation (G4int Z, G4double kineticEnergy) const
void	PrintData () const
Public Member Functions inherited from G4VEnergySpectrum
	G4VEnergySpectrum ()
virtual	~G4VEnergySpectrum ()
virtual G4double	Probability (G4int Z, G4double minKineticEnergy, G4double maxKineticEnergy, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const =0
virtual G4double	AverageEnergy (G4int Z, G4double minKineticEnergy, G4double maxKineticEnergy, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const =0
virtual G4double	SampleEnergy (G4int Z, G4double minKineticEnergy, G4double maxKineticEnergy, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const =0
virtual G4double	MaxEnergyOfSecondaries (G4double kineticEnergy, G4int Z=0, const G4ParticleDefinition *pd=0) const =0
virtual G4double	Excitation (G4int Z, G4double kineticEnergy) const =0
virtual void	PrintData () const =0

Detailed Description

Definition at line 64 of file G4eBremsstrahlungSpectrum.hh.

Constructor & Destructor Documentation

G4eBremsstrahlungSpectrum()

G4eBremsstrahlungSpectrum::G4eBremsstrahlungSpectrum	(	const G4DataVector &	bins,
		const G4String &	name
	)

Definition at line 54 of file G4eBremsstrahlungSpectrum.cc.

                       :G4VEnergySpectrum(),
  lowestE(0.1*eV),
  xp(bins)
{
  length = xp.size();
  theBRparam = new G4BremsstrahlungParameters(name,length+1);
 
  verbose = 0;
}

~G4eBremsstrahlungSpectrum()

G4eBremsstrahlungSpectrum::~G4eBremsstrahlungSpectrum

(

)

Definition at line 66 of file G4eBremsstrahlungSpectrum.cc.

{
  delete theBRparam;
}

Member Function Documentation

AverageEnergy()

G4double G4eBremsstrahlungSpectrum::AverageEnergy ( G4int  Z, G4double  tMin, G4double  tMax, G4double  kineticEnergy, G4int  shell = 0, const G4ParticleDefinition *  pd = 0  ) const

virtual

Implements G4VEnergySpectrum.

Definition at line 119 of file G4eBremsstrahlungSpectrum.cc.

{
  G4double tm = std::min(tmax, e);
  G4double t0 = std::max(tmin, lowestE);
  if(t0 >= tm) return 0.0;
 
  t0 /= e;
  tm /= e;
 
  G4double z0 = lowestE/e;
 
  G4DataVector p;
 
  // Access parameters
  for (size_t i=0; i<=length; i++) {
      p.push_back(theBRparam->Parameter(i, Z, e));
  }
 
  G4double x  = AverageValue(t0, tm, p);
  G4double y  = IntSpectrum(z0, 1.0, p);
 
  // Add integrant over lowest energies
  G4double zmin = tmin/e;
  if(zmin < t0) {
    G4double c  = std::sqrt(theBRparam->ParameterC(Z));
    x += p[0]*(t0 - zmin - c*(std::atan(t0/c) - std::atan(zmin/c)));
  }
  x *= e;
 
  if(1 < verbose) {
    G4cout << "tcut(MeV)= " << tmin/MeV
           << "; tMax(MeV)= " << tmax/MeV
           << "; e(MeV)= " << e/MeV
           << "; t0= " << t0
           << "; tm= " << tm
           << "; y= " << y
           << "; x= " << x
           << G4endl;
  }
  p.clear();
 
  if(y > 0.0) x /= y;
  else        x  = 0.0;
  //  if(x < 0.0) x  = 0.0;
 
  return x;
}

Excitation()

G4double G4eBremsstrahlungSpectrum::Excitation ( G4int  Z, G4double  kineticEnergy  ) const

virtual

Implements G4VEnergySpectrum.

Definition at line 301 of file G4eBremsstrahlungSpectrum.cc.

{
  return 0.0;
}

MaxEnergyOfSecondaries()

G4double G4eBremsstrahlungSpectrum::MaxEnergyOfSecondaries ( G4double  kineticEnergy, G4int  Z = 0, const G4ParticleDefinition *  pd = 0  ) const

virtual

Implements G4VEnergySpectrum.

Definition at line 306 of file G4eBremsstrahlungSpectrum.cc.

{
  return kineticEnergy;
}

PrintData()

void G4eBremsstrahlungSpectrum::PrintData ( ) const

virtual

Implements G4VEnergySpectrum.

Definition at line 298 of file G4eBremsstrahlungSpectrum.cc.

{ theBRparam->PrintData(); }

Probability()

G4double G4eBremsstrahlungSpectrum::Probability ( G4int  Z, G4double  tMin, G4double  tMax, G4double  kineticEnergy, G4int  shell = 0, const G4ParticleDefinition *  pd = 0  ) const

virtual

Implements G4VEnergySpectrum.

Definition at line 72 of file G4eBremsstrahlungSpectrum.cc.

{
  G4double tm = std::min(tmax, e);
  G4double t0 = std::max(tmin, lowestE);
  if(t0 >= tm) return 0.0;
 
  t0 /= e;
  tm /= e;
 
  G4double z0 = lowestE/e;
  G4DataVector p;
 
  // Access parameters
  for (size_t i=0; i<=length; i++) {
    p.push_back(theBRparam->Parameter(i, Z, e));
  }
 
  G4double x  = IntSpectrum(t0, tm, p);
  G4double y  = IntSpectrum(z0, 1.0, p);
 
 
  if(1 < verbose) {
    G4cout << "tcut(MeV)= " << tmin/MeV
           << "; tMax(MeV)= " << tmax/MeV
           << "; t0= " << t0
           << "; tm= " << tm
           << "; xp[0]= " << xp[0]
           << "; z= " << z0
           << "; val= " << x
           << "; nor= " << y
           << G4endl;
  }
  p.clear();
 
  if(y > 0.0) x /= y;
  else        x  = 0.0;
  //  if(x < 0.0) x  = 0.0;
 
  return x;
}

SampleEnergy()

G4double G4eBremsstrahlungSpectrum::SampleEnergy ( G4int  Z, G4double  tMin, G4double  tMax, G4double  kineticEnergy, G4int  shell = 0, const G4ParticleDefinition *  pd = 0  ) const

virtual

Implements G4VEnergySpectrum.

Definition at line 173 of file G4eBremsstrahlungSpectrum.cc.

{
  G4double tm = std::min(tmax, e);
  G4double t0 = std::max(tmin, lowestE);
  if(t0 >= tm) return 0.0;
 
  t0 /= e;
  tm /= e;
 
  G4DataVector p;
 
  for (size_t i=0; i<=length; i++) {
    p.push_back(theBRparam->Parameter(i, Z, e));
  }
  G4double amaj = std::max(p[length], 1. - (p[1] - p[0])*xp[0]/(xp[1] - xp[0]) );
 
  G4double amax = std::log(tm);
  G4double amin = std::log(t0);
  G4double tgam, q, fun;
 
  do {
    G4double x = amin + G4UniformRand()*(amax - amin);
    tgam = G4Exp(x);
    fun = Function(tgam, p);
 
    if(fun > amaj) {
          G4cout << "WARNING in G4eBremsstrahlungSpectrum::SampleEnergy:"
                 << " Majoranta " << amaj
                 << " < " << fun
                 << G4endl;
    }
 
    q = amaj * G4UniformRand();
  } while (q > fun);
 
  tgam *= e;
 
  p.clear();
 
  return tgam;
}

---

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

• G4eBremsstrahlungSpectrum.hh
• G4eBremsstrahlungSpectrum.cc