#include <G4PenelopeBremsstrahlungAngular.hh>

Inheritance diagram for G4PenelopeBremsstrahlungAngular:

Public Member Functions
	G4PenelopeBremsstrahlungAngular ()

const G4String &	GetName () const

G4int	GetVerbosityLevel ()

void	Initialize ()

void	PrepareTables (const G4Material *material, G4bool isMaster)
	Reserved for Master Model. More...

virtual void	PrintGeneratorInformation () const

G4ThreeVector &	SampleDirection (const G4DynamicParticle dp, G4double out_energy, G4int Z, const G4Material mat=nullptr) override
	Samples the direction of the outgoing photon (in global coordinates). More...

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)

void	SetVerbosityLevel (G4int vl)
	Set/Get Verbosity level. More...

	~G4PenelopeBremsstrahlungAngular ()

Protected Attributes
G4ThreeVector	fLocalDirection

G4bool	fPolarisation

Private Member Functions
G4double	CalculateEffectiveZ (const G4Material *material)

void	ClearTables ()

void	ReadDataFile ()

Private Attributes
G4bool	fDataRead

std::map< const G4Material , G4double >	fEffectiveZSq

std::map< G4double, G4PhysicsTable * > *	fLorentzTables1

std::map< G4double, G4PhysicsTable * > *	fLorentzTables2

G4String	fName

G4double	fQQ1 [fNumberofZPoints][fNumberofEPoints][fNumberofKPoints]

G4double	fQQ2 [fNumberofZPoints][fNumberofEPoints][fNumberofKPoints]

G4int	fVerbosityLevel

Static Private Attributes
static const G4int	fNumberofEPoints =6

static const G4int	fNumberofKPoints =4

static const G4int	fNumberofZPoints =6

Detailed Description

Definition at line 55 of file G4PenelopeBremsstrahlungAngular.hh.

Constructor & Destructor Documentation

◆ G4PenelopeBremsstrahlungAngular()

G4PenelopeBremsstrahlungAngular::G4PenelopeBremsstrahlungAngular ( )

explicit

Definition at line 60 of file G4PenelopeBremsstrahlungAngular.cc.

                                                                 :
  G4VEmAngularDistribution("Penelope"), fEffectiveZSq(nullptr),
  fLorentzTables1(nullptr),fLorentzTables2(nullptr)
 
{
  fDataRead = false;
  fVerbosityLevel = 0;
}

References fDataRead, and fVerbosityLevel.

◆ ~G4PenelopeBremsstrahlungAngular()

G4PenelopeBremsstrahlungAngular::~G4PenelopeBremsstrahlungAngular ( )

Definition at line 71 of file G4PenelopeBremsstrahlungAngular.cc.

{
  ClearTables();
}

References ClearTables().

Member Function Documentation

◆ CalculateEffectiveZ()

G4double G4PenelopeBremsstrahlungAngular::CalculateEffectiveZ ( const G4Material * material )

private

Definition at line 465 of file G4PenelopeBremsstrahlungAngular.cc.

{
  if (!fEffectiveZSq)
    fEffectiveZSq = new std::map<const G4Material*,G4double>;
 
  //Already exists: return it
  if (fEffectiveZSq->count(material))
    return fEffectiveZSq->find(material)->second;
 
  //Helper for the calculation
  std::vector<G4double> *StechiometricFactors = new std::vector<G4double>;
  G4int nElements = material->GetNumberOfElements();
  const G4ElementVector* elementVector = material->GetElementVector();
  const G4double* fractionVector = material->GetFractionVector();
  for (G4int i=0;i<nElements;i++)
    {
      G4double fraction = fractionVector[i];
      G4double atomicWeigth = (*elementVector)[i]->GetA()/(g/mole);
      StechiometricFactors->push_back(fraction/atomicWeigth);
    }
  //Find max
  G4double MaxStechiometricFactor = 0.;
  for (G4int i=0;i<nElements;i++)
    {
      if ((*StechiometricFactors)[i] > MaxStechiometricFactor)
        MaxStechiometricFactor = (*StechiometricFactors)[i];
    }
  //Normalize
  for (G4int i=0;i<nElements;i++)
    (*StechiometricFactors)[i] /=  MaxStechiometricFactor;
 
  G4double sumz2 = 0;
  G4double sums = 0;
  for (G4int i=0;i<nElements;i++)
    {
      G4double Z = (*elementVector)[i]->GetZ();
      sumz2 += (*StechiometricFactors)[i]*Z*Z;
      sums  += (*StechiometricFactors)[i];
    }
  delete StechiometricFactors;
 
  G4double ZBR = std::sqrt(sumz2/sums);
  fEffectiveZSq->insert(std::make_pair(material,ZBR));
 
  return ZBR;
}

References fEffectiveZSq, g, eplot::material, mole, and Z.

Referenced by PrepareTables().

◆ ClearTables()

void G4PenelopeBremsstrahlungAngular::ClearTables ( )

private

Definition at line 85 of file G4PenelopeBremsstrahlungAngular.cc.

{
  if (fLorentzTables1)
    {
      for (auto j=fLorentzTables1->begin(); j != fLorentzTables1->end(); ++j)
        {
      G4PhysicsTable* tab = j->second;
          tab->clearAndDestroy();
          delete tab;
        }
      fLorentzTables1->clear();
      delete fLorentzTables1;
      fLorentzTables1 = nullptr;
    }
 
  if (fLorentzTables2)
    {
      for (auto j=fLorentzTables2->begin(); j != fLorentzTables2->end(); ++j)
        {
      G4PhysicsTable* tab = j->second;
      tab->clearAndDestroy();
          delete tab;
        }
      fLorentzTables2->clear();
      delete fLorentzTables2;
      fLorentzTables2 = nullptr;
    }
  if (fEffectiveZSq)
    {
      delete fEffectiveZSq;
      fEffectiveZSq = nullptr;
    }
}

References G4PhysicsTable::clearAndDestroy(), fEffectiveZSq, fLorentzTables1, and fLorentzTables2.

Referenced by Initialize(), and ~G4PenelopeBremsstrahlungAngular().

◆ GetName()

const G4String & G4VEmAngularDistribution::GetName ( ) const

inlineinherited

Definition at line 111 of file G4VEmAngularDistribution.hh.

{
  return fName;
}

References G4VEmAngularDistribution::fName.

◆ GetVerbosityLevel()

G4int G4PenelopeBremsstrahlungAngular::GetVerbosityLevel ( )

inline

Definition at line 69 of file G4PenelopeBremsstrahlungAngular.hh.

69{return fVerbosityLevel;};

References fVerbosityLevel.

◆ Initialize()

void G4PenelopeBremsstrahlungAngular::Initialize ( )

Reserved for Master Model The Initialize() method forces the cleaning of tables

Definition at line 78 of file G4PenelopeBremsstrahlungAngular.cc.

{
  ClearTables();
}

References ClearTables().

Referenced by G4PenelopeBremsstrahlungModel::Initialise(), and G4PenelopeBremsstrahlungModel::InitialiseLocal().

◆ PrepareTables()

void G4PenelopeBremsstrahlungAngular::PrepareTables	(	const G4Material *	material,
		G4bool	isMaster
	)

Reserved for Master Model.

Definition at line 174 of file G4PenelopeBremsstrahlungAngular.cc.

{
  //Unused at the moment: the G4PenelopeBremsstrahlungAngular is thread-local, so each worker
  //builds its own version of the tables.
  
  //Check if data file has already been read
  if (!fDataRead)
    {
      ReadDataFile();
      if (!fDataRead)
    G4Exception("G4PenelopeBremsstrahlungAngular::PrepareInterpolationTables()",
            "em2001",FatalException,"Unable to build interpolation table");
    }
 
  if (!fLorentzTables1)
      fLorentzTables1 = new std::map<G4double,G4PhysicsTable*>;
  if (!fLorentzTables2)
    fLorentzTables2 = new std::map<G4double,G4PhysicsTable*>;
 
  G4double Zmat = CalculateEffectiveZ(material);
 
  const G4int reducedEnergyGrid=21;
  //Support arrays.
  G4double betas[fNumberofEPoints]; //betas for interpolation
  //tables for interpolation
  G4double Q1[fNumberofEPoints][fNumberofKPoints];
  G4double Q2[fNumberofEPoints][fNumberofKPoints];
  //expanded tables for interpolation
  G4double Q1E[fNumberofEPoints][reducedEnergyGrid];
  G4double Q2E[fNumberofEPoints][reducedEnergyGrid];
  G4double pZ[fNumberofZPoints] = {2.0,8.0,13.0,47.0,79.0,92.0};
 
  G4int i=0,j=0,k=0; // i=index for Z, j=index for E, k=index for K
  //Interpolation in Z
  for (i=0;i<fNumberofEPoints;i++)
    {
      for (j=0;j<fNumberofKPoints;j++)
    {
      G4PhysicsFreeVector* QQ1vector = 
        new G4PhysicsFreeVector(fNumberofZPoints,/*spline =*/true);
      G4PhysicsFreeVector* QQ2vector = 
        new G4PhysicsFreeVector(fNumberofZPoints,/*spline =*/true);
 
      //fill vectors
      for (k=0;k<fNumberofZPoints;k++)
        {
          QQ1vector->PutValues(k,pZ[k],G4Log(fQQ1[k][i][j]));
          QQ2vector->PutValues(k,pZ[k],fQQ2[k][i][j]);
        }
      //Filled: now calculate derivatives
      QQ1vector->FillSecondDerivatives();
      QQ2vector->FillSecondDerivatives();
 
      Q1[i][j]= G4Exp(QQ1vector->Value(Zmat));
      Q2[i][j]=QQ2vector->Value(Zmat);
      delete QQ1vector;
      delete QQ2vector;
    }
    }
  G4double pE[fNumberofEPoints] = {1.0e-03*MeV,5.0e-03*MeV,1.0e-02*MeV,5.0e-02*MeV,
                  1.0e-01*MeV,5.0e-01*MeV};
  G4double pK[fNumberofKPoints] = {0.0,0.6,0.8,0.95};
  G4double ppK[reducedEnergyGrid];
 
  for(i=0;i<reducedEnergyGrid;i++)
    ppK[i]=((G4double) i) * 0.05;
 
 
  for(i=0;i<fNumberofEPoints;i++)
    betas[i]=std::sqrt(pE[i]*(pE[i]+2*electron_mass_c2))/(pE[i]+electron_mass_c2);
 
 
  for (i=0;i<fNumberofEPoints;i++)
    {
      for (j=0;j<fNumberofKPoints;j++)
    Q1[i][j]=Q1[i][j]/Zmat;
    }
 
  //Expanded table of distribution parameters
  for (i=0;i<fNumberofEPoints;i++)
    {
      G4PhysicsFreeVector* Q1vector = new G4PhysicsFreeVector(fNumberofKPoints);
      G4PhysicsFreeVector* Q2vector = new G4PhysicsFreeVector(fNumberofKPoints);
 
      for (j=0;j<fNumberofKPoints;j++)
    {
      Q1vector->PutValues(j,pK[j],G4Log(Q1[i][j])); //logarithmic
      Q2vector->PutValues(j,pK[j],Q2[i][j]);
    }
 
      for (j=0;j<reducedEnergyGrid;j++)
    {
      Q1E[i][j]=Q1vector->Value(ppK[j]);
      Q2E[i][j]=Q2vector->Value(ppK[j]);
    }
      delete Q1vector;
      delete Q2vector;
    }
  //
  //TABLES to be stored
  //
  G4PhysicsTable* theTable1 = new G4PhysicsTable();
  G4PhysicsTable* theTable2 = new G4PhysicsTable();
  // the table will contain reducedEnergyGrid G4PhysicsFreeVectors with different
  // values of k,
  // Each of the G4PhysicsFreeVectors has a profile of
  // y vs. E
  //
  //reserve space of the vectors.
  for (j=0;j<reducedEnergyGrid;j++)
    {
      G4PhysicsFreeVector* thevec = new G4PhysicsFreeVector(fNumberofEPoints,/*spline=*/true);
      theTable1->push_back(thevec);
      G4PhysicsFreeVector* thevec2 = new G4PhysicsFreeVector(fNumberofEPoints,/*spline=*/true);
      theTable2->push_back(thevec2);
    }
 
  for (j=0;j<reducedEnergyGrid;j++)
    {
      G4PhysicsFreeVector* thevec = (G4PhysicsFreeVector*) (*theTable1)[j];
      G4PhysicsFreeVector* thevec2 = (G4PhysicsFreeVector*) (*theTable2)[j];
      for (i=0;i<fNumberofEPoints;i++)
    {
      thevec->PutValues(i,betas[i],Q1E[i][j]);
      thevec2->PutValues(i,betas[i],Q2E[i][j]);
    }
      //Vectors are filled: calculate derivatives
      thevec->FillSecondDerivatives();
      thevec2->FillSecondDerivatives();
    }
 
  if (fLorentzTables1 && fLorentzTables2)
    {
      fLorentzTables1->insert(std::make_pair(Zmat,theTable1));
      fLorentzTables2->insert(std::make_pair(Zmat,theTable2));
    }
  else
    {
      G4ExceptionDescription ed;
      ed << "Unable to create tables of Lorentz coefficients for " << G4endl;
      ed << "<Z>= "  << Zmat << " in G4PenelopeBremsstrahlungAngular" << G4endl;
      delete theTable1;
      delete theTable2;
      G4Exception("G4PenelopeBremsstrahlungAngular::PrepareInterpolationTables()",
          "em2005",FatalException,ed);
    }
 
  return;
}

References CalculateEffectiveZ(), source.hepunit::electron_mass_c2, FatalException, fDataRead, G4PhysicsVector::FillSecondDerivatives(), fLorentzTables1, fLorentzTables2, fNumberofEPoints, fNumberofKPoints, fNumberofZPoints, fQQ1, fQQ2, G4endl, G4Exception(), G4Exp(), G4Log(), eplot::material, MeV, G4PhysicsTable::push_back(), G4PhysicsFreeVector::PutValues(), ReadDataFile(), and G4PhysicsVector::Value().

Referenced by G4PenelopeBremsstrahlungModel::Initialise(), and G4PenelopeBremsstrahlungModel::InitialiseLocal().

◆ PrintGeneratorInformation()

void G4VEmAngularDistribution::PrintGeneratorInformation ( ) const

virtualinherited

Reimplemented in G4DeltaAngleFreeScat, G4DipBustGenerator, G4ModifiedTsai, G4DNABornAngle, G4DNARuddAngle, G4Generator2BN, G4Generator2BS, G4PhotoElectricAngularGeneratorPolarized, G4PhotoElectricAngularGeneratorSauterGavrila, G4ModifiedMephi, and G4SauterGavrilaAngularDistribution.

Definition at line 92 of file G4VEmAngularDistribution.cc.

93{}

◆ ReadDataFile()

void G4PenelopeBremsstrahlungAngular::ReadDataFile ( void )

private

Definition at line 121 of file G4PenelopeBremsstrahlungAngular.cc.

{
   //Read information from DataBase file
  char* path = std::getenv("G4LEDATA");
  if (!path)
    {
      G4String excep =
    "G4PenelopeBremsstrahlungAngular - G4LEDATA environment variable not set!";
      G4Exception("G4PenelopeBremsstrahlungAngular::ReadDataFile()",
          "em0006",FatalException,excep);
      return;
    }
  G4String pathString(path);
  G4String pathFile = pathString + "/penelope/bremsstrahlung/pdbrang.p08";
  std::ifstream file(pathFile);
 
  if (!file.is_open())
    {
      G4String excep = "G4PenelopeBremsstrahlungAngular - data file " + pathFile + " not found!";
      G4Exception("G4PenelopeBremsstrahlungAngular::ReadDataFile()",
          "em0003",FatalException,excep);
      return;
    }
  G4int i=0,j=0,k=0; // i=index for Z, j=index for E, k=index for K
 
  for (k=0;k<fNumberofKPoints;k++)
    for (i=0;i<fNumberofZPoints;i++)
      for (j=0;j<fNumberofEPoints;j++)
    {
      G4double a1,a2;
      G4int ik1,iz1,ie1;
      G4double zr,er,kr;
      file >> iz1 >> ie1 >> ik1 >> zr >> er >> kr >> a1 >> a2;
      //check the data are correct
      if ((iz1-1 == i) && (ik1-1 == k) && (ie1-1 == j))
        {
          fQQ1[i][j][k]=a1;
          fQQ2[i][j][k]=a2;
        }
      else
        {
          G4ExceptionDescription ed;
          ed << "Corrupted data file " << pathFile << "?" << G4endl;
          G4Exception("G4PenelopeBremsstrahlungAngular::ReadDataFile()",
          "em0005",FatalException,ed);
        }
    }
  file.close();
  fDataRead = true;
}

References FatalException, fDataRead, geant4_check_module_cycles::file, fNumberofEPoints, fNumberofKPoints, fNumberofZPoints, fQQ1, fQQ2, G4endl, and G4Exception().

Referenced by PrepareTables().

◆ SampleDirection()

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

overridevirtual

Samples the direction of the outgoing photon (in global coordinates).

Implements G4VEmAngularDistribution.

Definition at line 326 of file G4PenelopeBremsstrahlungAngular.cc.

{
  if (!material)
    {
      G4Exception("G4PenelopeBremsstrahlungAngular::SampleDirection()",
          "em2040",FatalException,"The pointer to G4Material* is nullptr");
      return fLocalDirection;
    }
 
  //Retrieve the effective Z
  G4double Zmat = 0;
 
  if (!fEffectiveZSq)
    {
      G4Exception("G4PenelopeBremsstrahlungAngular::SampleDirection()",
          "em2040",FatalException,"EffectiveZ table not available");
      return fLocalDirection;
    }
 
  //found in the table: return it
  if (fEffectiveZSq->count(material))
    Zmat = fEffectiveZSq->find(material)->second;
  else
    {
      G4Exception("G4PenelopeBremsstrahlungAngular::SampleDirection()",
          "em2040",FatalException,"Material not found in the effectiveZ table");
      return fLocalDirection;
    }
 
  if (fVerbosityLevel > 0)
    {
      G4cout << "Effective <Z> for material : " << material->GetName() <<
    " = " << Zmat << G4endl;
    }
 
  G4double ePrimary = dp->GetKineticEnergy();
 
  G4double beta = std::sqrt(ePrimary*(ePrimary+2*electron_mass_c2))/
    (ePrimary+electron_mass_c2);
  G4double cdt = 0;
  G4double sinTheta = 0;
  G4double phi = 0;
 
  //Use a pure dipole distribution for energy above 500 keV
  if (ePrimary > 500*keV)
    {
      cdt = 2.0*G4UniformRand() - 1.0;
      if (G4UniformRand() > 0.75)
    {
      if (cdt<0)
        cdt = -1.0*std::pow(-cdt,1./3.);
      else
        cdt = std::pow(cdt,1./3.);
    }
      cdt = (cdt+beta)/(1.0+beta*cdt);
      //Get primary kinematics
      sinTheta = std::sqrt(1. - cdt*cdt);
      phi  = twopi * G4UniformRand();
      fLocalDirection.set(sinTheta* std::cos(phi),
              sinTheta* std::sin(phi),
              cdt);
      //rotate
      fLocalDirection.rotateUz(dp->GetMomentumDirection());
      //return
      return fLocalDirection;
    }
 
  if (!(fLorentzTables1->count(Zmat)) || !(fLorentzTables2->count(Zmat)))
    {
      G4ExceptionDescription ed;
      ed << "Unable to retrieve Lorentz tables for Z= " << Zmat << G4endl;
      G4Exception("G4PenelopeBremsstrahlungAngular::SampleDirection()",
          "em2006",FatalException,ed);
    }
 
  //retrieve actual tables
  const G4PhysicsTable* theTable1 = fLorentzTables1->find(Zmat)->second;
  const G4PhysicsTable* theTable2 = fLorentzTables2->find(Zmat)->second;
 
  G4double RK=20.0*eGamma/ePrimary;
  G4int ik=std::min((G4int) RK,19);
 
  G4double P10=0,P11=0,P1=0;
  G4double P20=0,P21=0,P2=0;
 
  //First coefficient
  const G4PhysicsFreeVector* v1 = (G4PhysicsFreeVector*) (*theTable1)[ik];
  const G4PhysicsFreeVector* v2 = (G4PhysicsFreeVector*) (*theTable1)[ik+1];
  P10 = v1->Value(beta);
  P11 = v2->Value(beta);
  P1=P10+(RK-(G4double) ik)*(P11-P10);
 
  //Second coefficient
  const G4PhysicsFreeVector* v3 = (G4PhysicsFreeVector*) (*theTable2)[ik];
  const G4PhysicsFreeVector* v4 = (G4PhysicsFreeVector*) (*theTable2)[ik+1];
  P20=v3->Value(beta);
  P21=v4->Value(beta);
  P2=P20+(RK-(G4double) ik)*(P21-P20);
 
  //Sampling from the Lorenz-trasformed dipole distributions
  P1=std::min(G4Exp(P1)/beta,1.0);
  G4double betap = std::min(std::max(beta*(1.0+P2/beta),0.0),0.9999);
 
  G4double testf=0;
 
  if (G4UniformRand() < P1)
    {
      do{
    cdt = 2.0*G4UniformRand()-1.0;
    testf=2.0*G4UniformRand()-(1.0+cdt*cdt);
      }while(testf>0);
    }
  else
    {
      do{
    cdt = 2.0*G4UniformRand()-1.0;
    testf=G4UniformRand()-(1.0-cdt*cdt);
      }while(testf>0);
    }
  cdt = (cdt+betap)/(1.0+betap*cdt);
 
  //Get primary kinematics
  sinTheta = std::sqrt(1. - cdt*cdt);
  phi  = twopi * G4UniformRand();
  fLocalDirection.set(sinTheta* std::cos(phi),
              sinTheta* std::sin(phi),
              cdt);
  //rotate
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
  //return
  return fLocalDirection;
 
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, source.hepunit::electron_mass_c2, FatalException, fEffectiveZSq, G4VEmAngularDistribution::fLocalDirection, fLorentzTables1, fLorentzTables2, fVerbosityLevel, G4cout, G4endl, G4Exception(), G4Exp(), G4UniformRand, G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), keV, eplot::material, G4INCL::Math::max(), G4INCL::Math::min(), P1, P10, P11, P2, P20, P21, CLHEP::Hep3Vector::rotateUz(), CLHEP::Hep3Vector::set(), twopi, and G4PhysicsVector::Value().

Referenced by G4PenelopeBremsstrahlungModel::SampleSecondaries().

◆ 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().

◆ SetVerbosityLevel()

void G4PenelopeBremsstrahlungAngular::SetVerbosityLevel ( G4int vl )

inline

Set/Get Verbosity level.

Definition at line 68 of file G4PenelopeBremsstrahlungAngular.hh.

68{fVerbosityLevel = vl;};

References fVerbosityLevel.

Field Documentation

◆ fDataRead

G4bool G4PenelopeBremsstrahlungAngular::fDataRead

private

Definition at line 98 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by G4PenelopeBremsstrahlungAngular(), PrepareTables(), and ReadDataFile().

◆ fEffectiveZSq

std::map<const G4Material*,G4double>* G4PenelopeBremsstrahlungAngular::fEffectiveZSq

private

Definition at line 83 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by CalculateEffectiveZ(), ClearTables(), and SampleDirection().

◆ fLocalDirection

G4ThreeVector G4VEmAngularDistribution::fLocalDirection

protectedinherited

Definition at line 103 of file G4VEmAngularDistribution.hh.

◆ fLorentzTables1

std::map<G4double,G4PhysicsTable*>* G4PenelopeBremsstrahlungAngular::fLorentzTables1

private

Definition at line 87 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by ClearTables(), PrepareTables(), and SampleDirection().

◆ fLorentzTables2

std::map<G4double,G4PhysicsTable*>* G4PenelopeBremsstrahlungAngular::fLorentzTables2

private

Definition at line 88 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by ClearTables(), PrepareTables(), and SampleDirection().

◆ fName

G4String G4VEmAngularDistribution::fName

privateinherited

Definition at line 108 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::GetName().

◆ fNumberofEPoints

const G4int G4PenelopeBremsstrahlungAngular::fNumberofEPoints =6

staticprivate

Definition at line 91 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by PrepareTables(), and ReadDataFile().

◆ fNumberofKPoints

const G4int G4PenelopeBremsstrahlungAngular::fNumberofKPoints =4

staticprivate

Definition at line 92 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by PrepareTables(), and ReadDataFile().

◆ fNumberofZPoints

const G4int G4PenelopeBremsstrahlungAngular::fNumberofZPoints =6

staticprivate

Definition at line 90 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by PrepareTables(), and ReadDataFile().

◆ fPolarisation

G4bool G4VEmAngularDistribution::fPolarisation

protectedinherited

Definition at line 104 of file G4VEmAngularDistribution.hh.

Referenced by G4VEmAngularDistribution::G4VEmAngularDistribution().

◆ fQQ1

G4double G4PenelopeBremsstrahlungAngular::fQQ1[fNumberofZPoints][fNumberofEPoints][fNumberofKPoints]

private

Definition at line 94 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by PrepareTables(), and ReadDataFile().

◆ fQQ2

G4double G4PenelopeBremsstrahlungAngular::fQQ2[fNumberofZPoints][fNumberofEPoints][fNumberofKPoints]

private

Definition at line 95 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by PrepareTables(), and ReadDataFile().

◆ fVerbosityLevel

G4int G4PenelopeBremsstrahlungAngular::fVerbosityLevel

private

Definition at line 97 of file G4PenelopeBremsstrahlungAngular.hh.

Referenced by G4PenelopeBremsstrahlungAngular(), GetVerbosityLevel(), SampleDirection(), and SetVerbosityLevel().

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

source/processes/electromagnetic/lowenergy/include/G4PenelopeBremsstrahlungAngular.hh
source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungAngular.cc

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4PenelopeBremsstrahlungAngular()

◆ ~G4PenelopeBremsstrahlungAngular()

Member Function Documentation

◆ CalculateEffectiveZ()

◆ ClearTables()

◆ GetName()

◆ GetVerbosityLevel()

◆ Initialize()

◆ PrepareTables()

◆ PrintGeneratorInformation()

◆ ReadDataFile()

◆ SampleDirection()

◆ SampleDirectionForShell()

◆ SamplePairDirections()

◆ SetVerbosityLevel()

Field Documentation

◆ fDataRead

◆ fEffectiveZSq

◆ fLocalDirection

◆ fLorentzTables1

◆ fLorentzTables2

◆ fName

◆ fNumberofEPoints

◆ fNumberofKPoints

◆ fNumberofZPoints

◆ fPolarisation

◆ fQQ1

◆ fQQ2

◆ fVerbosityLevel