G4SauterGavrilaAngularDistribution Class Reference

#include <G4SauterGavrilaAngularDistribution.hh>

Inheritance diagram for G4SauterGavrilaAngularDistribution:

Public Member Functions
	G4SauterGavrilaAngularDistribution ()
	G4SauterGavrilaAngularDistribution (const G4SauterGavrilaAngularDistribution &)=delete
const G4String &	GetName () const
G4SauterGavrilaAngularDistribution &	operator= (const G4SauterGavrilaAngularDistribution &right)=delete
void	PrintGeneratorInformation () const override
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double e=0.0, G4int shellId=0, const G4Material *mat=0) final
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)
	~G4SauterGavrilaAngularDistribution () override

Protected Attributes
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Private Attributes
G4String	fName

Detailed Description

Definition at line 52 of file G4SauterGavrilaAngularDistribution.hh.

Constructor & Destructor Documentation

G4SauterGavrilaAngularDistribution() [1/2]

G4SauterGavrilaAngularDistribution::G4SauterGavrilaAngularDistribution ( )

explicit

Definition at line 48 of file G4SauterGavrilaAngularDistribution.cc.

  : G4VEmAngularDistribution("SauterGavrila")
{}    

~G4SauterGavrilaAngularDistribution()

G4SauterGavrilaAngularDistribution::~G4SauterGavrilaAngularDistribution ( )

override

Definition at line 52 of file G4SauterGavrilaAngularDistribution.cc.

{}

G4SauterGavrilaAngularDistribution() [2/2]

G4SauterGavrilaAngularDistribution::G4SauterGavrilaAngularDistribution ( const G4SauterGavrilaAngularDistribution &  )

delete

Member Function Documentation

GetName()

const G4String & G4VEmAngularDistribution::GetName ( ) const

inlineinherited

Definition at line 111 of file G4VEmAngularDistribution.hh.

{
  return fName;
}

References G4VEmAngularDistribution::fName.

operator=()

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

delete

PrintGeneratorInformation()

void G4SauterGavrilaAngularDistribution::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 110 of file G4SauterGavrilaAngularDistribution.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Non-polarized photoelectric effect angular generator." << G4endl;
  G4cout << "The Sauter-Gavrila distribution for the K-shell is used."<<G4endl;
  G4cout << "Originally developed by M.Maire for Geant3" 
         << G4endl;
} 

References G4cout, and G4endl.

SampleDirection()

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

finalvirtual

Implements G4VEmAngularDistribution.

Definition at line 55 of file G4SauterGavrilaAngularDistribution.cc.

{
  static const G4double emin = 1*CLHEP::eV;
  static const G4double emax = 100*CLHEP::MeV;
 
  G4double energy = std::max(dp->GetKineticEnergy(), emin);
  if (energy > emax) {
    fLocalDirection = dp->GetMomentumDirection();
  } else {
    // Initial algorithm according Penelope 2008 manual and 
    // F.Sauter Ann. Physik 9, 217(1931); 11, 454(1931).
    // Modified according Penelope 2014 manual 
    G4double costheta = 0.0;
 
    // 1) initialize energy-dependent variables
    // Variable naming according to Eq. (2.24) of Penelope Manual
    // (pag. 44)
    G4double tau = energy/electron_mass_c2;
    G4double gamma = 1.0 + tau;
    G4double beta = std::sqrt(tau*(tau + 2.0))/gamma;
    
    // ac corresponds to "A" of Eq. (2.31)
    //
    G4double ac = (1.0 - beta)/beta;
    G4double a1 = 0.5*beta*gamma*tau*(gamma-2.0);
    G4double a2 = ac + 2.0;
    // gtmax = maximum of the rejection function according to Eq. (2.28), 
    // obtained for tsam=0
    G4double gtmax = 2.0*(a1 + 1.0/ac);
    
    G4double tsam = 0.0;
    G4double gtr  = 0.0;
    
    //2) sampling. Eq. (2.31) of Penelope Manual
    // tsam = 1-std::cos(theta)
    // gtr = rejection function according to Eq. (2.28)
    do{
      G4double rand = G4UniformRand();
      tsam = 2.0*ac * (2.0*rand + a2*std::sqrt(rand)) / (a2*a2 - 4.0*rand);
      gtr = (2.0 - tsam) * (a1 + 1.0/(ac+tsam));
      // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
    } while(G4UniformRand()*gtmax > gtr);
 
    costheta = 1.0 - tsam;
        
    G4double sint = std::sqrt(tsam*(2.0 - tsam));
    G4double phi  = CLHEP::twopi*G4UniformRand();
        
    fLocalDirection.set(sint*std::cos(phi), sint*std::sin(phi), costheta);
    fLocalDirection.rotateUz(dp->GetMomentumDirection());
  }
  return fLocalDirection;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::electron_mass_c2, emax, G4INCL::KinematicsUtils::energy(), CLHEP::eV, G4VEmAngularDistribution::fLocalDirection, G4UniformRand, G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), G4INCL::Math::max(), CLHEP::MeV, CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), and CLHEP::twopi.

SampleDirectionForShell()

G4ThreeVector & G4VEmAngularDistribution::SampleDirectionForShell ( const G4DynamicParticle *  dp, G4double  finalTotalEnergy, G4int  Z, G4int  shellID, const G4Material *  mat  )

virtualinherited

Reimplemented in G4DeltaAngle, G4DNABornAngle, and G4DNARuddAngle.

Definition at line 69 of file G4VEmAngularDistribution.cc.

{
  return SampleDirection(dp, finalTotalEnergy, Z, mat); 
}

References G4VEmAngularDistribution::SampleDirection(), and Z.

Referenced by G4DNABornIonisationModel1::SampleSecondaries(), G4DNABornIonisationModel2::SampleSecondaries(), G4DNAEmfietzoglouIonisationModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), G4DNARuddIonisationModel::SampleSecondaries(), G4MicroElecInelasticModel::SampleSecondaries(), and G4MicroElecInelasticModel_new::SampleSecondaries().

SamplePairDirections()

void G4VEmAngularDistribution::SamplePairDirections ( const G4DynamicParticle *  dp, G4double  elecKinEnergy, G4double  posiKinEnergy, G4ThreeVector &  dirElectron, G4ThreeVector &  dirPositron, G4int  Z = 0, const G4Material *  mat = nullptr  )

virtualinherited

Reimplemented in G4ModifiedMephi, G4DipBustGenerator, and G4ModifiedTsai.

Definition at line 80 of file G4VEmAngularDistribution.cc.

{
  dirElectron = dp->GetMomentumDirection();
  dirPositron = dp->GetMomentumDirection();
}

References G4DynamicParticle::GetMomentumDirection().

Referenced by G4MuPairProductionModel::SampleSecondaries(), G4BetheHeitlerModel::SampleSecondaries(), and G4PairProductionRelModel::SampleSecondaries().

Field Documentation

fLocalDirection

G4ThreeVector G4VEmAngularDistribution::fLocalDirection

protectedinherited

Definition at line 103 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::G4VEmAngularDistribution(), SampleDirection(), G4PhotoElectricAngularGeneratorSauterGavrila::SampleDirection(), G4PhotoElectricAngularGeneratorPolarized::SampleDirection(), G4ModifiedMephi::SampleDirection(), G4DeltaAngle::SampleDirection(), G4DeltaAngleFreeScat::SampleDirection(), G4DipBustGenerator::SampleDirection(), G4ModifiedTsai::SampleDirection(), G4Generator2BN::SampleDirection(), G4Generator2BS::SampleDirection(), G4PenelopeBremsstrahlungAngular::SampleDirection(), G4RayleighAngularGenerator::SampleDirection(), G4AngleDirect::SampleDirection(), G4DNABornAngle::SampleDirectionForShell(), and G4DNARuddAngle::SampleDirectionForShell().

fName

G4String G4VEmAngularDistribution::fName

privateinherited

Definition at line 108 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::GetName().

fPolarisation

G4bool G4VEmAngularDistribution::fPolarisation

protectedinherited

Definition at line 104 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::G4VEmAngularDistribution().

---

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

• source/processes/electromagnetic/standard/include/G4SauterGavrilaAngularDistribution.hh
• source/processes/electromagnetic/standard/src/G4SauterGavrilaAngularDistribution.cc