G4AdjointPhotoElectricModel Class Reference

#include <G4AdjointPhotoElectricModel.hh>

Inheritance diagram for G4AdjointPhotoElectricModel:

Public Member Functions
G4double	AdjointCrossSection (const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool isScatProjToProj) override
G4double	AdjointCrossSectionPerAtom (const G4Element *anElement, G4double electronEnergy)
std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerAtomForScatProj (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)
std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerAtomForSecond (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)
std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerVolumeForScatProj (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)
std::vector< std::vector< double > * >	ComputeAdjointCrossSectionVectorPerVolumeForSecond (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)
void	DefineCurrentMaterial (const G4MaterialCutsCouple *couple)
void	DefineDirectEMModel (G4VEmModel *aModel)
virtual G4double	DiffCrossSectionPerAtomPrimToScatPrim (G4double kinEnergyProj, G4double kinEnergyScatProj, G4double Z, G4double A=0.)
virtual G4double	DiffCrossSectionPerAtomPrimToSecond (G4double kinEnergyProj, G4double kinEnergyProd, G4double Z, G4double A=0.)
virtual G4double	DiffCrossSectionPerVolumePrimToScatPrim (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyScatProj)
virtual G4double	DiffCrossSectionPerVolumePrimToSecond (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)
	G4AdjointPhotoElectricModel ()
	G4AdjointPhotoElectricModel (G4AdjointPhotoElectricModel &)=delete
G4ParticleDefinition *	GetAdjointEquivalentOfDirectPrimaryParticleDefinition ()
G4ParticleDefinition *	GetAdjointEquivalentOfDirectSecondaryParticleDefinition ()
G4bool	GetApplyCutInRange ()
G4double	GetHighEnergyLimit ()
G4double	GetLowEnergyLimit ()
G4String	GetName ()
virtual G4double	GetSecondAdjEnergyMaxForProdToProj (G4double primAdjEnergy)
virtual G4double	GetSecondAdjEnergyMaxForScatProjToProj (G4double primAdjEnergy)
virtual G4double	GetSecondAdjEnergyMinForProdToProj (G4double primAdjEnergy)
virtual G4double	GetSecondAdjEnergyMinForScatProjToProj (G4double primAdjEnergy, G4double tcut=0.)
G4bool	GetSecondPartOfSameType ()
G4bool	GetUseMatrix ()
G4bool	GetUseMatrixPerElement ()
G4bool	GetUseOnlyOneMatrixForAllElements ()
G4AdjointPhotoElectricModel &	operator= (const G4AdjointPhotoElectricModel &right)=delete
void	SampleSecondaries (const G4Track &aTrack, G4bool isScatProjToProj, G4ParticleChange *fParticleChange) override
void	SetAdditionalWeightCorrectionFactorForPostStepOutsideModel (G4double factor)
void	SetAdjointEquivalentOfDirectPrimaryParticleDefinition (G4ParticleDefinition *aPart)
void	SetAdjointEquivalentOfDirectSecondaryParticleDefinition (G4ParticleDefinition *aPart)
void	SetApplyCutInRange (G4bool aBool)
void	SetCorrectWeightForPostStepInModel (G4bool aBool)
virtual void	SetCSBiasingFactor (G4double aVal)
void	SetCSMatrices (std::vector< G4AdjointCSMatrix * > *Vec1CSMatrix, std::vector< G4AdjointCSMatrix * > *Vec2CSMatrix)
void	SetHighEnergyLimit (G4double aVal)
void	SetLowEnergyLimit (G4double aVal)
void	SetSecondPartOfSameType (G4bool aBool)
void	SetUseMatrix (G4bool aBool)
void	SetUseMatrixPerElement (G4bool aBool)
void	SetUseOnlyOneMatrixForAllElements (G4bool aBool)
	~G4AdjointPhotoElectricModel () override

Protected Member Functions
void	CorrectPostStepWeight (G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool isScatProjToProj) override
G4double	DiffCrossSectionFunction1 (G4double kinEnergyProj)
G4double	DiffCrossSectionFunction2 (G4double kinEnergyProj)
G4double	SampleAdjSecEnergyFromCSMatrix (G4double prim_energy, G4bool isScatProjToProj)
G4double	SampleAdjSecEnergyFromCSMatrix (size_t MatrixIndex, G4double prim_energy, G4bool isScatProjToProj)
virtual G4double	SampleAdjSecEnergyFromDiffCrossSectionPerAtom (G4double prim_energy, G4bool isScatProjToProj)
void	SelectCSMatrix (G4bool isScatProjToProj)

Protected Attributes
G4ParticleDefinition *	fAdjEquivDirectPrimPart = nullptr
G4ParticleDefinition *	fAdjEquivDirectSecondPart = nullptr
G4bool	fApplyCutInRange = true
G4int	fASelectedNucleus = 0
G4double	fCsBiasingFactor = 1.
G4AdjointCSManager *	fCSManager
std::vector< G4AdjointCSMatrix * > *	fCSMatrixProdToProjBackScat = nullptr
std::vector< G4AdjointCSMatrix * > *	fCSMatrixProjToProjBackScat = nullptr
size_t	fCSMatrixUsed = 0
G4MaterialCutsCouple *	fCurrentCouple = nullptr
G4Material *	fCurrentMaterial = nullptr
G4VEmModel *	fDirectModel = nullptr
G4ParticleDefinition *	fDirectPrimaryPart = nullptr
std::vector< G4double >	fElementCSProdToProj
std::vector< G4double >	fElementCSScatProjToProj
G4double	fHighEnergyLimit = 0.
G4bool	fInModelWeightCorr
G4double	fKinEnergyProdForIntegration = 0.
G4double	fKinEnergyScatProjForIntegration = 0.
G4double	fLastAdjointCSForProdToProj = 0.
G4double	fLastAdjointCSForScatProjToProj = 0.
G4double	fLastCS = 0.
G4double	fLowEnergyLimit = 0.
const G4String	fName
G4bool	fOneMatrixForAllElements = false
G4double	fOutsideWeightFactor = 1.
G4double	fPreStepEnergy = 0.
G4bool	fSecondPartSameType = false
G4Material *	fSelectedMaterial = nullptr
G4double	fTcutPrim = 0.
G4double	fTcutSecond = 0.
G4bool	fUseMatrix = false
G4bool	fUseMatrixPerElement = false
G4int	fZSelectedNucleus = 0

Private Member Functions
void	DefineCurrentMaterialAndElectronEnergy (const G4MaterialCutsCouple *aCouple, G4double eEnergy)

Private Attributes
G4double	fCurrenteEnergy = 0.
G4double	fFactorCSBiasing = 1.
size_t	fIndexElement = 0
G4double	fPostStepAdjointCS = 0.
G4double	fPreStepAdjointCS = 0.
G4double	fShellProb [40][40]
G4double	fTotAdjointCS = 0.
G4double	fXsec [40]

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

References G4AdjointElectron::AdjointElectron(), G4AdjointGamma::AdjointGamma(), G4VEmAdjointModel::fAdjEquivDirectPrimPart, G4VEmAdjointModel::fAdjEquivDirectSecondPart, G4VEmAdjointModel::fDirectModel, G4VEmAdjointModel::fDirectPrimaryPart, G4VEmAdjointModel::fSecondPartSameType, G4Gamma::Gamma(), G4VEmAdjointModel::SetApplyCutInRange(), and G4VEmAdjointModel::SetUseMatrix().

~G4AdjointPhotoElectricModel()

G4AdjointPhotoElectricModel::~G4AdjointPhotoElectricModel ( )

override

Definition at line 54 of file G4AdjointPhotoElectricModel.cc.

{}

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

References AdjointCrossSectionPerAtom(), DefineCurrentMaterialAndElectronEnergy(), G4VEmAdjointModel::fCurrentCouple, fCurrenteEnergy, G4VEmAdjointModel::fCurrentMaterial, fFactorCSBiasing, fIndexElement, fTotAdjointCS, fXsec, G4Material::GetElementVector(), G4Material::GetNumberOfElements(), G4Material::GetVecNbOfAtomsPerVolume(), and G4INCL::Math::min().

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

References G4VEmModel::ComputeCrossSectionPerAtom(), G4VEmAdjointModel::fDirectModel, fIndexElement, fShellProb, G4Gamma::Gamma(), G4Element::GetAtomicShell(), G4Element::GetNbOfAtomicShells(), G4Element::GetZ(), and Z.

Referenced by AdjointCrossSection().

ComputeAdjointCrossSectionVectorPerAtomForScatProj()

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

inherited

Definition at line 252 of file G4VEmAdjointModel.cc.

{
  G4Integrator<G4VEmAdjointModel, G4double (G4VEmAdjointModel::*)(G4double)>
    integral;
  fASelectedNucleus                = G4lrint(A);
  fZSelectedNucleus                = G4lrint(Z);
  fKinEnergyScatProjForIntegration = kinEnergyScatProj;
 
  // compute the vector of integrated cross sections
  G4double minEProj = GetSecondAdjEnergyMinForScatProjToProj(kinEnergyScatProj);
  G4double maxEProj = GetSecondAdjEnergyMaxForScatProjToProj(kinEnergyScatProj);
  G4double dEmax    = maxEProj - kinEnergyScatProj;
  G4double dEmin    = GetLowEnergyLimit();
  G4double dE1      = dEmin;
  G4double dE2      = dEmin;
 
  std::vector<G4double>* log_ESec_vector = new std::vector<G4double>();
  std::vector<G4double>* log_Prob_vector = new std::vector<G4double>();
  log_ESec_vector->push_back(std::log(dEmin));
  log_Prob_vector->push_back(-50.);
  G4int nbins = std::max(G4int(std::log10(dEmax / dEmin)) * nbin_pro_decade, 5);
  G4double fE = std::pow(dEmax / dEmin, 1. / nbins);
 
  G4double int_cross_section = 0.;
  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
  while(dE1 < dEmax * 0.9999999999999)
  {
    dE2 = dE1 * fE;
    int_cross_section +=
      integral.Simpson(this, &G4VEmAdjointModel::DiffCrossSectionFunction2,
                       minEProj + dE1, std::min(minEProj + dE2, maxEProj), 5);
    log_ESec_vector->push_back(std::log(std::min(dE2, maxEProj - minEProj)));
    log_Prob_vector->push_back(std::log(int_cross_section));
    dE1 = dE2;
  }
 
  std::vector<std::vector<G4double>*> res_mat;
  if(int_cross_section > 0.)
  {
    res_mat.push_back(log_ESec_vector);
    res_mat.push_back(log_Prob_vector);
  }
  else {
    delete  log_ESec_vector;
    delete  log_Prob_vector;
  }
 
  return res_mat;
}

References A, G4VEmAdjointModel::DiffCrossSectionFunction2(), G4VEmAdjointModel::fASelectedNucleus, G4VEmAdjointModel::fKinEnergyScatProjForIntegration, G4VEmAdjointModel::fZSelectedNucleus, G4lrint(), G4VEmAdjointModel::G4VEmAdjointModel(), G4VEmAdjointModel::GetLowEnergyLimit(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForScatProjToProj(), G4INCL::Math::max(), G4INCL::Math::min(), and Z.

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndElement().

ComputeAdjointCrossSectionVectorPerAtomForSecond()

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

inherited

Definition at line 198 of file G4VEmAdjointModel.cc.

{
  G4Integrator<G4VEmAdjointModel, G4double (G4VEmAdjointModel::*)(G4double)>
    integral;
  fASelectedNucleus            = G4lrint(A);
  fZSelectedNucleus            = G4lrint(Z);
  fKinEnergyProdForIntegration = kinEnergyProd;
 
  // compute the vector of integrated cross sections
  G4double minEProj = GetSecondAdjEnergyMinForProdToProj(kinEnergyProd);
  G4double maxEProj = GetSecondAdjEnergyMaxForProdToProj(kinEnergyProd);
  G4double E1       = minEProj;
  std::vector<G4double>* log_ESec_vector = new std::vector<G4double>();
  std::vector<G4double>* log_Prob_vector = new std::vector<G4double>();
  log_ESec_vector->push_back(std::log(E1));
  log_Prob_vector->push_back(-50.);
 
  G4double E2 =
    std::pow(10., G4double(G4int(std::log10(minEProj) * nbin_pro_decade) + 1) /
                    nbin_pro_decade);
  G4double fE = std::pow(10., 1. / nbin_pro_decade);
 
  if(std::pow(fE, 5.) > (maxEProj / minEProj))
    fE = std::pow(maxEProj / minEProj, 0.2);
 
  G4double int_cross_section = 0.;
  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
  while(E1 < maxEProj * 0.9999999)
  {
    int_cross_section +=
      integral.Simpson(this, &G4VEmAdjointModel::DiffCrossSectionFunction1, E1,
                       std::min(E2, maxEProj * 0.99999999), 5);
    log_ESec_vector->push_back(std::log(std::min(E2, maxEProj)));
    log_Prob_vector->push_back(std::log(int_cross_section));
    E1 = E2;
    E2 *= fE;
  }
  std::vector<std::vector<G4double>*> res_mat;
  if(int_cross_section > 0.)
  {
    res_mat.push_back(log_ESec_vector);
    res_mat.push_back(log_Prob_vector);
  }
  else {
    delete  log_ESec_vector;
    delete  log_Prob_vector;
  }
  return res_mat;
}

References A, G4VEmAdjointModel::DiffCrossSectionFunction1(), G4VEmAdjointModel::fASelectedNucleus, G4VEmAdjointModel::fKinEnergyProdForIntegration, G4VEmAdjointModel::fZSelectedNucleus, G4lrint(), G4VEmAdjointModel::G4VEmAdjointModel(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj(), G4INCL::Math::min(), and Z.

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndElement().

ComputeAdjointCrossSectionVectorPerVolumeForScatProj()

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

inherited

Definition at line 360 of file G4VEmAdjointModel.cc.

{
  G4Integrator<G4VEmAdjointModel, G4double (G4VEmAdjointModel::*)(G4double)>
    integral;
  fSelectedMaterial                = aMaterial;
  fKinEnergyScatProjForIntegration = kinEnergyScatProj;
 
  // compute the vector of integrated cross sections
  G4double minEProj = GetSecondAdjEnergyMinForScatProjToProj(kinEnergyScatProj);
  G4double maxEProj = GetSecondAdjEnergyMaxForScatProjToProj(kinEnergyScatProj);
 
  G4double dEmax = maxEProj - kinEnergyScatProj;
  G4double dEmin = GetLowEnergyLimit();
  G4double dE1   = dEmin;
  G4double dE2   = dEmin;
 
  std::vector<G4double>* log_ESec_vector = new std::vector<G4double>();
  std::vector<G4double>* log_Prob_vector = new std::vector<G4double>();
  log_ESec_vector->push_back(std::log(dEmin));
  log_Prob_vector->push_back(-50.);
  G4int nbins = std::max(int(std::log10(dEmax / dEmin)) * nbin_pro_decade, 5);
  G4double fE = std::pow(dEmax / dEmin, 1. / nbins);
 
  G4double int_cross_section = 0.;
  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
  while(dE1 < dEmax * 0.9999999999999)
  {
    dE2 = dE1 * fE;
    int_cross_section +=
      integral.Simpson(this, &G4VEmAdjointModel::DiffCrossSectionFunction2,
                       minEProj + dE1, std::min(minEProj + dE2, maxEProj), 5);
    log_ESec_vector->push_back(std::log(std::min(dE2, maxEProj - minEProj)));
    log_Prob_vector->push_back(std::log(int_cross_section));
    dE1 = dE2;
  }
 
  std::vector<std::vector<G4double>*> res_mat;
  if(int_cross_section > 0.)
  {
    res_mat.push_back(log_ESec_vector);
    res_mat.push_back(log_Prob_vector);
  }
  else {
    delete  log_ESec_vector;
    delete  log_Prob_vector;
  }
 
  return res_mat;
}

References G4VEmAdjointModel::DiffCrossSectionFunction2(), G4VEmAdjointModel::fKinEnergyScatProjForIntegration, G4VEmAdjointModel::fSelectedMaterial, G4VEmAdjointModel::G4VEmAdjointModel(), G4VEmAdjointModel::GetLowEnergyLimit(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForScatProjToProj(), G4INCL::Math::max(), and G4INCL::Math::min().

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndMaterial().

ComputeAdjointCrossSectionVectorPerVolumeForSecond()

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

inherited

Definition at line 306 of file G4VEmAdjointModel.cc.

{
  G4Integrator<G4VEmAdjointModel, G4double (G4VEmAdjointModel::*)(G4double)>
    integral;
  fSelectedMaterial            = aMaterial;
  fKinEnergyProdForIntegration = kinEnergyProd;
 
  // compute the vector of integrated cross sections
  G4double minEProj = GetSecondAdjEnergyMinForProdToProj(kinEnergyProd);
  G4double maxEProj = GetSecondAdjEnergyMaxForProdToProj(kinEnergyProd);
  G4double E1       = minEProj;
  std::vector<G4double>* log_ESec_vector = new std::vector<G4double>();
  std::vector<G4double>* log_Prob_vector = new std::vector<G4double>();
  log_ESec_vector->push_back(std::log(E1));
  log_Prob_vector->push_back(-50.);
 
  G4double E2 =
    std::pow(10., G4double(G4int(std::log10(minEProj) * nbin_pro_decade) + 1) /
                    nbin_pro_decade);
  G4double fE = std::pow(10., 1. / nbin_pro_decade);
 
  if(std::pow(fE, 5.) > (maxEProj / minEProj))
    fE = std::pow(maxEProj / minEProj, 0.2);
 
  G4double int_cross_section = 0.;
  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
  while(E1 < maxEProj * 0.9999999)
  {
    int_cross_section +=
      integral.Simpson(this, &G4VEmAdjointModel::DiffCrossSectionFunction1, E1,
                       std::min(E2, maxEProj * 0.99999999), 5);
    log_ESec_vector->push_back(std::log(std::min(E2, maxEProj)));
    log_Prob_vector->push_back(std::log(int_cross_section));
    E1 = E2;
    E2 *= fE;
  }
  std::vector<std::vector<G4double>*> res_mat;
 
  if(int_cross_section > 0.)
  {
    res_mat.push_back(log_ESec_vector);
    res_mat.push_back(log_Prob_vector);
  }
  else {
    delete  log_ESec_vector;
    delete  log_Prob_vector;
  }
  return res_mat;
}

References G4VEmAdjointModel::DiffCrossSectionFunction1(), G4VEmAdjointModel::fKinEnergyProdForIntegration, G4VEmAdjointModel::fSelectedMaterial, G4VEmAdjointModel::G4VEmAdjointModel(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj(), and G4INCL::Math::min().

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndMaterial().

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

References fFactorCSBiasing, fPostStepAdjointCS, fPreStepAdjointCS, G4AdjointCSManager::GetAdjointCSManager(), G4AdjointCSManager::GetPostStepWeightCorrection(), G4VParticleChange::ProposeParentWeight(), G4VParticleChange::SetParentWeightByProcess(), and G4VParticleChange::SetSecondaryWeightByProcess().

Referenced by SampleSecondaries().

DefineCurrentMaterial()

void G4VEmAdjointModel::DefineCurrentMaterial ( const G4MaterialCutsCouple *  couple )

inherited

Definition at line 684 of file G4VEmAdjointModel.cc.

{
  if(couple != fCurrentCouple)
  {
    fCurrentCouple   = const_cast<G4MaterialCutsCouple*>(couple);
    fCurrentMaterial = const_cast<G4Material*>(couple->GetMaterial());
    size_t idx       = 56;
    fTcutSecond      = 1.e-11;
    if(fAdjEquivDirectSecondPart)
    {
      if(fAdjEquivDirectSecondPart == G4AdjointGamma::AdjointGamma())
        idx = 0;
      else if(fAdjEquivDirectSecondPart == G4AdjointElectron::AdjointElectron())
        idx = 1;
      else if(fAdjEquivDirectSecondPart == G4AdjointPositron::AdjointPositron())
        idx = 2;
      if(idx < 56)
      {
        const std::vector<G4double>* aVec =
          G4ProductionCutsTable::GetProductionCutsTable()->GetEnergyCutsVector(
            idx);
        fTcutSecond = (*aVec)[couple->GetIndex()];
      }
    }
  }
}

References G4AdjointElectron::AdjointElectron(), G4AdjointGamma::AdjointGamma(), G4AdjointPositron::AdjointPositron(), G4VEmAdjointModel::fAdjEquivDirectSecondPart, G4VEmAdjointModel::fCurrentCouple, G4VEmAdjointModel::fCurrentMaterial, G4VEmAdjointModel::fTcutSecond, G4ProductionCutsTable::GetEnergyCutsVector(), G4MaterialCutsCouple::GetIndex(), G4MaterialCutsCouple::GetMaterial(), and G4ProductionCutsTable::GetProductionCutsTable().

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4AdjointComptonModel::AdjointCrossSection(), G4AdjointhIonisationModel::AdjointCrossSection(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4AdjointhIonisationModel::RapidSampleSecondaries(), and G4AdjointBremsstrahlungModel::SampleSecondaries().

DefineCurrentMaterialAndElectronEnergy()

void G4AdjointPhotoElectricModel::DefineCurrentMaterialAndElectronEnergy ( const G4MaterialCutsCouple *  aCouple, G4double  eEnergy  )

private

Definition at line 228 of file G4AdjointPhotoElectricModel.cc.

{
  fCurrentCouple   = const_cast<G4MaterialCutsCouple*>(couple);
  fCurrentMaterial = const_cast<G4Material*>(couple->GetMaterial());
  fCurrenteEnergy = anEnergy;
  fDirectModel->SetCurrentCouple(couple);
}

References G4VEmAdjointModel::fCurrentCouple, fCurrenteEnergy, G4VEmAdjointModel::fCurrentMaterial, G4VEmAdjointModel::fDirectModel, G4MaterialCutsCouple::GetMaterial(), and G4VEmModel::SetCurrentCouple().

Referenced by AdjointCrossSection().

DefineDirectEMModel()

void G4VEmAdjointModel::DefineDirectEMModel ( G4VEmModel *  aModel )

inlineinherited

Definition at line 160 of file G4VEmAdjointModel.hh.

{ fDirectModel = aModel; }

References G4VEmAdjointModel::fDirectModel.

DiffCrossSectionFunction1()

G4double G4VEmAdjointModel::DiffCrossSectionFunction1 ( G4double  kinEnergyProj )

protectedinherited

Definition at line 155 of file G4VEmAdjointModel.cc.

{
  G4double bias_factor =
    fCsBiasingFactor * fKinEnergyProdForIntegration / kinEnergyProj;
 
  if(fUseMatrixPerElement)
  {
    return DiffCrossSectionPerAtomPrimToSecond(
             kinEnergyProj, fKinEnergyProdForIntegration, fZSelectedNucleus,
             fASelectedNucleus) *
           bias_factor;
  }
  else
  {
    return DiffCrossSectionPerVolumePrimToSecond(
             fSelectedMaterial, kinEnergyProj, fKinEnergyProdForIntegration) *
           bias_factor;
  }
}

References G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond(), G4VEmAdjointModel::fASelectedNucleus, G4VEmAdjointModel::fCsBiasingFactor, G4VEmAdjointModel::fKinEnergyProdForIntegration, G4VEmAdjointModel::fSelectedMaterial, G4VEmAdjointModel::fUseMatrixPerElement, and G4VEmAdjointModel::fZSelectedNucleus.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForSecond(), and G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForSecond().

DiffCrossSectionFunction2()

G4double G4VEmAdjointModel::DiffCrossSectionFunction2 ( G4double  kinEnergyProj )

protectedinherited

Definition at line 176 of file G4VEmAdjointModel.cc.

{
  G4double bias_factor =
    fCsBiasingFactor * fKinEnergyScatProjForIntegration / kinEnergyProj;
  if(fUseMatrixPerElement)
  {
    return DiffCrossSectionPerAtomPrimToScatPrim(
             kinEnergyProj, fKinEnergyScatProjForIntegration, fZSelectedNucleus,
             fASelectedNucleus) *
           bias_factor;
  }
  else
  {
    return DiffCrossSectionPerVolumePrimToScatPrim(
             fSelectedMaterial, kinEnergyProj,
             fKinEnergyScatProjForIntegration) *
           bias_factor;
  }
}

References G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToScatPrim(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToScatPrim(), G4VEmAdjointModel::fASelectedNucleus, G4VEmAdjointModel::fCsBiasingFactor, G4VEmAdjointModel::fKinEnergyScatProjForIntegration, G4VEmAdjointModel::fSelectedMaterial, G4VEmAdjointModel::fUseMatrixPerElement, and G4VEmAdjointModel::fZSelectedNucleus.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), and G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForScatProj().

DiffCrossSectionPerAtomPrimToScatPrim()

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

virtualinherited

Reimplemented in G4AdjointComptonModel.

Definition at line 105 of file G4VEmAdjointModel.cc.

{
  G4double kinEnergyProd = kinEnergyProj - kinEnergyScatProj;
  G4double dSigmadEprod;
  if(kinEnergyProd <= 0.)
    dSigmadEprod = 0.;
  else
    dSigmadEprod =
      DiffCrossSectionPerAtomPrimToSecond(kinEnergyProj, kinEnergyProd, Z, A);
  return dSigmadEprod;
}

References A, G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), and Z.

Referenced by G4VEmAdjointModel::DiffCrossSectionFunction2(), and G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom().

DiffCrossSectionPerAtomPrimToSecond()

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

virtualinherited

Reimplemented in G4AdjointComptonModel, G4AdjointeIonisationModel, G4AdjointhIonisationModel, and G4AdjointIonIonisationModel.

Definition at line 82 of file G4VEmAdjointModel.cc.

{
  G4double dSigmadEprod = 0.;
  G4double Emax_proj    = GetSecondAdjEnergyMaxForProdToProj(kinEnergyProd);
  G4double Emin_proj    = GetSecondAdjEnergyMinForProdToProj(kinEnergyProd);
 
  // the produced particle should have a kinetic energy less than the projectile
  if(kinEnergyProj > Emin_proj && kinEnergyProj <= Emax_proj)
  {
    G4double E1     = kinEnergyProd;
    G4double E2     = kinEnergyProd * 1.000001;
    G4double sigma1 = fDirectModel->ComputeCrossSectionPerAtom(
      fDirectPrimaryPart, kinEnergyProj, Z, A, E1, 1.e20);
    G4double sigma2 = fDirectModel->ComputeCrossSectionPerAtom(
      fDirectPrimaryPart, kinEnergyProj, Z, A, E2, 1.e20);
 
    dSigmadEprod = (sigma1 - sigma2) / (E2 - E1);
  }
  return dSigmadEprod;
}

References A, G4VEmModel::ComputeCrossSectionPerAtom(), G4VEmAdjointModel::fDirectModel, G4VEmAdjointModel::fDirectPrimaryPart, G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj(), and Z.

Referenced by G4VEmAdjointModel::DiffCrossSectionFunction1(), G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToScatPrim(), and G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom().

DiffCrossSectionPerVolumePrimToScatPrim()

G4double G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToScatPrim ( const G4Material *  aMaterial, G4double  kinEnergyProj, G4double  kinEnergyScatProj  )

virtualinherited

Definition at line 140 of file G4VEmAdjointModel.cc.

{
  G4double kinEnergyProd = kinEnergyProj - kinEnergyScatProj;
  G4double dSigmadEprod;
  if(kinEnergyProd <= 0.)
    dSigmadEprod = 0.;
  else
    dSigmadEprod = DiffCrossSectionPerVolumePrimToSecond(
      aMaterial, kinEnergyProj, kinEnergyProd);
  return dSigmadEprod;
}

References G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond().

Referenced by G4VEmAdjointModel::DiffCrossSectionFunction2(), and G4AdjointeIonisationModel::SampleSecondaries().

DiffCrossSectionPerVolumePrimToSecond()

G4double G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond ( const G4Material *  aMaterial, G4double  kinEnergyProj, G4double  kinEnergyProd  )

virtualinherited

Reimplemented in G4AdjointBremsstrahlungModel.

Definition at line 119 of file G4VEmAdjointModel.cc.

{
  G4double dSigmadEprod = 0.;
  G4double Emax_proj    = GetSecondAdjEnergyMaxForProdToProj(kinEnergyProd);
  G4double Emin_proj    = GetSecondAdjEnergyMinForProdToProj(kinEnergyProd);
 
  if(kinEnergyProj > Emin_proj && kinEnergyProj <= Emax_proj)
  {
    G4double E1     = kinEnergyProd;
    G4double E2     = kinEnergyProd * 1.0001;
    G4double sigma1 = fDirectModel->CrossSectionPerVolume(
      aMaterial, fDirectPrimaryPart, kinEnergyProj, E1, 1.e20);
    G4double sigma2 = fDirectModel->CrossSectionPerVolume(
      aMaterial, fDirectPrimaryPart, kinEnergyProj, E2, 1.e20);
    dSigmadEprod = (sigma1 - sigma2) / (E2 - E1);
  }
  return dSigmadEprod;
}

References G4VEmModel::CrossSectionPerVolume(), G4VEmAdjointModel::fDirectModel, G4VEmAdjointModel::fDirectPrimaryPart, G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), and G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj().

Referenced by G4VEmAdjointModel::DiffCrossSectionFunction1(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToScatPrim(), G4AdjointBremsstrahlungModel::DiffCrossSectionPerVolumePrimToSecond(), and G4AdjointeIonisationModel::SampleSecondaries().

GetAdjointEquivalentOfDirectPrimaryParticleDefinition()

G4ParticleDefinition * G4VEmAdjointModel::GetAdjointEquivalentOfDirectPrimaryParticleDefinition ( )

inlineinherited

Definition at line 141 of file G4VEmAdjointModel.hh.

  {
    return fAdjEquivDirectPrimPart;
  }

References G4VEmAdjointModel::fAdjEquivDirectPrimPart.

Referenced by G4AdjointCSManager::ComputeTotalAdjointCS().

GetAdjointEquivalentOfDirectSecondaryParticleDefinition()

G4ParticleDefinition * G4VEmAdjointModel::GetAdjointEquivalentOfDirectSecondaryParticleDefinition ( )

inlineinherited

Definition at line 147 of file G4VEmAdjointModel.hh.

  {
    return fAdjEquivDirectSecondPart;
  }

References G4VEmAdjointModel::fAdjEquivDirectSecondPart.

Referenced by G4AdjointCSManager::ComputeTotalAdjointCS().

GetApplyCutInRange()

G4bool G4VEmAdjointModel::GetApplyCutInRange ( )

inlineinherited

Definition at line 201 of file G4VEmAdjointModel.hh.

{ return fApplyCutInRange; }

References G4VEmAdjointModel::fApplyCutInRange.

Referenced by G4AdjointCSManager::ComputeAdjointCS(), and G4AdjointCSManager::ComputeTotalAdjointCS().

GetHighEnergyLimit()

G4double G4VEmAdjointModel::GetHighEnergyLimit ( )

inlineinherited

Definition at line 152 of file G4VEmAdjointModel.hh.

{ return fHighEnergyLimit; }

References G4VEmAdjointModel::fHighEnergyLimit.

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndElement(), G4AdjointCSManager::BuildCrossSectionsModelAndMaterial(), G4AdjointhIonisationModel::GetSecondAdjEnergyMaxForProdToProj(), G4AdjointIonIonisationModel::GetSecondAdjEnergyMaxForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj(), G4AdjointComptonModel::GetSecondAdjEnergyMaxForScatProjToProj(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4AdjointhIonisationModel::RapidSampleSecondaries(), G4AdjointBremsstrahlungModel::SampleSecondaries(), G4AdjointComptonModel::SampleSecondaries(), G4AdjointeIonisationModel::SampleSecondaries(), G4AdjointhIonisationModel::SampleSecondaries(), and G4AdjointIonIonisationModel::SampleSecondaries().

GetLowEnergyLimit()

G4double G4VEmAdjointModel::GetLowEnergyLimit ( )

inlineinherited

Definition at line 154 of file G4VEmAdjointModel.hh.

{ return fLowEnergyLimit; }

References G4VEmAdjointModel::fLowEnergyLimit.

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndElement(), G4AdjointCSManager::BuildCrossSectionsModelAndMaterial(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForScatProj(), and G4AdjointCSManager::ComputeTotalAdjointCS().

GetName()

G4String G4VEmAdjointModel::GetName ( )

inlineinherited

Definition at line 203 of file G4VEmAdjointModel.hh.

{ return fName; }

References G4VEmAdjointModel::fName.

GetSecondAdjEnergyMaxForProdToProj()

G4double G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj ( G4double  primAdjEnergy )

virtualinherited

Reimplemented in G4AdjointhIonisationModel, and G4AdjointIonIonisationModel.

Definition at line 668 of file G4VEmAdjointModel.cc.

{
  return fHighEnergyLimit;
}

References G4VEmAdjointModel::fHighEnergyLimit.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4AdjointComptonModel::AdjointCrossSection(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForSecond(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForSecond(), G4AdjointCSManager::ComputeAdjointCS(), G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), G4AdjointeIonisationModel::DiffCrossSectionPerAtomPrimToSecond(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom(), and G4AdjointeIonisationModel::SampleSecondaries().

GetSecondAdjEnergyMaxForScatProjToProj()

G4double G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj ( G4double  primAdjEnergy )

virtualinherited

Reimplemented in G4AdjointComptonModel, G4AdjointhIonisationModel, and G4AdjointIonIonisationModel.

Definition at line 648 of file G4VEmAdjointModel.cc.

{
  G4double maxEProj = GetHighEnergyLimit();
  if(fSecondPartSameType)
    maxEProj = std::min(kinEnergyScatProj * 2., maxEProj);
  return maxEProj;
}

References G4VEmAdjointModel::fSecondPartSameType, G4VEmAdjointModel::GetHighEnergyLimit(), and G4INCL::Math::min().

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForScatProj(), G4AdjointCSManager::ComputeAdjointCS(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom(), and G4AdjointeIonisationModel::SampleSecondaries().

GetSecondAdjEnergyMinForProdToProj()

G4double G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj ( G4double  primAdjEnergy )

virtualinherited

Reimplemented in G4AdjointComptonModel, G4AdjointhIonisationModel, and G4AdjointIonIonisationModel.

Definition at line 674 of file G4VEmAdjointModel.cc.

{
  G4double minEProj = primAdjEnergy;
  if(fSecondPartSameType)
    minEProj = primAdjEnergy * 2.;
  return minEProj;
}

References G4VEmAdjointModel::fSecondPartSameType.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForSecond(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForSecond(), G4AdjointCSManager::ComputeAdjointCS(), G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), G4AdjointeIonisationModel::DiffCrossSectionPerAtomPrimToSecond(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom(), and G4AdjointeIonisationModel::SampleSecondaries().

GetSecondAdjEnergyMinForScatProjToProj()

G4double G4VEmAdjointModel::GetSecondAdjEnergyMinForScatProjToProj ( G4double  primAdjEnergy, G4double  tcut = 0.  )

virtualinherited

Reimplemented in G4AdjointhIonisationModel, and G4AdjointIonIonisationModel.

Definition at line 658 of file G4VEmAdjointModel.cc.

{
  G4double Emin = primAdjEnergy;
  if(fApplyCutInRange)
    Emin += tcut;
  return Emin;
}

References Emin, and G4VEmAdjointModel::fApplyCutInRange.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4AdjointComptonModel::AdjointCrossSection(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForScatProj(), G4AdjointCSManager::ComputeAdjointCS(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), and G4AdjointeIonisationModel::SampleSecondaries().

GetSecondPartOfSameType()

G4bool G4VEmAdjointModel::GetSecondPartOfSameType ( )

inlineinherited

Definition at line 176 of file G4VEmAdjointModel.hh.

{ return fSecondPartSameType; }

References G4VEmAdjointModel::fSecondPartSameType.

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndElement(), and G4AdjointCSManager::BuildCrossSectionsModelAndMaterial().

GetUseMatrix()

G4bool G4VEmAdjointModel::GetUseMatrix ( )

inlineinherited

Definition at line 192 of file G4VEmAdjointModel.hh.

{ return fUseMatrix; }

References G4VEmAdjointModel::fUseMatrix.

Referenced by G4AdjointCSManager::ComputeAdjointCS().

GetUseMatrixPerElement()

G4bool G4VEmAdjointModel::GetUseMatrixPerElement ( )

inlineinherited

Definition at line 194 of file G4VEmAdjointModel.hh.

{ return fUseMatrixPerElement; }

References G4VEmAdjointModel::fUseMatrixPerElement.

Referenced by G4AdjointCSManager::ComputeAdjointCS().

GetUseOnlyOneMatrixForAllElements()

G4bool G4VEmAdjointModel::GetUseOnlyOneMatrixForAllElements ( )

inlineinherited

Definition at line 196 of file G4VEmAdjointModel.hh.

  {
    return fOneMatrixForAllElements;
  }

References G4VEmAdjointModel::fOneMatrixForAllElements.

Referenced by G4AdjointCSManager::ComputeAdjointCS().

operator=()

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

delete

SampleAdjSecEnergyFromCSMatrix() [1/2]

G4double G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix ( G4double  prim_energy, G4bool  isScatProjToProj  )

protectedinherited

Definition at line 518 of file G4VEmAdjointModel.cc.

{
  SelectCSMatrix(isScatProjToProj);
  return SampleAdjSecEnergyFromCSMatrix(fCSMatrixUsed, aPrimEnergy,
                                        isScatProjToProj);
}

References G4VEmAdjointModel::fCSMatrixUsed, G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), and G4VEmAdjointModel::SelectCSMatrix().

SampleAdjSecEnergyFromCSMatrix() [2/2]

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

protectedinherited

Definition at line 413 of file G4VEmAdjointModel.cc.

{
  G4AdjointCSMatrix* theMatrix = (*fCSMatrixProdToProjBackScat)[MatrixIndex];
  if(isScatProjToProj)
    theMatrix = (*fCSMatrixProjToProjBackScat)[MatrixIndex];
  std::vector<G4double>* theLogPrimEnergyVector =
    theMatrix->GetLogPrimEnergyVector();
 
  if(theLogPrimEnergyVector->empty())
  {
    G4cout << "No data are contained in the given AdjointCSMatrix!" << G4endl;
    G4cout << "The sampling procedure will be stopped." << G4endl;
    return 0.;
  }
 
  G4AdjointInterpolator* theInterpolator = G4AdjointInterpolator::GetInstance();
  G4double aLogPrimEnergy                = std::log(aPrimEnergy);
  size_t ind = theInterpolator->FindPositionForLogVector(
    aLogPrimEnergy, *theLogPrimEnergyVector);
 
  G4double aLogPrimEnergy1, aLogPrimEnergy2;
  G4double aLogCS1, aLogCS2;
  G4double log01, log02;
  std::vector<double>* aLogSecondEnergyVector1 = nullptr;
  std::vector<double>* aLogSecondEnergyVector2 = nullptr;
  std::vector<double>* aLogProbVector1         = nullptr;
  std::vector<double>* aLogProbVector2         = nullptr;
  std::vector<size_t>* aLogProbVectorIndex1    = nullptr;
  std::vector<size_t>* aLogProbVectorIndex2    = nullptr;
 
  theMatrix->GetData(ind, aLogPrimEnergy1, aLogCS1, log01,
                     aLogSecondEnergyVector1, aLogProbVector1,
                     aLogProbVectorIndex1 );
  theMatrix->GetData(ind + 1, aLogPrimEnergy2, aLogCS2, log02,
                     aLogSecondEnergyVector2, aLogProbVector2,
                     aLogProbVectorIndex2);
 
  if (! (aLogProbVector1 && aLogProbVector2 &&
               aLogSecondEnergyVector1 && aLogSecondEnergyVector2)){
     return  0.;
  }
 
  G4double rand_var     = G4UniformRand();
  G4double log_rand_var = std::log(rand_var);
  G4double log_Tcut     = std::log(fTcutSecond);
  G4double Esec         = 0.;
  G4double log_dE1, log_dE2;
  G4double log_rand_var1, log_rand_var2;
  G4double log_E1, log_E2;
  log_rand_var1 = log_rand_var;
  log_rand_var2 = log_rand_var;
 
  G4double Emin = 0.;
  G4double Emax = 0.;
  if(theMatrix->IsScatProjToProj())
  {  // case where Tcut plays a role
    Emin = GetSecondAdjEnergyMinForScatProjToProj(aPrimEnergy, fTcutSecond);
    Emax = GetSecondAdjEnergyMaxForScatProjToProj(aPrimEnergy);
    G4double dE = 0.;
    if(Emin < Emax)
    {
      if(fApplyCutInRange)
      {
        if(fSecondPartSameType && fTcutSecond > aPrimEnergy)
          return aPrimEnergy;
 
        log_rand_var1 = log_rand_var +
                        theInterpolator->InterpolateForLogVector(
                          log_Tcut, *aLogSecondEnergyVector1, *aLogProbVector1);
        log_rand_var2 = log_rand_var +
                        theInterpolator->InterpolateForLogVector(
                          log_Tcut, *aLogSecondEnergyVector2, *aLogProbVector2);
      }
      log_dE1 = theInterpolator->Interpolate(log_rand_var1, *aLogProbVector1,
                                             *aLogSecondEnergyVector1, "Lin");
      log_dE2 = theInterpolator->Interpolate(log_rand_var2, *aLogProbVector2,
                                             *aLogSecondEnergyVector2, "Lin");
      dE      = std::exp(theInterpolator->LinearInterpolation(
        aLogPrimEnergy, aLogPrimEnergy1, aLogPrimEnergy2, log_dE1, log_dE2));
    }
 
    Esec = aPrimEnergy + dE;
    Esec = std::max(Esec, Emin);
    Esec = std::min(Esec, Emax);
  }
  else
  {  // Tcut condition is already full-filled
 
    log_E1 = theInterpolator->Interpolate(log_rand_var, *aLogProbVector1,
                                          *aLogSecondEnergyVector1, "Lin");
    log_E2 = theInterpolator->Interpolate(log_rand_var, *aLogProbVector2,
                                          *aLogSecondEnergyVector2, "Lin");
 
    Esec = std::exp(theInterpolator->LinearInterpolation(
      aLogPrimEnergy, aLogPrimEnergy1, aLogPrimEnergy2, log_E1, log_E2));
    Emin = GetSecondAdjEnergyMinForProdToProj(aPrimEnergy);
    Emax = GetSecondAdjEnergyMaxForProdToProj(aPrimEnergy);
    Esec = std::max(Esec, Emin);
    Esec = std::min(Esec, Emax);
  }
  return Esec;
}

References dE, Emax, Emin, G4VEmAdjointModel::fApplyCutInRange, G4AdjointInterpolator::FindPositionForLogVector(), G4VEmAdjointModel::fSecondPartSameType, G4VEmAdjointModel::fTcutSecond, G4cout, G4endl, G4UniformRand, G4AdjointCSMatrix::GetData(), G4AdjointInterpolator::GetInstance(), G4AdjointCSMatrix::GetLogPrimEnergyVector(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForScatProjToProj(), G4AdjointInterpolator::Interpolate(), G4AdjointInterpolator::InterpolateForLogVector(), G4AdjointCSMatrix::IsScatProjToProj(), G4AdjointInterpolator::LinearInterpolation(), G4INCL::Math::max(), and G4INCL::Math::min().

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), G4AdjointBremsstrahlungModel::SampleSecondaries(), G4AdjointComptonModel::SampleSecondaries(), G4AdjointeIonisationModel::SampleSecondaries(), G4AdjointhIonisationModel::SampleSecondaries(), and G4AdjointIonIonisationModel::SampleSecondaries().

SampleAdjSecEnergyFromDiffCrossSectionPerAtom()

G4double G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom ( G4double  prim_energy, G4bool  isScatProjToProj  )

protectedvirtualinherited

Definition at line 557 of file G4VEmAdjointModel.cc.

{
  // here we try to use the rejection method
  constexpr G4int iimax = 1000;
  G4double E            = 0.;
  G4double x, xmin, greject;
  if(isScatProjToProj)
  {
    G4double Emax = GetSecondAdjEnergyMaxForScatProjToProj(prim_energy);
    G4double Emin = prim_energy + fTcutSecond;
    xmin          = Emin / Emax;
    G4double grejmax =
      DiffCrossSectionPerAtomPrimToScatPrim(Emin, prim_energy, 1) * prim_energy;
 
    G4int ii = 0;
    do
    {
      // q = G4UniformRand();
      x = 1. / (G4UniformRand() * (1. / xmin - 1.) + 1.);
      E = x * Emax;
      greject =
        DiffCrossSectionPerAtomPrimToScatPrim(E, prim_energy, 1) * prim_energy;
      ++ii;
      if(ii >= iimax)
      {
        break;
      }
    }
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
    while(greject < G4UniformRand() * grejmax);
  }
  else
  {
    G4double Emax = GetSecondAdjEnergyMaxForProdToProj(prim_energy);
    G4double Emin = GetSecondAdjEnergyMinForProdToProj(prim_energy);
    xmin          = Emin / Emax;
    G4double grejmax =
      DiffCrossSectionPerAtomPrimToSecond(Emin, prim_energy, 1);
    G4int ii = 0;
    do
    {
      x       = std::pow(xmin, G4UniformRand());
      E       = x * Emax;
      greject = DiffCrossSectionPerAtomPrimToSecond(E, prim_energy, 1);
      ++ii;
      if(ii >= iimax)
      {
        break;
      }
    }
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
    while(greject < G4UniformRand() * grejmax);
  }
 
  return E;
}

References G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToScatPrim(), G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), Emax, Emin, G4VEmAdjointModel::fTcutSecond, G4UniformRand, G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj(), and G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj().

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

References G4ParticleChange::AddSecondary(), AdjointCrossSection(), G4AdjointGamma::AdjointGamma(), anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CorrectPostStepWeight(), source.hepunit::electron_mass_c2, G4VEmAdjointModel::fCurrentMaterial, fIndexElement, fPostStepAdjointCS, fPreStepAdjointCS, fShellProb, fStopAndKill, fTotAdjointCS, fXsec, G4UniformRand, G4Track::GetDynamicParticle(), G4Material::GetElementVector(), G4DynamicParticle::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), G4DynamicParticle::GetMomentumDirection(), G4Material::GetNumberOfElements(), G4Track::GetWeight(), G4VParticleChange::ProposeTrackStatus(), CLHEP::Hep3Vector::rotateUz(), and twopi.

SelectCSMatrix()

void G4VEmAdjointModel::SelectCSMatrix ( G4bool  isScatProjToProj )

protectedinherited

Definition at line 527 of file G4VEmAdjointModel.cc.

{
  fCSMatrixUsed = 0;
  if(!fUseMatrixPerElement)
    fCSMatrixUsed = fCurrentMaterial->GetIndex();
  else if(!fOneMatrixForAllElements)
  {  // Select Material
    std::vector<G4double>* CS_Vs_Element = &fElementCSScatProjToProj;
    fLastCS                              = fLastAdjointCSForScatProjToProj;
    if(!isScatProjToProj)
    {
      CS_Vs_Element = &fElementCSProdToProj;
      fLastCS       = fLastAdjointCSForProdToProj;
    }
    G4double SumCS = 0.;
    size_t ind     = 0;
    for(size_t i = 0; i < CS_Vs_Element->size(); ++i)
    {
      SumCS += (*CS_Vs_Element)[i];
      if(G4UniformRand() <= SumCS / fLastCS)
      {
        ind = i;
        break;
      }
    }
    fCSMatrixUsed = fCurrentMaterial->GetElement(ind)->GetIndex();
  }
}

References G4VEmAdjointModel::fCSMatrixUsed, G4VEmAdjointModel::fCurrentMaterial, G4VEmAdjointModel::fElementCSProdToProj, G4VEmAdjointModel::fElementCSScatProjToProj, G4VEmAdjointModel::fLastAdjointCSForProdToProj, G4VEmAdjointModel::fLastAdjointCSForScatProjToProj, G4VEmAdjointModel::fLastCS, G4VEmAdjointModel::fOneMatrixForAllElements, G4VEmAdjointModel::fUseMatrixPerElement, G4UniformRand, G4Material::GetElement(), G4Element::GetIndex(), and G4Material::GetIndex().

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

SetAdditionalWeightCorrectionFactorForPostStepOutsideModel()

void G4VEmAdjointModel::SetAdditionalWeightCorrectionFactorForPostStepOutsideModel ( G4double  factor )

inlineinherited

Definition at line 215 of file G4VEmAdjointModel.hh.

  {
    fOutsideWeightFactor = factor;
  }

References G4VEmAdjointModel::fOutsideWeightFactor.

Referenced by G4AdjointForcedInteractionForGamma::PostStepDoIt().

SetAdjointEquivalentOfDirectPrimaryParticleDefinition()

void G4VEmAdjointModel::SetAdjointEquivalentOfDirectPrimaryParticleDefinition ( G4ParticleDefinition *  aPart )

inherited

Definition at line 729 of file G4VEmAdjointModel.cc.

{
  fAdjEquivDirectPrimPart = aPart;
  if(fAdjEquivDirectPrimPart->GetParticleName() == "adj_e-")
    fDirectPrimaryPart = G4Electron::Electron();
  else if(fAdjEquivDirectPrimPart->GetParticleName() == "adj_gamma")
    fDirectPrimaryPart = G4Gamma::Gamma();
}

References G4Electron::Electron(), G4VEmAdjointModel::fAdjEquivDirectPrimPart, G4VEmAdjointModel::fDirectPrimaryPart, G4Gamma::Gamma(), and G4ParticleDefinition::GetParticleName().

SetAdjointEquivalentOfDirectSecondaryParticleDefinition()

void G4VEmAdjointModel::SetAdjointEquivalentOfDirectSecondaryParticleDefinition ( G4ParticleDefinition *  aPart )

inlineinherited

Definition at line 165 of file G4VEmAdjointModel.hh.

  {
    fAdjEquivDirectSecondPart = aPart;
  }

References G4VEmAdjointModel::fAdjEquivDirectSecondPart.

SetApplyCutInRange()

void G4VEmAdjointModel::SetApplyCutInRange ( G4bool  aBool )

inlineinherited

Definition at line 190 of file G4VEmAdjointModel.hh.

{ fApplyCutInRange = aBool; }

References G4VEmAdjointModel::fApplyCutInRange.

Referenced by G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointPhotoElectricModel(), and G4AdjointBremsstrahlungModel::Initialize().

SetCorrectWeightForPostStepInModel()

void G4VEmAdjointModel::SetCorrectWeightForPostStepInModel ( G4bool  aBool )

inlineinherited

Definition at line 210 of file G4VEmAdjointModel.hh.

  {
    fInModelWeightCorr = aBool;
  }

References G4VEmAdjointModel::fInModelWeightCorr.

Referenced by G4AdjointForcedInteractionForGamma::PostStepDoIt().

SetCSBiasingFactor()

virtual void G4VEmAdjointModel::SetCSBiasingFactor ( G4double  aVal )

inlinevirtualinherited

Definition at line 205 of file G4VEmAdjointModel.hh.

  {
    fCsBiasingFactor = aVal;
  }

References G4VEmAdjointModel::fCsBiasingFactor.

SetCSMatrices()

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

inlineinherited

Definition at line 133 of file G4VEmAdjointModel.hh.

  {
    fCSMatrixProdToProjBackScat = Vec1CSMatrix;
    fCSMatrixProjToProjBackScat = Vec2CSMatrix;
  };

References G4VEmAdjointModel::fCSMatrixProdToProjBackScat, and G4VEmAdjointModel::fCSMatrixProjToProjBackScat.

SetHighEnergyLimit()

void G4VEmAdjointModel::SetHighEnergyLimit ( G4double  aVal )

inherited

Definition at line 713 of file G4VEmAdjointModel.cc.

{
  fHighEnergyLimit = aVal;
  if(fDirectModel)
    fDirectModel->SetHighEnergyLimit(aVal);
}

References G4VEmAdjointModel::fDirectModel, G4VEmAdjointModel::fHighEnergyLimit, and G4VEmModel::SetHighEnergyLimit().

SetLowEnergyLimit()

void G4VEmAdjointModel::SetLowEnergyLimit ( G4double  aVal )

inherited

Definition at line 721 of file G4VEmAdjointModel.cc.

{
  fLowEnergyLimit = aVal;
  if(fDirectModel)
    fDirectModel->SetLowEnergyLimit(aVal);
}

References G4VEmAdjointModel::fDirectModel, G4VEmAdjointModel::fLowEnergyLimit, and G4VEmModel::SetLowEnergyLimit().

SetSecondPartOfSameType()

void G4VEmAdjointModel::SetSecondPartOfSameType ( G4bool  aBool )

inlineinherited

Definition at line 171 of file G4VEmAdjointModel.hh.

  {
    fSecondPartSameType = aBool;
  }

References G4VEmAdjointModel::fSecondPartSameType.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4eInverseCompton::G4eInverseCompton(), G4eInverseIonisation::G4eInverseIonisation(), G4hInverseIonisation::G4hInverseIonisation(), G4InversePEEffect::G4InversePEEffect(), and G4IonInverseIonisation::G4IonInverseIonisation().

SetUseMatrix()

void G4VEmAdjointModel::SetUseMatrix ( G4bool  aBool )

inlineinherited

Definition at line 178 of file G4VEmAdjointModel.hh.

{ fUseMatrix = aBool; }

References G4VEmAdjointModel::fUseMatrix.

Referenced by G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointPhotoElectricModel(), and G4AdjointBremsstrahlungModel::Initialize().

SetUseMatrixPerElement()

void G4VEmAdjointModel::SetUseMatrixPerElement ( G4bool  aBool )

inlineinherited

Definition at line 180 of file G4VEmAdjointModel.hh.

  {
    fUseMatrixPerElement = aBool;
  }

References G4VEmAdjointModel::fUseMatrixPerElement.

Referenced by G4AdjointComptonModel::G4AdjointComptonModel(), and G4AdjointBremsstrahlungModel::Initialize().

SetUseOnlyOneMatrixForAllElements()

void G4VEmAdjointModel::SetUseOnlyOneMatrixForAllElements ( G4bool  aBool )

inlineinherited

Definition at line 185 of file G4VEmAdjointModel.hh.

  {
    fOneMatrixForAllElements = aBool;
  }

References G4VEmAdjointModel::fOneMatrixForAllElements.

Referenced by G4AdjointComptonModel::G4AdjointComptonModel().

Field Documentation

fAdjEquivDirectPrimPart

G4ParticleDefinition* G4VEmAdjointModel::fAdjEquivDirectPrimPart = nullptr

protectedinherited

Definition at line 259 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointPhotoElectricModel(), G4VEmAdjointModel::GetAdjointEquivalentOfDirectPrimaryParticleDefinition(), G4AdjointBremsstrahlungModel::Initialize(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4AdjointhIonisationModel::RapidSampleSecondaries(), G4AdjointBremsstrahlungModel::SampleSecondaries(), G4AdjointComptonModel::SampleSecondaries(), G4AdjointeIonisationModel::SampleSecondaries(), G4AdjointhIonisationModel::SampleSecondaries(), G4AdjointIonIonisationModel::SampleSecondaries(), G4VEmAdjointModel::SetAdjointEquivalentOfDirectPrimaryParticleDefinition(), and G4AdjointIonIonisationModel::SetIon().

fAdjEquivDirectSecondPart

G4ParticleDefinition* G4VEmAdjointModel::fAdjEquivDirectSecondPart = nullptr

protectedinherited

Definition at line 260 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::DefineCurrentMaterial(), G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4AdjointPhotoElectricModel(), G4VEmAdjointModel::GetAdjointEquivalentOfDirectSecondaryParticleDefinition(), G4AdjointBremsstrahlungModel::Initialize(), G4AdjointhIonisationModel::RapidSampleSecondaries(), G4AdjointeIonisationModel::SampleSecondaries(), G4AdjointhIonisationModel::SampleSecondaries(), G4AdjointIonIonisationModel::SampleSecondaries(), and G4VEmAdjointModel::SetAdjointEquivalentOfDirectSecondaryParticleDefinition().

fApplyCutInRange

G4bool G4VEmAdjointModel::fApplyCutInRange = true

protectedinherited

Definition at line 304 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4VEmAdjointModel::GetApplyCutInRange(), G4VEmAdjointModel::GetSecondAdjEnergyMinForScatProjToProj(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), and G4VEmAdjointModel::SetApplyCutInRange().

fASelectedNucleus

G4int G4VEmAdjointModel::fASelectedNucleus = 0

protectedinherited

Definition at line 295 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForSecond(), G4VEmAdjointModel::DiffCrossSectionFunction1(), and G4VEmAdjointModel::DiffCrossSectionFunction2().

fCsBiasingFactor

G4double G4VEmAdjointModel::fCsBiasingFactor = 1.

protectedinherited

Definition at line 287 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4VEmAdjointModel::CorrectPostStepWeight(), G4AdjointIonIonisationModel::CorrectPostStepWeight(), G4VEmAdjointModel::DiffCrossSectionFunction1(), G4VEmAdjointModel::DiffCrossSectionFunction2(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), and G4VEmAdjointModel::SetCSBiasingFactor().

fCSManager

G4AdjointCSManager* G4VEmAdjointModel::fCSManager

protectedinherited

Definition at line 249 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4VEmAdjointModel::CorrectPostStepWeight(), and G4VEmAdjointModel::G4VEmAdjointModel().

fCSMatrixProdToProjBackScat

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

protectedinherited

Definition at line 264 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::SetCSMatrices(), and G4VEmAdjointModel::~G4VEmAdjointModel().

fCSMatrixProjToProjBackScat

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

protectedinherited

Definition at line 265 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::SetCSMatrices(), and G4VEmAdjointModel::~G4VEmAdjointModel().

fCSMatrixUsed

size_t G4VEmAdjointModel::fCSMatrixUsed = 0

protectedinherited

Definition at line 298 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), and G4VEmAdjointModel::SelectCSMatrix().

fCurrentCouple

G4MaterialCutsCouple* G4VEmAdjointModel::fCurrentCouple = nullptr

protectedinherited

Definition at line 256 of file G4VEmAdjointModel.hh.

Referenced by AdjointCrossSection(), G4VEmAdjointModel::CorrectPostStepWeight(), G4VEmAdjointModel::DefineCurrentMaterial(), DefineCurrentMaterialAndElectronEnergy(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), and G4AdjointComptonModel::RapidSampleSecondaries().

fCurrenteEnergy

G4double G4AdjointPhotoElectricModel::fCurrenteEnergy = 0.

private

Definition at line 89 of file G4AdjointPhotoElectricModel.hh.

Referenced by AdjointCrossSection(), and DefineCurrentMaterialAndElectronEnergy().

fCurrentMaterial

G4Material* G4VEmAdjointModel::fCurrentMaterial = nullptr

protectedinherited

Definition at line 255 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4AdjointComptonModel::AdjointCrossSection(), G4AdjointhIonisationModel::AdjointCrossSection(), AdjointCrossSection(), G4AdjointIonIonisationModel::CorrectPostStepWeight(), G4VEmAdjointModel::DefineCurrentMaterial(), DefineCurrentMaterialAndElectronEnergy(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4AdjointhIonisationModel::RapidSampleSecondaries(), G4AdjointeIonisationModel::SampleSecondaries(), SampleSecondaries(), and G4VEmAdjointModel::SelectCSMatrix().

fDirectModel

G4VEmModel* G4VEmAdjointModel::fDirectModel = nullptr

protectedinherited

Definition at line 250 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), AdjointCrossSectionPerAtom(), G4AdjointIonIonisationModel::CorrectPostStepWeight(), DefineCurrentMaterialAndElectronEnergy(), G4VEmAdjointModel::DefineDirectEMModel(), G4AdjointComptonModel::DiffCrossSectionPerAtomPrimToScatPrim(), G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), G4AdjointhIonisationModel::DiffCrossSectionPerAtomPrimToSecond(), G4AdjointIonIonisationModel::DiffCrossSectionPerAtomPrimToSecond(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond(), G4AdjointBremsstrahlungModel::G4AdjointBremsstrahlungModel(), G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointPhotoElectricModel(), G4AdjointBremsstrahlungModel::Initialize(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4VEmAdjointModel::SetHighEnergyLimit(), and G4VEmAdjointModel::SetLowEnergyLimit().

fDirectPrimaryPart

G4ParticleDefinition* G4VEmAdjointModel::fDirectPrimaryPart = nullptr

protectedinherited

Definition at line 261 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4AdjointIonIonisationModel::CorrectPostStepWeight(), G4AdjointhIonisationModel::DefineProjectileProperty(), G4AdjointIonIonisationModel::DefineProjectileProperty(), G4VEmAdjointModel::DiffCrossSectionPerAtomPrimToSecond(), G4AdjointhIonisationModel::DiffCrossSectionPerAtomPrimToSecond(), G4AdjointIonIonisationModel::DiffCrossSectionPerAtomPrimToSecond(), G4VEmAdjointModel::DiffCrossSectionPerVolumePrimToSecond(), G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4AdjointPhotoElectricModel(), G4AdjointBremsstrahlungModel::Initialize(), G4VEmAdjointModel::SetAdjointEquivalentOfDirectPrimaryParticleDefinition(), and G4AdjointIonIonisationModel::SetIon().

fElementCSProdToProj

std::vector<G4double> G4VEmAdjointModel::fElementCSProdToProj

protectedinherited

Definition at line 268 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), and G4VEmAdjointModel::SelectCSMatrix().

fElementCSScatProjToProj

std::vector<G4double> G4VEmAdjointModel::fElementCSScatProjToProj

protectedinherited

Definition at line 267 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), and G4VEmAdjointModel::SelectCSMatrix().

fFactorCSBiasing

G4double G4AdjointPhotoElectricModel::fFactorCSBiasing = 1.

private

Definition at line 86 of file G4AdjointPhotoElectricModel.hh.

Referenced by AdjointCrossSection(), and CorrectPostStepWeight().

fHighEnergyLimit

G4double G4VEmAdjointModel::fHighEnergyLimit = 0.

protectedinherited

Definition at line 283 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::GetHighEnergyLimit(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForProdToProj(), and G4VEmAdjointModel::SetHighEnergyLimit().

fIndexElement

size_t G4AdjointPhotoElectricModel::fIndexElement = 0

private

Definition at line 91 of file G4AdjointPhotoElectricModel.hh.

Referenced by AdjointCrossSection(), AdjointCrossSectionPerAtom(), and SampleSecondaries().

fInModelWeightCorr

G4bool G4VEmAdjointModel::fInModelWeightCorr

protectedinherited

Initial value:

=
    false

Definition at line 301 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), and G4VEmAdjointModel::SetCorrectWeightForPostStepInModel().

fKinEnergyProdForIntegration

G4double G4VEmAdjointModel::fKinEnergyProdForIntegration = 0.

protectedinherited

Definition at line 270 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForSecond(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForSecond(), and G4VEmAdjointModel::DiffCrossSectionFunction1().

fKinEnergyScatProjForIntegration

G4double G4VEmAdjointModel::fKinEnergyScatProjForIntegration = 0.

protectedinherited

Definition at line 271 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForScatProj(), and G4VEmAdjointModel::DiffCrossSectionFunction2().

fLastAdjointCSForProdToProj

G4double G4VEmAdjointModel::fLastAdjointCSForProdToProj = 0.

protectedinherited

Definition at line 275 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4VEmAdjointModel::CorrectPostStepWeight(), G4AdjointeIonisationModel::SampleSecondaries(), and G4VEmAdjointModel::SelectCSMatrix().

fLastAdjointCSForScatProjToProj

G4double G4VEmAdjointModel::fLastAdjointCSForScatProjToProj = 0.

protectedinherited

Definition at line 274 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4VEmAdjointModel::CorrectPostStepWeight(), G4AdjointeIonisationModel::SampleSecondaries(), and G4VEmAdjointModel::SelectCSMatrix().

fLastCS

G4double G4VEmAdjointModel::fLastCS = 0.

protectedinherited

Definition at line 273 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4AdjointComptonModel::AdjointCrossSection(), G4AdjointhIonisationModel::AdjointCrossSection(), G4VEmAdjointModel::CorrectPostStepWeight(), G4AdjointhIonisationModel::RapidSampleSecondaries(), G4AdjointeIonisationModel::SampleSecondaries(), and G4VEmAdjointModel::SelectCSMatrix().

fLowEnergyLimit

G4double G4VEmAdjointModel::fLowEnergyLimit = 0.

protectedinherited

Definition at line 284 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::GetLowEnergyLimit(), and G4VEmAdjointModel::SetLowEnergyLimit().

fName

const G4String G4VEmAdjointModel::fName

protectedinherited

Definition at line 252 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::GetName().

fOneMatrixForAllElements

G4bool G4VEmAdjointModel::fOneMatrixForAllElements = false

protectedinherited

Definition at line 309 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4VEmAdjointModel::GetUseOnlyOneMatrixForAllElements(), G4VEmAdjointModel::SelectCSMatrix(), and G4VEmAdjointModel::SetUseOnlyOneMatrixForAllElements().

fOutsideWeightFactor

G4double G4VEmAdjointModel::fOutsideWeightFactor = 1.

protectedinherited

Definition at line 292 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), and G4VEmAdjointModel::SetAdditionalWeightCorrectionFactorForPostStepOutsideModel().

fPostStepAdjointCS

G4double G4AdjointPhotoElectricModel::fPostStepAdjointCS = 0.

private

Definition at line 88 of file G4AdjointPhotoElectricModel.hh.

Referenced by CorrectPostStepWeight(), and SampleSecondaries().

fPreStepAdjointCS

G4double G4AdjointPhotoElectricModel::fPreStepAdjointCS = 0.

private

Definition at line 87 of file G4AdjointPhotoElectricModel.hh.

Referenced by CorrectPostStepWeight(), and SampleSecondaries().

fPreStepEnergy

G4double G4VEmAdjointModel::fPreStepEnergy = 0.

protectedinherited

Definition at line 277 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), and G4VEmAdjointModel::CorrectPostStepWeight().

fSecondPartSameType

G4bool G4VEmAdjointModel::fSecondPartSameType = false

protectedinherited

Definition at line 300 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointComptonModel::G4AdjointComptonModel(), G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4AdjointPhotoElectricModel(), G4VEmAdjointModel::GetSecondAdjEnergyMaxForScatProjToProj(), G4VEmAdjointModel::GetSecondAdjEnergyMinForProdToProj(), G4VEmAdjointModel::GetSecondPartOfSameType(), G4AdjointBremsstrahlungModel::Initialize(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), and G4VEmAdjointModel::SetSecondPartOfSameType().

fSelectedMaterial

G4Material* G4VEmAdjointModel::fSelectedMaterial = nullptr

protectedinherited

Definition at line 254 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerVolumeForSecond(), G4VEmAdjointModel::DiffCrossSectionFunction1(), and G4VEmAdjointModel::DiffCrossSectionFunction2().

fShellProb

G4double G4AdjointPhotoElectricModel::fShellProb[40][40]

private

Definition at line 83 of file G4AdjointPhotoElectricModel.hh.

Referenced by AdjointCrossSectionPerAtom(), and SampleSecondaries().

fTcutPrim

G4double G4VEmAdjointModel::fTcutPrim = 0.

protectedinherited

Definition at line 279 of file G4VEmAdjointModel.hh.

fTcutSecond

G4double G4VEmAdjointModel::fTcutSecond = 0.

protectedinherited

Definition at line 280 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::AdjointCrossSection(), G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4AdjointhIonisationModel::AdjointCrossSection(), G4AdjointIonIonisationModel::CorrectPostStepWeight(), G4VEmAdjointModel::DefineCurrentMaterial(), G4AdjointBremsstrahlungModel::RapidSampleSecondaries(), G4AdjointComptonModel::RapidSampleSecondaries(), G4AdjointhIonisationModel::RapidSampleSecondaries(), G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix(), G4VEmAdjointModel::SampleAdjSecEnergyFromDiffCrossSectionPerAtom(), and G4AdjointeIonisationModel::SampleSecondaries().

fTotAdjointCS

G4double G4AdjointPhotoElectricModel::fTotAdjointCS = 0.

private

Definition at line 85 of file G4AdjointPhotoElectricModel.hh.

Referenced by AdjointCrossSection(), and SampleSecondaries().

fUseMatrix

G4bool G4VEmAdjointModel::fUseMatrix = false

protectedinherited

Definition at line 307 of file G4VEmAdjointModel.hh.

Referenced by G4AdjointBremsstrahlungModel::AdjointCrossSection(), G4AdjointComptonModel::AdjointCrossSection(), G4AdjointhIonisationModel::AdjointCrossSection(), G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4VEmAdjointModel::GetUseMatrix(), G4AdjointBremsstrahlungModel::SampleSecondaries(), G4AdjointComptonModel::SampleSecondaries(), G4AdjointhIonisationModel::SampleSecondaries(), and G4VEmAdjointModel::SetUseMatrix().

fUseMatrixPerElement

G4bool G4VEmAdjointModel::fUseMatrixPerElement = false

protectedinherited

Definition at line 308 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::DiffCrossSectionFunction1(), G4VEmAdjointModel::DiffCrossSectionFunction2(), G4AdjointeIonisationModel::G4AdjointeIonisationModel(), G4AdjointhIonisationModel::G4AdjointhIonisationModel(), G4AdjointIonIonisationModel::G4AdjointIonIonisationModel(), G4VEmAdjointModel::GetUseMatrixPerElement(), G4VEmAdjointModel::SelectCSMatrix(), and G4VEmAdjointModel::SetUseMatrixPerElement().

fXsec

G4double G4AdjointPhotoElectricModel::fXsec[40]

private

Definition at line 84 of file G4AdjointPhotoElectricModel.hh.

Referenced by AdjointCrossSection(), and SampleSecondaries().

fZSelectedNucleus

G4int G4VEmAdjointModel::fZSelectedNucleus = 0

protectedinherited

Definition at line 296 of file G4VEmAdjointModel.hh.

Referenced by G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForScatProj(), G4VEmAdjointModel::ComputeAdjointCrossSectionVectorPerAtomForSecond(), G4VEmAdjointModel::DiffCrossSectionFunction1(), and G4VEmAdjointModel::DiffCrossSectionFunction2().

---

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

• source/processes/electromagnetic/adjoint/include/G4AdjointPhotoElectricModel.hh
• source/processes/electromagnetic/adjoint/src/G4AdjointPhotoElectricModel.cc