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)

---

Detailed Description

Definition at line 52 of file G4ProtonInelasticCrossSection.hh.

---

Constructor & Destructor Documentation

G4ProtonInelasticCrossSection::G4ProtonInelasticCrossSection

(

)

Definition at line 41 of file G4ProtonInelasticCrossSection.cc.

References G4NistManager::Instance(), and G4Proton::Proton().

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

G4ProtonInelasticCrossSection::~G4ProtonInelasticCrossSection

(

)

Definition at line 48 of file G4ProtonInelasticCrossSection.cc.

{}

---

Member Function Documentation

G4double G4ProtonInelasticCrossSection::GetElementCrossSection	(	const G4DynamicParticle *	,
		G4int	Z,
		const G4Material *
	)			[virtual]

Reimplemented from G4VCrossSectionDataSet.

Definition at line 59 of file G4ProtonInelasticCrossSection.cc.

References G4DynamicParticle::GetKineticEnergy(), and GetProtonCrossSection().

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

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

Definition at line 66 of file G4ProtonInelasticCrossSection.cc.

References G4lrint(), G4NistManager::GetAtomicMassAmu(), and G4INCL::Math::pi.

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

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

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

Reimplemented from G4VCrossSectionDataSet.

Definition at line 52 of file G4ProtonInelasticCrossSection.cc.

References G4DynamicParticle::GetDefinition().

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

---

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

• G4ProtonInelasticCrossSection.hh
• G4ProtonInelasticCrossSection.cc

---