G4RDeBremsstrahlungSpectrum Class Reference

#include <G4RDeBremsstrahlungSpectrum.hh>

Inheritance 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 ()

Detailed Description

Definition at line 66 of file G4RDeBremsstrahlungSpectrum.hh.

Constructor & Destructor Documentation

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

(

)

Definition at line 68 of file G4RDeBremsstrahlungSpectrum.cc.

{
  delete theBRparam;
}

Member Function Documentation

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.

References test::c, G4cout, G4endl, G4INCL::Math::max(), python.hepunit::MeV, G4INCL::Math::min(), G4RDBremsstrahlungParameters::Parameter(), G4RDBremsstrahlungParameters::ParameterC(), G4InuclParticleNames::tm, test::x, and G4InuclParticleNames::z0.

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

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

virtual

Implements G4RDVEnergySpectrum.

Definition at line 303 of file G4RDeBremsstrahlungSpectrum.cc.

{
  return 0.0;
}

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

void G4RDeBremsstrahlungSpectrum::PrintData ( void  ) const

virtual

Implements G4RDVEnergySpectrum.

Definition at line 300 of file G4RDeBremsstrahlungSpectrum.cc.

References G4RDBremsstrahlungParameters::PrintData().

{ theBRparam->PrintData(); }

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.

References G4cout, G4endl, G4INCL::Math::max(), python.hepunit::MeV, G4INCL::Math::min(), G4RDBremsstrahlungParameters::Parameter(), G4InuclParticleNames::tm, test::x, and G4InuclParticleNames::z0.

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

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.

References G4cout, G4endl, G4UniformRand, G4INCL::Math::max(), G4INCL::Math::min(), G4RDBremsstrahlungParameters::Parameter(), G4InuclParticleNames::tm, and test::x.

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

---

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

• G4RDeBremsstrahlungSpectrum.hh
• G4RDeBremsstrahlungSpectrum.cc