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;
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;
68 G4double Probability = gg*gg*gg *(E-Epauli)*(E-Epauli);
static constexpr double pi2
G4int GetNumberOfParticles() const
G4int GetNumberOfHoles() const
void SetNumberOfCharged(G4int value)
G4double GetExcitationEnergy() const
void SetNumberOfHoles(G4int valueTot, G4int valueP=0)
void SetNumberOfParticles(G4int value)
G4int GetNumberOfCharged() const
G4DeexPrecoParameters * theParameters
virtual G4double CalculateProbability(const G4Fragment &aFragment)
virtual ~G4GNASHTransitions()
virtual void PerformTransition(G4Fragment &aFragment)
G4double GetLevelDensity(G4int Z, G4int A, G4double U)
G4DeexPrecoParameters * GetParameters()
static G4NuclearLevelData * GetInstance()
static constexpr double h_Planck
static constexpr double MeV