G4HadronCrossSections Class Reference

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

Detailed Description

Definition at line 60 of file G4HadronCrossSections.hh.

Constructor & Destructor Documentation

G4HadronCrossSections::G4HadronCrossSections

(

)

Definition at line 1214 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;
  correctInelasticNearZero = false;
  g4pow = G4Pow::GetInstance();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4HadronCrossSections::~G4HadronCrossSections

(

)

Definition at line 1225 of file G4HadronCrossSections.cc.

{}

Member Function Documentation

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 = G4Pow::GetInstance()->powA(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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Here is the call graph for this function:

Here is the caller graph for this function:

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); /* Loop checking, 08.01.2016, W. Pokorski */
  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;
}

Here is the call graph for this function:

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4HadronCrossSections::GetVerboseLevel ( )

inline

Definition at line 93 of file G4HadronCrossSections.hh.

{return verboseLevel;}

G4HadronCrossSections * G4HadronCrossSections::Instance ( void  )

static

Definition at line 1228 of file G4HadronCrossSections.cc.

{
  if (!theInstance) {
    static G4ThreadLocal G4HadronCrossSections *xsection_G4MT_TLS_ = 0 ; if (!xsection_G4MT_TLS_) xsection_G4MT_TLS_ = new  G4HadronCrossSections  ;  G4HadronCrossSections &xsection = *xsection_G4MT_TLS_;
    theInstance = &xsection;
  }
  return theInstance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4HadronCrossSections::IsApplicable

(

const G4DynamicParticle *

aParticle

)

Definition at line 1237 of file G4HadronCrossSections.cc.

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

Here is the caller graph for this function:

void G4HadronCrossSections::SetVerboseLevel ( G4int  value )

inline

Definition at line 91 of file G4HadronCrossSections.hh.

{verboseLevel = value;}

---

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

• source/geant4.10.03.p02/source/processes/hadronic/cross_sections/include/G4HadronCrossSections.hh
• source/geant4.10.03.p02/source/processes/hadronic/cross_sections/src/G4HadronCrossSections.cc