#include <G4PhotoNuclearCrossSection.hh>

Inheritance diagram for G4PhotoNuclearCrossSection:

Collaboration diagram for G4PhotoNuclearCrossSection:

Public Member Functions
	G4PhotoNuclearCrossSection ()

virtual	~G4PhotoNuclearCrossSection ()

virtual void	CrossSectionDescription (std::ostream &) const

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

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

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

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

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

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

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

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

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

const G4String &	GetName () const

Static Public Member Functions
static const char *	Default_Name ()

Private Member Functions
G4int	GetFunctions (G4double a, G4double y, G4double z)

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

G4double	ThresholdEnergy (G4int Z, G4int N)

Private Attributes
G4int	lastZ

G4double	lastSig

G4double *	lastGDR

G4double *	lastHEN

G4double	lastE

G4double	lastTH

G4double	lastSP

std::vector< G4double * >	GDR

std::vector< G4double * >	HEN

std::vector< G4double >	spA

std::vector< G4double >	eTH

G4NistManager *	nistmngr

G4double	mNeut

G4double	mProt

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 42 of file G4PhotoNuclearCrossSection.hh.

Constructor & Destructor Documentation

◆ G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::G4PhotoNuclearCrossSection ( )

Definition at line 1492 of file G4PhotoNuclearCrossSection.cc.

 : G4VCrossSectionDataSet(Default_Name()), lastZ(0), lastSig(0), lastGDR(0), lastHEN(0), lastE(0), lastTH(0), lastSP(0),
 mNeut(G4NucleiProperties::GetNuclearMass(1,0)), mProt(G4NucleiProperties::GetNuclearMass(1,1))
 {
     nistmngr = G4NistManager::Instance();
     
     for (G4int i=0;i<120;i++)
     {
         GDR.push_back(0);
         HEN.push_back(0);
         spA.push_back(0);
         eTH.push_back(0);
     }
 }

Here is the call graph for this function:

◆ ~G4PhotoNuclearCrossSection()

G4PhotoNuclearCrossSection::~G4PhotoNuclearCrossSection ( )

virtual

Definition at line 1507 of file G4PhotoNuclearCrossSection.cc.

 {
     std::vector<G4double*>::iterator posi;
     for(posi=GDR.begin(); posi<GDR.end(); posi++)
     { delete [] *posi; }
     GDR.clear();
     for(posi=HEN.begin(); posi<HEN.end(); posi++)
     { delete [] *posi; }
     HEN.clear();
 }

Member Function Documentation

◆ CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1519 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";
 }

Here is the caller graph for this function:

◆ Default_Name()

static const char* G4PhotoNuclearCrossSection::Default_Name ( )

inlinestatic

Definition at line 49 of file G4PhotoNuclearCrossSection.hh.

49 {return "PhotoNuclearXS";}

Here is the call graph for this function:

Here is the caller graph for this function:

◆ EquLinearFit()

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

private

Definition at line 1646 of file G4PhotoNuclearCrossSection.cc.

 {
     if(DX<=0. || N<2)
     {
         G4cout<<"***G4PhotoNuclearCrossSection::EquLinearFit: DX="<<DX<<", N="<<N<<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;
 }

Here is the caller graph for this function:

◆ GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1537 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 = (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
     
     if(sigma<0.) return 0.;
     return sigma*millibarn;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetFunctions()

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

private

Definition at line 1669 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(!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(!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;
 }

Here is the caller graph for this function:

◆ IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1528 of file G4PhotoNuclearCrossSection.cc.

 {
     return true;
 }

Here is the caller graph for this function:

◆ ThresholdEnergy()

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

private

Definition at line 1610 of file G4PhotoNuclearCrossSection.cc.

 {
     // ---------
     
     G4int A=Z+N;
     if(A<1) return infEn;
     else if(A==1) return 134.9766; // Pi0 threshold for the nucleon
     
     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;
 }