G4GNASHTransitions Class Reference

#include <G4GNASHTransitions.hh>

Inheritance diagram for G4GNASHTransitions:

Public Member Functions
	G4GNASHTransitions ()
virtual	~G4GNASHTransitions ()
virtual G4double	CalculateProbability (const G4Fragment &aFragment)
virtual void	PerformTransition (G4Fragment &aFragment)
Public Member Functions inherited from G4VPreCompoundTransitions
	G4VPreCompoundTransitions ()
virtual	~G4VPreCompoundTransitions ()
G4double	GetTransitionProb1 () const
G4double	GetTransitionProb2 () const
G4double	GetTransitionProb3 () const
void	UseNGB (G4bool use)
void	UseCEMtr (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundTransitions
G4bool	useNGB
G4bool	useCEMtr
G4double	TransitionProb1
G4double	TransitionProb2
G4double	TransitionProb3

Detailed Description

Definition at line 36 of file G4GNASHTransitions.hh.

Constructor & Destructor Documentation

G4GNASHTransitions::G4GNASHTransitions

(

)

Definition at line 36 of file G4GNASHTransitions.cc.

{
  theParameters = new G4PreCompoundParameters();
}

G4GNASHTransitions::~G4GNASHTransitions ( )

virtual

Definition at line 41 of file G4GNASHTransitions.cc.

{
  delete theParameters;
}

Member Function Documentation

G4double G4GNASHTransitions::CalculateProbability ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundTransitions.

Definition at line 47 of file G4GNASHTransitions.cc.

References G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4PreCompoundParameters::GetLevelDensity(), G4Fragment::GetNumberOfHoles(), G4Fragment::GetNumberOfParticles(), python.hepunit::h_Planck, python.hepunit::MeV, N, python.hepunit::pi2, and test::x.

{
  static const G4double k = 135.0 * MeV*MeV*MeV;
  G4double E = aFragment.GetExcitationEnergy();
  G4double P = aFragment.GetNumberOfParticles();
  G4double H = aFragment.GetNumberOfHoles();
  G4double N = P + H;
  G4double A = aFragment.GetA_asInt();

  G4double theMatrixElement(k*N/(A*A*A*E));
  G4double x = E/N;
  if ( x < 2.0*MeV ) theMatrixElement *= x/std::sqrt(14.0*MeV*MeV);
  else if ( x < 7.0*MeV ) x *= std::sqrt(x/7.0*MeV);
  else if ( x < 15.0*MeV ) ;
  else x *= std::sqrt(15.0*MeV/x);

  // gg = (6.0/pi2)*a*A
  G4double gg =  (6.0/pi2)*theParameters->GetLevelDensity()*A;

  G4double Epauli = ((P+1.0)*(P+1.0) + (H+1.0)*(H+1.0) + (P+1.0) - 3.0*(H-1.0))/4.0;

  G4double Probability = gg*gg*gg *(E-Epauli)*(E-Epauli);
  Probability /= 2.0*(N+1.0)*h_Planck;
  Probability *= theMatrixElement;

  return Probability;
}

void G4GNASHTransitions::PerformTransition ( G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundTransitions.

Definition at line 75 of file G4GNASHTransitions.cc.

References G4UniformRand, G4Fragment::GetA_asInt(), G4Fragment::GetNumberOfCharged(), G4Fragment::GetNumberOfHoles(), G4Fragment::GetNumberOfParticles(), G4Fragment::GetZ_asInt(), G4Fragment::SetNumberOfCharged(), G4Fragment::SetNumberOfHoles(), and G4Fragment::SetNumberOfParticles().

{
  result.SetNumberOfParticles(result.GetNumberOfParticles()+1);
  result.SetNumberOfHoles(result.GetNumberOfHoles()+1);
  if (G4UniformRand()*result.GetA_asInt() <= G4double(result.GetZ_asInt()))
    {
      result.SetNumberOfCharged(result.GetNumberOfCharged()+1);
    }

  if (result.GetNumberOfParticles() < result.GetNumberOfCharged())
    {
      result.SetNumberOfCharged(result.GetNumberOfParticles());
    }
}

---

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

• G4GNASHTransitions.hh
• G4GNASHTransitions.cc