G4RayleighAngularGenerator Class Reference

#include <G4RayleighAngularGenerator.hh>

Inheritance diagram for G4RayleighAngularGenerator:

Public Member Functions
	G4RayleighAngularGenerator ()
virtual	~G4RayleighAngularGenerator ()
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double out_energy, G4int Z, const G4Material *mat=nullptr) override
G4RayleighAngularGenerator &	operator= (const G4RayleighAngularGenerator &right)=delete
	G4RayleighAngularGenerator (const G4RayleighAngularGenerator &)=delete
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
virtual G4ThreeVector &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material *)
virtual void	SamplePairDirections (const G4DynamicParticle *dp, G4double elecKinEnergy, G4double posiKinEnergy, G4ThreeVector &dirElectron, G4ThreeVector &dirPositron, G4int Z=0, const G4Material *mat=nullptr)
virtual void	PrintGeneratorInformation () const
const G4String &	GetName () const
G4VEmAngularDistribution &	operator= (const G4VEmAngularDistribution &right)=delete
	G4VEmAngularDistribution (const G4VEmAngularDistribution &)=delete

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Detailed Description

Definition at line 59 of file G4RayleighAngularGenerator.hh.

Constructor & Destructor Documentation

G4RayleighAngularGenerator() [1/2]

G4RayleighAngularGenerator::G4RayleighAngularGenerator ( )

explicit

Definition at line 63 of file G4RayleighAngularGenerator.cc.

  : G4VEmAngularDistribution("CullenGenerator")
{
  G4double x = cm/(h_Planck*c_light);
  fFactor = 0.5*x*x;
}

~G4RayleighAngularGenerator()

G4RayleighAngularGenerator::~G4RayleighAngularGenerator ( )

virtual

Definition at line 72 of file G4RayleighAngularGenerator.cc.

{}

G4RayleighAngularGenerator() [2/2]

G4RayleighAngularGenerator::G4RayleighAngularGenerator ( const G4RayleighAngularGenerator &  )

delete

Member Function Documentation

operator=()

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

delete

SampleDirection()

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

overridevirtual

Implements G4VEmAngularDistribution.

Definition at line 78 of file G4RayleighAngularGenerator.cc.

{
  G4double ekin = dp->GetKineticEnergy();
  G4double xx = fFactor*ekin*ekin;
 
  G4double n0 = PP6[Z] - 1.0;
  G4double n1 = PP7[Z] - 1.0;
  G4double n2 = PP8[Z] - 1.0;
  G4double b0 = PP3[Z];
  G4double b1 = PP4[Z];
  G4double b2 = PP5[Z];
 
  static const G4double numlim = 0.02;
  G4double x  = 2.*xx*b0;
  G4double w0 = (x < numlim) ? n0*x*(1.0 - 0.5*(n0 - 1.0)*x*(1.0 - (n0 - 2.0)*x/3.))
    : 1.0 - G4Exp(-n0*G4Log(1.0 + x)); 
 
  x  = 2.*xx*b1;
  G4double w1 = (x < numlim) ? n1*x*(1.0 - 0.5*(n1 - 1.0)*x*(1.0 - (n1 - 2.0)*x/3.))
    : 1.0 - G4Exp(-n1*G4Log(1.0 + x));
 
  x  = 2.*xx*b2;
  G4double w2 = (x < numlim) ? n2*x*(1.0 - 0.5*(n2 - 1.0)*x*(1.0 - (n2 - 2.0)*x/3.))
    : 1.0 - G4Exp(-n2*G4Log(1.0 + x));
 
  G4double x0= w0*PP0[Z]/(b0*n0);
  G4double x1= w1*PP1[Z]/(b1*n1);
  G4double x2= w2*PP2[Z]/(b2*n2);
 
  G4double cost;
  do {
    G4double w = w0;
    G4double n = n0;
    G4double b = b0;
  
    x = G4UniformRand()*(x0 + x1 + x2);
    if(x > x0) {
      x -= x0;
      if(x <= x1 ) {
    w = w1;
    n = n1;
    b = b1;
      } else { 
    w = w2;
    n = n2;
    b = b2;
      }
    }
    n = 1.0/n;
 
    // sampling of angle
    G4double y = w*G4UniformRand();
    if(y < numlim) { x = y*n*( 1. + 0.5*(n + 1.)*y*(1. - (n + 2.)*y/3.)); }   
    else           { x = G4Exp(-n*G4Log(1. - y)) - 1.0; }
    cost = 1.0 - x/(b*xx);
  } while (2*G4UniformRand() > 1.0 + cost*cost || cost < -1.0);
 
  G4double phi  = twopi*G4UniformRand();
  G4double sint = sqrt((1. - cost)*(1.0 + cost));
  fLocalDirection.set(sint*cos(phi),sint*sin(phi),cost);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
 
  return fLocalDirection;
}

---

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

• G4RayleighAngularGenerator.hh
• G4RayleighAngularGenerator.cc