55{
60
61 static const G4double sixoverpi2 = 6.0/CLHEP::pi2;
64
65 G4double A0 = (P*P+H*H+P-3*H)/(4.0*g0);
67
69 if (E0 <= 0.0) { return 0.0; }
70
72 if (E1 <= 0.0) { return 0.0; }
73
76
77 if (rj < 0.0 || xs < 0.0) { return 0.0; }
78
79 static const G4double fact = 2*CLHEP::millibarn
80 /(CLHEP::pi2*CLHEP::hbarc*CLHEP::hbarc*CLHEP::hbarc);
82 *
g4calc->
powN(g1*E1/(g0*E0),
N-2) * g1 / (E0*g0*g0);
83
84
85
86
88}
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
G4NuclearLevelData * fNucData
G4double theBindingEnergy
G4double Probability(const G4int A, const G4FermiFragment *f1, const G4FermiFragment *f2, const G4double mass, const G4double exc)