#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	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

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 ( )

Definition at line 459 of file G4NucleonNuclearCrossSection.cc.

References G4Neutron::Neutron(), and G4Proton::Proton().

  : 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 ( )

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

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

static const char* G4NucleonNuclearCrossSection::Default_Name ( )

inlinestatic

Definition at line 54 of file G4NucleonNuclearCrossSection.hh.

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

54 {return "G4NucleonNuclearCrossSection";}

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

inline

Definition at line 184 of file G4NucleonNuclearCrossSection.hh.

References GetElementCrossSection().

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

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

G4double G4NucleonNuclearCrossSection::GetElasticXsc ( )

inline

Definition at line 70 of file G4NucleonNuclearCrossSection.hh.

70 { return fElasticXsc; };

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 578 of file G4NucleonNuclearCrossSection.cc.

References G4DynamicParticle::GetDefinition(), and G4DynamicParticle::GetKineticEnergy().

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

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

G4double G4NucleonNuclearCrossSection::GetTotalXsc ( )

inline

Definition at line 69 of file G4NucleonNuclearCrossSection.hh.

69 { return fTotalXsc; };

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 562 of file G4NucleonNuclearCrossSection.cc.

References G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), and python.hepunit::GeV.

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

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation