G4IonsShenCrossSection Class Reference

#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

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

(

)

Definition at line 42 of file G4IonsShenCrossSection.cc.

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

G4IonsShenCrossSection::~G4IonsShenCrossSection ( )

virtual

Definition at line 49 of file G4IonsShenCrossSection.cc.

{}

Member Function Documentation

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

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:

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:

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:

---

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

• geant4.10.03.p01/source/processes/hadronic/cross_sections/include/G4IonsShenCrossSection.hh
• geant4.10.03.p01/source/processes/hadronic/cross_sections/src/G4IonsShenCrossSection.cc