40 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMicroPartition::copy_constructor meant to not be accessable");
48 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMicroPartition::operator= meant to not be accessable");
71 CoulombConstFactor = elm_coupling * (3./5.) *
78 if (anA == 0 || anA == 1)
82 else if (anA == 2 || anA == 3 || anA == 4)
104 ZA*ZA*std::pow(static_cast<G4double>(
_thePartition[i]),5./3.)/
107 return CoulombEnergy;
164 PartitionEnergy += elm_coupling*(3./5.)*
theZ*
theZ*CoulombFactor/
168 return PartitionEnergy;
178 if (std::fabs(U + FreeInternalE0 - PartitionEnergy) < 0.003)
return -1.0;
191 while (Da*Db > 0.0 && maxit < 1000)
200 for (
G4int i = 0; i < 1000; i++)
203 if (std::fabs(Ta-Tb) <=
eps)
return Tmid;
205 if (std::fabs(Dmid) < 0.003)
return Tmid;
218 G4cerr <<
"G4StatMFMicroPartition::CalcPartitionTemperature: I can't calculate the temperature"
279 ThermalWaveLenght3 = ThermalWaveLenght3*ThermalWaveLenght3*ThermalWaveLenght3;
282 G4double kappa = (1. + elm_coupling*(std::pow(static_cast<G4double>(
_thePartition.size()),1./3.)-1.0)
284 kappa = kappa*kappa*kappa;
290 (
_thePartition.size()-1.0)*std::log(FreeVolume/ThermalWaveLenght3) +
293 PartitionEntropy += std::log(ProbDegeneracy) + TranslationalS;
297 G4double exponent = PartitionEntropy-SCompound;
298 if (exponent > 700.0) exponent = 700.0;
307 if (A > 4) DegFactor = 1.0;
308 else if (A == 1) DegFactor = 4.0;
309 else if (A == 2) DegFactor = 3.0;
310 else if (A == 3) DegFactor = 4.0;
311 else if (A == 4) DegFactor = 1.0;
318 std::vector<G4int> FragmentsZ;
329 if (Af > 1.5 && Af < 4.5) ZMean = 0.5*Af;
330 else ZMean = Af*Z0/A0;
331 G4double ZDispersion = std::sqrt(Af * MeanT/CC);
337 while (Zf < 0 || Zf > Af);
338 FragmentsZ.push_back(Zf);
341 ZBalance = Z0 - SumZ;
343 while (std::abs(ZBalance) > 1);
344 FragmentsZ[0] += ZBalance;
G4double InvLevelDensity(G4double Af)
static G4double GetGamma0()
G4bool operator!=(const G4StatMFMicroPartition &right) const
ThreeVector shoot(const G4int Ap, const G4int Af)
G4StatMFMicroPartition & operator=(const G4StatMFMicroPartition &right)
static G4double GetKappaCoulomb()
G4double GetDegeneracyFactor(G4int A)
G4double CalcPartitionTemperature(G4double U, G4double FreeInternalE0)
static const G4double eps
G4bool operator==(const G4StatMFMicroPartition &right) const
void CreateFragment(G4int A, G4int Z)
G4StatMFChannel * ChooseZ(G4int A0, G4int Z0, G4double MeanT)
G4double GetPartitionEnergy(G4double T)
std::vector< G4double > _theCoulombFreeEnergy
static const G4double A[nN]
void CoulombFreeEnergy(G4int anA)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
std::vector< G4int > _thePartition
G4double CalcPartitionProbability(G4double U, G4double FreeInternalE0, G4double SCompound)
static G4double DBetaDT(G4double T)
G4double GetCoulombEnergy(void)
static G4double Beta(G4double T)
static const double fermi
G4GLOB_DLL std::ostream G4cerr