G4PhotoElectricAngularGeneratorSauterGavrila Class Reference

#include <G4PhotoElectricAngularGeneratorSauterGavrila.hh>

Inheritance diagram for G4PhotoElectricAngularGeneratorSauterGavrila:

Public Member Functions
	G4PhotoElectricAngularGeneratorSauterGavrila ()
	~G4PhotoElectricAngularGeneratorSauterGavrila ()
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double e=0.0, G4int shellId=0, const G4Material *mat=nullptr) override
void	PrintGeneratorInformation () const override
G4PhotoElectricAngularGeneratorSauterGavrila &	operator= (const G4PhotoElectricAngularGeneratorSauterGavrila &right)=delete
	G4PhotoElectricAngularGeneratorSauterGavrila (const G4PhotoElectricAngularGeneratorSauterGavrila &)=delete
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 *)
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() [1/2]

G4PhotoElectricAngularGeneratorSauterGavrila::G4PhotoElectricAngularGeneratorSauterGavrila ( )

explicit

Definition at line 55 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

                                                                                          :
  G4VEmAngularDistribution("AngularGenSauterGavrilaLowE")
{}

Referenced by G4PhotoElectricAngularGeneratorSauterGavrila(), and operator=().

~G4PhotoElectricAngularGeneratorSauterGavrila()

G4PhotoElectricAngularGeneratorSauterGavrila::~G4PhotoElectricAngularGeneratorSauterGavrila

(

)

Definition at line 61 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

{}

G4PhotoElectricAngularGeneratorSauterGavrila() [2/2]

G4PhotoElectricAngularGeneratorSauterGavrila::G4PhotoElectricAngularGeneratorSauterGavrila ( const G4PhotoElectricAngularGeneratorSauterGavrila & )

delete

Member Function Documentation

operator=()

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

delete

PrintGeneratorInformation()

void G4PhotoElectricAngularGeneratorSauterGavrila::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 110 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 = nullptr )

overridevirtual

Implements G4VEmAngularDistribution.

Definition at line 67 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:

• G4PhotoElectricAngularGeneratorSauterGavrila.hh
• G4PhotoElectricAngularGeneratorSauterGavrila.cc