G4IonsKoxCrossSection Class Reference

#include <G4IonsKoxCrossSection.hh>

Inheritance diagram for G4IonsKoxCrossSection:

Public Member Functions
	G4IonsKoxCrossSection ()
	~G4IonsKoxCrossSection ()
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 G4IonsKoxCrossSection.hh.

Constructor & Destructor Documentation

G4IonsKoxCrossSection::G4IonsKoxCrossSection

(

)

Definition at line 41 of file G4IonsKoxCrossSection.cc.

  : G4VCrossSectionDataSet("IonsKox"),
// lowerLimit ( 10*MeV ), 
    r0 ( 1.1*fermi ), rc ( 1.3*fermi )
{}

G4IonsKoxCrossSection::~G4IonsKoxCrossSection

(

)

Definition at line 47 of file G4IonsKoxCrossSection.cc.

{}

Member Function Documentation

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 51 of file G4IonsKoxCrossSection.cc.

{
  outFile << "G4IonsKoxCrossSection calculates the total reaction cross\n"
          << "section for nucleus-nucleus scattering using the Kox\n"
          << "parameterization.  It is valid for projectiles and targets\n"
          << "of all Z, at projectile energies up to 10 GeV/n.  If the\n"
          << "projectile energy is less than 10 MeV/n, a zero cross section\n"
          << "is returned.\n";
}

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 68 of file G4IonsKoxCrossSection.cc.

References test::a, test::c, python.hepunit::eplus, python.hepunit::fermi, G4lrint(), G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMass(), G4DynamicParticle::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGCharge(), G4NistManager::Instance(), CLHEP::Hep3Vector::mag(), python.hepunit::MeV, and python.hepunit::pi.

{
   G4double xsection = 0.0;

   G4int Ap = aParticle->GetDefinition()->GetBaryonNumber();
   G4int Zp = G4int(aParticle->GetDefinition()->GetPDGCharge() / eplus + 0.5); 
   G4double ke_per_N = aParticle->GetKineticEnergy() / Ap;

   // Apply energy check, if less than lower limit then 0 value is returned
   // if (  ke_per_N < lowerLimit ) return xsection;

   G4int At = G4lrint(G4NistManager::Instance()->GetAtomicMassAmu(ZZ));
   G4int Zt = ZZ;  

   G4double one_third = 1.0 / 3.0;

   G4double cubicrAt = std::pow ( G4double(At) , G4double(one_third) );  
   G4double cubicrAp = std::pow ( G4double(Ap) , G4double(one_third) );  

   // rc divide fermi
   G4double Bc = Zt * Zp / ( (rc/fermi) * (cubicrAp+cubicrAt) );

   G4double targ_mass = G4NucleiProperties::GetNuclearMass(At, Zt);
   G4double proj_mass = aParticle->GetMass();
   G4double proj_momentum = aParticle->GetMomentum().mag(); 

   G4double Ecm = calEcm ( proj_mass , targ_mass , proj_momentum ); 
   if( Ecm <= Bc) return xsection;

   G4double Rvol = r0 * (  cubicrAp + cubicrAt );

//   G4double ke_per_N = aParticle->GetKineticEnergy() / Ap;
   G4double c = calCeValue ( ke_per_N / MeV  );  

   G4double a = 1.85;
   G4double Rsurf = r0 * (a*cubicrAp * cubicrAt/(cubicrAp + cubicrAt) - c); 
   G4double D = 5.0 * ( At - 2 * Zt ) * Zp / ( Ap * At );
   Rsurf = Rsurf + D * fermi;  // multiply D by fermi 

   G4double Rint = Rvol + Rsurf;
   xsection = pi * Rint * Rint * ( 1 - Bc / ( Ecm / MeV ) );
  
   return xsection; 
}

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 61 of file G4IonsKoxCrossSection.cc.

References G4ParticleDefinition::GetBaryonNumber(), and G4DynamicParticle::GetDefinition().

{
  return (1 <= aDP->GetDefinition()->GetBaryonNumber());
}

---

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

• G4IonsKoxCrossSection.hh
• G4IonsKoxCrossSection.cc