G4RDPhotoElectricAngularGeneratorPolarized Class Reference

#include <G4RDPhotoElectricAngularGeneratorPolarized.hh>

Inheritance diagram for G4RDPhotoElectricAngularGeneratorPolarized:

Public Member Functions
	G4RDPhotoElectricAngularGeneratorPolarized (const G4String &name)
	~G4RDPhotoElectricAngularGeneratorPolarized ()
G4ThreeVector	GetPhotoElectronDirection (const G4ThreeVector &direction, const G4double kineticEnergy, const G4ThreeVector &polarization, const G4int shellId) const
void	PrintGeneratorInformation () const
Public Member Functions inherited from G4RDVPhotoElectricAngularDistribution
	G4RDVPhotoElectricAngularDistribution (const G4String &name)
virtual	~G4RDVPhotoElectricAngularDistribution ()

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

Detailed Description

Definition at line 54 of file G4RDPhotoElectricAngularGeneratorPolarized.hh.

Constructor & Destructor Documentation

G4RDPhotoElectricAngularGeneratorPolarized::G4RDPhotoElectricAngularGeneratorPolarized

(

const G4String &

name

)

Definition at line 71 of file G4RDPhotoElectricAngularGeneratorPolarized.cc.

References FatalException, and G4Exception().

                                                                                                          :G4RDVPhotoElectricAngularDistribution(name)
{
  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 = getenv("G4LEDATA");
    if (!path)
      {
        G4String excep = "G4LEDATA environment variable not set!";
        G4Exception("G4RDPhotoElectricAngularGeneratorPolarized()",
                    "InvalidSetup", FatalException, excep);
      }

    G4String pathString(path);
    G4String dirFile = pathString + "/photoelectric_angular/" + filename;
    FILE *infile;
    infile = fopen(dirFile,"r"); 
    if (infile == 0)
      {
    G4String excep = "Data file: " + dirFile + " not found";
    G4Exception("G4RDPhotoElectricAngularGeneratorPolarized()",
                    "DataNotFound", FatalException, excep);
      }

    // Read parameters into tables. The parameters are function of incident electron energy and shell level
    G4float aRead,cRead, beta;
    for(G4int i=0 ; i<arrayDim ;i++){
      fscanf(infile,"%f\t %e\t %e",&beta,&aRead,&cRead);
      aMajorantSurfaceParameterTable[i][level] = aRead;    
      cMajorantSurfaceParameterTable[i][level] = cRead;
    }
    fclose(infile);

  }
}

G4RDPhotoElectricAngularGeneratorPolarized::~G4RDPhotoElectricAngularGeneratorPolarized

(

)

Definition at line 125 of file G4RDPhotoElectricAngularGeneratorPolarized.cc.

{;}

Member Function Documentation

G4ThreeVector G4RDPhotoElectricAngularGeneratorPolarized::GetPhotoElectronDirection ( const G4ThreeVector &  direction, const G4double  kineticEnergy, const G4ThreeVector &  polarization, const G4int  shellId  ) const

virtual

Implements G4RDVPhotoElectricAngularDistribution.

Definition at line 130 of file G4RDPhotoElectricAngularGeneratorPolarized.cc.

References python.hepunit::electron_mass_c2.

{
  // Calculate Lorentz term (gamma) and beta parameters
  G4double gamma   = 1. + eKineticEnergy/electron_mass_c2;
  G4double beta  = std::sqrt(gamma*gamma-1.)/gamma;

  G4double theta, phi = 0;
  G4double aBeta = 0; // Majorant surface parameter (function of the outgoing electron kinetic energy) 
  G4double cBeta = 0; // Majorant surface parameter (function of the outgoing electron kinetic energy)

  G4int shellLevel = 0;
  if(shellId <  2) shellLevel = 0; // K-shell  // Polarized model for K-shell
  if(shellId >= 2) shellLevel = 1; // L1-shell // Polarized model for L1 and higher shells

  // For the outgoing kinetic energy find the current majorant surface parameters
  PhotoElectronGetMajorantSurfaceAandCParameters( shellLevel, beta, &aBeta, &cBeta);

  // Generate pho and theta according to the shell level and beta parameter of the electron
  PhotoElectronGeneratePhiAndTheta(shellLevel, beta, aBeta, cBeta, &phi, &theta);

  // Determine the rotation matrix
  G4RotationMatrix rotation = PhotoElectronRotationMatrix(direction, polarization);

  // Compute final direction of the outgoing electron
  G4ThreeVector final_direction = PhotoElectronComputeFinalDirection(rotation, theta, phi);

  return final_direction;
}

void G4RDPhotoElectricAngularGeneratorPolarized::PrintGeneratorInformation ( ) const

virtual

Implements G4RDVPhotoElectricAngularDistribution.

Definition at line 397 of file G4RDPhotoElectricAngularGeneratorPolarized.cc.

References G4cout, and G4endl.

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

G4ThreeVector G4RDPhotoElectricAngularGeneratorPolarized::SetPerpendicularVector ( const G4ThreeVector &  a ) const

protected

Definition at line 407 of file G4RDPhotoElectricAngularGeneratorPolarized.cc.

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

{
  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);
  }
}

---

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

• G4RDPhotoElectricAngularGeneratorPolarized.hh
• G4RDPhotoElectricAngularGeneratorPolarized.cc