Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
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 ()
 
virtual G4double CalculateProbability (const G4Fragment &aFragment)=0
 
virtual void PerformTransition (G4Fragment &aFragment)=0
 
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 34 of file G4GNASHTransitions.hh.

Constructor & Destructor Documentation

◆ G4GNASHTransitions()

G4GNASHTransitions::G4GNASHTransitions ( )

Definition at line 36 of file G4GNASHTransitions.cc.

37{}

◆ ~G4GNASHTransitions()

G4GNASHTransitions::~G4GNASHTransitions ( )
virtual

Definition at line 39 of file G4GNASHTransitions.cc.

40{}

Member Function Documentation

◆ CalculateProbability()

G4double G4GNASHTransitions::CalculateProbability ( const G4Fragment aFragment)
virtual

Implements G4VPreCompoundTransitions.

Definition at line 42 of file G4GNASHTransitions.cc.

44{
45 const G4double k = 135.0 * MeV*MeV*MeV;
46 G4double E = aFragment.GetExcitationEnergy();
47 G4double P = aFragment.GetNumberOfParticles();
48 G4double H = aFragment.GetNumberOfHoles();
49 G4double N = P + H;
50 G4double A = aFragment.GetA();
51
52 G4double theMatrixElement(k*N/(A*A*A*E));
53 G4double x = E/N;
54 if ( x < 2.0*MeV ) theMatrixElement *= x/std::sqrt(14.0*MeV*MeV);
55 else if ( x < 7.0*MeV ) x *= std::sqrt(x/7.0*MeV);
56 else if ( x < 15.0*MeV ) ;
57 else x *= std::sqrt(15.0*MeV/x);
58
59 // gg = (6.0/pi2)*a*A
61
62 G4double Epauli = ((P+1.0)*(P+1.0) + (H+1.0)*(H+1.0) + (P+1.0) - 3.0*(H-1.0))/4.0;
63
64 G4double Probability = gg*gg*gg *(E-Epauli)*(E-Epauli);
65 Probability /= 2.0*(N+1.0)*h_Planck;
66 Probability *= theMatrixElement;
67
68 return Probability;
69}
double G4double
Definition: G4Types.hh:64
G4int GetNumberOfParticles() const
Definition: G4Fragment.hh:305
G4int GetNumberOfHoles() const
Definition: G4Fragment.hh:325
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:235
G4double GetA() const
Definition: G4Fragment.hh:283
static G4PreCompoundParameters * GetAddress()

◆ PerformTransition()

void G4GNASHTransitions::PerformTransition ( G4Fragment aFragment)
virtual

Implements G4VPreCompoundTransitions.

Definition at line 71 of file G4GNASHTransitions.cc.

72{
73 result.SetNumberOfParticles(result.GetNumberOfParticles()+1);
74 result.SetNumberOfHoles(result.GetNumberOfHoles()+1);
75 if (G4UniformRand() <= result.GetZ()/result.GetA())
76 {
77 result.SetNumberOfCharged(result.GetNumberOfCharged()+1);
78 }
79
80 if (result.GetNumberOfParticles() < result.GetNumberOfCharged())
81 {
82 result.SetNumberOfCharged(result.GetNumberOfParticles());
83 }
84}
#define G4UniformRand()
Definition: Randomize.hh:53

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