#include <G4AdjointPhotoElectricModel.hh>

Inheritance diagram for G4AdjointPhotoElectricModel:

Public Member Functions
	G4AdjointPhotoElectricModel ()

	~G4AdjointPhotoElectricModel () override

void	SampleSecondaries (const G4Track &aTrack, G4bool isScatProjToProj, G4ParticleChange *fParticleChange) override

G4double	AdjointCrossSection (const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool isScatProjToProj) override

G4double	AdjointCrossSectionPerAtom (const G4Element *anElement, G4double electronEnergy)

	G4AdjointPhotoElectricModel (G4AdjointPhotoElectricModel &)=delete

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

Public Member Functions inherited from G4VEmAdjointModel
	G4VEmAdjointModel (const G4String &nam)

virtual	~G4VEmAdjointModel ()

virtual G4double	DiffCrossSectionPerAtomPrimToSecond (G4double kinEnergyProj, G4double kinEnergyProd, G4double Z, G4double A=0.)

virtual G4double	DiffCrossSectionPerAtomPrimToScatPrim (G4double kinEnergyProj, G4double kinEnergyScatProj, G4double Z, G4double A=0.)

virtual G4double	DiffCrossSectionPerVolumePrimToSecond (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)

virtual G4double	DiffCrossSectionPerVolumePrimToScatPrim (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyScatProj)

virtual G4double	GetSecondAdjEnergyMaxForScatProjToProj (G4double primAdjEnergy)

virtual G4double	GetSecondAdjEnergyMinForScatProjToProj (G4double primAdjEnergy, G4double tcut=0.)

virtual G4double	GetSecondAdjEnergyMaxForProdToProj (G4double primAdjEnergy)

virtual G4double	GetSecondAdjEnergyMinForProdToProj (G4double primAdjEnergy)

void	DefineCurrentMaterial (const G4MaterialCutsCouple *couple)

std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerAtomForSecond (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)

std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerAtomForScatProj (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)

std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerVolumeForSecond (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)

std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerVolumeForScatProj (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)

void	SetCSMatrices (std::vector< G4AdjointCSMatrix * > Vec1CSMatrix, std::vector< G4AdjointCSMatrix > *Vec2CSMatrix)

G4ParticleDefinition *	GetAdjointEquivalentOfDirectPrimaryParticleDefinition () const

G4ParticleDefinition *	GetAdjointEquivalentOfDirectSecondaryParticleDefinition () const

G4double	GetHighEnergyLimit () const

G4double	GetLowEnergyLimit () const

void	SetHighEnergyLimit (G4double aVal)

void	SetLowEnergyLimit (G4double aVal)

void	DefineDirectEMModel (G4VEmModel *aModel)

void	SetAdjointEquivalentOfDirectPrimaryParticleDefinition (G4ParticleDefinition *aPart)

void	SetAdjointEquivalentOfDirectSecondaryParticleDefinition (G4ParticleDefinition *aPart)

void	SetSecondPartOfSameType (G4bool aBool)

G4bool	GetSecondPartOfSameType () const

void	SetUseMatrix (G4bool aBool)

void	SetUseMatrixPerElement (G4bool aBool)

void	SetUseOnlyOneMatrixForAllElements (G4bool aBool)

void	SetApplyCutInRange (G4bool aBool)

G4bool	GetUseMatrix () const

G4bool	GetUseMatrixPerElement () const

G4bool	GetUseOnlyOneMatrixForAllElements () const

G4bool	GetApplyCutInRange () const

const G4String &	GetName () const

virtual void	SetCSBiasingFactor (G4double aVal)

void	SetCorrectWeightForPostStepInModel (G4bool aBool)

void	SetAdditionalWeightCorrectionFactorForPostStepOutsideModel (G4double factor)

	G4VEmAdjointModel (G4VEmAdjointModel &)=delete

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

Protected Member Functions
void	CorrectPostStepWeight (G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool isScatProjToProj) override

Protected Member Functions inherited from G4VEmAdjointModel
G4double	DiffCrossSectionFunction1 (G4double kinEnergyProj)

G4double	DiffCrossSectionFunction2 (G4double kinEnergyProj)

G4double	SampleAdjSecEnergyFromCSMatrix (std::size_t MatrixIndex, G4double prim_energy, G4bool isScatProjToProj)

G4double	SampleAdjSecEnergyFromCSMatrix (G4double prim_energy, G4bool isScatProjToProj)

void	SelectCSMatrix (G4bool isScatProjToProj)

virtual G4double	SampleAdjSecEnergyFromDiffCrossSectionPerAtom (G4double prim_energy, G4bool isScatProjToProj)

Additional Inherited Members
Protected Attributes inherited from G4VEmAdjointModel
G4AdjointCSManager *	fCSManager

G4VEmModel *	fDirectModel = nullptr

const G4String	fName

G4Material *	fSelectedMaterial = nullptr

G4Material *	fCurrentMaterial = nullptr

G4MaterialCutsCouple *	fCurrentCouple = nullptr

G4ParticleDefinition *	fAdjEquivDirectPrimPart = nullptr

G4ParticleDefinition *	fAdjEquivDirectSecondPart = nullptr

G4ParticleDefinition *	fDirectPrimaryPart = nullptr

std::vector< G4AdjointCSMatrix * > *	fCSMatrixProdToProjBackScat = nullptr

std::vector< G4AdjointCSMatrix * > *	fCSMatrixProjToProjBackScat = nullptr

std::vector< G4double >	fElementCSScatProjToProj

std::vector< G4double >	fElementCSProdToProj

G4double	fKinEnergyProdForIntegration = 0.

G4double	fKinEnergyScatProjForIntegration = 0.

G4double	fLastCS = 0.

G4double	fLastAdjointCSForScatProjToProj = 0.

G4double	fLastAdjointCSForProdToProj = 0.

G4double	fPreStepEnergy = 0.

G4double	fTcutPrim = 0.

G4double	fTcutSecond = 0.

G4double	fHighEnergyLimit = 0.

G4double	fLowEnergyLimit = 0.

G4double	fCsBiasingFactor = 1.

G4double	fOutsideWeightFactor = 1.

G4int	fASelectedNucleus = 0

G4int	fZSelectedNucleus = 0

std::size_t	fCSMatrixUsed = 0

G4bool	fSecondPartSameType = false

G4bool	fInModelWeightCorr

G4bool	fApplyCutInRange = true

G4bool	fUseMatrix = false

G4bool	fUseMatrixPerElement = false

G4bool	fOneMatrixForAllElements = false

Detailed Description

Definition at line 53 of file G4AdjointPhotoElectricModel.hh.

Constructor & Destructor Documentation

◆ G4AdjointPhotoElectricModel() [1/2]

G4AdjointPhotoElectricModel::G4AdjointPhotoElectricModel ( )

Definition at line 39 of file G4AdjointPhotoElectricModel.cc.

  : G4VEmAdjointModel("AdjointPEEffect")
 
{
  SetUseMatrix(false);
  SetApplyCutInRange(false);
 
  fAdjEquivDirectPrimPart   = G4AdjointGamma::AdjointGamma();
  fAdjEquivDirectSecondPart = G4AdjointElectron::AdjointElectron();
  fDirectPrimaryPart        = G4Gamma::Gamma();
  fSecondPartSameType       = false;
  fDirectModel              = new G4PEEffectFluoModel();
}

Referenced by G4AdjointPhotoElectricModel(), and operator=().

◆ ~G4AdjointPhotoElectricModel()

G4AdjointPhotoElectricModel::~G4AdjointPhotoElectricModel ( )

override

Definition at line 54 of file G4AdjointPhotoElectricModel.cc.

54{}

◆ G4AdjointPhotoElectricModel() [2/2]

G4AdjointPhotoElectricModel::G4AdjointPhotoElectricModel ( G4AdjointPhotoElectricModel & )

delete

Member Function Documentation

◆ AdjointCrossSection()

G4double G4AdjointPhotoElectricModel::AdjointCrossSection	(	const G4MaterialCutsCouple *	aCouple,
		G4double	primEnergy,
		G4bool	isScatProjToProj )

overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 165 of file G4AdjointPhotoElectricModel.cc.

{
  if(isScatProjToProj)
    return 0.;
 
  G4double totBiasedAdjointCS = 0.;
  if(aCouple != fCurrentCouple || fCurrenteEnergy != electronEnergy)
  {
    fTotAdjointCS = 0.;
    DefineCurrentMaterialAndElectronEnergy(aCouple, electronEnergy);
    const G4ElementVector* theElementVector =
      fCurrentMaterial->GetElementVector();
    const G4double* theAtomNumDensityVector =
      fCurrentMaterial->GetVecNbOfAtomsPerVolume();
    size_t nelm = fCurrentMaterial->GetNumberOfElements();
    for(fIndexElement = 0; fIndexElement < nelm; ++fIndexElement)
    {
      fTotAdjointCS += AdjointCrossSectionPerAtom(
                         (*theElementVector)[fIndexElement], electronEnergy) *
                       theAtomNumDensityVector[fIndexElement];
      fXsec[fIndexElement] = fTotAdjointCS;
    }
 
    totBiasedAdjointCS = std::min(fTotAdjointCS, 0.01);
    fFactorCSBiasing   = totBiasedAdjointCS / fTotAdjointCS;
  }
  return totBiasedAdjointCS;
}

Referenced by SampleSecondaries().

◆ AdjointCrossSectionPerAtom()

G4double G4AdjointPhotoElectricModel::AdjointCrossSectionPerAtom	(	const G4Element *	anElement,
		G4double	electronEnergy )

Definition at line 197 of file G4AdjointPhotoElectricModel.cc.

{
  G4int nShells        = anElement->GetNbOfAtomicShells();
  G4double Z           = anElement->GetZ();
  G4double gammaEnergy = electronEnergy + anElement->GetAtomicShell(0);
  G4double CS          = fDirectModel->ComputeCrossSectionPerAtom(
    G4Gamma::Gamma(), gammaEnergy, Z, 0., 0., 0.);
  G4double adjointCS = 0.;
  if(CS > 0.)
    adjointCS += CS / gammaEnergy;
  fShellProb[fIndexElement][0] = adjointCS;
  for(G4int i = 1; i < nShells; ++i)
  {
    G4double Bi1 = anElement->GetAtomicShell(i - 1);
    G4double Bi  = anElement->GetAtomicShell(i);
    if(electronEnergy < Bi1 - Bi)
    {
      gammaEnergy = electronEnergy + Bi;
      CS          = fDirectModel->ComputeCrossSectionPerAtom(G4Gamma::Gamma(),
                                                    gammaEnergy, Z, 0., 0., 0.);
      if(CS > 0.)
        adjointCS += CS / gammaEnergy;
    }
    fShellProb[fIndexElement][i] = adjointCS;
  }
  adjointCS *= electronEnergy;
  return adjointCS;
}

Referenced by AdjointCrossSection().

◆ CorrectPostStepWeight()

void G4AdjointPhotoElectricModel::CorrectPostStepWeight	(	G4ParticleChange *	fParticleChange,
		G4double	old_weight,
		G4double	adjointPrimKinEnergy,
		G4double	projectileKinEnergy,
		G4bool	isScatProjToProj )

overrideprotectedvirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 147 of file G4AdjointPhotoElectricModel.cc.

{
  G4double new_weight = old_weight;
 
  G4double w_corr =
    G4AdjointCSManager::GetAdjointCSManager()->GetPostStepWeightCorrection() /
    fFactorCSBiasing;
  w_corr *= fPostStepAdjointCS / fPreStepAdjointCS;
 
  new_weight *= w_corr * projectileKinEnergy / adjointPrimKinEnergy;
  fParticleChange->SetParentWeightByProcess(false);
  fParticleChange->SetSecondaryWeightByProcess(false);
  fParticleChange->ProposeParentWeight(new_weight);
}

Referenced by SampleSecondaries().

◆ operator=()

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

delete

◆ SampleSecondaries()

void G4AdjointPhotoElectricModel::SampleSecondaries	(	const G4Track &	aTrack,
		G4bool	isScatProjToProj,
		G4ParticleChange *	fParticleChange )

overridevirtual

Implements G4VEmAdjointModel.

Definition at line 57 of file G4AdjointPhotoElectricModel.cc.

{
  if(isScatProjToProj)
    return;
 
  // Compute the fTotAdjointCS vectors if not already done for the current
  // couple and electron energy
  const G4DynamicParticle* aDynPart = aTrack.GetDynamicParticle();
  G4double electronEnergy           = aDynPart->GetKineticEnergy();
  G4ThreeVector electronDirection   = aDynPart->GetMomentumDirection();
  fPreStepAdjointCS =
    fTotAdjointCS;  // The last computed CS was at pre step point
  AdjointCrossSection(aTrack.GetMaterialCutsCouple(), electronEnergy,
                      isScatProjToProj);
  fPostStepAdjointCS = fTotAdjointCS;
 
  // Sample element
  const G4ElementVector* theElementVector =
    fCurrentMaterial->GetElementVector();
  size_t nelm      = fCurrentMaterial->GetNumberOfElements();
  G4double rand_CS = G4UniformRand() * fXsec[nelm - 1];
  for(fIndexElement = 0; fIndexElement < nelm - 1; ++fIndexElement)
  {
    if(rand_CS < fXsec[fIndexElement])
      break;
  }
 
  // Sample shell and binding energy
  G4int nShells = (*theElementVector)[fIndexElement]->GetNbOfAtomicShells();
  rand_CS       = fShellProb[fIndexElement][nShells - 1] * G4UniformRand();
  G4int i;
  for(i = 0; i < nShells - 1; ++i)
  {
    if(rand_CS < fShellProb[fIndexElement][i])
      break;
  }
  G4double gammaEnergy =
    electronEnergy + (*theElementVector)[fIndexElement]->GetAtomicShell(i);
 
  // Sample cos theta
  // Copy of the G4PEEfectFluoModel cos theta sampling method
  // ElecCosThetaDistribution. This method cannot be used directly from
  // G4PEEffectFluoModel because it is a friend method.
  G4double cos_theta = 1.;
  G4double gamma     = 1. + electronEnergy / electron_mass_c2;
  if(gamma <= 5.)
  {
    G4double beta = std::sqrt(gamma * gamma - 1.) / gamma;
    G4double b    = 0.5 * gamma * (gamma - 1.) * (gamma - 2.);
 
    G4double rndm, term, greject, grejsup;
    if(gamma < 2.)
      grejsup = gamma * gamma * (1. + b - beta * b);
    else
      grejsup = gamma * gamma * (1. + b + beta * b);
 
    do
    {
      rndm      = 1. - 2. * G4UniformRand();
      cos_theta = (rndm + beta) / (rndm * beta + 1.);
      term      = 1. - beta * cos_theta;
      greject = (1. - cos_theta * cos_theta) * (1. + b * term) / (term * term);
      // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
    } while(greject < G4UniformRand() * grejsup);
  }
 
  // direction of the adjoint gamma electron
  G4double sin_theta = std::sqrt(1. - cos_theta * cos_theta);
  G4double phi       = twopi * G4UniformRand();
  G4double dirx      = sin_theta * std::cos(phi);
  G4double diry      = sin_theta * std::sin(phi);
  G4double dirz      = cos_theta;
  G4ThreeVector adjoint_gammaDirection(dirx, diry, dirz);
  adjoint_gammaDirection.rotateUz(electronDirection);
 
  // Weight correction
  CorrectPostStepWeight(fParticleChange, aTrack.GetWeight(), electronEnergy,
                        gammaEnergy, isScatProjToProj);
 
  // Create secondary and modify fParticleChange
  G4DynamicParticle* anAdjointGamma = new G4DynamicParticle(
    G4AdjointGamma::AdjointGamma(), adjoint_gammaDirection, gammaEnergy);
 
  fParticleChange->ProposeTrackStatus(fStopAndKill);
  fParticleChange->AddSecondary(anAdjointGamma);
}

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

Public Member Functions

Protected Member Functions

Additional Inherited Members