128 std::vector<G4int> ANumbers(A);
132 std::vector<G4int> FragmentsA;
135 for (i = 0; i <
A; i++)
137 for (
G4int j = 0; j < ANumbers[i]; j++) FragmentsA.push_back(i+1);
142 for (
G4int j = 0; j < Multiplicity; j++)
144 G4int FragmentsAMax = 0;
146 for (i = j; i < Multiplicity; i++)
148 if (FragmentsA[i] <= FragmentsAMax) {
continue; }
152 FragmentsAMax = FragmentsA[im];
157 FragmentsA[im] = FragmentsA[j];
158 FragmentsA[j] = FragmentsAMax;
170 std::vector<G4double> AcumMultiplicity;
171 AcumMultiplicity.reserve(A);
173 AcumMultiplicity.push_back((*(
_theClusters.begin()))->GetMeanMultiplicity());
174 for (std::vector<G4VStatMFMacroCluster*>::iterator
it =
_theClusters.begin()+1;
177 AcumMultiplicity.push_back((*it)->GetMeanMultiplicity()+AcumMultiplicity.back());
186 for (i = 0; i <
A; i++) ANumbers[i] = 0;
189 for (i = 0; i <
A; i++) {
190 if (RandNumber < AcumMultiplicity[i]) {
196 ANumbers[ThisOne] = ANumbers[ThisOne]+1;
201 }
while (CheckA > 0);
210 std::vector<G4int> & FragmentsA)
214 std::vector<G4int> FragmentsZ;
218 G4int multiplicity = FragmentsA.size();
223 for (
G4int i = 0; i < multiplicity; i++)
231 FragmentsZ.push_back(1);
232 SumZ += FragmentsZ[i];
234 else FragmentsZ.push_back(0);
240 + 2*CP*g4calc->
Z23(FragmentsA[i]);
242 if (FragmentsA[i] > 1 && FragmentsA[i] < 5) { ZMean = 0.5*FragmentsA[i]; }
254 }
while (z < 0 || z > A);
255 FragmentsZ.push_back(z);
267 while (FragmentsZ[idx] < 1) { ++
idx; }
269 FragmentsZ[
idx] += DeltaZ;
272 for (
G4int i = multiplicity-1; i >= 0; i--)