G4ProtonInelasticCrossSection Class Reference

#include <G4ProtonInelasticCrossSection.hh>

Inheritance diagram for G4ProtonInelasticCrossSection:

Public Member Functions
	G4ProtonInelasticCrossSection ()
	~G4ProtonInelasticCrossSection ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *aPart, G4int Z, const G4Material *)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *)
G4double	GetProtonCrossSection (G4double kineticEnergy, G4int Z)
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
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	GetElementCrossSection (const G4DynamicParticle *, G4int Z, 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 void	CrossSectionDescription (std::ostream &) const
void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

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 52 of file G4ProtonInelasticCrossSection.hh.

Constructor & Destructor Documentation

G4ProtonInelasticCrossSection()

G4ProtonInelasticCrossSection::G4ProtonInelasticCrossSection

(

)

Definition at line 41 of file G4ProtonInelasticCrossSection.cc.

  : G4VCrossSectionDataSet("Axen-Wellisch"),thEnergy(19.8*CLHEP::GeV)
{
  nist = G4NistManager::Instance();
  theProton = G4Proton::Proton();
}

~G4ProtonInelasticCrossSection()

G4ProtonInelasticCrossSection::~G4ProtonInelasticCrossSection

(

)

Definition at line 48 of file G4ProtonInelasticCrossSection.cc.

{}

Member Function Documentation

GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 59 of file G4ProtonInelasticCrossSection.cc.

{
  return GetProtonCrossSection(aPart->GetKineticEnergy(), Z);
}

GetProtonCrossSection()

G4double G4ProtonInelasticCrossSection::GetProtonCrossSection	(	G4double	kineticEnergy,
		G4int	Z
	)

Definition at line 66 of file G4ProtonInelasticCrossSection.cc.

{
  if(kineticEnergy <= 0.0) { return 0.0; }
 
  // constant cross section above ~20GeV
  if (kineticEnergy > thEnergy) { kineticEnergy = thEnergy; }
 
  G4double a = nist->GetAtomicMassAmu(Z);
  G4double a13 = std::pow(a,-0.3333333333);
  G4int nOfNeutrons = G4lrint(a) - Z;
  kineticEnergy /=GeV;
  G4double alog10E = std::log10(kineticEnergy);
 
  const G4double nuleonRadius=1.36E-15;
  const G4double fac=CLHEP::pi*nuleonRadius*nuleonRadius;
 
  G4double b0   = 2.247-0.915*(1 - a13);
  G4double fac1 = b0*(1 - a13);
  G4double fac2 = 1.;
  if(nOfNeutrons > 1) { fac2=std::log((G4double(nOfNeutrons))); }
  G4double crossSection = 1.0E31*fac*fac2*(1. + 1./a13 - fac1);
 
  // high energy correction
  crossSection *= (1 - 0.15*std::exp(-kineticEnergy))/(1.0 - 0.0007*a);
 
  // first try on low energies: rise
  G4double ff1= 0.70-0.002*a;  // slope of the drop at medium energies.
  G4double ff2= 1.00+1/a;  // start of the slope.
  G4double ff3= 0.8+18/a-0.002*a; // stephight
 
  G4double ff4= 1.0 - (1.0/(1+std::exp(-8*ff1*(alog10E + 1.37*ff2))));
 
  crossSection *= (1 + ff3*ff4);
 
  // low energy return to zero
 
  ff1=1. - 1./a - 0.001*a;       // slope of the rise
  ff2=1.17 - 2.7/a - 0.0014*a;   // start of the rise
 
  ff4=-8.*ff1*(alog10E + 2.0*ff2);
 
  crossSection *= millibarn/(1. + std::exp(ff4));
  return crossSection;
}

Referenced by GetElementCrossSection(), and G4IonProtonCrossSection::GetElementCrossSection().

IsElementApplicable()

G4bool G4ProtonInelasticCrossSection::IsElementApplicable ( const G4DynamicParticle *  aPart, G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 52 of file G4ProtonInelasticCrossSection.cc.

{
  return ((1 < Z) && (aPart->GetDefinition() == theProton));
}

---

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

• G4ProtonInelasticCrossSection.hh
• G4ProtonInelasticCrossSection.cc