62 theCoulombBarrierPtr(nullptr)
91 for (
size_t i = 0; i <
nn; ++i) {
100 probability +=
width;
136 LevelDensityParameter(ResidualA,ResidualZ,MaximalKineticEnergy+V-delta0);
139 G4double T = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
142 - 1.25*
G4Log(Ux/MeV) + 2.0*std::sqrt(a*Ux));
151 G4double TCN = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
156 G4double t = MaximalKineticEnergy/T;
157 if ( MaximalKineticEnergy < Ex ) {
159 Width = (
I1(t,t)*T + (Beta+V)*
I0(t))/
G4Exp(E0/T);
166 static const G4double sqrt2 = std::sqrt(2.0);
169 G4double s0 = 2.0*std::sqrt(a*(MaximalKineticEnergy-delta0));
170 G4double sx = 2.0*std::sqrt(a*(Ex-delta0));
172 if(s0 > 350.) { s0 = 350.; }
173 Width =
I1(t,tx)*T/expE0T +
I3(s0,sx)*
G4Exp(s0)/(sqrt2*
a);
194 Rb = 1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85;
201 Rb=1.5*(Aj+Ad)*
fermi;
220 - 1.25*
G4Log(UxCN/MeV)
221 + 2.0*std::sqrt(aCN*UxCN));
223 InitialLevelDensity = (
pi/12.0)*
G4Exp((U-E0CN)/TCN)/TCN;
231 InitialLevelDensity = (
pi/12.0)*
G4Exp(2*x1)/(x*std::sqrt(x1));
236 Width *=
pi*gg*GeometricalXS*Alpha/(12.0*InitialLevelDensity);
250 G4double p1 = S *(2.0 + S2 *( 4.0 + S2 *( 13.5 + S2 *( 60.0 + S2 * 325.125 ))));
253 (1.5*s2+0.5*sx2) + Sx2 *(
254 (3.75*s2+0.25*sx2) + Sx2 *(
255 (12.875*s2+0.625*sx2) + Sx2 *(
256 (59.0625*s2+0.9375*sx2) + Sx2 *(324.8*s2+3.28*sx2))))));
272 G4cout <<
"GEM: List of Excited States for Isotope Z= "
273 <<
theZ <<
" A= " <<
theA <<
" Nlevels= " << nlev
274 <<
" Efermi(MeV)= " << efermi
276 for(
G4int i=0; i< nlev; ++i) {
278 <<
" Mass(GeV)= " << mass/
GeV
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
std::vector< G4double > ExcitSpins
G4double logZ(G4int Z) const
std::vector< ExP01TrackerHit * > a
G4double Z13(G4int Z) const
static constexpr double hbar_Planck
static constexpr double MeV
std::vector< G4double > ExcitEnergies
virtual G4double LevelDensityParameter(G4int A, G4int Z, G4double U) const =0
Float_t x1[n_points_granero]
static G4NuclearLevelData * GetInstance()
static constexpr double hbarc
G4double CalcProbability(const G4Fragment &fragment, G4double MaximalKineticEnergy, G4double V)
G4PairingCorrection * fPairCorr
virtual ~G4GEMProbability()
G4double CalcBetaParam(const G4Fragment &) const
G4double G4Log(G4double x)
G4VLevelDensityParameter * theEvapLDPptr
static constexpr double neutron_mass_c2
G4double I3(G4double s0, G4double sx)
std::vector< G4double > ExcitLifetimes
static constexpr double fermi
G4double CalcAlphaParam(const G4Fragment &) const
static G4Pow * GetInstance()
virtual G4double EmissionProbability(const G4Fragment &fragment, G4double maxKineticEnergy)
G4double ComputeGroundStateMass(G4int Z, G4int A) const
G4double GetCoulombBarrier(const G4Fragment &fragment) const
double A(double temperature)
static G4double GetNuclearMass(const G4double A, const G4double Z)
static constexpr double pi2
G4GLOB_DLL std::ostream G4cout
G4double I1(G4double t, G4double tx)
G4double GetPairingCorrection(G4int A, G4int Z) const
static constexpr double pi
G4PairingCorrection * GetPairingCorrection()
G4double GetExcitationEnergy() const
static constexpr double GeV