G4PhotoElectricAngularGeneratorSauterGavrila Class Reference

#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 &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material *)
const G4String &	GetName () const

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection

Detailed Description

Definition at line 57 of file G4PhotoElectricAngularGeneratorSauterGavrila.hh.

Constructor & Destructor Documentation

G4PhotoElectricAngularGeneratorSauterGavrila::G4PhotoElectricAngularGeneratorSauterGavrila

(

)

Definition at line 54 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

                                                                                          :
  G4VEmAngularDistribution("AngularGenSauterGavrilaLowE")
{}

G4PhotoElectricAngularGeneratorSauterGavrila::~G4PhotoElectricAngularGeneratorSauterGavrila

(

)

Definition at line 58 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

{}

Member Function Documentation

void G4PhotoElectricAngularGeneratorSauterGavrila::PrintGeneratorInformation

(

)

const

Definition at line 104 of file G4PhotoElectricAngularGeneratorSauterGavrila.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;
} 

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.

References test::b, python.hepunit::electron_mass_c2, G4VEmAngularDistribution::fLocalDirection, G4UniformRand, G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), rndm(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), 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. 

  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