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=nullptr)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4int	GetVerboseLevel () const
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
bool	ForAllAtomsAndEnergies () const
void	SetForAllAtomsAndEnergies (G4bool val)
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 = G4Pow::GetInstance()->powA(a,-0.3333333333);
  G4int nOfNeutrons = G4lrint(a) - Z;
  kineticEnergy /=GeV;
  G4double alog10E = std::log10(kineticEnergy);
 
  static const G4double nuleonRadius=1.36e-15;
  static 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=G4Log((G4double(nOfNeutrons))); }
  G4double crossSection = 1.0E31*fac*fac2*(1. + 1./a13 - fac1);
 
  // high energy correction
  crossSection *= (1 - 0.15*G4Exp(-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+G4Exp(-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. + G4Exp(ff4));
  return crossSection;
}

Referenced by 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