G4RDPhotoElectricAngularGeneratorSauterGavrila Class Reference

#include <G4RDPhotoElectricAngularGeneratorSauterGavrila.hh>

Inheritance diagram for G4RDPhotoElectricAngularGeneratorSauterGavrila:

Public Member Functions
	G4RDPhotoElectricAngularGeneratorSauterGavrila (const G4String &name)
	~G4RDPhotoElectricAngularGeneratorSauterGavrila ()
G4ThreeVector	GetPhotoElectronDirection (const G4ThreeVector &direction, const G4double kineticEnergy, const G4ThreeVector &polarization, const G4int shellId) const
void	PrintGeneratorInformation () const
Public Member Functions inherited from G4RDVPhotoElectricAngularDistribution
	G4RDVPhotoElectricAngularDistribution (const G4String &name)
virtual	~G4RDVPhotoElectricAngularDistribution ()

Detailed Description

Definition at line 57 of file G4RDPhotoElectricAngularGeneratorSauterGavrila.hh.

Constructor & Destructor Documentation

G4RDPhotoElectricAngularGeneratorSauterGavrila::G4RDPhotoElectricAngularGeneratorSauterGavrila

(

const G4String &

name

)

Definition at line 57 of file G4RDPhotoElectricAngularGeneratorSauterGavrila.cc.

                                                                                                                  :G4RDVPhotoElectricAngularDistribution(name)
{;}

G4RDPhotoElectricAngularGeneratorSauterGavrila::~G4RDPhotoElectricAngularGeneratorSauterGavrila

(

)

Definition at line 62 of file G4RDPhotoElectricAngularGeneratorSauterGavrila.cc.

{;}

Member Function Documentation

G4ThreeVector G4RDPhotoElectricAngularGeneratorSauterGavrila::GetPhotoElectronDirection ( const G4ThreeVector &  direction, const G4double  kineticEnergy, const G4ThreeVector &  polarization, const G4int  shellId  ) const

virtual

Implements G4RDVPhotoElectricAngularDistribution.

Definition at line 67 of file G4RDPhotoElectricAngularGeneratorSauterGavrila.cc.

References test::b, python.hepunit::electron_mass_c2, G4UniformRand, rndm(), CLHEP::Hep3Vector::rotateUz(), and python.hepunit::twopi.

{

  // Compute Theta distribution of the emitted electron, with respect to the
  // incident Gamma.
  // The Sauter-Gavrila distribution for the K-shell is used. (adapted from G4PhotoElectricEffect)

  G4double costeta = 1.;
  G4double Phi     = twopi * G4UniformRand();
  G4double cosphi = std::cos(Phi);
  G4double sinphi = std::sin(Phi);
  G4double sinteta = 0;
  G4double gamma   = 1. + eKineticEnergy/electron_mass_c2;

  if (gamma > 5.) {
    G4ThreeVector direction (sinteta*cosphi, sinteta*sinphi, costeta);
    return direction;
  }

  G4double beta  = std::sqrt(gamma*gamma-1.)/gamma;
  G4double b     = 0.5*gamma*(gamma-1.)*(gamma-2);
    
  G4double rndm,term,greject,grejsup;
  if (gamma < 2.) grejsup = gamma*gamma*(1.+b-beta*b);
  else            grejsup = gamma*gamma*(1.+b+beta*b);
  
  do { rndm = 1.-2*G4UniformRand();
       costeta = (rndm+beta)/(rndm*beta+1.);
       term = 1.-beta*costeta;
       greject = (1.-costeta*costeta)*(1.+b*term)/(term*term);
  } while(greject < G4UniformRand()*grejsup);
       

  sinteta = std::sqrt(1.-costeta*costeta);
  G4ThreeVector photoelectrondirection (sinteta*cosphi, sinteta*sinphi, costeta);
  photoelectrondirection.rotateUz(direction);
  return photoelectrondirection;
}

void G4RDPhotoElectricAngularGeneratorSauterGavrila::PrintGeneratorInformation ( ) const

virtual

Implements G4RDVPhotoElectricAngularDistribution.

Definition at line 108 of file G4RDPhotoElectricAngularGeneratorSauterGavrila.cc.

References G4cout, and G4endl.

{
  G4cout << "\n" << G4endl;
  G4cout << "" << G4endl;
  G4cout << "Re-implementation of the photolectric angular distribution" << G4endl;
  G4cout << "developed my M. Maire for the Standard EM Physics G4PhotoElectricEffect" << G4endl;
  G4cout << "It computes the theta distribution of the emitted electron, with respect to the" << G4endl;
  G4cout << "incident Gamma, using the Sauter-Gavrila distribution for the K-shell\n" << G4endl;
} 

---

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

• G4RDPhotoElectricAngularGeneratorSauterGavrila.hh
• G4RDPhotoElectricAngularGeneratorSauterGavrila.cc