G4DipBustGenerator Class Reference

#include <G4DipBustGenerator.hh>

Inheritance diagram for G4DipBustGenerator:

Public Member Functions
	G4DipBustGenerator (const G4String &name="")
virtual	~G4DipBustGenerator ()
virtual G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double out_energy, G4int Z, const G4Material *mat=0)
G4double	PolarAngle (const G4double initial_energy, const G4double final_energy, const G4int Z)
void	PrintGeneratorInformation () const
Public Member Functions inherited from G4VEmAngularDistribution
	G4VEmAngularDistribution (const G4String &name)
virtual	~G4VEmAngularDistribution ()
virtual G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
const G4String &	GetName () const

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection

Detailed Description

Definition at line 55 of file G4DipBustGenerator.hh.

Constructor & Destructor Documentation

G4DipBustGenerator()

G4DipBustGenerator::G4DipBustGenerator

(

const G4String &

name = ""

)

Definition at line 57 of file G4DipBustGenerator.cc.

  : G4VEmAngularDistribution("DipBustGen")
{}    

~G4DipBustGenerator()

G4DipBustGenerator::~G4DipBustGenerator ( )

virtual

Definition at line 61 of file G4DipBustGenerator.cc.

{}

Member Function Documentation

PolarAngle()

G4double G4DipBustGenerator::PolarAngle	(	const G4double	initial_energy,
		const G4double	final_energy,
		const G4int	Z
	)

Definition at line 103 of file G4DipBustGenerator.cc.

{
  G4double c, cosTheta, delta, cofA, signc = 1., a, power = 1./3.;
  G4double gamma, beta, theta;
 
  c = 4. - 8.*G4UniformRand();
  a = c;
 
  if( c < 0. )
  {
    signc = -1.;
    a     = -c;
  }
  delta  = std::sqrt(a*a+4.);
  delta += a;
  delta *= 0.5; 
 
  cofA = -signc*std::pow(delta, power);
 
  cosTheta = cofA - 1./cofA;
 
  gamma = 1. + eTkin/electron_mass_c2;
  beta = std::sqrt(1. - 1./gamma/gamma);
 
  cosTheta = (cosTheta + beta)/(1 + cosTheta*beta);
 
  theta = std::acos(cosTheta);
 
  if( theta < 0. )  theta = 0.;
  if( theta > pi )  theta = pi;
  // G4cout <<"theta = "<<theta<<"; ";
 
  return theta;
}

PrintGeneratorInformation()

void G4DipBustGenerator::PrintGeneratorInformation

(

)

const

Definition at line 140 of file G4DipBustGenerator.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Angular Generator based on classical formula from" << G4endl;
  G4cout << "J.D. Jackson, Classical Electrodynamics, Wiley, New York 1975" 
     << G4endl;
} 

SampleDirection()

G4ThreeVector & G4DipBustGenerator::SampleDirection ( const G4DynamicParticle *  dp, G4double  out_energy, G4int  Z, const G4Material *  mat = 0  )

virtual

Implements G4VEmAngularDistribution.

Definition at line 65 of file G4DipBustGenerator.cc.

{
  G4double c, cosTheta, delta, cofA, signc = 1., a, power = 1./3.;
 
  G4double eTkin = dp->GetKineticEnergy();
 
  c = 4. - 8.*G4UniformRand();
  a = c;
 
  if( c < 0. )
  {
    signc = -1.;
    a     = -c;
  }
  delta  = std::sqrt(a*a+4.);
  delta += a;
  delta *= 0.5; 
 
  cofA = -signc*std::pow(delta, power);
 
  cosTheta = cofA - 1./cofA;
 
  G4double tau = eTkin/electron_mass_c2;
  G4double beta = std::sqrt(tau*(tau + 2.))/(tau + 1.);
 
  cosTheta = (cosTheta + beta)/(1 + cosTheta*beta);
 
  G4double sinTheta = std::sqrt((1 - cosTheta)*(1 + cosTheta));
  G4double phi  = twopi*G4UniformRand(); 
 
  fLocalDirection.set(sinTheta*std::cos(phi), sinTheta*std::sin(phi),cosTheta);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
 
  return fLocalDirection;
 
}

---

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

• G4DipBustGenerator.hh
• G4DipBustGenerator.cc