41 #ifndef G4GEMProbability_h
42 #define G4GEMProbability_h 1
158 }
else if (aZ > 20) {
159 C=(0.123482-0.00534691*aZ-0.0000610624*aZ*aZ+5.93719*1
e-7*aZ*aZ*aZ+
196 return G4Exp(t) - 1.0;
201 return (t - tx + 1.0)*
G4Exp(tx) - t - 1.0;
210 G4double p1 = S*S*S*( 1.0 + S*S*( 1.5 + 3.75*S*
S) );
211 G4double p2 = Sx*Sx*Sx*( 1.0 + Sx*Sx*( 1.5 + 3.75*Sx*Sx) )*
G4Exp(sx-s0);
G4bool operator==(const G4GEMProbability &right) const
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
std::vector< G4double > ExcitSpins
G4double Z13(G4int Z) const
std::vector< G4double > ExcitEnergies
G4double CCoeficient(G4int) const
G4double I2(G4double s0, G4double sx)
G4double CalcProbability(const G4Fragment &fragment, G4double MaximalKineticEnergy, G4double V)
G4PairingCorrection * fPairCorr
virtual ~G4GEMProbability()
G4double CalcBetaParam(const G4Fragment &) const
G4VLevelDensityParameter * theEvapLDPptr
G4double I3(G4double s0, G4double sx)
std::vector< G4double > ExcitLifetimes
G4double CalcAlphaParam(const G4Fragment &) const
virtual G4double EmissionProbability(const G4Fragment &fragment, G4double maxKineticEnergy)
const G4GEMProbability & operator=(const G4GEMProbability &right)
G4bool operator!=(const G4GEMProbability &right) const
G4double GetCoulombBarrier(const G4Fragment &fragment) const
static constexpr double MeV
G4double GetSpin(void) const
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const =0
G4double I1(G4double t, G4double tx)
G4double GetPairingCorrection(G4int A, G4int Z) const
void SetCoulomBarrier(const G4VCoulombBarrier *aCoulombBarrierStrategy)
G4double GetExcitationEnergy() const
const G4VCoulombBarrier * theCoulombBarrierPtr
G4double Z23(G4int Z) const