92 if(U < 10*
eV || 0==N) {
return 0.0; }
104 G4bool ChargedNucleon(
false);
106 ChargedNucleon =
true;
112 if (ChargedNucleon) {
113 RelativeVelocitySqr = 2*RelativeEnergy/CLHEP::proton_mass_c2;
115 RelativeVelocitySqr = 2*RelativeEnergy/CLHEP::neutron_mass_c2;
118 G4double RelativeVelocity = std::sqrt(RelativeVelocitySqr);
122 (10.63/RelativeVelocitySqr - 29.92/RelativeVelocity + 42.9)
126 (34.10/RelativeVelocitySqr - 82.20/RelativeVelocity + 82.2)
134 AveragedXSection = ((Z-1)*ppXSection + (A-Z)*npXSection)/
G4double(A-1);
138 AveragedXSection = ((A-Z-1)*ppXSection + Z*npXSection)/
G4double(A-1);
145 G4double PauliFactor = 1.0 - 1.4*FermiRelRatio;
146 if (FermiRelRatio > 0.5) {
148 PauliFactor += 0.4*FermiRelRatio*x*x*std::sqrt(x);
151 G4double xx = 2*
r0 + CLHEP::hbarc/(CLHEP::proton_mass_c2*RelativeVelocity);
157 *std::sqrt(2.0*RelativeEnergy/CLHEP::proton_mass_c2)/Vint);
196 /(16*CLHEP::c_light);
242 if( (Ncharged == Npart) ||
248 }
else if ( deltaN > 0 ) {
259 if ( Npart < Ncharged )
G4PreCompoundTransitions()
G4double GetTransitionsr0()
void SetNumberOfHoles(G4int valueTot, G4int valueP=0)
G4int GetNumberOfParticles() const
static const G4double sixdpi2
G4int GetNumberOfHoles() const
virtual G4double CalculateProbability(const G4Fragment &aFragment)
double A(double temperature)
static G4Proton * Proton()
G4double GetFermiEnergy()
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()
const G4double x[NPOINTSGL]
static const double millibarn
void SetNumberOfCharged(G4int value)
G4int GetNumberOfCharged() const
G4double GetLevelDensity()
G4double GetExcitationEnergy() const