43 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMicroPartition::copy_constructor meant to not be accessible");
51 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMicroPartition::operator= meant to not be accessible");
78 if (anA == 0 || anA == 1)
82 else if (anA == 2 || anA == 3 || anA == 4)
109 return CoulombEnergy;
170 return PartitionEnergy;
180 if (std::fabs(U + FreeInternalE0 - PartitionEnergy) < 0.003)
return -1.0;
194 while (Da*Db > 0.0 && maxit < 1000)
203 for (
G4int i = 0; i < 1000; i++)
206 if (std::fabs(Ta-Tb) <=
eps)
return Tmid;
208 if (std::fabs(Dmid) < 0.003)
return Tmid;
221 G4cout <<
"G4StatMFMicroPartition::CalcPartitionTemperature: I can't calculate the temperature"
281 ThermalWaveLenght3 = ThermalWaveLenght3*ThermalWaveLenght3*ThermalWaveLenght3;
286 kappa = kappa*kappa*kappa;
295 PartitionEntropy +=
G4Log(ProbDegeneracy) + TranslationalS;
299 G4double exponent = PartitionEntropy-SCompound;
300 if (exponent > 300.0) exponent = 300.0;
309 if (
A > 4) DegFactor = 1.0;
310 else if (
A == 1) DegFactor = 4.0;
311 else if (
A == 2) DegFactor = 3.0;
312 else if (
A == 3) DegFactor = 4.0;
313 else if (
A == 4) DegFactor = 1.0;
320 std::vector<G4int> FragmentsZ;
331 if (Af > 1.5 && Af < 4.5) ZMean = 0.5*Af;
332 else ZMean = Af*Z0/A0;
333 G4double ZDispersion = std::sqrt(Af * MeanT/CC);
340 while (Zf < 0 || Zf > Af);
341 FragmentsZ.push_back(Zf);
344 ZBalance = Z0 - SumZ;
347 while (std::abs(ZBalance) > 1);
348 FragmentsZ[0] += ZBalance;
static const G4double eps
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double G4Log(G4double x)
static constexpr double fermi
static constexpr double MeV
static constexpr double pi
G4GLOB_DLL std::ostream G4cout
static G4Pow * GetInstance()
G4double logZ(G4int Z) const
G4double A13(G4double A) const
G4double Z13(G4int Z) const
G4double Z23(G4int Z) const
void CreateFragment(G4int A, G4int Z)
G4double GetCoulombEnergy(void)
G4double GetDegeneracyFactor(G4int A)
G4bool operator==(const G4StatMFMicroPartition &right) const
G4double CalcPartitionProbability(G4double U, G4double FreeInternalE0, G4double SCompound)
std::vector< G4int > _thePartition
void CoulombFreeEnergy(G4int anA)
G4StatMFMicroPartition & operator=(const G4StatMFMicroPartition &right)
std::vector< G4double > _theCoulombFreeEnergy
G4double InvLevelDensity(G4double Af)
G4double CalcPartitionTemperature(G4double U, G4double FreeInternalE0)
G4bool operator!=(const G4StatMFMicroPartition &right) const
G4double GetPartitionEnergy(G4double T)
G4StatMFChannel * ChooseZ(G4int A0, G4int Z0, G4double MeanT)
static G4double DBetaDT(G4double T)
static G4double GetGamma0()
static G4double Beta(G4double T)
static G4double GetCoulomb()
static G4double GetKappaCoulomb()
ThreeVector shoot(const G4int Ap, const G4int Af)
T max(const T t1, const T t2)
brief Return the largest of the two arguments