G4IonsKoxCrossSection Class Reference

#include <G4IonsKoxCrossSection.hh>

Inheritance diagram for G4IonsKoxCrossSection:

Collaboration 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

Private Member Functions
G4double	calEcm (G4double, G4double, G4double)
G4double	calCeValue (G4double)

Private Attributes
G4double	r0
G4double	rc

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::G4IonsKoxCrossSection

(

)

Definition at line 41 of file G4IonsKoxCrossSection.cc.

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

~G4IonsKoxCrossSection()

G4IonsKoxCrossSection::~G4IonsKoxCrossSection

(

)

Definition at line 47 of file G4IonsKoxCrossSection.cc.

{}

Member Function Documentation

calCeValue()

G4double G4IonsKoxCrossSection::calCeValue ( G4double  ke )

private

Definition at line 125 of file G4IonsKoxCrossSection.cc.

{
   // Calculate c value 
   // This value is indepenent from projectile and target particle 
   // ke is projectile kinetic energy per nucleon in the Lab system with MeV unit 
   // fitting function is made by T. Koi 
   // There are no data below 30 MeV/n in Kox et al., 

   G4double Ce; 
   G4double log10_ke = std::log10 ( ke );   
   if (log10_ke > 1.5) 
   {
      Ce = - 10.0 / G4Pow::GetInstance()->powA ( G4double(log10_ke) , G4double(5) ) + 2.0;
   }
   else
   {
      Ce = (-10.0/G4Pow::GetInstance()->powA(G4double(1.5), G4double(5) ) + 2.0) /
           G4Pow::GetInstance()->powA(G4double(1.5), G4double(3)) * G4Pow::GetInstance()->powA(G4double(log10_ke), G4double(3) );

   }
   return Ce;
}

Here is the call graph for this function:

Here is the caller graph for this function:

calEcm()

G4double G4IonsKoxCrossSection::calEcm ( G4double  mp, G4double  mt, G4double  Plab  )

private

Definition at line 115 of file G4IonsKoxCrossSection.cc.

{
   G4double Elab = std::sqrt ( mp * mp + Plab * Plab );
   G4double Ecm = std::sqrt ( mp * mp + mt * mt + 2 * Elab * mt );
   G4double Pcm = Plab * mt / Ecm;
   G4double KEcm = std::sqrt ( Pcm * Pcm + mp * mp ) - mp;
   return KEcm;
}

Here is the caller graph for this function:

CrossSectionDescription()

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

GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 68 of file G4IonsKoxCrossSection.cc.

{
   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 = G4Pow::GetInstance()->powA ( G4double(At) , G4double(one_third) );  
   G4double cubicrAp = G4Pow::GetInstance()->powA ( 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; 
}

Here is the call graph for this function:

IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 61 of file G4IonsKoxCrossSection.cc.

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

Here is the call graph for this function:

Member Data Documentation

r0

G4double G4IonsKoxCrossSection::r0

private

Definition at line 70 of file G4IonsKoxCrossSection.hh.

rc

G4double G4IonsKoxCrossSection::rc

private

Definition at line 71 of file G4IonsKoxCrossSection.hh.

---

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

• G4IonsKoxCrossSection.hh
• G4IonsKoxCrossSection.cc