G4Generator2BS Class Reference

#include <G4Generator2BS.hh>

Inheritance diagram for G4Generator2BS:

Public Member Functions
	G4Generator2BS (const G4Generator2BS &)=delete
	G4Generator2BS (const G4String &name="")
const G4String &	GetName () const
G4Generator2BS &	operator= (const G4Generator2BS &right)=delete
void	PrintGeneratorInformation () const override
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double out_energy, G4int Z, 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)
virtual	~G4Generator2BS ()

Protected Member Functions
G4double	RejectionFunction (G4double value) const

Protected Attributes
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Private Attributes
G4double	delta
G4String	fName
G4double	fz
G4Pow *	g4pow
G4int	nwarn
G4double	ratio
G4double	ratio1
G4double	ratio2

Detailed Description

Definition at line 64 of file G4Generator2BS.hh.

Constructor & Destructor Documentation

G4Generator2BS() [1/2]

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

explicit

Definition at line 63 of file G4Generator2BS.cc.

  : G4VEmAngularDistribution("AngularGen2BS"),fz(1.),ratio(1.),
    ratio1(1.),ratio2(1.),delta(0.)
{
  g4pow = G4Pow::GetInstance();
  nwarn = 0;
}

References g4pow, G4Pow::GetInstance(), and nwarn.

~G4Generator2BS()

G4Generator2BS::~G4Generator2BS ( )

virtual

Definition at line 71 of file G4Generator2BS.cc.

{}

G4Generator2BS() [2/2]

G4Generator2BS::G4Generator2BS ( const G4Generator2BS &  )

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

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

delete

PrintGeneratorInformation()

void G4Generator2BS::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 134 of file G4Generator2BS.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Bremsstrahlung Angular Generator is 2BS Generator "
     << "from 2BS Koch & Motz distribution (Rev Mod Phys 31(4), 920 (1959))" << G4endl;
  G4cout << "Sampling algorithm adapted from PIRS-0203" << G4endl;
  G4cout << "\n" << G4endl;
} 

References G4cout, and G4endl.

RejectionFunction()

G4double G4Generator2BS::RejectionFunction ( G4double  value ) const

inlineprotected

Definition at line 99 of file G4Generator2BS.hh.

{
  G4double y2 = (1 + y)*(1 + y);
  G4double x  = 4*y*ratio/y2;
  return 4*x - ratio1 - (ratio2 - x)*G4Log(delta + fz/y2); 
}

References delta, fz, G4Log(), ratio, ratio1, and ratio2.

Referenced by SampleDirection().

SampleDirection()

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

overridevirtual

Implements G4VEmAngularDistribution.

Definition at line 74 of file G4Generator2BS.cc.

{
  // Adapted from "Improved bremsstrahlung photon angular sampling in the EGS4 code system"
  // by Alex F. Bielajew, Rahde Mohan anc Chen-Shou Chui, PIRS-0203
  // Ionizing Radiation Standards
  // Institute for National Measurement Standards 
  // National Research Council of Canada
  // Departement of Medical Physics, Memorial Sloan-Kettering Cancer Center, New York
 
  G4double energy = dp->GetTotalEnergy();
  ratio = final_energy/energy;
  ratio1 = (1 + ratio)*(1 + ratio);
  ratio2 = 1 + ratio*ratio;
 
  G4double gamma = energy/electron_mass_c2;
  G4double beta  = std::sqrt((gamma - 1)*(gamma + 1))/gamma;
 
  // VI speadup
  fz = 0.00008116224*g4pow->Z13(Z)*g4pow->Z13(Z+1);
 
  // majoranta
  G4double ymax = 2*beta*(1 + beta)*gamma*gamma;
  G4double gMax = RejectionFunction(0.0);
  gMax = std::max(gMax,RejectionFunction(ymax));
 
  G4double y, gfun;
 
  do{
    G4double q = G4UniformRand();
    y = q*ymax/(1 + ymax*(1 - q));
    gfun = RejectionFunction(y);
 
    // violation point
    if(gfun > gMax && nwarn >= 20) {
      ++nwarn;
      G4cout << "### WARNING in G4Generator2BS: Etot(MeV)= " << energy/MeV 
         << "  Egamma(MeV)" << (energy - final_energy)/MeV
         << " gMax= " << gMax << "  < " << gfun
         << "  results are not reliable!" 
         << G4endl;
      if(20 == nwarn) { 
    G4cout << "   WARNING in G4Generator2BS is closed" << G4endl; 
      }
    }
 
  } while(G4UniformRand()*gMax > gfun || y > ymax);
  
  G4double cost = 1 - 2*y/ymax;
  G4double sint = std::sqrt((1 - cost)*(1 + cost));
  G4double phi  = twopi*G4UniformRand(); 
 
  fLocalDirection.set(sint*std::cos(phi), sint*std::sin(phi),cost);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
 
  return fLocalDirection;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::electron_mass_c2, G4INCL::KinematicsUtils::energy(), G4VEmAngularDistribution::fLocalDirection, fz, G4cout, G4endl, g4pow, G4UniformRand, G4DynamicParticle::GetMomentumDirection(), G4DynamicParticle::GetTotalEnergy(), G4INCL::Math::max(), MeV, nwarn, ratio, ratio1, ratio2, RejectionFunction(), CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), twopi, Z, and G4Pow::Z13().

Referenced by G4Generator2BN::SampleDirection().

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

delta

G4double G4Generator2BS::delta

private

Definition at line 93 of file G4Generator2BS.hh.

Referenced by RejectionFunction().

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(), G4DipBustGenerator::SampleDirection(), G4ModifiedTsai::SampleDirection(), G4Generator2BN::SampleDirection(), 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().

fz

G4double G4Generator2BS::fz

private

Definition at line 89 of file G4Generator2BS.hh.

Referenced by RejectionFunction(), and SampleDirection().

g4pow

G4Pow* G4Generator2BS::g4pow

private

Definition at line 87 of file G4Generator2BS.hh.

Referenced by G4Generator2BS(), and SampleDirection().

nwarn

G4int G4Generator2BS::nwarn

private

Definition at line 95 of file G4Generator2BS.hh.

Referenced by G4Generator2BS(), and SampleDirection().

ratio

G4double G4Generator2BS::ratio

private

Definition at line 90 of file G4Generator2BS.hh.

Referenced by RejectionFunction(), and SampleDirection().

ratio1

G4double G4Generator2BS::ratio1

private

Definition at line 91 of file G4Generator2BS.hh.

Referenced by RejectionFunction(), and SampleDirection().

ratio2

G4double G4Generator2BS::ratio2

private

Definition at line 92 of file G4Generator2BS.hh.

Referenced by RejectionFunction(), and SampleDirection().

---

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

• source/processes/electromagnetic/lowenergy/include/G4Generator2BS.hh
• source/processes/electromagnetic/lowenergy/src/G4Generator2BS.cc