Geant4-11
Public Member Functions | Protected Attributes | Private Attributes
G4PhotoElectricAngularGeneratorSauterGavrila Class Reference

#include <G4PhotoElectricAngularGeneratorSauterGavrila.hh>

Inheritance diagram for G4PhotoElectricAngularGeneratorSauterGavrila:
G4VEmAngularDistribution

Public Member Functions

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

Protected Attributes

G4ThreeVector fLocalDirection
 
G4bool fPolarisation
 

Private Attributes

G4String fName
 

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.

55 :
56 G4VEmAngularDistribution("AngularGenSauterGavrilaLowE")
57{}
G4VEmAngularDistribution::G4VEmAngularDistribution
G4VEmAngularDistribution(const G4String &name)
Definition: G4VEmAngularDistribution.cc:55

◆ ~G4PhotoElectricAngularGeneratorSauterGavrila()

G4PhotoElectricAngularGeneratorSauterGavrila::~G4PhotoElectricAngularGeneratorSauterGavrila ( )

Definition at line 61 of file G4PhotoElectricAngularGeneratorSauterGavrila.cc.

62{}

◆ G4PhotoElectricAngularGeneratorSauterGavrila() [2/2]

G4PhotoElectricAngularGeneratorSauterGavrila::G4PhotoElectricAngularGeneratorSauterGavrila ( const G4PhotoElectricAngularGeneratorSauterGavrila &  )
delete

Member Function Documentation

◆ GetName()

const G4String & G4VEmAngularDistribution::GetName ( ) const
inlineinherited

Definition at line 111 of file G4VEmAngularDistribution.hh.

112{
113 return fName;
114}

References G4VEmAngularDistribution::fName.

◆ 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.

111{
112 G4cout << "\n" << G4endl;
113 G4cout << "" << G4endl;
114 G4cout << "Re-implementation of the photolectric angular distribution" << G4endl;
115 G4cout << "developed my M. Maire for the Standard EM Physics G4PhotoElectricEffect" << G4endl;
116 G4cout << "It computes the theta distribution of the emitted electron, with respect to the" << G4endl;
117 G4cout << "incident Gamma, using the Sauter-Gavrila distribution for the K-shell\n" << G4endl;
118}
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout

References G4cout, and 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.

70{
71
72 // Compute Theta distribution of the emitted electron, with respect to the
73 // incident Gamma.
74 // The Sauter-Gavrila distribution for the K-shell is used.
75 G4double costeta = 1.;
76 G4double Phi = twopi * G4UniformRand();
77 G4double cosphi = std::cos(Phi);
78 G4double sinphi = std::sin(Phi);
79 G4double sinteta = 0;
80 G4double gamma = 1. + dp->GetKineticEnergy()/electron_mass_c2;
81
82 if (gamma > 5.) {
83 fLocalDirection = dp->GetMomentumDirection();
84 return fLocalDirection;
85 // Bugzilla 1120
86 // SI on 05/09/2010 as suggested by JG 04/09/10
87 }
88
89 G4double beta = std::sqrt((gamma - 1)*(gamma + 1))/gamma;
90 G4double b = 0.5*gamma*(gamma - 1)*(gamma - 2);
91
92 G4double rndm,term,greject,grejsup;
93 if (gamma < 2.) grejsup = gamma*gamma*(1.+b-beta*b);
94 else grejsup = gamma*gamma*(1.+b+beta*b);
95
96 do { rndm = 1.-2*G4UniformRand();
97 costeta = (rndm+beta)/(rndm*beta+1.);
98 term = 1.-beta*costeta;
99 greject = (1.-costeta*costeta)*(1.+b*term)/(term*term);
100 } while(greject < G4UniformRand()*grejsup);
101
102 sinteta = std::sqrt((1 - costeta)*(1 + costeta));
103 fLocalDirection.set(sinteta*cosphi, sinteta*sinphi, costeta);
104 fLocalDirection.rotateUz(dp->GetMomentumDirection());
105 return fLocalDirection;
106}
twopi
static constexpr double twopi
Definition: G4SIunits.hh:56
G4double
double G4double
Definition: G4Types.hh:83
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
CLHEP::Hep3Vector::set
void set(double x, double y, double z)
CLHEP::Hep3Vector::rotateUz
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:33
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta
const G4double beta
Definition: G4PionRadiativeDecayChannel.cc:44
source.hepunit.electron_mass_c2
float electron_mass_c2
Definition: hepunit.py:273

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::electron_mass_c2, G4VEmAngularDistribution::fLocalDirection, G4UniformRand, G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), and 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.

74{
75 return SampleDirection(dp, finalTotalEnergy, Z, mat);
76}
Z
const G4int Z[17]
Definition: G4WaterStopping.cc:51
G4VEmAngularDistribution::SampleDirection
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0

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.

85{
86 dirElectron = dp->GetMomentumDirection();
87 dirPositron = dp->GetMomentumDirection();
88}

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(), G4SauterGavrilaAngularDistribution::SampleDirection(), 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:
Generated by doxygen 1.9.3