G4NucleonNuclearCrossSection Class Reference

#include <G4NucleonNuclearCrossSection.hh>

Inheritance diagram for G4NucleonNuclearCrossSection:

Public Member Functions
	G4NucleonNuclearCrossSection ()
virtual	~G4NucleonNuclearCrossSection ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *aParticle, G4int Z, const G4Material *mat=0)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *aParticle, G4int Z, const G4Material *mat=0)
virtual void	CrossSectionDescription (std::ostream &) const
G4double	GetElasticCrossSection (const G4DynamicParticle *aParticle, G4int Z)
G4double	GetTotalXsc ()
G4double	GetElasticXsc ()
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
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	GetElementCrossSection (const G4DynamicParticle *, G4int Z, 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 void	CrossSectionDescription (std::ostream &) const
void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

Static Public Member Functions
static const char *	Default_Name ()

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 47 of file G4NucleonNuclearCrossSection.hh.

Constructor & Destructor Documentation

G4NucleonNuclearCrossSection()

G4NucleonNuclearCrossSection::G4NucleonNuclearCrossSection

(

)

Definition at line 459 of file G4NucleonNuclearCrossSection.cc.

 : G4VCrossSectionDataSet(Default_Name()),
   fTotalXsc(0.0), fElasticXsc(0.0)
{
  theNeutron = G4Neutron::Neutron();
  theProton  = G4Proton::Proton();
  
  // He, Be, C
   
   thePimData.push_back(new G4PiData(he_m_t, he_m_in, e1, 44));
   thePipData.push_back(new G4PiData(he_m_t, he_p_in, e1, 44));
 
   thePimData.push_back(new G4PiData(be_m_t, be_m_in, e1, 44));
   thePipData.push_back(new G4PiData(be_m_t, be_p_in, e1, 44));
 
   thePimData.push_back(new G4PiData(c_m_t,  c_m_in,  e1, 44));
   thePipData.push_back(new G4PiData(c_m_t,  c_p_in,  e1, 44));
 
   // N, O, Na
 
   thePimData.push_back(new G4PiData(n_m_t,  n_m_in,  e2, 44));
   thePipData.push_back(new G4PiData(n_m_t,  n_p_in,  e2, 44));
 
   thePimData.push_back(new G4PiData(o_m_t,  o_m_in,  e2, 44));
   thePipData.push_back(new G4PiData(o_m_t,  o_p_in,  e2, 44));
 
   thePimData.push_back(new G4PiData(na_m_t, na_m_in, e2, 44));
   thePipData.push_back(new G4PiData(na_m_t, na_p_in, e2, 44));
 
   // Al, Si, Ca
 
   thePimData.push_back(new G4PiData(al_m_t, al_m_in, e3, 45));
   thePipData.push_back(new G4PiData(al_m_t, al_p_in, e3, 45));
 
   thePimData.push_back(new G4PiData(si_m_t, si_m_in, e3, 45));
   thePipData.push_back(new G4PiData(si_m_t, si_p_in, e3, 45));
 
   thePimData.push_back(new G4PiData(ca_m_t, ca_m_in, e3, 45));
   thePipData.push_back(new G4PiData(ca_m_t, ca_p_in, e3, 45));
 
   // Fe, Cu, Mo
 
   thePimData.push_back(new G4PiData(fe_m_t, fe_m_in, e4, 47));
   thePipData.push_back(new G4PiData(fe_m_t, fe_p_in, e4, 47));
 
   thePimData.push_back(new G4PiData(cu_m_t, cu_m_in, e4, 47));
   thePipData.push_back(new G4PiData(cu_m_t, cu_p_in, e4, 47));
 
   thePimData.push_back(new G4PiData(mo_m_t, mo_m_in, e4, 47));
   thePipData.push_back(new G4PiData(mo_m_t, mo_p_in, e4, 47));
 
   // Cd, Sn, W
 
   thePimData.push_back(new G4PiData(cd_m_t, cd_m_in, e5, 48));
   thePipData.push_back(new G4PiData(cd_m_t, cd_p_in, e5, 48));
 
   thePimData.push_back(new G4PiData(sn_m_t, sn_m_in, e5, 48));
   thePipData.push_back(new G4PiData(sn_m_t, sn_p_in, e5, 48));
 
   thePimData.push_back(new G4PiData(w_m_t,  w_m_in,  e5, 48));
   thePipData.push_back(new G4PiData(w_m_t,  w_p_in,  e5, 48));
 
   // Pb, U
 
   thePimData.push_back(new G4PiData(pb_m_t, pb_m_in, e6, 46));
   thePipData.push_back(new G4PiData(pb_m_t, pb_p_in, e6, 46));
 
   thePimData.push_back(new G4PiData(u_m_t,  u_m_in,  e6, 46));
   thePipData.push_back(new G4PiData(u_m_t,  u_p_in,  e6, 46));
 
   theZ.push_back(2); // He
   theZ.push_back(4); // Be
   theZ.push_back(6); // C
   theZ.push_back(7); // N
   theZ.push_back(8); // O
   theZ.push_back(11); // Na
   theZ.push_back(13); // Al
   theZ.push_back(14); // Si
   theZ.push_back(20); // Ca
   theZ.push_back(26); // Fe
   theZ.push_back(29); // Cu
   theZ.push_back(42); // Mo
   theZ.push_back(48); // Cd
   theZ.push_back(50); // Sn
   theZ.push_back(74); // W
   theZ.push_back(82); // Pb
   theZ.push_back(92); // U
 
}

~G4NucleonNuclearCrossSection()

G4NucleonNuclearCrossSection::~G4NucleonNuclearCrossSection ( )

virtual

Definition at line 552 of file G4NucleonNuclearCrossSection.cc.

{
   std::for_each(thePimData.begin(), thePimData.end(), G4PiData::Delete());
   std::for_each(thePipData.begin(), thePipData.end(), G4PiData::Delete());
}

Member Function Documentation

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 674 of file G4NucleonNuclearCrossSection.cc.

{
  outFile << "G4NucleonNuclearCrossSection is a variant of the Barashenkov\n"
          << "cross section parameterization to be used of protons and\n"
          << "nucleons on targets heavier than hydrogen.  It is intended for\n"
          << "use as a cross section component and is currently used by\n"
          << "G4BGGNucleonInelasticXS.  It is valid for incident energies up\n"
          << "to 1 TeV.\n"; 
}

Default_Name()

static const char * G4NucleonNuclearCrossSection::Default_Name ( )

inlinestatic

Definition at line 54 of file G4NucleonNuclearCrossSection.hh.

{return "G4NucleonNuclearCrossSection";}

Referenced by G4BGGNucleonElasticXS::BuildPhysicsTable(), and G4BGGNucleonInelasticXS::BuildPhysicsTable().

GetElasticCrossSection()

G4double G4NucleonNuclearCrossSection::GetElasticCrossSection ( const G4DynamicParticle *  aParticle, G4int  Z  )

inline

Definition at line 184 of file G4NucleonNuclearCrossSection.hh.

{
  GetElementCrossSection(dp, Z);
  return fElasticXsc;
}

Referenced by G4BGGNucleonElasticXS::BuildPhysicsTable(), and G4BGGNucleonElasticXS::GetElementCrossSection().

GetElasticXsc()

G4double G4NucleonNuclearCrossSection::GetElasticXsc ( )

inline

Definition at line 70 of file G4NucleonNuclearCrossSection.hh.

{ return fElasticXsc; };

GetElementCrossSection()

G4double G4NucleonNuclearCrossSection::GetElementCrossSection ( const G4DynamicParticle *  aParticle, G4int  Z, const G4Material *  mat = 0  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 578 of file G4NucleonNuclearCrossSection.cc.

{
   G4double kineticEnergy = aParticle->GetKineticEnergy();
  
   G4double result = 0;
   // G4cout<<"Z = "<<Z<<G4endl;
 
   size_t it = 0;
   size_t itmax = theZ.size() - 1;
   for(; it <= itmax; ++it) { if(Z <= theZ[it]) { break; } }
   if( it > itmax ) { it = itmax; }
   G4int Z1, Z2;
   G4double x1, x2, xt1, xt2;
 
   std::vector<G4PiData *> * theData = &thePimData;
   if(aParticle->GetDefinition() == theProton) { theData = &thePipData; }
 
   if( theZ[it] == Z )
     {
       result = (*theData)[it]->ReactionXSection(kineticEnergy);
       fTotalXsc = (*theData)[it]->TotalXSection(kineticEnergy);
     }
   else
     {
       if(0 == it) { it = 1; }
       x1  = (*theData)[it-1]->ReactionXSection(kineticEnergy);
       xt1 = (*theData)[it-1]->TotalXSection(kineticEnergy);
       Z1  = theZ[it-1];
       x2  = (*theData)[it]->ReactionXSection(kineticEnergy);
       xt2 = (*theData)[it]->TotalXSection(kineticEnergy);
       Z2  = theZ[it];
 
       result = Interpolate(Z1, Z2, Z, x1, x2);
       fTotalXsc = Interpolate(Z1, Z2, Z, xt1, xt2);
     }
 
   fElasticXsc = fTotalXsc - result;
   if( fElasticXsc < 0.) { fElasticXsc = 0.; }
 
   return result;
}

Referenced by G4BGGNucleonInelasticXS::BuildPhysicsTable(), GetElasticCrossSection(), and G4BGGNucleonInelasticXS::GetElementCrossSection().

GetTotalXsc()

G4double G4NucleonNuclearCrossSection::GetTotalXsc ( )

inline

Definition at line 69 of file G4NucleonNuclearCrossSection.hh.

{ return fTotalXsc;   };

IsElementApplicable()

G4bool G4NucleonNuclearCrossSection::IsElementApplicable ( const G4DynamicParticle *  aParticle, G4int  Z, const G4Material *  mat = 0  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 562 of file G4NucleonNuclearCrossSection.cc.

{
  G4bool result = false;
  if(aParticle->GetDefinition() == theNeutron ) result = true;
  if(aParticle->GetDefinition() == theProton)   result = true;
  if(Z < 2)                                     result = false;
  if(aParticle->GetKineticEnergy() > 999.9*GeV) result = false;
  return result;
}

---

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

• G4NucleonNuclearCrossSection.hh
• G4NucleonNuclearCrossSection.cc