G4ComponentGGHadronNucleusXsc Class Reference

#include <G4ComponentGGHadronNucleusXsc.hh>

Inheritance diagram for G4ComponentGGHadronNucleusXsc:

Public Member Functions
	G4ComponentGGHadronNucleusXsc ()
	~G4ComponentGGHadronNucleusXsc () final
G4double	GetTotalElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A) final
G4double	GetTotalIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A) final
G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A) final
G4double	GetInelasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A) final
G4double	GetElasticElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A) final
G4double	GetElasticIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A) final
G4double	ComputeQuasiElasticRatio (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A) final
void	ComputeCrossSections (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A, G4int nL=0)
G4double	GetProductionElementCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4double A)
G4double	GetProductionIsotopeCrossSection (const G4ParticleDefinition *aParticle, G4double kinEnergy, G4int Z, G4int A)
G4double	GetRatioSD (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetRatioQE (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetHadronNucleonXsc (const G4DynamicParticle *, const G4Element *)
G4double	GetHadronNucleonXsc (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetHadronNucleonXscPDG (const G4DynamicParticle *, const G4Element *)
G4double	GetHadronNucleonXscPDG (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetHadronNucleonXscNS (const G4DynamicParticle *, const G4Element *)
G4double	GetHadronNucleonXscNS (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetHNinelasticXsc (const G4DynamicParticle *, const G4Element *)
G4double	GetHNinelasticXsc (const G4DynamicParticle *, G4int At, G4int Zt)
G4double	GetHNinelasticXscVU (const G4DynamicParticle *, G4int At, G4int Zt)
void	Description (std::ostream &) const final
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
G4double	GetElasticGlauberGribov (const G4DynamicParticle *, G4int Z, G4int A)
G4double	GetInelasticGlauberGribov (const G4DynamicParticle *, G4int Z, G4int A)
G4double	GetAxsc2piR2 () const
G4double	GetModelInLog () const
G4double	GetTotalGlauberGribovXsc () const
G4double	GetElasticGlauberGribovXsc () const
G4double	GetInelasticGlauberGribovXsc () const
G4double	GetProductionGlauberGribovXsc () const
G4double	GetDiffractionGlauberGribovXsc () const
G4double	GetParticleBarCorTot (const G4ParticleDefinition *theParticle, G4int Z)
G4double	GetParticleBarCorIn (const G4ParticleDefinition *theParticle, G4int Z)
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 void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
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

Static Public Member Functions
static const char *	Default_Name ()

Detailed Description

Definition at line 51 of file G4ComponentGGHadronNucleusXsc.hh.

Constructor & Destructor Documentation

G4ComponentGGHadronNucleusXsc()

G4ComponentGGHadronNucleusXsc::G4ComponentGGHadronNucleusXsc ( )

explicit

Definition at line 52 of file G4ComponentGGHadronNucleusXsc.cc.

 : G4VComponentCrossSection(Default_Name()),
   fTotalXsc(0.0),fElasticXsc(0.0),fInelasticXsc(0.0),fProductionXsc(0.0),
   fDiffractionXsc(0.0),fAxsc2piR2(0.0),fModelInLog(0.0),fEnergy(0.0),
   fParticle(nullptr),fZ(0),fA(0), fL(0)
{
  theGamma    = G4Gamma::Gamma();
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  theAProton  = G4AntiProton::AntiProton();
  theANeutron = G4AntiNeutron::AntiNeutron();
  thePiPlus   = G4PionPlus::PionPlus();
  thePiMinus  = G4PionMinus::PionMinus();
  theKPlus    = G4KaonPlus::KaonPlus();
  theKMinus   = G4KaonMinus::KaonMinus();
  theK0S      = G4KaonZeroShort::KaonZeroShort();
  theK0L      = G4KaonZeroLong::KaonZeroLong();
  theLambda = G4Lambda::Lambda();
  
  hnXsc = new G4HadronNucleonXsc();
}

~G4ComponentGGHadronNucleusXsc()

G4ComponentGGHadronNucleusXsc::~G4ComponentGGHadronNucleusXsc ( )

final

Definition at line 76 of file G4ComponentGGHadronNucleusXsc.cc.

{
  delete hnXsc;
}

Member Function Documentation

ComputeCrossSections()

void G4ComponentGGHadronNucleusXsc::ComputeCrossSections	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4int	A,
		G4int	nL = 0 )

Definition at line 182 of file G4ComponentGGHadronNucleusXsc.cc.

{
  // check cache
  if(aParticle == fParticle && fZ == Z && fA == A  && fL == nL && kinEnergy == fEnergy)
    { return; }
  fParticle = aParticle;
  fZ = Z;
  fA = A;
  fL = nL;
  fEnergy = kinEnergy;
  G4Pow* pG4Pow=G4Pow::GetInstance();
  //
  G4double cofInelastic = 2.4;
  static const G4double cofTotal = 2.0;
  G4double sigma(0.0), hpInXsc(0.0), hnInXsc(0.0), R(0.0); 
  
  G4int N = std::max(A - Z, 0);  // number of neutrons
 
  if( aParticle == theKPlus || aParticle == theKMinus  || 
      aParticle == theK0S   || aParticle == theK0L) 
  {
    sigma = (1 == Z) 
      ? hnXsc->KaonNucleonXscNS(aParticle, theProton, kinEnergy)
      : Z*hnXsc->KaonNucleonXscGG(aParticle, theProton, kinEnergy);
    hpInXsc = hnXsc->GetInelasticHadronNucleonXsc();
 
    if(N > 0) { 
      sigma += N*hnXsc->KaonNucleonXscGG(aParticle, theNeutron, kinEnergy);
      hnInXsc = hnXsc->GetInelasticHadronNucleonXsc();
    }
    R = G4NuclearRadii::RadiusKNGG(A);
    cofInelastic = 2.2;
  }
  else
  {
    sigma = Z*hnXsc->HadronNucleonXsc(aParticle, theProton, kinEnergy);
    hpInXsc = hnXsc->GetInelasticHadronNucleonXsc();
 
    if(N > 0) { 
      sigma += N*hnXsc->HadronNucleonXsc(aParticle, theNeutron, kinEnergy);
      hnInXsc = hnXsc->GetInelasticHadronNucleonXsc();
    }
    R = G4NuclearRadii::RadiusHNGG(A);
    if( nL > 0 ) {
      G4double mp = theProton->GetPDGMass();
      G4double ml = theLambda->GetPDGMass();
      G4double kinCof = ml/mp; // moving hyperon - rest nucleon
      G4double cHN(0.88);
      sigma += nL*hnXsc->HadronNucleonXsc(theLambda, theProton, kinEnergy*kinCof);
      R *= std::sqrt( pG4Pow->Z23( A - nL ) + cHN*pG4Pow->Z23( nL ) )/pG4Pow->Z13(A);
    }
  }
 
  G4double nucleusSquare = cofTotal*pi*R*R;   // basically 2piRR
  G4double ratio = sigma/nucleusSquare;
  G4double difratio = ratio/(1.+ratio);
  fDiffractionXsc = 0.5*nucleusSquare*( difratio - G4Log( 1. + difratio ) );
 
  if( A > 1 )
  { 
    fTotalXsc = nucleusSquare*G4Log( 1. + ratio )
      *GetParticleBarCorTot(aParticle, Z);
 
    // inelastic xsc
    fAxsc2piR2 = cofInelastic*ratio;
    fModelInLog = G4Log( 1. + fAxsc2piR2 );
    fInelasticXsc = nucleusSquare*fModelInLog/cofInelastic;
    G4double barcorr = GetParticleBarCorIn(aParticle, Z);
    fInelasticXsc *= barcorr;
    fElasticXsc = std::max(fTotalXsc - fInelasticXsc, 0.);
 
    G4double xratio = (Z*hpInXsc + N*hnInXsc)/nucleusSquare;
    fProductionXsc = 
      nucleusSquare*G4Log(1. + cofInelastic*xratio)*barcorr/cofInelastic;
    fProductionXsc = std::min(fProductionXsc, fInelasticXsc);
  }
  else // H
  {
    fTotalXsc = sigma;
    fInelasticXsc = hpInXsc;
    fElasticXsc   = std::max(fTotalXsc - fInelasticXsc, 0.);
    fProductionXsc = fInelasticXsc;
    fDiffractionXsc = 0.2*fInelasticXsc;
    // G4double xratio = hpInXsc/nucleusSquare;
    // fProductionXsc = nucleusSquare*G4Log(1. + cofInelastic*xratio)/cofInelastic;
    // fProductionXsc = std::min(fProductionXsc, fInelasticXsc);
  }
  /*  
  G4cout << "GGXsc: Z= " << Z << " A= " << A << " E= " << kinEnergy 
     << " xtot(b)= " << fTotalXsc/barn
     << " xel(b)= " <<  fElasticXsc/barn << " xinel(b)= " << fInelasticXsc/barn
     << G4endl;
  */
}

Referenced by ComputeQuasiElasticRatio(), GetElasticElementCrossSection(), GetElasticGlauberGribov(), GetElasticIsotopeCrossSection(), GetHadronNucleonXscPDG(), GetInelasticElementCrossSection(), GetInelasticGlauberGribov(), GetInelasticIsotopeCrossSection(), GetIsoCrossSection(), GetProductionElementCrossSection(), GetProductionIsotopeCrossSection(), GetRatioQE(), GetRatioSD(), GetTotalElementCrossSection(), and GetTotalIsotopeCrossSection().

ComputeQuasiElasticRatio()

G4double G4ComponentGGHadronNucleusXsc::ComputeQuasiElasticRatio ( const G4ParticleDefinition * aParticle, G4double kinEnergy, G4int Z, G4int A )

finalvirtual

Reimplemented from G4VComponentCrossSection.

Definition at line 143 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, A);
  G4double ratio = (fInelasticXsc > 0.) 
    ? (fInelasticXsc - fProductionXsc)/fInelasticXsc : 0.;
  ratio = std::max(ratio, 0.);
  return ratio;
}

Default_Name()

static const char * G4ComponentGGHadronNucleusXsc::Default_Name ( )

inlinestatic

Definition at line 58 of file G4ComponentGGHadronNucleusXsc.hh.

{ return "Glauber-Gribov"; }

Referenced by G4ComponentGGHadronNucleusXsc().

Description()

void G4ComponentGGHadronNucleusXsc::Description ( std::ostream & outFile ) const

finalvirtual

Reimplemented from G4VComponentCrossSection.

Definition at line 458 of file G4ComponentGGHadronNucleusXsc.cc.

{
  outFile << "G4ComponentGGHadronNucleusXsc calculates total, inelastic and\n"
          << "elastic cross sections for hadron-nucleus cross sections using\n"
          << "the Glauber model with Gribov corrections.  It is valid for all\n"
          << "targets except hydrogen, and for incident p, pbar, n, sigma-,\n"
          << "pi+, pi-, K+, K- and gammas with energies above 3 GeV.  This is\n"
          << "a cross section component which is to be used to build a cross\n"
          << "data set.\n";
}

GetAxsc2piR2()

G4double G4ComponentGGHadronNucleusXsc::GetAxsc2piR2 ( ) const

inline

Definition at line 127 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fAxsc2piR2; };

GetDiffractionGlauberGribovXsc()

G4double G4ComponentGGHadronNucleusXsc::GetDiffractionGlauberGribovXsc ( ) const

inline

Definition at line 133 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fDiffractionXsc; }; 

GetElasticElementCrossSection()

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

finalvirtual

Implements G4VComponentCrossSection.

Definition at line 123 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, G4lrint(A));
  return fElasticXsc;
}

GetElasticGlauberGribov()

G4double G4ComponentGGHadronNucleusXsc::GetElasticGlauberGribov ( const G4DynamicParticle * dp, G4int Z, G4int A )

inline

Definition at line 193 of file G4ComponentGGHadronNucleusXsc.hh.

{
  ComputeCrossSections(dp->GetDefinition(), dp->GetKineticEnergy(), Z, A);
  return fElasticXsc;
}

GetElasticGlauberGribovXsc()

G4double G4ComponentGGHadronNucleusXsc::GetElasticGlauberGribovXsc ( ) const

inline

Definition at line 130 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fElasticXsc; }; 

GetElasticIsotopeCrossSection()

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

finalvirtual

Implements G4VComponentCrossSection.

Definition at line 133 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, A);
  return fElasticXsc;
}

GetHadronNucleonXsc() [1/2]

G4double G4ComponentGGHadronNucleusXsc::GetHadronNucleonXsc	(	const G4DynamicParticle *	aParticle,
		const G4Element *	anElement )

Definition at line 310 of file G4ComponentGGHadronNucleusXsc.cc.

{
  G4int At = G4lrint(anElement->GetN());  // number of nucleons 
  G4int Zt = anElement->GetZasInt();  // number of protons
 
  return GetHadronNucleonXsc(aParticle, At, Zt);
}

Referenced by GetHadronNucleonXsc().

GetHadronNucleonXsc() [2/2]

G4double G4ComponentGGHadronNucleusXsc::GetHadronNucleonXsc	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 323 of file G4ComponentGGHadronNucleusXsc.cc.

{
  return hnXsc->HadronNucleonXscEL(aParticle->GetDefinition(), theProton, 
                   aParticle->GetKineticEnergy());
}

GetHadronNucleonXscNS() [1/2]

G4double G4ComponentGGHadronNucleusXsc::GetHadronNucleonXscNS	(	const G4DynamicParticle *	aParticle,
		const G4Element *	anElement )

Definition at line 372 of file G4ComponentGGHadronNucleusXsc.cc.

{
  G4int At = G4lrint(anElement->GetN());  // number of nucleons 
  G4int Zt = anElement->GetZasInt();      // number of protons
 
  return GetHadronNucleonXscNS(aParticle, At, Zt);
}

Referenced by GetHadronNucleonXscNS().

GetHadronNucleonXscNS() [2/2]

G4double G4ComponentGGHadronNucleusXsc::GetHadronNucleonXscNS	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 386 of file G4ComponentGGHadronNucleusXsc.cc.

{
  G4double res = 0.0;
  if(1 == At && 1 == Zt) {
    res = hnXsc->HadronNucleonXscNS(aParticle->GetDefinition(), theProton, 
                    aParticle->GetKineticEnergy());
  } else if(1 == At && 0 == Zt) {
    res = hnXsc->HadronNucleonXscNS(aParticle->GetDefinition(), theNeutron, 
                    aParticle->GetKineticEnergy());
  }
  return res;
}

GetHadronNucleonXscPDG() [1/2]

G4double G4ComponentGGHadronNucleusXsc::GetHadronNucleonXscPDG	(	const G4DynamicParticle *	aParticle,
		const G4Element *	anElement )

Definition at line 335 of file G4ComponentGGHadronNucleusXsc.cc.

{
  G4int At = G4lrint(anElement->GetN());  // number of nucleons 
  G4int Zt = anElement->GetZasInt();      // number of protons
 
  return GetHadronNucleonXscPDG(aParticle, At, Zt);
}

Referenced by GetHadronNucleonXscPDG().

GetHadronNucleonXscPDG() [2/2]

G4double G4ComponentGGHadronNucleusXsc::GetHadronNucleonXscPDG	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 349 of file G4ComponentGGHadronNucleusXsc.cc.

{
  G4double res = 0.0;
  if(1 == At && 1 == Zt) {
    res = hnXsc->HadronNucleonXscPDG(aParticle->GetDefinition(), theProton, 
                     aParticle->GetKineticEnergy());
  } else if(1 == At && 0 == Zt) {
    res = hnXsc->HadronNucleonXscPDG(aParticle->GetDefinition(), theNeutron, 
                     aParticle->GetKineticEnergy());
  } else {
    ComputeCrossSections(aParticle->GetDefinition(),
             aParticle->GetKineticEnergy(), Zt, At);
    res = fTotalXsc;
  }
  return res;
}

GetHNinelasticXsc() [1/2]

G4double G4ComponentGGHadronNucleusXsc::GetHNinelasticXsc	(	const G4DynamicParticle *	aParticle,
		const G4Element *	anElement )

Definition at line 405 of file G4ComponentGGHadronNucleusXsc.cc.

{
  G4int At = G4lrint(anElement->GetN());  // number of nucleons 
  G4int Zt = anElement->GetZasInt();      // number of protons
 
  return GetHNinelasticXsc(aParticle, At, Zt);
}

Referenced by GetHNinelasticXsc().

GetHNinelasticXsc() [2/2]

G4double G4ComponentGGHadronNucleusXsc::GetHNinelasticXsc	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 418 of file G4ComponentGGHadronNucleusXsc.cc.

{
  const G4ParticleDefinition* hadron = aParticle->GetDefinition();
  G4double e = aParticle->GetKineticEnergy();
  G4int Nt = std::max(At - Zt, 0);
 
  hnXsc->HadronNucleonXscNS(hadron, theProton, e);
  G4double sumInelastic = Zt*hnXsc->GetInelasticHadronNucleonXsc();
  if(Nt > 0) {
    hnXsc->HadronNucleonXscNS(hadron, theNeutron, e);
    sumInelastic += Nt*hnXsc->GetInelasticHadronNucleonXsc();
  }
  return sumInelastic;
}

GetHNinelasticXscVU()

G4double G4ComponentGGHadronNucleusXsc::GetHNinelasticXscVU	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 438 of file G4ComponentGGHadronNucleusXsc.cc.

{
  const G4ParticleDefinition* hadron = aParticle->GetDefinition();
  G4double e = aParticle->GetKineticEnergy();
  G4int Nt = std::max(At - Zt, 0);
 
  hnXsc->HadronNucleonXscVU(hadron, theProton, e);
  G4double sumInelastic = Zt*hnXsc->GetInelasticHadronNucleonXsc();
  if(Nt > 0) {
    hnXsc->HadronNucleonXscVU(hadron, theNeutron, e);
    sumInelastic += Nt*hnXsc->GetInelasticHadronNucleonXsc();
  }
  return sumInelastic;
}

GetInelasticElementCrossSection()

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

finalvirtual

Implements G4VComponentCrossSection.

Definition at line 103 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, G4lrint(A));
  return fInelasticXsc;
}

GetInelasticGlauberGribov()

G4double G4ComponentGGHadronNucleusXsc::GetInelasticGlauberGribov ( const G4DynamicParticle * dp, G4int Z, G4int A )

inline

Definition at line 203 of file G4ComponentGGHadronNucleusXsc.hh.

{
  ComputeCrossSections(dp->GetDefinition(), dp->GetKineticEnergy(), Z, A);
  return fInelasticXsc;
}

GetInelasticGlauberGribovXsc()

G4double G4ComponentGGHadronNucleusXsc::GetInelasticGlauberGribovXsc ( ) const

inline

Definition at line 131 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fInelasticXsc; }; 

GetInelasticIsotopeCrossSection()

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

finalvirtual

Implements G4VComponentCrossSection.

Definition at line 113 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, A);
  return fInelasticXsc;
}

GetIsoCrossSection()

G4double G4ComponentGGHadronNucleusXsc::GetIsoCrossSection ( const G4DynamicParticle * dp, G4int Z, G4int A, const G4Isotope * iso = nullptr, const G4Element * elm = nullptr, const G4Material * mat = nullptr )

inline

Definition at line 182 of file G4ComponentGGHadronNucleusXsc.hh.

{
  ComputeCrossSections(dp->GetDefinition(), dp->GetKineticEnergy(), Z, A);
  return fTotalXsc;
}

GetModelInLog()

G4double G4ComponentGGHadronNucleusXsc::GetModelInLog ( ) const

inline

Definition at line 128 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fModelInLog; };

GetParticleBarCorIn()

G4double G4ComponentGGHadronNucleusXsc::GetParticleBarCorIn ( const G4ParticleDefinition * theParticle, G4int Z )

inline

Definition at line 233 of file G4ComponentGGHadronNucleusXsc.hh.

{
  G4double cor = 1.0;
  G4int z = std::min(92, std::max(ZZ, 1));
  if(      theParticle == theProton ) cor = fProtonBarCorrectionIn[z]; 
  else if( theParticle == theNeutron) cor = fNeutronBarCorrectionIn[z]; 
  else if( theParticle == thePiPlus ) cor = fPionPlusBarCorrectionIn[z];
  else if( theParticle == thePiMinus) cor = fPionMinusBarCorrectionIn[z];
  return cor;
}

Referenced by ComputeCrossSections().

GetParticleBarCorTot()

G4double G4ComponentGGHadronNucleusXsc::GetParticleBarCorTot ( const G4ParticleDefinition * theParticle, G4int Z )

inline

Definition at line 215 of file G4ComponentGGHadronNucleusXsc.hh.

{
  G4double cor = 1.0;
  G4int z = std::min(92, std::max(ZZ, 1));
  if(      theParticle == theProton ) cor = fProtonBarCorrectionTot[z]; 
  else if( theParticle == theNeutron) cor = fNeutronBarCorrectionTot[z]; 
  else if( theParticle == thePiPlus ) cor = fPionPlusBarCorrectionTot[z];
  else if( theParticle == thePiMinus) cor = fPionMinusBarCorrectionTot[z];
  return cor;
}

Referenced by ComputeCrossSections().

GetProductionElementCrossSection()

G4double G4ComponentGGHadronNucleusXsc::GetProductionElementCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4double	A )

Definition at line 156 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, G4lrint(A));
  return fProductionXsc;
}

GetProductionGlauberGribovXsc()

G4double G4ComponentGGHadronNucleusXsc::GetProductionGlauberGribovXsc ( ) const

inline

Definition at line 132 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fProductionXsc; }; 

GetProductionIsotopeCrossSection()

G4double G4ComponentGGHadronNucleusXsc::GetProductionIsotopeCrossSection	(	const G4ParticleDefinition *	aParticle,
		G4double	kinEnergy,
		G4int	Z,
		G4int	A )

Definition at line 166 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, A);
  return fProductionXsc;
}

GetRatioQE()

G4double G4ComponentGGHadronNucleusXsc::GetRatioQE	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 296 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle->GetDefinition(), 
               aParticle->GetKineticEnergy(), Z, A);
 
  return (fInelasticXsc > std::max(fProductionXsc, 0.)) 
    ? 1.0 - fProductionXsc/fInelasticXsc : 0.0;
}

GetRatioSD()

G4double G4ComponentGGHadronNucleusXsc::GetRatioSD	(	const G4DynamicParticle *	aParticle,
		G4int	At,
		G4int	Zt )

Definition at line 283 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle->GetDefinition(), 
               aParticle->GetKineticEnergy(), Z, A);
 
  return (fInelasticXsc > 0.0) ? fDiffractionXsc/fInelasticXsc : 0.0;
}

GetTotalElementCrossSection()

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

finalvirtual

Implements G4VComponentCrossSection.

Definition at line 83 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, G4lrint(A));
  return fTotalXsc;
}

GetTotalGlauberGribovXsc()

G4double G4ComponentGGHadronNucleusXsc::GetTotalGlauberGribovXsc ( ) const

inline

Definition at line 129 of file G4ComponentGGHadronNucleusXsc.hh.

{ return fTotalXsc; }; 

GetTotalIsotopeCrossSection()

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

finalvirtual

Implements G4VComponentCrossSection.

Definition at line 93 of file G4ComponentGGHadronNucleusXsc.cc.

{
  ComputeCrossSections(aParticle, kinEnergy, Z, A);
  return fTotalXsc;
}

---

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

• G4ComponentGGHadronNucleusXsc.hh
• G4ComponentGGHadronNucleusXsc.cc