#include <G4PhotoNuclearCrossSection.hh>

Inheritance 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.

References G4NistManager::Instance().

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

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.

Referenced by G4ElectroVDNuclearModel::G4ElectroVDNuclearModel().

49 {return "PhotoNuclearXS";}

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 1537 of file G4PhotoNuclearCrossSection.cc.

References G4cerr, G4endl, G4NistManager::GetAtomicMassAmu(), G4DynamicParticle::GetKineticEnergy(), python.hepunit::MeV, python.hepunit::millibarn, and N.

 {
     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=std::log(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=std::log(Energy);
         sigma=EquLinearFit(lE,nH,milE,dlE,lastHEN);
     }
     else               // UHE region (calculation, but not so frequent)
     {
         G4double lE=std::log(Energy);
         sigma=lastSP*(poc*(lE-pos)+shd*std::exp(-reg*lE));
     }
     // End of "sigma" calculation
     
     if(sigma<0.) return 0.;
     return sigma*millibarn;
 }

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:

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation