G4DipBustGenerator Class Reference

#include <G4DipBustGenerator.hh>

Inheritance diagram for G4DipBustGenerator:

Public Member Functions
	G4DipBustGenerator (const G4DipBustGenerator &)=delete
	G4DipBustGenerator (const G4String &name="")
const G4String &	GetName () const
G4DipBustGenerator &	operator= (const G4DipBustGenerator &right)=delete
G4double	PolarAngle (G4double initial_energy, G4double final_energy, G4int Z)
void	PrintGeneratorInformation () const final
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double out_energy, G4int Z, const G4Material *mat=nullptr) final
virtual G4ThreeVector &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material *)
void	SamplePairDirections (const G4DynamicParticle *dp, G4double elecKinEnergy, G4double posiKinEnergy, G4ThreeVector &dirElectron, G4ThreeVector &dirPositron, G4int Z=0, const G4Material *mat=nullptr) final
	~G4DipBustGenerator () override

Protected Attributes
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Private Member Functions
G4double	SampleCosTheta (G4double kinEnergy)

Private Attributes
G4String	fName

Detailed Description

Definition at line 54 of file G4DipBustGenerator.hh.

Constructor & Destructor Documentation

G4DipBustGenerator() [1/2]

G4DipBustGenerator::G4DipBustGenerator ( const G4String &  name = "" )

explicit

Definition at line 62 of file G4DipBustGenerator.cc.

  : G4VEmAngularDistribution("DipBustGen")
{}

~G4DipBustGenerator()

G4DipBustGenerator::~G4DipBustGenerator ( )

override

Definition at line 68 of file G4DipBustGenerator.cc.

{}

G4DipBustGenerator() [2/2]

G4DipBustGenerator::G4DipBustGenerator ( const G4DipBustGenerator &  )

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=()

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

delete

PolarAngle()

G4double G4DipBustGenerator::PolarAngle	(	G4double	initial_energy,
		G4double	final_energy,
		G4int	Z
	)

Definition at line 106 of file G4DipBustGenerator.cc.

{
  const G4double cosTheta = SampleCosTheta(eTkin);
  G4double theta = std::acos(cosTheta);
  theta = std::min(std::max(theta, 0.), CLHEP::pi);
  return theta;
}

References G4INCL::Math::max(), G4INCL::Math::min(), CLHEP::pi, and SampleCosTheta().

PrintGeneratorInformation()

void G4DipBustGenerator::PrintGeneratorInformation ( ) const

finalvirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 144 of file G4DipBustGenerator.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Angular Generator based on classical formula from" << G4endl;
  G4cout << "J.D. Jackson, Classical Electrodynamics, Wiley, New York 1975"
     << G4endl;
}

References G4cout, and G4endl.

SampleCosTheta()

G4double G4DipBustGenerator::SampleCosTheta ( G4double  kinEnergy )

private

Definition at line 73 of file G4DipBustGenerator.cc.

{
  const G4double c        = 4. - 8.*G4UniformRand();
  const G4double delta    = 0.5*(std::sqrt(c*c+4.) + std::abs(c));
  const G4double signc    = (c < 0.) ? -1.0 : 1.0;
  // const G4double cofA = -signc*G4Exp(G4Log(delta)/3.0);
  const G4double cofA     = -signc*G4Pow::GetInstance()->A13(delta);
  const G4double cosTheta = std::min(1.,std::max(-1.,cofA - 1./cofA));
  const G4double tau      = kinEnergy/CLHEP::electron_mass_c2;
  const G4double beta     = std::sqrt(tau*(tau + 2.))/(tau + 1.);
 
  return (cosTheta + beta)/(1. + cosTheta*beta);
}

References G4Pow::A13(), anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CLHEP::electron_mass_c2, G4UniformRand, G4Pow::GetInstance(), G4INCL::Math::max(), and G4INCL::Math::min().

Referenced by PolarAngle(), SampleDirection(), and SamplePairDirections().

SampleDirection()

G4ThreeVector & G4DipBustGenerator::SampleDirection ( const G4DynamicParticle *  dp, G4double  out_energy, G4int  Z, const G4Material *  mat = nullptr  )

finalvirtual

Implements G4VEmAngularDistribution.

Definition at line 90 of file G4DipBustGenerator.cc.

{
  const G4double cosTheta = SampleCosTheta(dp->GetKineticEnergy());
 
  const G4double sinTheta = std::sqrt((1. - cosTheta)*(1. + cosTheta));
  const G4double phi      = CLHEP::twopi*G4UniformRand();
 
  fLocalDirection.set(sinTheta*std::cos(phi), sinTheta*std::sin(phi),cosTheta);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
 
  return fLocalDirection;
}

References G4VEmAngularDistribution::fLocalDirection, G4UniformRand, G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), CLHEP::Hep3Vector::rotateUz(), SampleCosTheta(), 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 G4DipBustGenerator::SamplePairDirections ( const G4DynamicParticle *  dp, G4double  elecKinEnergy, G4double  posiKinEnergy, G4ThreeVector &  dirElectron, G4ThreeVector &  dirPositron, G4int  Z = 0, const G4Material *  mat = nullptr  )

finalvirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 118 of file G4DipBustGenerator.cc.

{
  const G4double phi  = CLHEP::twopi * G4UniformRand();
  const G4double sinp = std::sin(phi);
  const G4double cosp = std::cos(phi);
 
  G4double cost = SampleCosTheta(elecKinEnergy);
  G4double sint = std::sqrt((1. - cost)*(1. + cost));
 
  dirElectron.set(sint*cosp, sint*sinp, cost);
  dirElectron.rotateUz(dp->GetMomentumDirection());
 
  cost = SampleCosTheta(posiKinEnergy);
  sint = std::sqrt((1. - cost)*(1. + cost));
 
  dirPositron.set(-sint*cosp, -sint*sinp, cost);
  dirPositron.rotateUz(dp->GetMomentumDirection());
}

References G4UniformRand, G4DynamicParticle::GetMomentumDirection(), CLHEP::Hep3Vector::rotateUz(), SampleCosTheta(), CLHEP::Hep3Vector::set(), and CLHEP::twopi.

Field Documentation

fLocalDirection

G4ThreeVector G4VEmAngularDistribution::fLocalDirection

protectedinherited

Definition at line 103 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::G4VEmAngularDistribution(), G4SauterGavrilaAngularDistribution::SampleDirection(), G4PhotoElectricAngularGeneratorSauterGavrila::SampleDirection(), G4PhotoElectricAngularGeneratorPolarized::SampleDirection(), G4ModifiedMephi::SampleDirection(), G4DeltaAngle::SampleDirection(), G4DeltaAngleFreeScat::SampleDirection(), 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/G4DipBustGenerator.hh
• source/processes/electromagnetic/standard/src/G4DipBustGenerator.cc