G4HETCTriton Class Reference

#include <G4HETCTriton.hh>

Inheritance diagram for G4HETCTriton:

G4HETCChargedFragment G4HETCFragment G4VPreCompoundFragment

Public Member Functions

 G4HETCTriton ()
 ~G4HETCTriton ()

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 G4HETCTriton.hh.


Constructor & Destructor Documentation

G4HETCTriton::G4HETCTriton (  ) 

Definition at line 37 of file G4HETCTriton.cc.

00038   : G4HETCChargedFragment(G4Triton::Triton(), &theTritonCoulombBarrier)
00039 {}

G4HETCTriton::~G4HETCTriton (  ) 

Definition at line 41 of file G4HETCTriton.cc.

00042 {}


Member Function Documentation

G4double G4HETCTriton::GetAlpha (  )  [protected, virtual]

Implements G4HETCFragment.

Definition at line 44 of file G4HETCTriton.cc.

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

00045 {
00046   G4double C = 0.0;
00047   G4int aZ = GetZ() + GetRestZ();
00048   if (aZ >= 70) 
00049     {
00050       C = 0.10;
00051     } 
00052   else 
00053     {
00054       C = ((((0.15417e-06*aZ) - 0.29875e-04)*aZ + 0.21071e-02)*aZ - 0.66612e-01)*aZ + 0.98375; 
00055     }
00056   
00057   return 1.0 + C/3.0;
00058 }

G4double G4HETCTriton::GetBeta (  )  [protected, virtual]

Implements G4HETCFragment.

Definition at line 60 of file G4HETCTriton.cc.

References G4VPreCompoundFragment::GetCoulombBarrier().

00061 {
00062   return -GetCoulombBarrier();
00063 }

G4double G4HETCTriton::GetSpinFactor (  )  [protected, virtual]

Implements G4HETCFragment.

Definition at line 65 of file G4HETCTriton.cc.

00066 {
00067   // 2s+1
00068   return 2.0;
00069 }

G4double G4HETCTriton::K ( const G4Fragment aFragment  )  [protected, virtual]

Implements G4HETCFragment.

Definition at line 71 of file G4HETCTriton.cc.

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

00072 {
00073   // Number of protons in emitted fragment
00074   G4int Pa = GetZ();
00075   // Number of neutrons in emitted fragment 
00076   G4int Na = GetA() - Pa;
00077 
00078   G4int TargetZ = GetRestZ();
00079   G4int TargetA = GetRestA();
00080   G4double r = G4double(TargetZ)/G4double(TargetA);
00081 
00082   G4int P = aFragment.GetNumberOfParticles();
00083   G4int H = aFragment.GetNumberOfHoles();
00084 
00085   G4double result = 0.0;
00086   if (P > 2)
00087     {
00088       result = 3.0/(P*(P-1.0)*(P-2.0)) * 
00089         (H*(H-1.0)*(H-2.0)*r*(r-1.0)*(r-1.0) +
00090          H*(H-1.0)*(2.0*Na*r*(1.0-r)+Pa*(1.0-r)*(1.0-r)) +
00091          H*(Na*(Na-1.0)*r+2.0*Na*Pa*(1.0-r)) +
00092          Pa*Na*(Na-1.0));
00093 
00094       result /= 3.0*r*(1.0 - r)*(1.0 - r);
00095     }
00096   return std::max(0.0,result);
00097 }


The documentation for this class was generated from the following files:
Generated on Mon May 27 17:52:13 2013 for Geant4 by  doxygen 1.4.7