Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4GNASHTransitions Class Reference

#include <G4GNASHTransitions.hh>

Inheritance diagram for G4GNASHTransitions:
Collaboration 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 ( )
explicit

Definition at line 38 of file G4GNASHTransitions.cc.

39 {
41 }
G4DeexPrecoParameters * GetParameters()
static G4NuclearLevelData * GetInstance()

Here is the call graph for this function:

G4GNASHTransitions::~G4GNASHTransitions ( )
virtual

Definition at line 43 of file G4GNASHTransitions.cc.

44 {}

Member Function Documentation

G4double G4GNASHTransitions::CalculateProbability ( const G4Fragment aFragment)
virtual

Implements G4VPreCompoundTransitions.

Definition at line 47 of file G4GNASHTransitions.cc.

48 {
49  static const G4double k = 135.0 *CLHEP::MeV*CLHEP::MeV*CLHEP::MeV;
50  G4double E = aFragment.GetExcitationEnergy();
51  G4double P = aFragment.GetNumberOfParticles();
52  G4double H = aFragment.GetNumberOfHoles();
53  G4double N = P + H;
54  G4double A = aFragment.GetA_asInt();
55 
56  G4double theMatrixElement(k*N/(A*A*A*E));
57  G4double x = E/(N*CLHEP::MeV);
58  static const G4double xf = std::sqrt(2.0/7.0);
59  if ( x < 2.0) { x *= xf; }
60  else if ( x < 7.0) { x *= std::sqrt(x/7.0); }
61  else if ( x > 15.0){ x *= std::sqrt(15.0/x); }
62  theMatrixElement *= x;
63 
64  G4double gg = (6.0/pi2)*theParameters->GetLevelDensity()*A;
65 
66  G4double Epauli = ((P+1.0)*(P+1.0) + (H+1.0)*(H+1.0) + (P+1.0) - 3.0*(H-1.0))*0.25;
67 
68  G4double Probability = gg*gg*gg *(E-Epauli)*(E-Epauli);
69  Probability *= theMatrixElement/(2.0*(N+1.0)*CLHEP::h_Planck);
70 
71  return Probability;
72 }
const int N
Definition: mixmax.h:43
static constexpr double h_Planck
G4int GetNumberOfParticles() const
Definition: G4Fragment.hh:345
G4double GetLevelDensity() const
static double P[]
G4int GetNumberOfHoles() const
Definition: G4Fragment.hh:365
double A(double temperature)
G4int GetA_asInt() const
Definition: G4Fragment.hh:266
static constexpr double MeV
double G4double
Definition: G4Types.hh:76
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:283
static constexpr double pi2
Definition: G4SIunits.hh:78

Here is the call graph for this function:

void G4GNASHTransitions::PerformTransition ( G4Fragment aFragment)
virtual

Implements G4VPreCompoundTransitions.

Definition at line 74 of file G4GNASHTransitions.cc.

75 {
76  result.SetNumberOfParticles(result.GetNumberOfParticles()+1);
77  result.SetNumberOfHoles(result.GetNumberOfHoles()+1);
78  if (G4UniformRand()*result.GetA_asInt() <= G4double(result.GetZ_asInt()))
79  {
80  result.SetNumberOfCharged(result.GetNumberOfCharged()+1);
81  }
82 
83  if (result.GetNumberOfParticles() < result.GetNumberOfCharged())
84  {
85  result.SetNumberOfCharged(result.GetNumberOfParticles());
86  }
87 }
G4double G4ParticleHPJENDLHEData::G4double result
#define G4UniformRand()
Definition: Randomize.hh:97
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:


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