G4DNABornAngle Class Reference

#include <G4DNABornAngle.hh>

Inheritance diagram for G4DNABornAngle:

Public Member Functions
	G4DNABornAngle (const G4String &name="")
	~G4DNABornAngle () override
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, const G4Material *mat=nullptr) override
G4ThreeVector &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double kinEnergyFinal, G4int Z, G4int shellIdx, const G4Material *mat=nullptr) override
void	PrintGeneratorInformation () const override
G4DNABornAngle &	operator= (const G4DNABornAngle &right)=delete
	G4DNABornAngle (const G4DNABornAngle &)=delete
Public Member Functions inherited from G4VEmAngularDistribution
	G4VEmAngularDistribution (const G4String &name)
virtual	~G4VEmAngularDistribution ()
virtual G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
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 void	PrintGeneratorInformation () const
const G4String &	GetName () const
G4VEmAngularDistribution &	operator= (const G4VEmAngularDistribution &right)=delete
	G4VEmAngularDistribution (const G4VEmAngularDistribution &)=delete

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Detailed Description

Definition at line 56 of file G4DNABornAngle.hh.

Constructor & Destructor Documentation

G4DNABornAngle() [1/2]

G4DNABornAngle::G4DNABornAngle

(

const G4String &

name = ""

)

Definition at line 58 of file G4DNABornAngle.cc.

  : G4VEmAngularDistribution("deltaBorn")
{
  fElectron = G4Electron::Electron();
}    

~G4DNABornAngle()

G4DNABornAngle::~G4DNABornAngle ( )

override

Definition at line 64 of file G4DNABornAngle.cc.

{}

G4DNABornAngle() [2/2]

G4DNABornAngle::G4DNABornAngle ( const G4DNABornAngle &  )

delete

Member Function Documentation

operator=()

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

delete

PrintGeneratorInformation()

void G4DNABornAngle::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 120 of file G4DNABornAngle.cc.

{} 

SampleDirection()

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

overridevirtual

Implements G4VEmAngularDistribution.

Definition at line 113 of file G4DNABornAngle.cc.

{
  return SampleDirectionForShell(dp, secEkin, Z, 0, mat);
}

SampleDirectionForShell()

G4ThreeVector & G4DNABornAngle::SampleDirectionForShell ( const G4DynamicParticle *  dp, G4double  kinEnergyFinal, G4int  Z, G4int  shellIdx, const G4Material *  mat = nullptr  )

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 68 of file G4DNABornAngle.cc.

{
  G4double k = dp->GetKineticEnergy();
  G4double cosTheta = 1.0;
  if (dp->GetDefinition() == fElectron)
    {
      if (secKinetic < 50.*eV) cosTheta = (2.*G4UniformRand())-1.;
      else if (secKinetic <= 200.*eV)
    {
      if (G4UniformRand() <= 0.1) cosTheta = (2.*G4UniformRand())-1.;
      else cosTheta = G4UniformRand()*(std::sqrt(2.)/2);
    }
      else
    {
      G4double sin2O = 
        (1.-secKinetic/k)/(1.+secKinetic/(2.*electron_mass_c2));
      cosTheta = std::sqrt(1.-sin2O);
    }
    }
  else 
    {
      G4double mass = dp->GetDefinition()->GetPDGMass();
      G4double maxSecKinetic = 4.* (electron_mass_c2 / mass) * k;
 
      // cosTheta = std::sqrt(secKinetic / maxSecKinetic);
 
      // Restriction below 100 eV from Emfietzoglou (2000)
 
      if (secKinetic>100*eV) cosTheta = std::sqrt(secKinetic / maxSecKinetic);
      else cosTheta = (2.*G4UniformRand())-1.;
    }
 
  G4double sint = sqrt((1.0 - cosTheta)*(1.0 + cosTheta));
  G4double phi  = twopi*G4UniformRand(); 
 
  fLocalDirection.set(sint*cos(phi), sint*sin(phi), cosTheta);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
 
  return fLocalDirection;
}

Referenced by SampleDirection().

---

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

• G4DNABornAngle.hh
• G4DNABornAngle.cc