#include <G4HadronCrossSections.hh>

Inheritance diagram for G4HadronCrossSections:

Public Member Functions
	G4HadronCrossSections ()

	~G4HadronCrossSections ()

G4bool	IsApplicable (const G4DynamicParticle *aParticle)

G4double	GetElasticCrossSection (const G4DynamicParticle *, G4int, G4int)

G4double	GetInelasticCrossSection (const G4DynamicParticle *, G4int, G4int)

G4double	GetCaptureCrossSection (const G4DynamicParticle *, G4int)

G4double	GetFissionCrossSection (const G4DynamicParticle *, G4int, G4int)

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel ()

Static Public Member Functions
static G4HadronCrossSections *	Instance ()

static void	SetCorrectInelasticNearZero (G4bool value)

static G4bool	GetCorrectInelasticNearZero ()

Detailed Description

Definition at line 63 of file G4HadronCrossSections.hh.

Constructor & Destructor Documentation

◆ G4HadronCrossSections()

G4HadronCrossSections::G4HadronCrossSections ( )

Definition at line 1217 of file G4HadronCrossSections.cc.

  : prevParticleDefinition(0), prevElement(0), prevZZ(0), prevAA(0), 
    prevKineticEnergy(DBL_MAX), lastEkx(0.), lastEkxPower(0.), verboseLevel(0)
{
  siginelastic = sigelastic = 0.0;
  g4pow = G4Pow::GetInstance();
}

◆ ~G4HadronCrossSections()

G4HadronCrossSections::~G4HadronCrossSections ( )

Definition at line 1225 of file G4HadronCrossSections.cc.

1226{}

Member Function Documentation

◆ GetCaptureCrossSection()

G4double G4HadronCrossSections::GetCaptureCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	ZZ
	)

Definition at line 1556 of file G4HadronCrossSections.cc.

{
  if (GetParticleCode(aParticle) != 16) return 0.;
  G4double ek = aParticle->GetKineticEnergy()/GeV;
  if (ek > 0.0327)  { return 0.; }
 
  G4double ekx = std::max(ek, 1.e-9);
  if (ekx != lastEkx) {
    lastEkx = ekx;
    lastEkxPower = std::pow(ekx*1.e6, 0.577);
  }
 
  G4int izno = ZZ;
  if (izno > 100) izno = 100;      // Not in GHESIG
  izno = izno - 1;      // For array indexing
  G4double sigcap = 11.12*cscap[izno]/lastEkxPower;
 
  // Convert cross section from mb to default units
  sigcap = sigcap*millibarn;
  return sigcap;
}

Referenced by G4HadronCaptureDataSet::GetElementCrossSection().

◆ GetCorrectInelasticNearZero()

static G4bool G4HadronCrossSections::GetCorrectInelasticNearZero ( )

inlinestatic

Definition at line 89 of file G4HadronCrossSections.hh.

90 {return correctInelasticNearZero;}

◆ GetElasticCrossSection()

G4double G4HadronCrossSections::GetElasticCrossSection	(	const G4DynamicParticle *	particle,
		G4int	ZZ,
		G4int	AA
	)

Definition at line 1260 of file G4HadronCrossSections.cc.

{
  if (particle->GetDefinition() != prevParticleDefinition ||
      particle->GetKineticEnergy() != prevKineticEnergy ||
      ZZ != prevZZ || AA != prevAA)
    { CalcScatteringCrossSections(particle, ZZ, AA); }
 
  if (verboseLevel > 1) {
    G4cout << "        x-section elastic(mb)= " 
       << sigelastic/millibarn << G4endl; 
  }
  return sigelastic;
}

Referenced by G4HadronElasticDataSet::GetElementCrossSection().

◆ GetFissionCrossSection()

G4double G4HadronCrossSections::GetFissionCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	ZZ,
		G4int	AA
	)

Definition at line 1581 of file G4HadronCrossSections.cc.

{
  if (AA < 230) return 0;
 
  G4double ek = aParticle->GetKineticEnergy()/GeV;      
 
  G4int ie1 = 0;
  G4int ie2 = NFISS - 1;
  do {
    G4int midBin = (ie1 + ie2)/2;
    if (ek < ekfiss[midBin])
      ie2 = midBin;
    else
      ie1 = midBin;
  } while (ie2 - ie1 > 1); 
  G4int i = ie2;
  if (ek < ekfiss[0]) i = 0;
 
  G4int j = 4;
  if (ek <= 0.01) {
    if (ZZ == 92 && AA == 233) j = 1;
    else if (ZZ == 92 && AA == 235) j = 2;
    else if (ZZ == 94 && AA == 239) j = 3;
  }
 
  G4double z43ba;
  if (j == 4) {
    z43ba = ZZ*g4pow->Z13(ZZ)/G4double(AA);
    z43ba = std::max(-67. + 38.7*z43ba, 0.);
  } else {
    z43ba = 1.;
  }
  j = j - 1;      // For array indexing
 
  G4double sigfiss = csfiss[j][i]*z43ba;
  sigfiss = sigfiss*millibarn;
  return sigfiss;
}

Referenced by G4HadronFissionDataSet::GetElementCrossSection().

◆ GetInelasticCrossSection()

G4double G4HadronCrossSections::GetInelasticCrossSection	(	const G4DynamicParticle *	particle,
		G4int	ZZ,
		G4int	AA
	)

Definition at line 1243 of file G4HadronCrossSections.cc.

{
  if (particle->GetDefinition() != prevParticleDefinition ||
      particle->GetKineticEnergy() != prevKineticEnergy ||
      ZZ != prevZZ || AA != prevAA)
    { CalcScatteringCrossSections(particle, ZZ, AA); }
  if (verboseLevel > 1) {
    G4cout << "        x-section inelastic(mb)= " 
       << siginelastic/millibarn << G4endl; 
  }
 
  return siginelastic;
}

Referenced by G4HadronInelasticDataSet::GetElementCrossSection().

◆ GetVerboseLevel()

G4int G4HadronCrossSections::GetVerboseLevel ( )

inline

Definition at line 94 of file G4HadronCrossSections.hh.

94{return verboseLevel;}

◆ Instance()

G4HadronCrossSections * G4HadronCrossSections::Instance ( )

static

Definition at line 1228 of file G4HadronCrossSections.cc.

{
  if (!theInstance) {
    static G4HadronCrossSections xsection;
    theInstance = &xsection;
  }
  return theInstance;
}

Referenced by G4HadronCaptureDataSet::G4HadronCaptureDataSet(), G4HadronElasticDataSet::G4HadronElasticDataSet(), G4HadronFissionDataSet::G4HadronFissionDataSet(), and G4HadronInelasticDataSet::G4HadronInelasticDataSet().

◆ IsApplicable()

G4bool G4HadronCrossSections::IsApplicable ( const G4DynamicParticle * aParticle )

Definition at line 1237 of file G4HadronCrossSections.cc.

{
  return (GetParticleCode(aParticle) > 0);
}

Referenced by G4HadronElasticDataSet::IsElementApplicable(), G4HadronFissionDataSet::IsElementApplicable(), and G4HadronInelasticDataSet::IsElementApplicable().

◆ SetCorrectInelasticNearZero()

static void G4HadronCrossSections::SetCorrectInelasticNearZero ( G4bool value )

inlinestatic

Definition at line 86 of file G4HadronCrossSections.hh.

87 {correctInelasticNearZero = value;}

◆ SetVerboseLevel()

void G4HadronCrossSections::SetVerboseLevel ( G4int value )