40 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMicroManager::copy_constructor meant to not be accessable");
48 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMicroManager::operator= meant to not be accessable");
70 Initialize(theFragment,multiplicity,FreeIntE,SCompNuc);
77 if (!_Partition.empty())
79 std::for_each(_Partition.begin(),_Partition.end(),
88 void G4StatMFMicroManager::Initialize(
const G4Fragment & theFragment,
G4int im,
102 _MeanMultiplicity = 0.0;
105 _MeanTemperature = 0.0;
113 G4int FragmentAtomicNumbers[4];
119 FragmentAtomicNumbers[im-1] = A;
120 for (i = 0; i < (im - 1); i++) FragmentAtomicNumbers[i] = 0;
125 while (MakePartition(im,FragmentAtomicNumbers)) {
129 G4double PartitionProbability = 0.0;
133 _Partition.push_back(aPartition);
135 _WW += PartitionProbability;
136 _MeanMultiplicity += im*PartitionProbability;
137 _MeanTemperature += aPartition->
GetTemperature() * PartitionProbability;
138 if (PartitionProbability > 0.0)
139 _MeanEntropy += PartitionProbability * aPartition->
GetEntropy();
151 G4int tmp = ANumbers[l-1] + ANumbers[k-1];
154 if (ANumbers[l-1] > ANumbers[l] || ANumbers[k-2] > ANumbers[k-1]) {
156 ANumbers[k-1] = tmp - 1;
165 _Normalization = Norm;
167 _MeanMultiplicity /= Norm;
168 _MeanTemperature /= Norm;
169 _MeanEntropy /= Norm;
180 for (std::vector<G4StatMFMicroPartition*>::iterator i = _Partition.begin();
181 i != _Partition.end(); ++i)
183 AccumWeight += (*i)->GetProbability();
184 if (RandNumber < AccumWeight)
185 return (*i)->ChooseZ(A0,Z0,MeanT);
189 "G4StatMFMicroCanonical::ChooseChannel: Couldn't find a channel.");
void Normalize(G4double Norm)
G4StatMFChannel * ChooseChannel(G4int A0, G4int Z0, G4double MeanT)
G4double GetEntropy(void)
G4bool operator==(const G4StatMFMicroManager &right) const
G4StatMFMicroManager(const G4Fragment &theFragment, G4int multiplicity, G4double FreeIntE, G4double SCompNuc)
G4double CalcPartitionProbability(G4double U, G4double FreeInternalE0, G4double SCompound)
G4double GetTemperature(void)
G4double GetExcitationEnergy() const
void SetPartitionFragment(G4int anA)
G4bool operator!=(const G4StatMFMicroManager &right) const