G4PhotoElectricAngularGeneratorPolarized Class Reference

#include <G4PhotoElectricAngularGeneratorPolarized.hh>

Inheritance diagram for G4PhotoElectricAngularGeneratorPolarized:

Public Member Functions
	G4PhotoElectricAngularGeneratorPolarized ()
	G4PhotoElectricAngularGeneratorPolarized (const G4PhotoElectricAngularGeneratorPolarized &)=delete
const G4String &	GetName () const
G4PhotoElectricAngularGeneratorPolarized &	operator= (const G4PhotoElectricAngularGeneratorPolarized &right)=delete
void	PrintGeneratorInformation () const override
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double eKinEnergy, G4int shellId, 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)
	~G4PhotoElectricAngularGeneratorPolarized ()

Protected Member Functions
G4ThreeVector	PerpendicularVector (const G4ThreeVector &a) const

Protected Attributes
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

Private Member Functions
G4double	CrossSectionMajorantFunction (G4double theta, G4double cBeta) const
G4double	DSigmaKshellGavrila1959 (G4double beta, G4double theta, G4double phi) const
G4double	DSigmaL1shellGavrila (G4double beta, G4double theta, G4double phi) const
G4ThreeVector	PhotoElectronComputeFinalDirection (const G4RotationMatrix &rotation, G4double theta, G4double phi) const
void	PhotoElectronGeneratePhiAndTheta (G4int shellId, G4double beta, G4double aBeta, G4double cBeta, G4double *pphi, G4double *ptheta) const
void	PhotoElectronGetMajorantSurfaceAandCParameters (G4int shellId, G4double beta, G4double *majorantSurfaceParameterA, G4double *majorantSurfaceParameterC) const
G4RotationMatrix	PhotoElectronRotationMatrix (const G4ThreeVector &direction, const G4ThreeVector &polarization)

Private Attributes
G4double	aMajorantSurfaceParameterTable [980][2]
G4double	betaArray [3]
G4double	cMajorantSurfaceParameterTable [980][2]
G4String	fName

Detailed Description

Definition at line 52 of file G4PhotoElectricAngularGeneratorPolarized.hh.

Constructor & Destructor Documentation

G4PhotoElectricAngularGeneratorPolarized() [1/2]

G4PhotoElectricAngularGeneratorPolarized::G4PhotoElectricAngularGeneratorPolarized ( )

explicit

Definition at line 71 of file G4PhotoElectricAngularGeneratorPolarized.cc.

  :G4VEmAngularDistribution("AngularGenSauterGavrilaPolarized")
{
  const G4int arrayDim = 980;
 
  //minimum electron beta parameter allowed
  betaArray[0] = 0.02;
  //beta step
  betaArray[1] = 0.001;
  //maximum index array for a and c tables
  betaArray[2] = arrayDim - 1;
 
  // read Majorant Surface Parameters. This are required in order to generate 
  //Gavrila angular photoelectron distribution
  for(G4int level = 0; level < 2; level++){
    char nameChar0[100] = "ftab0.dat"; // K-shell Majorant Surface Parameters
    char nameChar1[100] = "ftab1.dat"; // L-shell Majorant Surface Parameters
 
    G4String filename;
    if(level == 0) filename = nameChar0;
    if(level == 1) filename = nameChar1;
 
    char* path = std::getenv("G4LEDATA");
    if (!path)
      {
        G4String excep = "G4EMDataSet - G4LEDATA environment variable not set";
        G4Exception("G4PhotoElectricAngularGeneratorPolarized::G4PhotoElectricAngularGeneratorPolarized",
            "em0006",FatalException,"G4LEDATA environment variable not set");
    return;
      }
 
    G4String pathString(path);
    G4String dirFile = pathString + "/photoelectric_angular/" + filename;
    std::ifstream infile(dirFile);
    if (!infile.is_open())
      {
    G4String excep = "data file: " + dirFile + " not found";
        G4Exception("G4PhotoElectricAngularGeneratorPolarized::G4PhotoElectricAngularGeneratorPolarized",
            "em0003",FatalException,excep);
    return;
      }
 
    // Read parameters into tables. The parameters are function of incident electron 
    // energy and shell level
    G4float aRead=0,cRead=0, beta=0;
    for(G4int i=0 ; i<arrayDim ;++i){     
      infile >> beta >> aRead >> cRead;
      aMajorantSurfaceParameterTable[i][level] = aRead;
      cMajorantSurfaceParameterTable[i][level] = cRead;
    }
    infile.close();
  }
}

References aMajorantSurfaceParameterTable, anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, betaArray, cMajorantSurfaceParameterTable, FatalException, and G4Exception().

~G4PhotoElectricAngularGeneratorPolarized()

G4PhotoElectricAngularGeneratorPolarized::~G4PhotoElectricAngularGeneratorPolarized

(

)

Definition at line 127 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{}

G4PhotoElectricAngularGeneratorPolarized() [2/2]

G4PhotoElectricAngularGeneratorPolarized::G4PhotoElectricAngularGeneratorPolarized ( const G4PhotoElectricAngularGeneratorPolarized &  )

delete

Member Function Documentation

CrossSectionMajorantFunction()

G4double G4PhotoElectricAngularGeneratorPolarized::CrossSectionMajorantFunction ( G4double  theta, G4double  cBeta  ) const

private

Definition at line 221 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  // Compute Majorant Function
  G4double crossSectionMajorantFunctionValue = 0;
  crossSectionMajorantFunctionValue = theta/(1+cBeta*theta*theta);
  return crossSectionMajorantFunctionValue;
}

Referenced by PhotoElectronGeneratePhiAndTheta().

DSigmaKshellGavrila1959()

G4double G4PhotoElectricAngularGeneratorPolarized::DSigmaKshellGavrila1959 ( G4double  beta, G4double  theta, G4double  phi  ) const

private

Definition at line 233 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  //Double differential K shell cross-section (Gavrila 1959)
 
  G4double beta2 = beta*beta;
  G4double oneBeta2 = 1 - beta2;
  G4double sqrtOneBeta2 = std::sqrt(oneBeta2);
  G4double oneBeta2_to_3_2 = std::pow(oneBeta2,1.5);
  G4double cosTheta = std::cos(theta);
  G4double sinTheta2 = std::sin(theta)*std::sin(theta);
  G4double cosPhi2 = std::cos(phi)*std::cos(phi);
  G4double oneBetaCosTheta = 1-beta*cosTheta;
  G4double dsigma = 0;
  G4double firstTerm = 0;
  G4double secondTerm = 0;
 
  firstTerm = sinTheta2*cosPhi2/std::pow(oneBetaCosTheta,4)-(1 - sqrtOneBeta2)/(2*oneBeta2) *
              (sinTheta2 * cosPhi2)/std::pow(oneBetaCosTheta,3) + (1-sqrtOneBeta2)*
              (1-sqrtOneBeta2)/(4*oneBeta2_to_3_2) * sinTheta2/std::pow(oneBetaCosTheta,3);
 
  secondTerm = std::sqrt(1 - sqrtOneBeta2)/(std::pow(2.,3.5)*beta2*std::pow(oneBetaCosTheta,2.5)) *
               (4*beta2/sqrtOneBeta2 * sinTheta2*cosPhi2/oneBetaCosTheta + 4*beta/oneBeta2 * cosTheta * cosPhi2
               - 4*(1-sqrtOneBeta2)/oneBeta2 *(1+cosPhi2) - beta2 * (1-sqrtOneBeta2)/oneBeta2 * sinTheta2/oneBetaCosTheta
               + 4*beta2*(1-sqrtOneBeta2)/oneBeta2_to_3_2 - 4*beta*(1-sqrtOneBeta2)*(1-sqrtOneBeta2)/oneBeta2_to_3_2 * cosTheta)
               + (1-sqrtOneBeta2)/(4*beta2*oneBetaCosTheta*oneBetaCosTheta) * (beta/oneBeta2 - 2/oneBeta2 * cosTheta * cosPhi2 +
               (1-sqrtOneBeta2)/oneBeta2_to_3_2 * cosTheta - beta * (1-sqrtOneBeta2)/oneBeta2_to_3_2);
 
  dsigma = ( firstTerm*(1-pi*fine_structure_const/beta) + secondTerm*(pi*fine_structure_const) )*std::sin(theta);
 
  return dsigma;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::fine_structure_const, and pi.

Referenced by PhotoElectronGeneratePhiAndTheta().

DSigmaL1shellGavrila()

G4double G4PhotoElectricAngularGeneratorPolarized::DSigmaL1shellGavrila ( G4double  beta, G4double  theta, G4double  phi  ) const

private

Definition at line 269 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  //Double differential L1 shell cross-section (Gavrila 1961)
  G4double beta2 = beta*beta;
  G4double oneBeta2 = 1-beta2;
  G4double sqrtOneBeta2 = std::sqrt(oneBeta2);
  G4double oneBeta2_to_3_2=std::pow(oneBeta2,1.5);
  G4double cosTheta = std::cos(theta);
  G4double sinTheta2 =std::sin(theta)*std::sin(theta);
  G4double cosPhi2 = std::cos(phi)*std::cos(phi);
  G4double oneBetaCosTheta = 1-beta*cosTheta;
 
  G4double dsigma = 0;
  G4double firstTerm = 0;
  G4double secondTerm = 0;
 
  firstTerm = sinTheta2*cosPhi2/std::pow(oneBetaCosTheta,4)-(1 - sqrtOneBeta2)/(2*oneBeta2)
              *  (sinTheta2 * cosPhi2)/std::pow(oneBetaCosTheta,3) + (1-sqrtOneBeta2)*
              (1-sqrtOneBeta2)/(4*oneBeta2_to_3_2) * sinTheta2/std::pow(oneBetaCosTheta,3);
 
  secondTerm = std::sqrt(1 - sqrtOneBeta2)/(std::pow(2.,3.5)*beta2*std::pow(oneBetaCosTheta,2.5)) *
               (4*beta2/sqrtOneBeta2 * sinTheta2*cosPhi2/oneBetaCosTheta + 4*beta/oneBeta2 * cosTheta * cosPhi2
               - 4*(1-sqrtOneBeta2)/oneBeta2 *(1+cosPhi2) - beta2 * (1-sqrtOneBeta2)/oneBeta2 * sinTheta2/oneBetaCosTheta
               + 4*beta2*(1-sqrtOneBeta2)/oneBeta2_to_3_2 - 4*beta*(1-sqrtOneBeta2)*(1-sqrtOneBeta2)/oneBeta2_to_3_2 * cosTheta)
               + (1-sqrtOneBeta2)/(4*beta2*oneBetaCosTheta*oneBetaCosTheta) * (beta/oneBeta2 - 2/oneBeta2 * cosTheta * cosPhi2 +
               (1-sqrtOneBeta2)/oneBeta2_to_3_2*cosTheta - beta*(1-sqrtOneBeta2)/oneBeta2_to_3_2);
 
  dsigma = ( firstTerm*(1-pi*fine_structure_const/beta) + secondTerm*(pi*fine_structure_const) )*std::sin(theta);
 
  return dsigma;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::fine_structure_const, and pi.

Referenced by PhotoElectronGeneratePhiAndTheta().

GetName()

const G4String & G4VEmAngularDistribution::GetName ( ) const

inlineinherited

Definition at line 111 of file G4VEmAngularDistribution.hh.

{
  return fName;
}

References G4VEmAngularDistribution::fName.

operator=()

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

delete

PerpendicularVector()

G4ThreeVector G4PhotoElectricAngularGeneratorPolarized::PerpendicularVector ( const G4ThreeVector &  a ) const

protected

Definition at line 423 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  G4double dx = a.x();
  G4double dy = a.y();
  G4double dz = a.z();
  G4double x = dx < 0.0 ? -dx : dx;
  G4double y = dy < 0.0 ? -dy : dy;
  G4double z = dz < 0.0 ? -dz : dz;
  if (x < y) {
    return x < z ? G4ThreeVector(-dy,dx,0) : G4ThreeVector(0,-dz,dy);
  }else{
    return y < z ? G4ThreeVector(dz,0,-dx) : G4ThreeVector(-dy,dx,0);
  }
}

References CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by PhotoElectronRotationMatrix().

PhotoElectronComputeFinalDirection()

G4ThreeVector G4PhotoElectricAngularGeneratorPolarized::PhotoElectronComputeFinalDirection ( const G4RotationMatrix &  rotation, G4double  theta, G4double  phi  ) const

private

Definition at line 394 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  //computes the photoelectron momentum unitary vector
  G4double sint = std::sin(theta);
  G4double px = std::cos(phi)*sint;
  G4double py = std::sin(phi)*sint;
  G4double pz = std::cos(theta);
 
  G4ThreeVector samplingDirection(px,py,pz);
 
  G4ThreeVector outgoingDirection = rotation*samplingDirection;
  return outgoingDirection;
}

Referenced by SampleDirection().

PhotoElectronGeneratePhiAndTheta()

void G4PhotoElectricAngularGeneratorPolarized::PhotoElectronGeneratePhiAndTheta ( G4int  shellId, G4double  beta, G4double  aBeta, G4double  cBeta, G4double *  pphi, G4double *  ptheta  ) const

private

Definition at line 176 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  G4double rand1, rand2, rand3 = 0;
  G4double phi = 0;
  G4double theta = 0;
  G4double crossSectionValue = 0;
  G4double crossSectionMajorantFunctionValue = 0;
  G4double maxBeta = 0;
 
  do {
 
    rand1 = G4UniformRand();
    rand2 = G4UniformRand();
    rand3 = G4UniformRand();
 
    phi=2*pi*rand1;
 
    if(shellLevel == 0){
      // Polarized Gavrila Cross-Section for K-shell (1959)
      theta=std::sqrt(((G4Exp(rand2*std::log(1+cBeta*pi*pi)))-1)/cBeta);
      crossSectionMajorantFunctionValue =
    CrossSectionMajorantFunction(theta, cBeta);
      crossSectionValue = DSigmaKshellGavrila1959(beta, theta, phi);
    } else {
      //  Polarized Gavrila Cross-Section for other shells (L1-shell) (1961)
      theta = std::sqrt(((G4Exp(rand2*std::log(1+cBeta*pi*pi)))-1)/cBeta);
      crossSectionMajorantFunctionValue =
    CrossSectionMajorantFunction(theta, cBeta);
      crossSectionValue = DSigmaL1shellGavrila(beta, theta, phi);
    }
 
    maxBeta=rand3*aBeta*crossSectionMajorantFunctionValue;
    if(crossSectionValue < 0.0) { crossSectionValue = maxBeta; }
 
  } while(maxBeta > crossSectionValue || theta > CLHEP::pi);
 
  *pphi = phi;
  *ptheta = theta;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CrossSectionMajorantFunction(), DSigmaKshellGavrila1959(), DSigmaL1shellGavrila(), G4Exp(), G4UniformRand, CLHEP::pi, and pi.

Referenced by SampleDirection().

PhotoElectronGetMajorantSurfaceAandCParameters()

void G4PhotoElectricAngularGeneratorPolarized::PhotoElectronGetMajorantSurfaceAandCParameters ( G4int  shellId, G4double  beta, G4double *  majorantSurfaceParameterA, G4double *  majorantSurfaceParameterC  ) const

private

Definition at line 348 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  // This member function finds for a given shell and beta value
  // of the outgoing electron the correct Majorant Surface parameters
  G4double aBeta,cBeta;
  G4double bMin,bStep;
  G4int indexMax;
  G4int level = 0;
  if(shellId > 0) { level = 1; }
 
  bMin = betaArray[0];
  bStep = betaArray[1];
  indexMax = (G4int)betaArray[2];
  const G4double kBias = 1e-9;
 
  G4int k = (G4int)((beta-bMin+kBias)/bStep);
 
  if(k < 0)
    k = 0;
  if(k > indexMax)
    k = indexMax;
 
  if(k == 0)
    aBeta = std::max(aMajorantSurfaceParameterTable[k][level],aMajorantSurfaceParameterTable[k+1][level]);
  else if(k==indexMax)
    aBeta = std::max(aMajorantSurfaceParameterTable[k-1][level],aMajorantSurfaceParameterTable[k][level]);
  else{
    aBeta = std::max(aMajorantSurfaceParameterTable[k-1][level],aMajorantSurfaceParameterTable[k][level]);
    aBeta = std::max(aBeta,aMajorantSurfaceParameterTable[k+1][level]);
  }
 
  if(k == 0)
    cBeta = std::max(cMajorantSurfaceParameterTable[k][level],cMajorantSurfaceParameterTable[k+1][level]);
  else if(k == indexMax)
    cBeta = std::max(cMajorantSurfaceParameterTable[k-1][level],cMajorantSurfaceParameterTable[k][level]);
  else{
    cBeta = std::max(cMajorantSurfaceParameterTable[k-1][level],cMajorantSurfaceParameterTable[k][level]);
    cBeta = std::max(cBeta,cMajorantSurfaceParameterTable[k+1][level]);
  }
 
  *majorantSurfaceParameterA = aBeta;
  *majorantSurfaceParameterC = cBeta;
}

References aMajorantSurfaceParameterTable, anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, betaArray, cMajorantSurfaceParameterTable, and G4INCL::Math::max().

Referenced by SampleDirection().

PhotoElectronRotationMatrix()

G4RotationMatrix G4PhotoElectricAngularGeneratorPolarized::PhotoElectronRotationMatrix ( const G4ThreeVector &  direction, const G4ThreeVector &  polarization  )

private

Definition at line 305 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  G4double mK = direction.mag();
  G4double mS = polarization.mag();
  G4ThreeVector polarization2 = polarization;
  const G4double kTolerance = 1e-6;
 
  if(!(polarization.isOrthogonal(direction,kTolerance)) || mS == 0){
    G4ThreeVector d0 = direction.unit();
    G4ThreeVector a1 = PerpendicularVector(d0);
    G4ThreeVector a0 = a1.unit();
    G4double rand1 = G4UniformRand();
    G4double angle = twopi*rand1;
    G4ThreeVector b0 = d0.cross(a0);
    G4ThreeVector c;
    c.setX(std::cos(angle)*(a0.x())+std::sin(angle)*b0.x());
    c.setY(std::cos(angle)*(a0.y())+std::sin(angle)*b0.y());
    c.setZ(std::cos(angle)*(a0.z())+std::sin(angle)*b0.z());
    polarization2 = c.unit();
    mS = polarization2.mag();
  }else
    {
      if ( polarization.howOrthogonal(direction) != 0)
    {
      polarization2 = polarization
        - polarization.dot(direction)/direction.dot(direction) * direction;
    }
    }
 
  G4ThreeVector direction2 = direction/mK;
  polarization2 = polarization2/mS;
 
  G4ThreeVector y = direction2.cross(polarization2);
 
  G4RotationMatrix R(polarization2,y,direction2);
  return R;
}

References a0, angle, CLHEP::Hep3Vector::cross(), CLHEP::Hep3Vector::dot(), G4UniformRand, CLHEP::Hep3Vector::howOrthogonal(), CLHEP::Hep3Vector::isOrthogonal(), CLHEP::Hep3Vector::mag(), anonymous_namespace{G4ChipsKaonMinusElasticXS.cc}::mK, PerpendicularVector(), CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), CLHEP::Hep3Vector::setZ(), twopi, CLHEP::Hep3Vector::unit(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by SampleDirection().

PrintGeneratorInformation()

void G4PhotoElectricAngularGeneratorPolarized::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 410 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Polarized Photoelectric Angular Generator" << G4endl;
  G4cout << "PhotoElectric Electron Angular Generator based on the general Gavrila photoelectron angular distribution" << G4endl;
  G4cout << "Includes polarization effects for K and L1 atomic shells, according to Gavrilla (1959, 1961)." << G4endl;
  G4cout << "For higher shells the L1 cross-section is used." << G4endl;
  G4cout << "(see Physics Reference Manual) \n" << G4endl;
}

References G4cout, and G4endl.

SampleDirection()

G4ThreeVector & G4PhotoElectricAngularGeneratorPolarized::SampleDirection ( const G4DynamicParticle *  dp, G4double  eKinEnergy, G4int  shellId, const G4Material *  mat = nullptr  )

overridevirtual

Implements G4VEmAngularDistribution.

Definition at line 133 of file G4PhotoElectricAngularGeneratorPolarized.cc.

{
  // (shellId == 0) - K-shell  - Polarized model for K-shell
  // (shellId > 0)  - L1-shell - Polarized model for L1 and higher shells
 
  // Calculate Lorentz term (gamma) and beta parameters
  G4double gamma = 1 + eKinEnergy/electron_mass_c2;
  G4double beta  = std::sqrt((gamma - 1)*(gamma + 1))/gamma;
 
  const G4ThreeVector& direction = dp->GetMomentumDirection();
  const G4ThreeVector& polarization = dp->GetPolarization();
 
  G4double theta, phi = 0;
  // Majorant surface parameters
  // function of the outgoing electron kinetic energy
  G4double aBeta = 0;
  G4double cBeta = 0;
 
  // For the outgoing kinetic energy
  // find the current majorant surface parameters
  PhotoElectronGetMajorantSurfaceAandCParameters(shellId, beta, &aBeta, &cBeta);
 
  // Generate pho and theta according to the shell level
  // and beta parameter of the electron
  PhotoElectronGeneratePhiAndTheta(shellId, beta, aBeta, cBeta, &phi, &theta);
 
  // Determine the rotation matrix
  const G4RotationMatrix rotation =
    PhotoElectronRotationMatrix(direction, polarization);
 
  // Compute final direction of the outgoing electron
  fLocalDirection = PhotoElectronComputeFinalDirection(rotation, theta, phi);
 
  return fLocalDirection;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::electron_mass_c2, G4VEmAngularDistribution::fLocalDirection, G4DynamicParticle::GetMomentumDirection(), G4DynamicParticle::GetPolarization(), PhotoElectronComputeFinalDirection(), PhotoElectronGeneratePhiAndTheta(), PhotoElectronGetMajorantSurfaceAandCParameters(), and PhotoElectronRotationMatrix().

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

aMajorantSurfaceParameterTable

G4double G4PhotoElectricAngularGeneratorPolarized::aMajorantSurfaceParameterTable[980][2]

private

Definition at line 94 of file G4PhotoElectricAngularGeneratorPolarized.hh.

Referenced by G4PhotoElectricAngularGeneratorPolarized(), and PhotoElectronGetMajorantSurfaceAandCParameters().

betaArray

G4double G4PhotoElectricAngularGeneratorPolarized::betaArray[3]

private

Definition at line 93 of file G4PhotoElectricAngularGeneratorPolarized.hh.

Referenced by G4PhotoElectricAngularGeneratorPolarized(), and PhotoElectronGetMajorantSurfaceAandCParameters().

cMajorantSurfaceParameterTable

G4double G4PhotoElectricAngularGeneratorPolarized::cMajorantSurfaceParameterTable[980][2]

private

Definition at line 95 of file G4PhotoElectricAngularGeneratorPolarized.hh.

Referenced by G4PhotoElectricAngularGeneratorPolarized(), and PhotoElectronGetMajorantSurfaceAandCParameters().

fLocalDirection

G4ThreeVector G4VEmAngularDistribution::fLocalDirection

protectedinherited

Definition at line 103 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::G4VEmAngularDistribution(), G4SauterGavrilaAngularDistribution::SampleDirection(), G4PhotoElectricAngularGeneratorSauterGavrila::SampleDirection(), 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:

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