72 if (E0 <= 0.0) {
return 0.0; }
75 if (E1 <= 0.0) {
return 0.0; }
80 if (rj <0.0 || xs < 0.0) {
return 0.0; }
85 *
g4calc->
powN(g1*E1/(g0*E0),N-2) * g1 / (E0*g0*g0);
G4int GetNumberOfParticles() const
G4int GetNumberOfHoles() const
G4double GetExcitationEnergy() const
G4int GetNumberOfCharged() const
G4double GetLevelDensity(G4int Z, G4int A, G4double U)
G4double powN(G4double x, G4int n) const
G4double CrossSection(G4double ekin) const
virtual G4double GetRj(G4int NumberParticles, G4int NumberCharged) const =0
G4double ProbabilityDistributionFunction(G4double eKin, const G4Fragment &) override
virtual ~G4PreCompoundNucleon()
G4PreCompoundNucleon(const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
G4NuclearLevelData * fNucData
G4double theBindingEnergy
static constexpr double millibarn
static constexpr double pi2
static constexpr double hbarc