#include <G4PhotoElectricAngularGeneratorSauterGavrila.hh>

Inheritance diagram for G4PhotoElectricAngularGeneratorSauterGavrila:

Public Member Functions
	G4PhotoElectricAngularGeneratorSauterGavrila ()

	~G4PhotoElectricAngularGeneratorSauterGavrila ()

virtual G4ThreeVector &	SampleDirection (const G4DynamicParticle dp, G4double e=0.0, G4int shellId=0, const G4Material mat=0)

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

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)

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 57 of file G4PhotoElectricAngularGeneratorSauterGavrila.hh.

Constructor & Destructor Documentation

◆ G4PhotoElectricAngularGeneratorSauterGavrila()

G4PhotoElectricAngularGeneratorSauterGavrila::G4PhotoElectricAngularGeneratorSauterGavrila ( )

Definition at line 54 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

                                                                                          :
  G4VEmAngularDistribution("AngularGenSauterGavrilaLowE")
{}

◆ ~G4PhotoElectricAngularGeneratorSauterGavrila()

G4PhotoElectricAngularGeneratorSauterGavrila::~G4PhotoElectricAngularGeneratorSauterGavrila ( )

Definition at line 58 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

59{}

Member Function Documentation

◆ PrintGeneratorInformation()

void G4PhotoElectricAngularGeneratorSauterGavrila::PrintGeneratorInformation ( ) const

Definition at line 104 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

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

◆ SampleDirection()

G4ThreeVector & G4PhotoElectricAngularGeneratorSauterGavrila::SampleDirection	(	const G4DynamicParticle *	dp,
		G4double	e = `0.0`,
		G4int	shellId = `0`,
		const G4Material *	mat = `0`
	)

virtual

Implements G4VEmAngularDistribution.

Definition at line 62 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

{
 
  // Compute Theta distribution of the emitted electron, with respect to the
  // incident Gamma.
  // The Sauter-Gavrila distribution for the K-shell is used. 
 
  G4double costeta = 1.;
  G4double Phi     = twopi * G4UniformRand();
  G4double cosphi = std::cos(Phi);
  G4double sinphi = std::sin(Phi);
  G4double sinteta = 0;
  G4double gamma   = 1. + dp->GetKineticEnergy()/electron_mass_c2;
 
  if (gamma > 5.) {
    fLocalDirection = dp->GetMomentumDirection(); 
    return fLocalDirection; 
    // Bugzilla 1120
    // SI on 05/09/2010 as suggested by JG 04/09/10 
  }
 
  G4double beta  = std::sqrt((gamma - 1)*(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)*(1 + costeta));
  fLocalDirection.set(sinteta*cosphi, sinteta*sinphi, costeta);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
  return fLocalDirection;
}

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

Public Member Functions