#include <G4IonsShenCrossSection.hh>

Inheritance diagram for G4IonsShenCrossSection:

Collaboration 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

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 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	calEcmValue (const G4double, const G4double, const G4double)

G4double	calCeValue (const G4double)

Private Attributes
const G4double	upperLimit

const G4double	r0

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 51 of file G4IonsShenCrossSection.hh.

Constructor & Destructor Documentation

◆ G4IonsShenCrossSection()

G4IonsShenCrossSection::G4IonsShenCrossSection ( )

Definition at line 42 of file G4IonsShenCrossSection.cc.

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

◆ ~G4IonsShenCrossSection()

G4IonsShenCrossSection::~G4IonsShenCrossSection ( )

virtual

Definition at line 49 of file G4IonsShenCrossSection.cc.

50 {}

Member Function Documentation

◆ calCeValue()

G4double G4IonsShenCrossSection::calCeValue ( const G4double ke )

private

Definition at line 144 of file G4IonsShenCrossSection.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/std::pow(G4double(log10_ke), G4double(5)) + 2.0;
    }
    else
    {
      Ce = (-10.0/std::pow(G4double(1.5), G4double(5) ) + 2.0) /
          std::pow(G4double(1.5) , G4double(3)) * std::pow(G4double(log10_ke), G4double(3));
    }
    return Ce;
 }

Here is the caller graph for this function:

◆ calEcmValue()

G4double G4IonsShenCrossSection::calEcmValue	(	const G4double	mp,
		const G4double	mt,
		const G4double	Plab
	)

private

Definition at line 133 of file G4IonsShenCrossSection.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 G4IonsShenCrossSection::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 53 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";
 }

◆ GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 70 of file G4IonsShenCrossSection.cc.

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

Here is the call graph for this function:

◆ GetIsoCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 78 of file G4IonsShenCrossSection.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

◆ IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 63 of file G4IonsShenCrossSection.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

◆ r0

const G4double G4IonsShenCrossSection::r0

private

Definition at line 78 of file G4IonsShenCrossSection.hh.

◆ upperLimit

const G4double G4IonsShenCrossSection::upperLimit

private

Definition at line 76 of file G4IonsShenCrossSection.hh.

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4IonsShenCrossSection()

◆ ~G4IonsShenCrossSection()

Member Function Documentation

◆ calCeValue()

◆ calEcmValue()

◆ CrossSectionDescription()

◆ GetElementCrossSection()

◆ GetIsoCrossSection()

◆ IsElementApplicable()

Member Data Documentation

◆ r0

◆ upperLimit