G4HETCAlpha Class Reference

#include <G4HETCAlpha.hh>

Inheritance diagram for G4HETCAlpha:

Public Member Functions
	G4HETCAlpha ()
	~G4HETCAlpha ()
Public Member Functions inherited from G4HETCChargedFragment
	G4HETCChargedFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4HETCChargedFragment ()
virtual G4double	GetKineticEnergy (const G4Fragment &aFragment)
Public Member Functions inherited from G4HETCFragment
	G4HETCFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4HETCFragment ()
G4double	CalcEmissionProbability (const G4Fragment &aFragment)
Public Member Functions inherited from G4VPreCompoundFragment
	G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
virtual	~G4VPreCompoundFragment ()
void	Initialize (const G4Fragment &aFragment)
G4ReactionProduct *	GetReactionProduct () const
G4int	GetA () const
G4int	GetZ () const
G4int	GetRestA () const
G4int	GetRestZ () const
G4double	ResidualA13 () const
G4double	GetCoulombBarrier () const
G4double	GetBindingEnergy () const
G4double	GetMaximalKineticEnergy () const
G4double	GetEnergyThreshold () const
G4double	GetEmissionProbability () const
G4double	GetNuclearMass () const
G4double	GetRestNuclearMass () const
G4double	GetReducedMass () const
const G4LorentzVector &	GetMomentum () const
void	SetMomentum (const G4LorentzVector &value)
const G4String	GetName () const
void	SetOPTxs (G4int)
void	UseSICB (G4bool)

Protected Member Functions
virtual G4double	GetAlpha ()
virtual G4double	GetBeta ()
virtual G4double	GetSpinFactor ()
virtual G4double	K (const G4Fragment &aFragment)
Protected Member Functions inherited from G4HETCFragment
G4double	BetaRand (const G4int N, const G4int L) const
Protected Member Functions inherited from G4VPreCompoundFragment
G4bool	IsItPossible (const G4Fragment &aFragment) const

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
G4PreCompoundParameters *	theParameters
G4Pow *	g4pow
G4double	theEmissionProbability
G4double	theCoulombBarrier
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 41 of file G4HETCAlpha.hh.

Constructor & Destructor Documentation

G4HETCAlpha::G4HETCAlpha

(

)

Definition at line 37 of file G4HETCAlpha.cc.

  : G4HETCChargedFragment(G4Alpha::Alpha(), &theAlphaCoulombBarrier)
{}

G4HETCAlpha::~G4HETCAlpha

(

)

Definition at line 41 of file G4HETCAlpha.cc.

{}

Member Function Documentation

G4double G4HETCAlpha::GetAlpha ( )

protectedvirtual

Implements G4HETCFragment.

Definition at line 44 of file G4HETCAlpha.cc.

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

{
  G4double C = 0.0;
  G4int aZ = GetZ() + GetRestZ();
  if (aZ <= 30) 
    {
      C = 0.10;
    } 
  else if (aZ <= 50) 
    {
      C = 0.1 + -((aZ-50.)/20.)*0.02;
    }
  else if (aZ < 70) 
    {
      C = 0.08 + -((aZ-70.)/20.)*0.02;
    }
  else 
    {
      C = 0.06;
    }
  return 1.0+C;
}

G4double G4HETCAlpha::GetBeta ( )

protectedvirtual

Implements G4HETCFragment.

Definition at line 67 of file G4HETCAlpha.cc.

References G4VPreCompoundFragment::GetCoulombBarrier().

{
  return -GetCoulombBarrier();
}

G4double G4HETCAlpha::GetSpinFactor ( )

protectedvirtual

Implements G4HETCFragment.

Definition at line 72 of file G4HETCAlpha.cc.

{
  return 1.0;
}

G4double G4HETCAlpha::K ( const G4Fragment &  aFragment )

protectedvirtual

Implements G4HETCFragment.

Definition at line 77 of file G4HETCAlpha.cc.

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

{
  // 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 > 3)
    {
      result = 3.0/(P*(P-1.0)*(P-2.0)*(P-3.0)) * 
    (H*(H-1.0)*(H-2.0)*(H-3.0)*r*r*(r-1.0)*(r-1.0) +
     2.0*H*(H-1.0)*(H-2.0)*(Pa*r*(1.0-r)*(1.0-r)+Na*r*r*(1.0-r)) +
     H*(H-1.0)*(Pa*(Pa-1.0)*(1.0-r)*(1.0-r)+4.0*Na*Pa*r*(1.0-r)+Na*(Na-1.0)*r*r) +
     2*H*(Pa*Na*(Na-1.0)*r+Pa*(Pa-1.0)*Na*(1.0-r)) +
     Pa*(Pa-1.0)*Na*(Na-1.0));

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

---

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

• G4HETCAlpha.hh
• G4HETCAlpha.cc