#include <G4BGGNucleonInelasticXS.hh>

Inheritance diagram for G4BGGNucleonInelasticXS:

Public Member Functions
	G4BGGNucleonInelasticXS (const G4ParticleDefinition *)

	~G4BGGNucleonInelasticXS () override

G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material mat) override

G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element elm, const G4Material *mat) override

G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat) override

G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=nullptr, const G4Element elm=nullptr, const G4Material mat=nullptr) override

void	BuildPhysicsTable (const G4ParticleDefinition &) override

void	CrossSectionDescription (std::ostream &) const override

G4BGGNucleonInelasticXS &	operator= (const G4BGGNucleonInelasticXS &right)=delete

	G4BGGNucleonInelasticXS (const G4BGGNucleonInelasticXS &)=delete

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

virtual	~G4VCrossSectionDataSet ()

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

virtual G4double	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition , const G4Element , const G4Material *mat=nullptr)

virtual G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition , G4int Z, G4int A, const G4Isotope iso=nullptr, const G4Element elm=nullptr, const G4Material mat=nullptr)

virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

bool	ForAllAtomsAndEnergies () const

void	SetForAllAtomsAndEnergies (G4bool val)

const G4String &	GetName () const

void	SetName (const G4String &nam)

G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete

	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Additional Inherited Members
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel {0}

G4String	name

Detailed Description

Definition at line 62 of file G4BGGNucleonInelasticXS.hh.

Constructor & Destructor Documentation

◆ G4BGGNucleonInelasticXS() [1/2]

G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS ( const G4ParticleDefinition * p )

explicit

Definition at line 71 of file G4BGGNucleonInelasticXS.cc.

 : G4VCrossSectionDataSet("BarashenkovGlauberGribov")
{
  verboseLevel = 0;
  fGlauberEnergy = 91.*CLHEP::GeV;
  fLowEnergy = 14.*CLHEP::MeV;
 
  fNucleon = new G4NucleonNuclearCrossSection();
  fGlauber = new G4ComponentGGHadronNucleusXsc();
  fHadron  = new G4HadronNucleonXsc();
 
  theProton= G4Proton::Proton();
  isProton = (theProton == p);
  SetForAllAtomsAndEnergies(true);
 
  if (0 == theA[0]) { Initialise(); }
}

Referenced by G4BGGNucleonInelasticXS(), G4ParticleHPBGGNucleonInelasticXS::G4ParticleHPBGGNucleonInelasticXS(), and operator=().

◆ ~G4BGGNucleonInelasticXS()

G4BGGNucleonInelasticXS::~G4BGGNucleonInelasticXS ( )

override

Definition at line 91 of file G4BGGNucleonInelasticXS.cc.

{
  delete fHadron;
}

◆ G4BGGNucleonInelasticXS() [2/2]

G4BGGNucleonInelasticXS::G4BGGNucleonInelasticXS ( const G4BGGNucleonInelasticXS & )

delete

Member Function Documentation

◆ BuildPhysicsTable()

void G4BGGNucleonInelasticXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 179 of file G4BGGNucleonInelasticXS.cc.

{
  if(&p == theProton || &p == G4Neutron::Neutron()) {
    isProton = (theProton == &p);
  } else {
    G4ExceptionDescription ed;
    ed << "This BGG cross section is applicable only to nucleons and not to " 
       << p.GetParticleName() << G4endl; 
    G4Exception("G4BGGNucleonInelasticXS::BuildPhysicsTable", "had001", 
              FatalException, ed);
    return;
  }
}

◆ CrossSectionDescription()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 290 of file G4BGGNucleonInelasticXS.cc.

{
  outFile << "The Barashenkov-Glauber-Gribov cross section calculates inelastic\n"
          << "scattering of protons and neutrons from nuclei using the\n"
          << "Barashenkov parameterization below 91 GeV and the Glauber-Gribov\n"
          << "parameterization above 91 GeV.  It uses the G4HadronNucleonXsc\n"
          << "cross section component for hydrogen targets, and the\n"
          << "G4ComponentGGHadronNucleusXsc component for other targets.\n";
}

◆ GetElementCrossSection()

G4double G4BGGNucleonInelasticXS::GetElementCrossSection	(	const G4DynamicParticle *	dp,
		G4int	Z,
		const G4Material *	mat )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 117 of file G4BGGNucleonInelasticXS.cc.

{
  G4double cross = 0.0;
  G4double ekin = dp->GetKineticEnergy();
  G4int Z = std::min(ZZ, 92);
  if (1 == Z) {
    cross = 1.0115*GetIsoCrossSection(dp,1,1);
  } else {
    if (ekin <= fLowEnergy) {
      cross = CoulombFactor(ekin, Z);
      cross *= (isProton) ? theCoulombFacP[Z] : theCoulombFacN[Z];
    } else if (ekin > fGlauberEnergy) {
      cross = fGlauber->GetInelasticGlauberGribov(dp, Z, theA[Z]);
      cross *= (isProton) ? theGlauberFacP[Z] : theGlauberFacN[Z];
    } else {
      cross = fNucleon->GetElementCrossSection(dp, Z);
    }
  }
 
#ifdef G4VERBOSE
  if (verboseLevel > 1) {
    G4cout << "G4BGGNucleonInelasticXS::GetCrossSection  for "
           << dp->GetDefinition()->GetParticleName()
           << "  Ekin(GeV)= " << dp->GetKineticEnergy()/CLHEP::GeV
           << " in nucleus Z= " << Z << "  A= " << theA[Z]
           << " XS(b)= " << cross/barn 
           << G4endl;
  }
#endif
  return cross;
}

◆ GetIsoCrossSection()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 153 of file G4BGGNucleonInelasticXS.cc.

{
  // this method should be called only for Z = 1
  fHadron->HadronNucleonXscNS(dp->GetDefinition(), theProton, 
                           dp->GetKineticEnergy());
  G4double cross = A*fHadron->GetInelasticHadronNucleonXsc();
 
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout << "G4BGGNucleonInelasticXS::GetIsoCrossSection  for "
           << dp->GetDefinition()->GetParticleName()
           << "  Ekin(GeV)= " << dp->GetKineticEnergy()/CLHEP::GeV
           << " in nucleus Z=1  A=" << A
           << " XS(b)= " << cross/barn 
           << G4endl;
  }
#endif
  return cross;
}