#include <G4ComponentBarNucleonNucleusXsc.hh>

Inheritance diagram for G4ComponentBarNucleonNucleusXsc:

Collaboration diagram for G4ComponentBarNucleonNucleusXsc:

Public Member Functions
	G4ComponentBarNucleonNucleusXsc ()

virtual	~G4ComponentBarNucleonNucleusXsc ()

virtual G4double	GetTotalIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int)

virtual G4double	GetTotalElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double)

virtual G4double	GetInelasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int)

virtual G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double)

virtual G4double	GetElasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double)

virtual G4double	GetElasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int)

G4bool	IsElementApplicable (const G4DynamicParticle *aParticle, G4int Z)

G4double	GetElementCrossSection (const G4DynamicParticle *aParticle, G4int Z)

void	CrossSectionDescription (std::ostream &) const

G4double	GetElasticCrossSection (const G4DynamicParticle *aParticle, G4int Z)

G4double	GetTotalXsc ()

G4double	GetElasticXsc ()

Public Member Functions inherited from G4VComponentCrossSection
	G4VComponentCrossSection (const G4String &nam="")

virtual	~G4VComponentCrossSection ()

G4double	GetTotalElementCrossSection (const G4ParticleDefinition , G4double kinEnergy, const G4Element )

G4double	GetInelasticElementCrossSection (const G4ParticleDefinition , G4double kinEnergy, const G4Element )

G4double	GetElasticElementCrossSection (const G4ParticleDefinition , G4double kinEnergy, const G4Element )

virtual G4double	ComputeQuasiElasticRatio (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	Description () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

const G4String &	GetName () const

Detailed Description

Definition at line 47 of file G4ComponentBarNucleonNucleusXsc.hh.

Constructor & Destructor Documentation

G4ComponentBarNucleonNucleusXsc::G4ComponentBarNucleonNucleusXsc ( )

Definition at line 451 of file G4ComponentBarNucleonNucleusXsc.cc.

  : G4VComponentCrossSection("G4ComponentBarNucleonNucleusXsc"),
    fTotalXsc(0.0), fInelasticXsc(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
 
 }

Here is the call graph for this function:

G4ComponentBarNucleonNucleusXsc::~G4ComponentBarNucleonNucleusXsc ( )

virtual

Definition at line 544 of file G4ComponentBarNucleonNucleusXsc.cc.

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

Member Function Documentation

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

Definition at line 752 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   outFile << "G4ComponentBarNucleonNucleusXsc 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"; 
 }

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

inline

Definition at line 221 of file G4ComponentBarNucleonNucleusXsc.hh.

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

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetElasticElementCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4double
	)

virtual

Implements G4VComponentCrossSection.

Definition at line 608 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   G4DynamicParticle* aDP = new G4DynamicParticle(aParticle,G4ParticleMomentum(1.,0.,0.), 
                                                 kinEnergy);
   fInelasticXsc = GetElementCrossSection(aDP, Z);
   delete aDP;
 
   return fElasticXsc;
 }

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetElasticIsotopeCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4int
	)

virtual

Implements G4VComponentCrossSection.

Definition at line 622 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   G4DynamicParticle* aDP = new G4DynamicParticle(aParticle,G4ParticleMomentum(1.,0.,0.), 
                          kinEnergy);
   fInelasticXsc = GetElementCrossSection(aDP, Z);
   delete aDP;
 
   return fElasticXsc;
 }

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetElasticXsc ( )

inline

Definition at line 106 of file G4ComponentBarNucleonNucleusXsc.hh.

106 { return fElasticXsc; };

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

Definition at line 656 of file G4ComponentBarNucleonNucleusXsc.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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetInelasticElementCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4double
	)

virtual

Implements G4VComponentCrossSection.

Definition at line 594 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   G4DynamicParticle* aDP = new G4DynamicParticle(aParticle,G4ParticleMomentum(1.,0.,0.), 
                                                 kinEnergy);
   fInelasticXsc = GetElementCrossSection(aDP, Z);
   delete aDP;
 
   return fInelasticXsc;
 }

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetInelasticIsotopeCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4int
	)

virtual

Implements G4VComponentCrossSection.

Definition at line 580 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   G4DynamicParticle* aDP = new G4DynamicParticle(aParticle,G4ParticleMomentum(1.,0.,0.), 
                                                 kinEnergy);
   fInelasticXsc = GetElementCrossSection(aDP, Z);
   delete aDP;
 
   return fInelasticXsc;
 }

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetTotalElementCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4double
	)

virtual

Implements G4VComponentCrossSection.

Definition at line 566 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   G4DynamicParticle* aDP = new G4DynamicParticle(aParticle,G4ParticleMomentum(1.,0.,0.), 
                                                 kinEnergy);
   fInelasticXsc = GetElementCrossSection(aDP, Z);
   delete aDP;
 
   return fTotalXsc;
 }

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetTotalIsotopeCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4int
	)

virtual

Implements G4VComponentCrossSection.

Definition at line 552 of file G4ComponentBarNucleonNucleusXsc.cc.

 {
   G4DynamicParticle* aDP = new G4DynamicParticle(aParticle,G4ParticleMomentum(1.,0.,0.), 
                                                 kinEnergy);
   fInelasticXsc = GetElementCrossSection(aDP, Z);
   delete aDP;
 
   return fTotalXsc;
 }

Here is the call graph for this function:

G4double G4ComponentBarNucleonNucleusXsc::GetTotalXsc ( )

inline

Definition at line 105 of file G4ComponentBarNucleonNucleusXsc.hh.

105 { return fTotalXsc; };

G4bool G4ComponentBarNucleonNucleusXsc::IsElementApplicable	(	const G4DynamicParticle *	aParticle,
		G4int	Z
	)

Definition at line 640 of file G4ComponentBarNucleonNucleusXsc.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;
 }

Here is the call graph for this function:

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

source/geant4.10.03.p03/source/processes/hadronic/cross_sections/include/G4ComponentBarNucleonNucleusXsc.hh
source/geant4.10.03.p03/source/processes/hadronic/cross_sections/src/G4ComponentBarNucleonNucleusXsc.cc

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation