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 &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material *)
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

(

const G4String &

name = ""

)

Definition at line 59 of file G4DipBustGenerator.cc.

  : G4VEmAngularDistribution("DipBustGen")
{}    

G4DipBustGenerator::~G4DipBustGenerator ( )

virtual

Definition at line 63 of file G4DipBustGenerator.cc.

{}

Member Function Documentation

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

Definition at line 105 of file G4DipBustGenerator.cc.

References test::a, test::c, python.hepunit::electron_mass_c2, G4Exp(), G4Log(), G4UniformRand, and python.hepunit::pi.

{
  G4double c, cosTheta, delta, cofA, signc = 1., a;
  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*G4Exp(G4Log(delta)/3.0);

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

void G4DipBustGenerator::PrintGeneratorInformation

(

)

const

Definition at line 142 of file G4DipBustGenerator.cc.

References G4cout, and G4endl.

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

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

virtual

Implements G4VEmAngularDistribution.

Definition at line 67 of file G4DipBustGenerator.cc.

References test::a, test::c, python.hepunit::electron_mass_c2, G4VEmAngularDistribution::fLocalDirection, G4Exp(), G4Log(), G4UniformRand, G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), and python.hepunit::twopi.

{
  G4double a, c, cosTheta, delta, cofA, signc = 1.;

  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*G4Exp(G4Log(delta)/3.0);

  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