G4HETCNeutron Class Reference

#include <G4HETCNeutron.hh>

Inheritance diagram for G4HETCNeutron:

Public Member Functions
	G4HETCNeutron ()
	~G4HETCNeutron ()
virtual G4double	GetKineticEnergy (const G4Fragment &aFragment)
Protected Member Functions
virtual G4double	GetAlpha ()
virtual G4double	GetBeta ()
virtual G4double	GetSpinFactor ()
virtual G4double	K (const G4Fragment &aFragment)

---

Detailed Description

Definition at line 41 of file G4HETCNeutron.hh.

---

Constructor & Destructor Documentation

G4HETCNeutron::G4HETCNeutron

(

)

Definition at line 39 of file G4HETCNeutron.cc.

  : G4HETCFragment(G4Neutron::Neutron(), &theNeutronCoulombBarrier)
{}

G4HETCNeutron::~G4HETCNeutron

(

)

Definition at line 43 of file G4HETCNeutron.cc.

{}

---

Member Function Documentation

G4double G4HETCNeutron::GetAlpha

(

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 46 of file G4HETCNeutron.cc.

References G4VPreCompoundFragment::g4pow, G4VPreCompoundFragment::GetRestA(), and G4Pow::Z13().

Referenced by GetBeta().

{
  return 0.76+2.2/g4pow->Z13(GetRestA());
}

G4double G4HETCNeutron::GetBeta

(

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 51 of file G4HETCNeutron.cc.

References G4VPreCompoundFragment::g4pow, GetAlpha(), G4VPreCompoundFragment::GetRestA(), and G4Pow::Z23().

Referenced by GetKineticEnergy().

{
  return (2.12/g4pow->Z23(GetRestA())-0.05)*MeV/GetAlpha();
}

G4double G4HETCNeutron::GetKineticEnergy

(

const G4Fragment &

aFragment

)

[virtual]

Implements G4VPreCompoundFragment.

Definition at line 85 of file G4HETCNeutron.cc.

References G4HETCFragment::BetaRand(), G4UniformRand, G4Fragment::GetA_asInt(), GetBeta(), G4PreCompoundParameters::GetLevelDensity(), G4VPreCompoundFragment::GetMaximalKineticEnergy(), G4Fragment::GetNumberOfHoles(), G4Fragment::GetNumberOfParticles(), and G4VPreCompoundFragment::theParameters.

{
  G4int H = aFragment.GetNumberOfHoles();
  G4int Pb = aFragment.GetNumberOfParticles();
  G4int Nb = Pb + H;
  G4double g0 = (6.0/pi2)*aFragment.GetA_asInt()*theParameters->GetLevelDensity();
  
  G4double Ab = std::max(0.0,G4double(Pb*Pb+H*H+Pb-3*H)/(4.0*g0));
  G4double Emax = GetMaximalKineticEnergy() - Ab;
  
  G4double cut = GetBeta() / (GetBeta()+Emax/G4double(Nb+1));
  G4double x(0.0);
  if (G4UniformRand() <= cut)
    {
      x = BetaRand(Nb,1);
    }
  else 
    {
      x = BetaRand(Nb,2);
    }

  return Emax * (1.0 - x);
}

G4double G4HETCNeutron::GetSpinFactor

(

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 56 of file G4HETCNeutron.cc.

{
  // (2s+1)
  return 2.0;
}

G4double G4HETCNeutron::K

(

const G4Fragment &

aFragment

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 62 of file G4HETCNeutron.cc.

References G4VPreCompoundFragment::GetA(), G4Fragment::GetNumberOfHoles(), G4Fragment::GetNumberOfParticles(), G4VPreCompoundFragment::GetRestA(), G4VPreCompoundFragment::GetRestZ(), and G4VPreCompoundFragment::GetZ().

{
  // Number of protons in emitted fragment
  G4int Pa = GetZ();
  // Number of neutrons in emitted fragment 
  G4int Na = GetA() - Pa;

  G4int TargetZ = GetRestZ();
  G4int TargetA = GetRestA();
  G4double r = G4double(TargetZ)/G4double(TargetA);
  
  G4int P = aFragment.GetNumberOfParticles();
  G4int H = aFragment.GetNumberOfHoles();
  
  G4double result = 0.0;
  if (P > 0)
    {
      result = (H + Na/(1.0-r))/P;
    }
  
  return std::max(0.0,result);
}

---

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

• G4HETCNeutron.hh
• G4HETCNeutron.cc

---