41 #ifndef G4GEMProbability_h
42 #define G4GEMProbability_h 1
163 return Normalization;
169 theCoulombBarrierPtr = aCoulombBarrierStrategy;
176 if (theCoulombBarrierPtr) {
186 inline G4double G4GEMProbability::CCoeficient(
G4int aZ)
const
194 }
else if (aZ > 20) {
195 C=(0.123482-0.00534691*aZ-0.0000610624*aZ*aZ+5.93719*1e-7*aZ*aZ*aZ+
196 1.95687*1e-8*aZ*aZ*aZ*aZ)/
G4double(theA);
232 return std::exp(t) - 1.0;
237 return (t - tx + 1.0)*std::exp(tx) - t - 1.0;
246 G4double p1 = S*S*S*( 1.0 + S*S*( 1.5 + 3.75*S*S) );
247 G4double p2 = Sx*Sx*Sx*( 1.0 + Sx*Sx*( 1.5 + 3.75*Sx*Sx) )*std::exp(sx-
s0);
G4double EmissionProbability(const G4Fragment &fragment, G4double anEnergy)
void SetCoulomBarrier(const G4VCoulombBarrier *aCoulombBarrierStrategy)
G4double CalcAlphaParam(const G4Fragment &) const
G4double GetCoulombBarrier(const G4Fragment &fragment) const
G4int GetZ_asInt(void) const
std::vector< G4double > ExcitLifetimes
G4int GetA_asInt(void) const
G4double Z13(G4int Z) const
std::vector< G4double > ExcitSpins
G4double GetSpin(void) const
G4double GetNormalization(void) const
G4double GetPairingCorrection(G4int A, G4int Z) const
G4double GetZ(void) const
virtual ~G4GEMProbability()
G4double CalcBetaParam(const G4Fragment &) const
std::vector< G4double > ExcitEnergies
G4double GetA(void) const
G4double Z23(G4int Z) const
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const =0
G4double GetExcitationEnergy() const