#include <G4RDModifiedTsai.hh>

Inheritance diagram for G4RDModifiedTsai:

Collaboration diagram for G4RDModifiedTsai:

Public Member Functions
	G4RDModifiedTsai (const G4String &name)

	~G4RDModifiedTsai ()

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

void	PrintGeneratorInformation () const

Public Member Functions inherited from G4RDVBremAngularDistribution
	G4RDVBremAngularDistribution (const G4String &name)

virtual	~G4RDVBremAngularDistribution ()

Private Member Functions
G4RDModifiedTsai &	operator= (const G4RDModifiedTsai &right)

	G4RDModifiedTsai (const G4RDModifiedTsai &)

Detailed Description

Definition at line 60 of file G4RDModifiedTsai.hh.

Constructor & Destructor Documentation

◆ G4RDModifiedTsai() [1/2]

G4RDModifiedTsai::G4RDModifiedTsai ( const G4String & name )

Definition at line 59 of file G4RDModifiedTsai.cc.

59 :G4RDVBremAngularDistribution(name)

60 {;}

G4RDVBremAngularDistribution::G4RDVBremAngularDistribution

G4RDVBremAngularDistribution(const G4String &name)

Definition: G4RDVBremAngularDistribution.cc:57

◆ ~G4RDModifiedTsai()

G4RDModifiedTsai::~G4RDModifiedTsai ( )

Definition at line 64 of file G4RDModifiedTsai.cc.

65 {;}

◆ G4RDModifiedTsai() [2/2]

G4RDModifiedTsai::G4RDModifiedTsai ( const G4RDModifiedTsai & )

private

Member Function Documentation

◆ operator=()

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

private

◆ PolarAngle()

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

virtual

Implements G4RDVBremAngularDistribution.

Definition at line 69 of file G4RDModifiedTsai.cc.

 {
 
   // Sample gamma angle (Z - axis along the parent particle).
   // Universal distribution suggested by L. Urban (Geant3 manual (1993) 
   // Phys211) derived from Tsai distribution (Rev Mod Phys 49,421(1977))
 
   G4double totalEnergy = initial_energy + electron_mass_c2;   
 
   const G4double a1 = 0.625, a2 = 3.*a1, d = 27.;
   G4double u, theta = 0;
 
   do{
   u = - std::log(G4UniformRand()*G4UniformRand());
 
   if (9./(9.+d) > G4UniformRand()) u /= a1;
   else                             u /= a2;
     
   theta = u*electron_mass_c2/totalEnergy;
   }while(u > (totalEnergy*pi/electron_mass_c2));
 
   return theta;
 }

◆ PrintGeneratorInformation()

void G4RDModifiedTsai::PrintGeneratorInformation ( ) const