94 if(U < 10*
eV || 0==N) {
return 0.0; }
107 G4bool ChargedNucleon(
false);
109 ChargedNucleon =
true;
115 if (ChargedNucleon) {
116 RelativeVelocitySqr = 2*RelativeEnergy/CLHEP::proton_mass_c2;
118 RelativeVelocitySqr = 2*RelativeEnergy/CLHEP::neutron_mass_c2;
121 G4double RelativeVelocity = std::sqrt(RelativeVelocitySqr);
125 (10.63/RelativeVelocitySqr - 29.92/RelativeVelocity + 42.9)
129 (34.10/RelativeVelocitySqr - 82.20/RelativeVelocity + 82.2)
137 AveragedXSection = ((Z-1)*ppXSection + (A-Z)*npXSection)/
G4double(A-1);
141 AveragedXSection = ((A-Z-1)*ppXSection + Z*npXSection)/
G4double(A-1);
148 G4double PauliFactor = 1.0 - 1.4*FermiRelRatio;
149 if (FermiRelRatio > 0.5) {
150 G4double x = 2.0 - 1.0/FermiRelRatio;
151 PauliFactor += 0.4*FermiRelRatio*x*x*std::sqrt(x);
154 G4double xx = 2*
r0 + CLHEP::hbarc/(CLHEP::proton_mass_c2*RelativeVelocity);
160 *std::sqrt(2.0*RelativeEnergy/CLHEP::proton_mass_c2)/Vint);
199 /(16*CLHEP::c_light);
245 if( (Ncharged == Npart) ||
251 }
else if ( deltaN > 0 ) {
262 if ( Npart < Ncharged )
G4PreCompoundTransitions()
void SetNumberOfHoles(G4int valueTot, G4int valueP=0)
G4int GetNumberOfParticles() const
G4double GetLevelDensity() const
G4int GetNumberOfHoles() const
virtual G4double CalculateProbability(const G4Fragment &aFragment)
double A(double temperature)
G4double GetFermiEnergy() const
G4DeexPrecoParameters * GetParameters()
static G4Proton * Proton()
static constexpr double eV
const G4ParticleDefinition * proton
virtual void PerformTransition(G4Fragment &aFragment)
void SetNumberOfParticles(G4int value)
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
virtual ~G4PreCompoundTransitions()
static constexpr double pi
void SetNumberOfCharged(G4int value)
G4double GetTransitionsR0() const
G4int GetNumberOfCharged() const
static G4NuclearLevelData * GetInstance()
static constexpr double millibarn
G4double GetExcitationEnergy() const
static constexpr double pi2