G4ModifiedTsai Class Reference

#include <G4ModifiedTsai.hh>

Inheritance diagram for G4ModifiedTsai:

Public Member Functions
	G4ModifiedTsai (const G4ModifiedTsai &)=delete
	G4ModifiedTsai (const G4String &name="")
const G4String &	GetName () const
G4ModifiedTsai &	operator= (const G4ModifiedTsai &right)=delete
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
	~G4ModifiedTsai () override

Protected Attributes
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Private Member Functions
G4double	SampleCosTheta (G4double kinEnergy)

Private Attributes
G4String	fName

Detailed Description

Definition at line 62 of file G4ModifiedTsai.hh.

Constructor & Destructor Documentation

G4ModifiedTsai() [1/2]

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

explicit

Definition at line 65 of file G4ModifiedTsai.cc.

  : G4VEmAngularDistribution("ModifiedTsai")
{}    

~G4ModifiedTsai()

G4ModifiedTsai::~G4ModifiedTsai ( )

override

Definition at line 71 of file G4ModifiedTsai.cc.

{}

G4ModifiedTsai() [2/2]

G4ModifiedTsai::G4ModifiedTsai ( const G4ModifiedTsai &  )

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

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

delete

PrintGeneratorInformation()

void G4ModifiedTsai::PrintGeneratorInformation ( ) const

finalvirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 144 of file G4ModifiedTsai.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Bremsstrahlung Angular Generator is Modified Tsai" << G4endl;
  G4cout << "Distribution suggested by L.Urban (Geant3 manual (1993) Phys211)" 
         << G4endl;
  G4cout << "Derived from Tsai distribution (Rev Mod Phys 49,421(1977)) \n" 
         << G4endl;
} 

References G4cout, and G4endl.

SampleCosTheta()

G4double G4ModifiedTsai::SampleCosTheta ( G4double  kinEnergy )

private

Definition at line 93 of file G4ModifiedTsai.cc.

{
  // Universal distribution suggested by L. Urban (Geant3 manual (1993) 
  // Phys211) derived from Tsai distribution (Rev Mod Phys 49,421(1977))
 
  G4double uMax = 2*(1. + kinEnergy/CLHEP::electron_mass_c2);   
 
  static const G4double a1     = 1.6;
  static const G4double a2     = a1/3.;
  static const G4double border = 0.25;
  G4double u;
  CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
 
  do {
    G4double uu = -G4Log(rndmEngine->flat()*rndmEngine->flat());
    u = (border > rndmEngine->flat()) ? uu*a1 : uu*a2;
    
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  } while(u > uMax);
 
  return 1.0 - 2.0*u*u/(uMax*uMax);
}

References CLHEP::electron_mass_c2, CLHEP::HepRandomEngine::flat(), and G4Log().

Referenced by SampleDirection(), and SamplePairDirections().

SampleDirection()

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

finalvirtual

Implements G4VEmAngularDistribution.

Definition at line 77 of file G4ModifiedTsai.cc.

{
  // Sample gamma angle (Z - axis along the parent particle).
  G4double cost = SampleCosTheta(dp->GetKineticEnergy());
  G4double sint = std::sqrt((1 - cost)*(1 + cost));
  G4double phi  = CLHEP::twopi*G4UniformRand(); 
 
  fLocalDirection.set(sint*std::cos(phi), sint*std::sin(phi), cost);
  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 G4ModifiedTsai::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 G4ModifiedTsai.cc.

{
  G4double phi  = CLHEP::twopi * G4UniformRand();
  G4double sinp = std::sin(phi);
  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(), G4DipBustGenerator::SampleDirection(), 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/G4ModifiedTsai.hh
• source/processes/electromagnetic/standard/src/G4ModifiedTsai.cc