G4HETCDeuteron Class Reference

#include <G4HETCDeuteron.hh>

Inheritance diagram for G4HETCDeuteron:

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

---

Detailed Description

Definition at line 42 of file G4HETCDeuteron.hh.

---

Constructor & Destructor Documentation

G4HETCDeuteron::G4HETCDeuteron

(

)

Definition at line 37 of file G4HETCDeuteron.cc.

  : G4HETCChargedFragment(G4Deuteron::Deuteron(), &theDeuteronCoulombBarrier)
{}

G4HETCDeuteron::~G4HETCDeuteron

(

)

Definition at line 41 of file G4HETCDeuteron.cc.

{}

---

Member Function Documentation

G4double G4HETCDeuteron::GetAlpha

(

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 44 of file G4HETCDeuteron.cc.

References G4VPreCompoundFragment::GetRestZ(), and G4VPreCompoundFragment::GetZ().

{
  G4double C = 0.0;
  G4int aZ = GetZ() + GetRestZ();
  if (aZ >= 70) 
    {
      C = 0.10;
    } 
  else 
    {
      C = ((((0.15417e-06*aZ) - 0.29875e-04)*aZ + 0.21071e-02)*aZ - 0.66612e-01)*aZ + 0.98375; 
    }
  return 1.0 + C/2.0;
}

G4double G4HETCDeuteron::GetBeta

(

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 59 of file G4HETCDeuteron.cc.

References G4VPreCompoundFragment::GetCoulombBarrier().

{
  return -GetCoulombBarrier();
}

G4double G4HETCDeuteron::GetSpinFactor

(

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 64 of file G4HETCDeuteron.cc.

{
  // 2s+1
  return 3.0;
}

G4double G4HETCDeuteron::K

(

const G4Fragment &

aFragment

)

[protected, virtual]

Implements G4HETCFragment.

Definition at line 70 of file G4HETCDeuteron.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 > 1)
    {
      result = 2.0* (H*(H-1.0)*r*(r-1.0)+H*(Na*r+Pa*(1.0-r)) + Pa*Na)/(P*(P-1.0));

      result /= r*(1.0 - r);
    }
  return std::max(0.0,result);
}

---

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

• G4HETCDeuteron.hh
• G4HETCDeuteron.cc

---