G4IonsShenCrossSection Class Reference

#include <G4IonsShenCrossSection.hh>

Inheritance diagram for G4IonsShenCrossSection:

Public Member Functions
	G4IonsShenCrossSection ()
virtual	~G4IonsShenCrossSection ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *aDP, G4int Z, const G4Material *)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual void	CrossSectionDescription (std::ostream &) const

---

Detailed Description

Definition at line 51 of file G4IonsShenCrossSection.hh.

---

Constructor & Destructor Documentation

G4IonsShenCrossSection::G4IonsShenCrossSection

(

)

Definition at line 43 of file G4IonsShenCrossSection.cc.

  : G4VCrossSectionDataSet("IonsShen"),
    upperLimit( 10*GeV ), lowerLimit( 10*MeV ), r0 ( 1.1 )
{}

G4IonsShenCrossSection::~G4IonsShenCrossSection

(

)

[virtual]

Definition at line 48 of file G4IonsShenCrossSection.cc.

{}

---

Member Function Documentation

void G4IonsShenCrossSection::CrossSectionDescription

(

std::ostream &

)

const [virtual]

Reimplemented from G4VCrossSectionDataSet.

Definition at line 52 of file G4IonsShenCrossSection.cc.

{
  outFile << "G4IonsShenCrossSection calculates the total reaction cross\n"
          << "section for nucleus-nucleus scattering using the Shen\n"
          << "parameterization.  It is valid for projectiles and targets of\n"
          << "all Z, and projectile energies up to 1 TeV/n.  Above 10 GeV/n"
          << "the cross section is constant.  Below 10 MeV/n zero cross\n"
          << "is returned.\n";
}

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

Reimplemented from G4VCrossSectionDataSet.

Definition at line 69 of file G4IonsShenCrossSection.cc.

References G4lrint(), GetIsoCrossSection(), and G4NistManager::Instance().

{
  G4int A = G4lrint(G4NistManager::Instance()->GetAtomicMassAmu(Z));
  return GetIsoCrossSection(aParticle, Z, A);
}

G4double G4IonsShenCrossSection::GetIsoCrossSection	(	const G4DynamicParticle *	,
		G4int	Z,
		G4int	A,
		const G4Isotope *	iso = 0,
		const G4Element *	elm = 0,
		const G4Material *	mat = 0
	)			[virtual]

Reimplemented from G4VCrossSectionDataSet.

Definition at line 77 of file G4IonsShenCrossSection.cc.

References G4Pow::A13(), G4lrint(), G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), G4Pow::GetInstance(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMass(), G4DynamicParticle::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGCharge(), G4INCL::Math::pi, and G4Pow::Z13().

Referenced by GetElementCrossSection().

{
   G4double xsection = 0.0;

   G4int Ap = aParticle->GetDefinition()->GetBaryonNumber();
   G4int Zp = G4lrint(aParticle->GetDefinition()->GetPDGCharge()/eplus); 
   G4double ke_per_N = aParticle->GetKineticEnergy() / Ap; 
   if ( ke_per_N > upperLimit ) { ke_per_N = upperLimit; }

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

   G4Pow* g4pow = G4Pow::GetInstance();
  
   G4double cubicrAt = g4pow->Z13(At);
   G4double cubicrAp = g4pow->Z13(Ap);
 
   G4double Rt = 1.12 * cubicrAt - 0.94 * ( 1.0 / cubicrAt );
   G4double Rp = 1.12 * cubicrAp - 0.94 * ( 1.0 / cubicrAp );

   G4double r = Rt + Rp + 3.2;   // in fm
   G4double b = 1.0;   // in MeV/fm
   G4double targ_mass = G4NucleiProperties::GetNuclearMass(At, Zt);

   G4double proj_mass = aParticle->GetMass();
   G4double proj_momentum = aParticle->GetMomentum().mag();

   G4double Ecm = calEcmValue (proj_mass, targ_mass, proj_momentum); 

   G4double B = 1.44 * Zt * Zp / r - b * Rt * Rp / ( Rt + Rp ); 
   if(Ecm <= B) { return xsection; }

   G4double c = calCeValue ( ke_per_N / MeV  );  

   G4double R1 = r0 * (cubicrAt + cubicrAp + 1.85*cubicrAt*cubicrAp/(cubicrAt + cubicrAp) - c); 

   G4double R2 = 1.0 * ( At - 2 * Zt ) * Zp / ( Ap * At );


   G4double R3 = (0.176 / g4pow->A13(Ecm)) * cubicrAt * cubicrAp /(cubicrAt + cubicrAp);

   G4double R = R1 + R2 + R3;

   xsection = 10 * pi * R * R * ( 1 - B / Ecm );   
   xsection = xsection * millibarn;   // mulitply xsection by millibarn

   return xsection; 
}

G4bool G4IonsShenCrossSection::IsElementApplicable	(	const G4DynamicParticle *	aDP,
		G4int	Z,
		const G4Material *
	)			[virtual]

Reimplemented from G4VCrossSectionDataSet.

Definition at line 62 of file G4IonsShenCrossSection.cc.

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

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

---

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

• G4IonsShenCrossSection.hh
• G4IonsShenCrossSection.cc

---