G4HadronCrossSections Class Reference

#include <G4HadronCrossSections.hh>

Public Member Functions
	~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 ()

Friends
class	G4ThreadLocalSingleton< G4HadronCrossSections >

Detailed Description

Definition at line 58 of file G4HadronCrossSections.hh.

Constructor & Destructor Documentation

~G4HadronCrossSections()

G4HadronCrossSections::~G4HadronCrossSections

(

)

Definition at line 1208 of file G4HadronCrossSections.cc.

{}

Member Function Documentation

GetCaptureCrossSection()

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

Definition at line 1549 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;
}

Referenced by G4HadronCaptureDataSet::GetElementCrossSection().

GetElasticCrossSection()

G4double G4HadronCrossSections::GetElasticCrossSection	(	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 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 1574 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;
}

GetInelasticCrossSection()

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

Definition at line 1226 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 83 of file G4HadronCrossSections.hh.

{return verboseLevel;}

Instance()

G4HadronCrossSections * G4HadronCrossSections::Instance ( )

static

Definition at line 1211 of file G4HadronCrossSections.cc.

{
  if (!instance) {
    static G4ThreadLocalSingleton<G4HadronCrossSections> inst;
    instance = inst.Instance();
  }
  return instance;
}

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

IsApplicable()

G4bool G4HadronCrossSections::IsApplicable

(

const G4DynamicParticle *

aParticle

)

Definition at line 1220 of file G4HadronCrossSections.cc.

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

SetVerboseLevel()

void G4HadronCrossSections::SetVerboseLevel ( G4int  value )

inline

Definition at line 81 of file G4HadronCrossSections.hh.

{verboseLevel = value;}

Friends And Related Function Documentation

G4ThreadLocalSingleton< G4HadronCrossSections >

friend class G4ThreadLocalSingleton< G4HadronCrossSections >

friend

Definition at line 1 of file G4HadronCrossSections.hh.

---

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

• G4HadronCrossSections.hh
• G4HadronCrossSections.cc