G4RDeBremsstrahlungSpectrum Class Reference

#include <G4RDeBremsstrahlungSpectrum.hh>

Inheritance diagram for G4RDeBremsstrahlungSpectrum:

Collaboration diagram for G4RDeBremsstrahlungSpectrum:

Public Member Functions
	G4RDeBremsstrahlungSpectrum (const G4DataVector &bins, const G4String &name)
	~G4RDeBremsstrahlungSpectrum ()
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 G4RDVEnergySpectrum
	G4RDVEnergySpectrum ()
virtual	~G4RDVEnergySpectrum ()

Private Member Functions
G4double	IntSpectrum (G4double xMin, G4double xMax, const G4DataVector &p) const
G4double	AverageValue (G4double xMin, G4double xMax, const G4DataVector &p) const
G4double	Function (G4double x, const G4DataVector &p) const
	G4RDeBremsstrahlungSpectrum (const G4RDeBremsstrahlungSpectrum &)
G4RDeBremsstrahlungSpectrum &	operator= (const G4RDeBremsstrahlungSpectrum &right)

Private Attributes
G4RDBremsstrahlungParameters *	theBRparam
G4double	lowestE
size_t	length
G4int	verbose
const G4DataVector	xp

Detailed Description

Definition at line 66 of file G4RDeBremsstrahlungSpectrum.hh.

Constructor & Destructor Documentation

G4RDeBremsstrahlungSpectrum() [1/2]

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

Definition at line 56 of file G4RDeBremsstrahlungSpectrum.cc.

                       :G4RDVEnergySpectrum(),
  lowestE(0.1*eV),
  xp(bins)
{
  length = xp.size();
  theBRparam = new G4RDBremsstrahlungParameters(name,length+1);

  verbose = 0;
}

~G4RDeBremsstrahlungSpectrum()

G4RDeBremsstrahlungSpectrum::~G4RDeBremsstrahlungSpectrum

(

)

Definition at line 68 of file G4RDeBremsstrahlungSpectrum.cc.

{
  delete theBRparam;
}

G4RDeBremsstrahlungSpectrum() [2/2]

G4RDeBremsstrahlungSpectrum::G4RDeBremsstrahlungSpectrum ( const G4RDeBremsstrahlungSpectrum &  )

private

Member Function Documentation

AverageEnergy()

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

virtual

Implements G4RDVEnergySpectrum.

Definition at line 121 of file G4RDeBremsstrahlungSpectrum.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;
}

Here is the call graph for this function:

AverageValue()

G4double G4RDeBremsstrahlungSpectrum::AverageValue ( G4double  xMin, G4double  xMax, const G4DataVector &  p  ) const

private

Definition at line 248 of file G4RDeBremsstrahlungSpectrum.cc.

{
  G4double x1 = std::min(xMin, xp[0]);
  G4double x2 = std::min(xMax, xp[0]);
  G4double z1 = x1;
  G4double z2 = x2;
  G4double sum = 0.0;

  if(x1 < x2) {
    G4double k = (p[1] - p[0])/(xp[1] - xp[0]);
    sum += (z2 - z1)*(1. - k*xp[0]);
    z1 *= x1;
    z2 *= x2;
    sum += 0.5*k*(z2 - z1);
  }

  for (size_t i=0; i<length-1; i++) {
    x1 = std::max(xMin, xp[i]);
    x2 = std::min(xMax, xp[i+1]);
    if(x1 < x2) {
      z1 = p[i];
      z2 = p[i+1];
      sum += 0.5*(z2 - z1)*(x2 + x1) + z1*x2 - z2*x1;
    }
  }
  if(sum < 0.0) sum = 0.0;
  return sum;
}

Here is the caller graph for this function:

Excitation()

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

virtual

Implements G4RDVEnergySpectrum.

Definition at line 303 of file G4RDeBremsstrahlungSpectrum.cc.

{
  return 0.0;
}

Function()

G4double G4RDeBremsstrahlungSpectrum::Function ( G4double  x, const G4DataVector &  p  ) const

private

Definition at line 279 of file G4RDeBremsstrahlungSpectrum.cc.

{
  G4double f = 0.0;

  if(x <= xp[0]) {
    f = p[0] + (p[1] - p[0])*(x - xp[0])/(xp[1] - xp[0]);

  } else {

    for (size_t i=0; i<length-1; i++) {

      if(x <= xp[i+1]) {
        f = p[i] + (p[i+1] - p[i])*(x - xp[i])/(xp[i+1] - xp[i]);
        break;
      }
    }
  }
  return f;
}

Here is the caller graph for this function:

IntSpectrum()

G4double G4RDeBremsstrahlungSpectrum::IntSpectrum ( G4double  xMin, G4double  xMax, const G4DataVector &  p  ) const

private

Definition at line 222 of file G4RDeBremsstrahlungSpectrum.cc.

{
  G4double x1 = std::min(xMin, xp[0]);
  G4double x2 = std::min(xMax, xp[0]);
  G4double sum = 0.0;

  if(x1 < x2) {
    G4double k = (p[1] - p[0])/(xp[1] - xp[0]);
    sum += (1. - k*xp[0])*std::log(x2/x1) + k*(x2 - x1);
  }

  for (size_t i=0; i<length-1; i++) {
    x1 = std::max(xMin, xp[i]);
    x2 = std::min(xMax, xp[i+1]);
    if(x1 < x2) {
      G4double z1 = p[i];
      G4double z2 = p[i+1];
      sum += z2 - z1 + std::log(x2/x1)*(z1*x2 - z2*x1)/(x2 - x1);
    }
  }
  if(sum < 0.0) sum = 0.0;
  return sum;
}

Here is the caller graph for this function:

MaxEnergyOfSecondaries()

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

virtual

Implements G4RDVEnergySpectrum.

Definition at line 308 of file G4RDeBremsstrahlungSpectrum.cc.

{
  return kineticEnergy;
}

operator=()

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

private

PrintData()

void G4RDeBremsstrahlungSpectrum::PrintData ( void  ) const

virtual

Implements G4RDVEnergySpectrum.

Definition at line 300 of file G4RDeBremsstrahlungSpectrum.cc.

{ theBRparam->PrintData(); }

Here is the call graph for this function:

Probability()

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

virtual

Implements G4RDVEnergySpectrum.

Definition at line 74 of file G4RDeBremsstrahlungSpectrum.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;
}

Here is the call graph for this function:

SampleEnergy()

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

virtual

Implements G4RDVEnergySpectrum.

Definition at line 175 of file G4RDeBremsstrahlungSpectrum.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 = std::exp(x);
    fun = Function(tgam, p);

    if(fun > amaj) {
          G4cout << "WARNING in G4RDeBremsstrahlungSpectrum::SampleEnergy:"
                 << " Majoranta " << amaj
                 << " < " << fun
                 << G4endl;
    }

    q = amaj * G4UniformRand();
  } while (q > fun);

  tgam *= e;

  p.clear();

  return tgam;
}

Here is the call graph for this function:

Member Data Documentation

length

size_t G4RDeBremsstrahlungSpectrum::length

private

Definition at line 120 of file G4RDeBremsstrahlungSpectrum.hh.

lowestE

G4double G4RDeBremsstrahlungSpectrum::lowestE

private

Definition at line 119 of file G4RDeBremsstrahlungSpectrum.hh.

theBRparam

G4RDBremsstrahlungParameters* G4RDeBremsstrahlungSpectrum::theBRparam

private

Definition at line 118 of file G4RDeBremsstrahlungSpectrum.hh.

verbose

G4int G4RDeBremsstrahlungSpectrum::verbose

private

Definition at line 121 of file G4RDeBremsstrahlungSpectrum.hh.

xp

const G4DataVector G4RDeBremsstrahlungSpectrum::xp

private

Definition at line 123 of file G4RDeBremsstrahlungSpectrum.hh.

---

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

• G4RDeBremsstrahlungSpectrum.hh
• G4RDeBremsstrahlungSpectrum.cc