#include <G4PhotoElectricAngularGeneratorPolarized.hh>

Inheritance diagram for G4PhotoElectricAngularGeneratorPolarized:

Public Member Functions
	G4PhotoElectricAngularGeneratorPolarized ()

	~G4PhotoElectricAngularGeneratorPolarized ()

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

void	PrintGeneratorInformation () const override

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

	G4PhotoElectricAngularGeneratorPolarized (const G4PhotoElectricAngularGeneratorPolarized &)=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

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

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection

G4bool	fPolarisation

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;
 
    const char* path = G4FindDataDir("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();
  }
}

◆ ~G4PhotoElectricAngularGeneratorPolarized()

G4PhotoElectricAngularGeneratorPolarized::~G4PhotoElectricAngularGeneratorPolarized ( )

Definition at line 127 of file G4PhotoElectricAngularGeneratorPolarized.cc.

128{}

◆ G4PhotoElectricAngularGeneratorPolarized() [2/2]

G4PhotoElectricAngularGeneratorPolarized::G4PhotoElectricAngularGeneratorPolarized ( const G4PhotoElectricAngularGeneratorPolarized & )

delete

Member Function Documentation

◆ operator=()

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

delete

◆ PerpendicularVector()

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

protected

Definition at line 424 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);
  }
}

◆ PrintGeneratorInformation()

void G4PhotoElectricAngularGeneratorPolarized::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 411 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;
}

◆ 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;
}

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

Public Member Functions

Protected Member Functions

Additional Inherited Members