G4PhotoNuclearCrossSection Class Reference

#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 ()

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

(

)

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

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

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

static const char* G4PhotoNuclearCrossSection::Default_Name ( )

inlinestatic

Definition at line 49 of file G4PhotoNuclearCrossSection.hh.

{return "PhotoNuclearXS";}

Here is the caller graph for this function:

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:

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1528 of file G4PhotoNuclearCrossSection.cc.

{
    return true;
}

---

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

• geant4.10.03.p01/source/processes/hadronic/cross_sections/include/G4PhotoNuclearCrossSection.hh
• geant4.10.03.p01/source/processes/hadronic/cross_sections/src/G4PhotoNuclearCrossSection.cc