72 :
A(theA),
Z(theZ), levelDensity(0.1)
124 if(
maxExc < 0.0) {
return prob; }
133 }
else if(
resZ > 20) {
134 C = (0.123482-0.00534691*Z-0.0000610624*(Z*
Z)+5.93719*1
e-7*(Z*Z*Z)+
189 for(
G4int i=0; i<100; ++i) {
193 }
while (exc + resExc >
maxExc);
198 G4cout <<
"### G4GEMChannelVI::EmittedFragment WARNING: majoranta "
199 <<
maxProb <<
" is exceeded " << prob <<
"\n"
201 <<
" Z= " <<
Z <<
" A= " <<
A
202 <<
" resZ= " <<
resZ <<
" resA= " <<
resA <<
"\n"
203 <<
" exc(MeV)= " << exc <<
" resExc(MeV)= " << resExc
226 p1 = std::sqrt((e1 - mass1)*(e1 + mass1));
234 lv1.
boost(boostVector);
244 lv0.
boost(boostVector);
284 if(s0 > 350.) { s0 = 350.; }
288 static const G4double sqrt2 = std::sqrt(2.0);
315 ild =
G4Exp((exc-E0CN)/TCN)/TCN;
319 ild =
G4Exp(2*x1)/(x*std::sqrt(x1));
322 Width *= (Rb*Rb/ild);
338 pr = (exc - e1 <= e2 - exc) ? 1.0 - (1.0 - pr)*2*(exc - e1)/(e2 - e1) :
339 2*pr*(e2 - exc)/(e2 - e1);
357 G4double p1 = S *(2.0 + S2 *( 4.0 + S2 *( 13.5 + S2 *( 60.0 + S2 * 325.125 ))));
359 + Sx2 *((1.5*s2+0.5*sx2)
360 + Sx2 *((3.75*s2+0.25*sx2)
361 + Sx2 *((12.875*s2+0.625*sx2)
362 + Sx2 *((59.0625*s2+0.9375*sx2)
363 + Sx2 *(324.8*s2+3.28*sx2))))));
G4double FindLevel(const G4LevelManager *, G4double exc, G4double exclim)
const G4LevelManager * GetLevelManager(G4int Z, G4int A)
G4GEMChannelVI(G4int theA, G4int theZ)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
const G4LorentzVector & GetMomentum() const
G4double I1(G4double t, G4double tx)
G4double MaxLevelEnergy() const
G4double Z13(G4int Z) const
G4DeexPrecoParameters * GetParameters()
G4int GetMinA(G4int Z) const
Float_t x1[n_points_granero]
void SetZandA_asInt(G4int Znew, G4int Anew)
static G4NuclearLevelData * GetInstance()
static constexpr double hbarc
virtual ~G4GEMChannelVI()
G4ThreeVector G4RandomDirection()
virtual void Initialise() final
G4double GetGroundStateMass() const
size_t NearestLowEdgeLevelIndex(G4double energy) const
G4double G4Log(G4double x)
const G4LevelManager * levelManager
G4double ProbabilityDistributionFunction(G4double exc, G4double resExc)
const G4PairingCorrection * pairingCorrection
const G4VCoulombBarrier * cBarrier
static constexpr double fermi
static const G4double ws[NPOINTSGEM]
static G4Pow * GetInstance()
G4NuclearLevelData * nData
static constexpr double MeV
virtual void Dump() const
double A(double temperature)
void SetMomentum(const G4LorentzVector &value)
static G4double GetNuclearMass(const G4double A, const G4double Z)
void set(double x, double y, double z, double t)
G4int GetMaxA(G4int Z) const
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus) final
G4double GetLevelDensity() const
G4double IntegratedProbability(G4double exc)
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const =0
G4int SpinTwo(size_t i) const
G4double GetMaxLevelEnergy(G4int Z, G4int A) const
static const G4double xs[NPOINTSGEM]
size_t NumberOfTransitions() const
virtual G4double GetEmissionProbability(G4Fragment *theNucleus) final
G4GLOB_DLL std::ostream G4cout
G4double GetPairingCorrection(G4int A, G4int Z) const
Hep3Vector boostVector() const
G4double I2(G4double s0, G4double sx)
G4PairingCorrection * GetPairingCorrection()
G4double I3(G4double s0, G4double sx)
G4double LevelEnergy(size_t i) const
G4double GetExcitationEnergy() const
static constexpr double pi
HepLorentzVector & boost(double, double, double)