G4IonsSihverCrossSection Class Reference

#include <G4IonsSihverCrossSection.hh>

Inheritance diagram for G4IonsSihverCrossSection:

Public Member Functions
	G4IonsSihverCrossSection ()
virtual	~G4IonsSihverCrossSection ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *aDP, G4int Z, const G4Material *)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *)
virtual void	CrossSectionDescription (std::ostream &) const
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

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 46 of file G4IonsSihverCrossSection.hh.

Constructor & Destructor Documentation

G4IonsSihverCrossSection::G4IonsSihverCrossSection

(

)

Definition at line 40 of file G4IonsSihverCrossSection.cc.

  : G4VCrossSectionDataSet("IonsSihver"), square_r0 ( (1.36*fermi) * (1.36*fermi) )
{}

G4IonsSihverCrossSection::~G4IonsSihverCrossSection ( )

virtual

Definition at line 44 of file G4IonsSihverCrossSection.cc.

{}

Member Function Documentation

void G4IonsSihverCrossSection::CrossSectionDescription ( std::ostream &  outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 48 of file G4IonsSihverCrossSection.cc.

{
  outFile << "G4IonsSihverCrossSection calculates the total reaction cross\n"
          << "section for nucleus-nucleus scattering using the Sihver\n"
          << "parameterization.  It is valid for projectiles and targets of\n"
          << "all Z, and for all projectile energies above 100 MeV/n.\n"; 
}

G4double G4IonsSihverCrossSection::GetElementCrossSection ( const G4DynamicParticle *  aParticle, G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 67 of file G4IonsSihverCrossSection.cc.

References G4lrint(), G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), G4Pow::GetInstance(), G4NistManager::Instance(), python.hepunit::pi, and G4Pow::Z13().

{
  G4double xsection = 0.0;
  G4int At = G4lrint(G4NistManager::Instance()->GetAtomicMassAmu(Z));

  G4int Ap = aParticle->GetDefinition()->GetBaryonNumber();
 
  G4Pow* g4pow = G4Pow::GetInstance();
  
  G4double cubicrAt = g4pow->Z13(At);
  G4double cubicrAp = g4pow->Z13(Ap);
 
  G4double b0 = 1.581 - 0.876 * (1.0/cubicrAp + 1.0/cubicrAt);

  G4double xr = cubicrAp + cubicrAt - b0 * (1.0/cubicrAp + 1.0/cubicrAt);
  xsection = pi * square_r0 * xr * xr;
  
  return xsection; 
}

G4bool G4IonsSihverCrossSection::IsElementApplicable ( const G4DynamicParticle *  aDP, G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 57 of file G4IonsSihverCrossSection.cc.

References G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), and python.hepunit::MeV.

{
  G4int BaryonNumber = aDP->GetDefinition()->GetBaryonNumber();
  G4double KineticEnergy = aDP->GetKineticEnergy(); 
  if ( KineticEnergy / BaryonNumber >= 100*MeV && BaryonNumber > 1 ) { return true; }
  return false;
}

---

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

• G4IonsSihverCrossSection.hh
• G4IonsSihverCrossSection.cc