#include <G4PhotoNuclearCrossSection.hh>

Inheritance diagram for G4PhotoNuclearCrossSection:

Public Member Functions
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

virtual void	CrossSectionDescription (std::ostream &) const

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

bool	ForAllAtomsAndEnergies () const

	G4PhotoNuclearCrossSection ()

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material )

virtual G4double	GetIsoCrossSection (const G4DynamicParticle dynPart, G4int Z, G4int A, const G4Isotope , const G4Element elm, const G4Material mat)

G4double	GetMaxKinEnergy () const

G4double	GetMinKinEnergy () const

const G4String &	GetName () const

virtual G4int	GetVerboseLevel () const

virtual G4bool	IsElementApplicable (const G4DynamicParticle particle, G4int Z, const G4Material )

virtual G4bool	IsIsoApplicable (const G4DynamicParticle particle, G4int Z, G4int A, const G4Element elm=0, const G4Material *mat=0)

virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)

void	SetForAllAtomsAndEnergies (G4bool val)

void	SetMaxKinEnergy (G4double value)

void	SetMinKinEnergy (G4double value)

void	SetName (const G4String &nam)

virtual void	SetVerboseLevel (G4int value)

virtual	~G4PhotoNuclearCrossSection ()

Static Public Member Functions
static const char *	Default_Name ()

Protected Attributes
G4String	name

G4int	verboseLevel

Private Member Functions
G4double	EquLinearFit (G4double X, G4int N, const G4double X0, const G4double XD, const G4double *Y)

G4int	GetFunctions (G4double a, G4double y, G4double z)

G4double	ThresholdEnergy (G4int Z, G4int N)

Private Attributes
G4double *	deuteron_GDR = nullptr

G4double *	deuteron_HR = nullptr

G4double	deuteron_SP = 0.0

G4double	deuteron_TH = 0.0

std::vector< G4double >	eTH

std::vector< G4double * >	GDR

G4double *	he3_GDR = nullptr

G4double *	he3_HR = nullptr

G4double	he3_SP = 0.0

G4double	he3_TH = 0.0

std::vector< G4double * >	HEN

G4bool	isForAllAtomsAndEnergies

G4double	lastE = 0.0

G4double *	lastGDR = nullptr

G4double *	lastHEN = nullptr

G4double	lastSig = 0.0

G4double	lastSP = 0.0

G4double	lastTH = 0.0

G4int	lastZ = 0

G4double	maxKinEnergy

G4double	minKinEnergy

G4double	mNeut

G4double	mProt

G4NistManager *	nistmngr

G4CrossSectionDataSetRegistry *	registry

std::vector< G4double >	spA

G4double *	triton_GDR = nullptr

G4double *	triton_HR = nullptr

G4double	triton_SP = 0.0

G4double	triton_TH = 0.0

Detailed Description

Definition at line 42 of file G4PhotoNuclearCrossSection.hh.

Constructor & Destructor Documentation

◆ G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::G4PhotoNuclearCrossSection ( )

Definition at line 1486 of file G4PhotoNuclearCrossSection.cc.

 : G4VCrossSectionDataSet(Default_Name()), 
   mNeut(G4NucleiProperties::GetNuclearMass(1,0)),
   mProt(G4NucleiProperties::GetNuclearMass(1,1))
{
  SetForAllAtomsAndEnergies(true);
  nistmngr = G4NistManager::Instance();
    
  GDR.resize(120, nullptr);
  HEN.resize(120, nullptr);
  spA.resize(120, 0.0);
  eTH.resize(120, 0.0);
}

References eTH, GDR, HEN, G4NistManager::Instance(), nistmngr, G4VCrossSectionDataSet::SetForAllAtomsAndEnergies(), and spA.

◆ ~G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::~G4PhotoNuclearCrossSection ( )

virtual

Definition at line 1501 of file G4PhotoNuclearCrossSection.cc.

{
  for (auto & ptr : GDR) { delete [] ptr; }
  for (auto & ptr : HEN) { delete [] ptr; }
}

References GDR, and HEN.

Member Function Documentation

◆ BuildPhysicsTable()

void G4VCrossSectionDataSet::BuildPhysicsTable ( const G4ParticleDefinition & )

virtualinherited

◆ ComputeCrossSection()

G4double G4VCrossSectionDataSet::ComputeCrossSection	(	const G4DynamicParticle *	part,
		const G4Element *	elm,
		const G4Material *	mat = `nullptr`
	)

inherited

Definition at line 81 of file G4VCrossSectionDataSet.cc.

{
  G4int Z = elm->GetZasInt();
 
  if (IsElementApplicable(part, Z, mat)) { 
    return GetElementCrossSection(part, Z, mat);
  }
 
  // isotope-wise cross section making sum over available
  // isotope cross sections, which may be incomplete, so
  // the result is corrected 
  size_t nIso = elm->GetNumberOfIsotopes();    
  G4double fact = 0.0;
  G4double xsec = 0.0;
 
  // user-defined isotope abundances        
  const G4IsotopeVector* isoVector = elm->GetIsotopeVector();
  const G4double* abundVector = elm->GetRelativeAbundanceVector();
 
  for (size_t j=0; j<nIso; ++j) {
    const G4Isotope* iso = (*isoVector)[j];
    G4int A = iso->GetN();
    if(abundVector[j] > 0.0 && IsIsoApplicable(part, Z, A, elm, mat)) {
      fact += abundVector[j];
      xsec += abundVector[j]*GetIsoCrossSection(part, Z, A, iso, elm, mat);
    }
  }
  return (fact > 0.0) ? xsec/fact : 0.0; 
}

References A, G4VCrossSectionDataSet::GetElementCrossSection(), G4VCrossSectionDataSet::GetIsoCrossSection(), G4Element::GetIsotopeVector(), G4Isotope::GetN(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZasInt(), G4VCrossSectionDataSet::IsElementApplicable(), G4VCrossSectionDataSet::IsIsoApplicable(), and Z.

Referenced by G4VCrossSectionDataSet::GetCrossSection().

◆ CrossSectionDescription()

void G4PhotoNuclearCrossSection::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1509 of file G4PhotoNuclearCrossSection.cc.

{
  outFile << "G4PhotoNuclearCrossSection provides the total inelastic\n"
          << "cross section for photon interactions with nuclei.  The\n"
          << "cross section is a parameterization of data which covers\n"
          << "all incident gamma energies.\n";
}

◆ Default_Name()

static const char * G4PhotoNuclearCrossSection::Default_Name ( )

inlinestatic

Definition at line 49 of file G4PhotoNuclearCrossSection.hh.

49{return "PhotoNuclearXS";}

Referenced by G4ElectroVDNuclearModel::G4ElectroVDNuclearModel().

◆ DumpPhysicsTable()

void G4VCrossSectionDataSet::DumpPhysicsTable ( const G4ParticleDefinition & )

virtualinherited

Reimplemented in G4CrossSectionElastic, G4CrossSectionInelastic, G4LENDCrossSection, G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, G4ParticleHPInelasticData, G4ParticleHPJENDLHEData, G4ParticleHPThermalScatteringData, and G4UPiNuclearCrossSection.

Definition at line 174 of file G4VCrossSectionDataSet.cc.

175{}

◆ EquLinearFit()

G4double G4PhotoNuclearCrossSection::EquLinearFit	(	G4double	X,
		G4int	N,
		const G4double	X0,
		const G4double	XD,
		const G4double *	Y
	)

private

Definition at line 1740 of file G4PhotoNuclearCrossSection.cc.

{
  if(DX<=0. || N<2)
    {
      G4cout<<"***G4PhotoNuclearCrossSection::EquLinearFit: DX="
        <<DX<<", N="<<N<<", X0="<<X0<<", X="<<X<<", Y[0]="<<Y[0]<<G4endl;
      return Y[0];
    }
  G4int    N2=N-2;
  G4double d=(X-X0)/DX;
  G4int         j=static_cast<int>(d);
  if     (j<0)  j=0;
  else if(j>N2) j=N2;
  d-=j; // excess
  G4double yi=Y[j];
  G4double sigma=yi+(Y[j+1]-yi)*d;
  return sigma;
}

References G4cout, G4endl, and Y().

Referenced by GetElementCrossSection(), and GetIsoCrossSection().

◆ ForAllAtomsAndEnergies()

bool G4VCrossSectionDataSet::ForAllAtomsAndEnergies ( ) const

inlineinherited

Definition at line 230 of file G4VCrossSectionDataSet.hh.

{
  return isForAllAtomsAndEnergies;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4CrossSectionDataStore::AddDataSet().

◆ GetCrossSection()

G4double G4VCrossSectionDataSet::GetCrossSection	(	const G4DynamicParticle *	dp,
		const G4Element *	elm,
		const G4Material *	mat = `nullptr`
	)

inlineinherited

Definition at line 187 of file G4VCrossSectionDataSet.hh.

{
  return ComputeCrossSection(dp, elm, mat);
}

References G4VCrossSectionDataSet::ComputeCrossSection().

◆ GetElementCrossSection()

G4double G4PhotoNuclearCrossSection::GetElementCrossSection	(	const G4DynamicParticle *	aPart,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1629 of file G4PhotoNuclearCrossSection.cc.

{
    const G4double Energy = aPart->GetKineticEnergy()/MeV;
    
    if (Energy<THmin) return 0.;
    G4double sigma=0.;
    
    if(ZZ!=lastZ)   // Otherwise the set of parameters is ready
    {
        lastZ = ZZ;             // The last Z of calculated nucleus
        if(GDR[ZZ]) // Calculated nuclei in DB
        {
            lastGDR=GDR[ZZ];   // Pointer to prepared GDR cross sections
            lastHEN=HEN[ZZ];   // Pointer to prepared High Energy cross sections
            lastTH =eTH[ZZ];   // Energy Threshold
            lastSP =spA[ZZ];   // Shadowing coefficient for UHE
        }
        else
        {
            G4double Aa = nistmngr->GetAtomicMassAmu(ZZ); // average A
            G4int N = static_cast<G4int>(Aa) - ZZ;
            
            G4double lnA=G4Log(Aa); // The nucleus is not found in DB. It is new.
            if(Aa==1.) lastSP=1.;      // The Reggeon shadowing (A=1)
            else lastSP=Aa*(1.-shc*lnA);      // The Reggeon shadowing
            lastTH=ThresholdEnergy(ZZ, N); // Energy Threshold
            lastGDR = new G4double[nL];  // Allocate memory for the new
            // GDR cross sections
            lastHEN = new G4double[nH];  // Allocate memory for the new
            // HEN cross sections
            G4int er=GetFunctions(Aa,lastGDR,lastHEN); // new ZeroPosition and
            // filling of the functions
            if(er<1) G4cerr << "***G4PhotoNucCrossSection::GetCrossSection: A="
                << Aa << " failed" << G4endl;
            
            GDR[ZZ]=lastGDR;     // added GDR, found by AH 10/7/02
            HEN[ZZ]=lastHEN;     // added HEN, found by AH 10/7/02
            eTH[ZZ]=lastTH;      // Threshold Energy
            spA[ZZ]=lastSP;      // Pomeron Shadowing
        } // End of creation of the new set of parameters
    }// End of parameters udate
    
    //
    // =================== NOW the Magic Formula ===================
    //
    if (Energy<lastTH)
    {
        lastE=Energy;
        lastSig=0.;
        return 0.;
    }
    else if (Energy<Emin)   // GDR region (approximated in E, not in lnE)
    {
        sigma=EquLinearFit(Energy,nL,THmin,dE,lastGDR);
    }
    else if (Energy<Emax)                     // High Energy region
    {
        G4double lE=G4Log(Energy);
        sigma=EquLinearFit(lE,nH,milE,dlE,lastHEN);
    }
    else               // UHE region (calculation, but not so frequent)
    {
        G4double lE=G4Log(Energy);
        sigma=lastSP*(poc*(lE-pos)+shd*G4Exp(-reg*lE));
    }
    // End of "sigma" calculation
    
    sigma = std::max(sigma, 0.);
    return sigma*millibarn;
}

References dE, dlE, Emax, Emin, EquLinearFit(), eTH, G4cerr, G4endl, G4Exp(), G4Log(), GDR, G4NistManager::GetAtomicMassAmu(), GetFunctions(), G4DynamicParticle::GetKineticEnergy(), HEN, lastE, lastGDR, lastHEN, lastSig, lastSP, lastTH, lastZ, G4INCL::Math::max(), MeV, milE, millibarn, nH, nistmngr, nL, poc, pos, reg, shc, shd, spA, THmin, and ThresholdEnergy().

Referenced by G4ElectroVDNuclearModel::CalculateEMVertex(), and GetIsoCrossSection().

◆ GetFunctions()

G4int G4PhotoNuclearCrossSection::GetFunctions	(	G4double	a,
		G4double *	y,
		G4double *	z
	)

private

Definition at line 1764 of file G4PhotoNuclearCrossSection.cc.

{    
  if(a<=.9)
    {
        G4cout << "***G4PhotoNuclearCS::GetFunctions: A=" << a
               << "(?). No CS returned!" << G4endl;
        return -1;
    }
  G4int r=0;                    // Low channel for GDR (filling-flag for GDR)
  for(G4int i=0; i<nLA; ++i) if(std::abs(a-LA[i])<.0005)
    {
        for(G4int k=0; k<nL; ++k) y[k]=SL[i][k];
        r=1;                          // Flag of filled GDR part
    }
  G4int h=0;
  for(G4int j=0; j<nHA; ++j) if(std::abs(a-HA[j])<.0005)
    {
        for(G4int k=0; k<nH; ++k) z[k]=SH[j][k];
        h=1;                          // Flag of filled GDR part
    }
  if(0 == r)           // GDR part is not filled
    {
        G4int k=0;   // !! To be good for different compilers !!
        for(k=1; k<nLA; ++k) if(a<LA[k]) break;
        if(k<1) k=1;                  // Extrapolation from the first bin (D/He)
        if(k>=nLA) k=nLA-1;           // Extrapolation from the last bin (U)
        G4int     k1=k-1;
        G4double  xi=LA[k1];
        G4double   b=(a-xi)/(LA[k]-xi);
        for(G4int q=0; q<nL; ++q)
        {
            if(a>1.5)
            {
                G4double yi=SL[k1][q];
                y[q]=yi+(SL[k][q]-yi)*b;
            }
            else y[q]=0.;
        }
        r=1;
    }
  if(0 == h)    // High Energy part is not filled
    {
        G4int k=0;
        for(k=1; k<nHA; ++k) if(a<HA[k]) break;
        if(k<1) k=1;         // Extrapolation from the first bin (D/He)
        if(k>=nHA) k=nHA-1;  // Extrapolation from the last bin (Pu)
        G4int     k1=k-1;
        G4double  xi=HA[k1];
        G4double   b=(a-xi)/(HA[k]-xi);
        for(G4int q=0; q<nH; ++q)
        {
            G4double zi=SH[k1][q];
            z[q]=zi+(SH[k][q]-zi)*b;
        }
        h=1;
    }
  return r*h;
}

References G4cout, G4endl, HA, LA, nH, nHA, nL, nLA, SH, and SL.

Referenced by GetElementCrossSection().

◆ GetIsoCrossSection()

G4double G4PhotoNuclearCrossSection::GetIsoCrossSection	(	const G4DynamicParticle *	dynPart,
		G4int	Z,
		G4int	A,
		const G4Isotope *	,
		const G4Element *	elm,
		const G4Material *	mat
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1542 of file G4PhotoNuclearCrossSection.cc.

{
  G4double Energy = aPart->GetKineticEnergy()/MeV;
  if (Energy < THmin) return 0.;
 
  G4double sigma;
  G4double lE;
  if (Z == 1 && A == 2) { 
    // init the XS table if need be
    if (deuteron_GDR == nullptr) {
      deuteron_TH = ThresholdEnergy(1,1); // threshold calculation is correct
      deuteron_GDR = new G4double[nL];    // direct copies
      for (G4int i = 0; i < nL; ++i) deuteron_GDR[i] = SL[0][i];   // A = 2 -> SL0
      deuteron_HR = new G4double[nH];
      for (G4int i = 0; i < nH; ++i) deuteron_HR[i] = SH[1][i];    // A = 2 -> SH1
      deuteron_SP = 1;                    // as would be assigned for proton
    }
    if (Energy < deuteron_TH) {
      sigma = 0.;
    } else if (Energy < Emin) {      // GDR region (approximated in E, not in lnE)
      sigma = EquLinearFit(Energy,nL,THmin,dE,deuteron_GDR);
    } else if (Energy < Emax) {      // High Energy region
      lE = G4Log(Energy);
      sigma = EquLinearFit(lE,nH,milE,dlE,deuteron_HR);
    } else {                         // Very high energy region
      lE = G4Log(Energy);
      sigma = deuteron_SP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
    }
 
  } else if (Z == 1 && A == 3) {
    if (triton_GDR == nullptr) {
      triton_TH = ThresholdEnergy(1,2);
      triton_GDR = new G4double[nL];      // same as for deuteron since no A = 3 entry
      for (G4int i = 0; i < nL; ++i) triton_GDR[i] = SL[0][i];     // A = 3 -> SL0
      triton_HR = new G4double[nH];
      for (G4int i = 0; i < nH; ++i) triton_HR[i] = SH[2][i];      // A = 3 -> SH2
      triton_SP = 1;
    }
    if (Energy < triton_TH) {
      sigma = 0.;
    } else if (Energy < Emin) {      // GDR region
      sigma = EquLinearFit(Energy,nL,THmin,dE,triton_GDR);
    } else if (Energy < Emax) {      // High Energy region
      lE = G4Log(Energy);
      sigma = EquLinearFit(lE,nH,milE,dlE,triton_HR);
    } else {
      lE = G4Log(Energy);
      sigma = triton_SP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
    }
 
  } else if (Z == 2 && A == 3) {
    if (he3_GDR == nullptr) {
      he3_TH = ThresholdEnergy(2,1);
      he3_GDR = new G4double[nL]; // same as for deuteron since no A = 3 entry
      for (G4int i = 0 ; i < nL ; ++i) he3_GDR[i] = SL[0][i]; // A = 3 -> SL0
      he3_HR = new G4double[nH];
      for (G4int i = 0 ; i < nH ; ++i) he3_HR[i] = SH[2][i];  // A = 3 -> SH2
      he3_SP = 2;
    }
    if (Energy < he3_TH) {
      sigma = 0.;
    } else if (Energy < Emin) {      // GDR region
      sigma = EquLinearFit(Energy,nL,THmin,dE,he3_GDR);
    } else if (Energy < Emax) {      // High Energy region
      lE = G4Log(Energy);
      sigma = EquLinearFit(lE,nH,milE,dlE,he3_HR);
    } else {
      lE = G4Log(Energy);
      sigma = he3_SP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
    }
 
  } else {
    return GetElementCrossSection(aPart, Z, mat);
  }
 
  if(sigma < 0.) sigma = 0.;
  return sigma*millibarn;
}

References A, dE, deuteron_GDR, deuteron_HR, deuteron_SP, deuteron_TH, dlE, Emax, Emin, EquLinearFit(), G4Log(), GetElementCrossSection(), G4DynamicParticle::GetKineticEnergy(), he3_GDR, he3_HR, he3_SP, he3_TH, MeV, milE, millibarn, nH, nL, poc, pos, reg, SH, shd, SL, THmin, ThresholdEnergy(), triton_GDR, triton_HR, triton_SP, triton_TH, and Z.

◆ GetMaxKinEnergy()

G4double G4VCrossSectionDataSet::GetMaxKinEnergy ( ) const

inlineinherited

Definition at line 220 of file G4VCrossSectionDataSet.hh.

{
  return maxKinEnergy;
}

References G4VCrossSectionDataSet::maxKinEnergy.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4LENDCrossSection::IsIsoApplicable(), G4ParticleHPCaptureData::IsIsoApplicable(), G4ParticleHPElasticData::IsIsoApplicable(), G4ParticleHPFissionData::IsIsoApplicable(), and G4ParticleHPInelasticData::IsIsoApplicable().

◆ GetMinKinEnergy()

G4double G4VCrossSectionDataSet::GetMinKinEnergy ( ) const

inlineinherited

Definition at line 210 of file G4VCrossSectionDataSet.hh.

{
  return minKinEnergy;
}

References G4VCrossSectionDataSet::minKinEnergy.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4LENDCrossSection::IsIsoApplicable(), G4ParticleHPCaptureData::IsIsoApplicable(), G4ParticleHPElasticData::IsIsoApplicable(), G4ParticleHPFissionData::IsIsoApplicable(), and G4ParticleHPInelasticData::IsIsoApplicable().

◆ GetName()

const G4String & G4VCrossSectionDataSet::GetName ( ) const

inlineinherited

Definition at line 225 of file G4VCrossSectionDataSet.hh.

{
  return name;
}

References G4VCrossSectionDataSet::name.

Referenced by G4LENDCrossSection::DumpLENDTargetInfo(), G4LENDCrossSection::DumpPhysicsTable(), G4LENDCrossSection::GetIsoCrossSection(), and G4CrossSectionDataStore::PrintCrossSectionHtml().

◆ GetVerboseLevel()

G4int G4VCrossSectionDataSet::GetVerboseLevel ( ) const

inlinevirtualinherited

Reimplemented in G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, and G4ParticleHPInelasticData.

Definition at line 195 of file G4VCrossSectionDataSet.hh.

{
  return verboseLevel;
}

References G4VCrossSectionDataSet::verboseLevel.

◆ IsElementApplicable()

G4bool G4PhotoNuclearCrossSection::IsElementApplicable	(	const G4DynamicParticle *	particle,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1533 of file G4PhotoNuclearCrossSection.cc.

{
  return true;
}

◆ IsIsoApplicable()

G4bool G4PhotoNuclearCrossSection::IsIsoApplicable	(	const G4DynamicParticle *	particle,
		G4int	Z,
		G4int	A,
		const G4Element *	elm = `0`,
		const G4Material *	mat = `0`
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1520 of file G4PhotoNuclearCrossSection.cc.

{
  // explicitly allow deuterium and tritium
  if ((Z == 1 && A == 2) || (Z == 1 && A == 3) ||
      (Z == 2 && A == 3) ) return true; 
  return false;
}

References A, and Z.

◆ SelectIsotope()

const G4Isotope * G4VCrossSectionDataSet::SelectIsotope	(	const G4Element *	anElement,
		G4double	kinEnergy,
		G4double	logE
	)

virtualinherited

Reimplemented in G4GammaNuclearXS, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticXS, and G4ParticleInelasticXS.

Definition at line 149 of file G4VCrossSectionDataSet.cc.

{
  size_t nIso = anElement->GetNumberOfIsotopes();
  const G4Isotope* iso = anElement->GetIsotope(0);
 
  // more than 1 isotope
  if(1 < nIso) {
    const G4double* abundVector = anElement->GetRelativeAbundanceVector();
    G4double sum = 0.0;
    G4double q = G4UniformRand();
    for (size_t j=0; j<nIso; ++j) {
      sum += abundVector[j];
      if(q <= sum) {
    iso = anElement->GetIsotope(j);
    break;
      }
    }
  }
  return iso;
}

References G4UniformRand, G4Element::GetIsotope(), G4Element::GetNumberOfIsotopes(), and G4Element::GetRelativeAbundanceVector().

◆ SetForAllAtomsAndEnergies()

void G4VCrossSectionDataSet::SetForAllAtomsAndEnergies ( G4bool val )

inlineinherited

Definition at line 235 of file G4VCrossSectionDataSet.hh.

{
  isForAllAtomsAndEnergies = val;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4BGGNucleonElasticXS::G4BGGNucleonElasticXS(), G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4ElectroNuclearCrossSection::G4ElectroNuclearCrossSection(), G4GammaNuclearXS::G4GammaNuclearXS(), G4NeutronElasticXS::G4NeutronElasticXS(), G4NeutronInelasticXS::G4NeutronInelasticXS(), G4ParticleInelasticXS::G4ParticleInelasticXS(), and G4PhotoNuclearCrossSection().

◆ SetMaxKinEnergy()

void G4VCrossSectionDataSet::SetMaxKinEnergy ( G4double value )

inlineinherited

Definition at line 215 of file G4VCrossSectionDataSet.hh.

{
  maxKinEnergy = value;
}

References G4VCrossSectionDataSet::maxKinEnergy.

Referenced by G4AlphaPHPBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4He3PHPBuilder::Build(), G4ProtonPHPBuilder::Build(), G4TritonPHPBuilder::Build(), G4CrossSectionElastic::BuildPhysicsTable(), G4CrossSectionInelastic::BuildPhysicsTable(), G4ParticleInelasticXS::BuildPhysicsTable(), G4IonPhysicsPHP::ConstructProcess(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionElastic::G4CrossSectionElastic(), G4CrossSectionInelastic::G4CrossSectionInelastic(), G4LENDCrossSection::G4LENDCrossSection(), G4ParticleHPCaptureData::G4ParticleHPCaptureData(), G4ParticleHPElasticData::G4ParticleHPElasticData(), G4ParticleHPFissionData::G4ParticleHPFissionData(), G4ParticleHPInelasticData::G4ParticleHPInelasticData(), and G4ParticleHPThermalScatteringData::G4ParticleHPThermalScatteringData().

◆ SetMinKinEnergy()

void G4VCrossSectionDataSet::SetMinKinEnergy ( G4double value )

inlineinherited

Definition at line 205 of file G4VCrossSectionDataSet.hh.

{
  minKinEnergy = value;
}

References G4VCrossSectionDataSet::minKinEnergy.

Referenced by G4AlphaPHPBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4He3PHPBuilder::Build(), G4ProtonPHPBuilder::Build(), G4TritonPHPBuilder::Build(), G4IonElasticPhysics::ConstructProcess(), G4IonPhysicsPHP::ConstructProcess(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4CrossSectionElastic::G4CrossSectionElastic(), G4CrossSectionInelastic::G4CrossSectionInelastic(), G4LENDCrossSection::G4LENDCrossSection(), G4ParticleHPCaptureData::G4ParticleHPCaptureData(), G4ParticleHPElasticData::G4ParticleHPElasticData(), G4ParticleHPFissionData::G4ParticleHPFissionData(), G4ParticleHPInelasticData::G4ParticleHPInelasticData(), and G4ParticleHPThermalScatteringData::G4ParticleHPThermalScatteringData().

◆ SetName()

void G4VCrossSectionDataSet::SetName ( const G4String & nam )

inlineinherited

Definition at line 240 of file G4VCrossSectionDataSet.hh.

{
  name = nam;
}

References G4VCrossSectionDataSet::name.

Referenced by G4ParticleHPInelasticData::G4ParticleHPInelasticData().

◆ SetVerboseLevel()

void G4VCrossSectionDataSet::SetVerboseLevel ( G4int value )

inlinevirtualinherited

Reimplemented in G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, and G4ParticleHPInelasticData.

Definition at line 200 of file G4VCrossSectionDataSet.hh.

{
  verboseLevel = value;
}

References G4VCrossSectionDataSet::verboseLevel.

Referenced by G4VCrossSectionRatio::G4VCrossSectionRatio().

◆ ThresholdEnergy()

G4double G4PhotoNuclearCrossSection::ThresholdEnergy	(	G4int	Z,
		G4int	N
	)

private

Definition at line 1703 of file G4PhotoNuclearCrossSection.cc.

{
  G4int A = Z + N;
  if (A < 1) return infEn;
 
  // Thresholds for p, d, t, 3He in lab frame
  else if (A == 1) return 144.6821;           // pi0 production
  else if (Z == 1 && N == 1) return 2.2263;   // disintegration
  else if (Z == 1 && N == 2) return 6.2650;   // n separation
  else if (Z == 2 && N == 1) return 5.4994;   // p separation
 
  G4double mT = 0.;
  if (G4NucleiProperties::IsInStableTable(A,Z) ) {
    mT = G4NucleiProperties::GetNuclearMass(A,Z);
  } else {
    return infEn;
  }
 
  G4double mP = infEn;    
  if (Z && G4NucleiProperties::IsInStableTable(A-1,Z-1) ) {
    mP = G4NucleiProperties::GetNuclearMass(A-1,Z-1);
  }
  G4double mN = infEn;
  if (N && G4NucleiProperties::IsInStableTable(A-1,Z) ) {
    mN = G4NucleiProperties::GetNuclearMass(A-1,Z);
  }
 
  G4double dP = mP + mProt - mT;
  G4double dN = mN + mNeut - mT;
  if (dP < dN) dN = dP;
  return dN;
}