imagefilledrectangle($img[$imagenum], 0, 0, $plotwidth - 1, $plotheight - 1, $color[11]);
// LEGEND Parameters
$legendx = $plotwidth - 60;
$legendy = $plotmarginother + 20;
echo "<div class=\"t2\">\n";
echo $A;
echo $Element[$Z];
echo "</div><div class=\"t3\">\n";
echo $i;
echo "</div><div class=\"t4\">\n";
$BeamEnergy = 1.4;
$X = 3.86 * (sqrt($BeamEnergy * $A / $A) / pow($Z, 0.45));
$meanQ = $Z * pow(1 + pow($X, -1 / 0.6), -0.6);
$s = 0.5 * pow($meanQ * (1 - pow($meanQ / $Z, 1 / 0.6)), 0.5);
// for Z>20
printf("%.1f", charge_fraction($i, $s, $meanQ));
echo "%";
echo "<BR>This A/Q = ";
if ($EnergyLoss > 0.0) {
$AQWindow2 = ($beammass - $i * $emass) / $i * sqrt(1.0 - $EnergyLoss / 100);
printf("%.3f", $AQWindow2);
} else {
$AQWindow2 = ($beammass - $i * $emass) / $i;
printf("%.3f", $AQWindow2);
}
echo "<BR>First A/Q = ";
$AQWindow1 = ($beammass - $Q * $emass) / $Q;
printf("%.3f", $AQWindow1);
echo "<BR>";
echo $A, $Element[$Z], $Q;
echo "+";
$thismass = getthismass($cocktailA[$j], $cocktailZ[$j]);
if (($thismass - $k * $emass) / $k > ($beammass - $Q * $emass) / $Q - ($beammass - $Q * $emass) / $Q * (0.5 / $Resolve) && ($thismass - $k * $emass) / $k < ($beammass - $Q * $emass) / $Q + ($beammass - $Q * $emass) / $Q * (0.5 / $Resolve)) {
for ($ii = 0; $ii < sizeof($BadZ[$Liner]); $ii++) {
if ($cocktailZ[$j] == $BadZ[$Liner][$ii]) {
$thismass = getthismass($cocktailA[$j], $cocktailZ[$j]);
$X = 3.86 * (sqrt($BeamEnergy * $cocktailA[$j] / $cocktailA[$j]) / pow($cocktailZ[$j], 0.45));
$meanQ = $cocktailZ[$j] * pow(1 + pow($X, -1 / 0.6), -0.6);
$s = 0.5 * pow($meanQ * (1 - pow($meanQ / $cocktailZ[$j], 1 / 0.6)), 0.5);
// for Z>20
$c = 2 * ($cocktailZ[$j] - $meanQ + 2) / ($s * $s);
$nu = 2 * ($cocktailZ[$j] - $meanQ + 2) * ($cocktailZ[$j] - $meanQ + 2) / ($s * $s);
printf("<BR><b>%d%s:</b><BR>", $cocktailA[$j], $Element[$cocktailZ[$j]]);
for ($ii = 1; $ii <= $cocktailZ[$j] + 2; $ii++) {
$t = $c * ($Z - $ii + 2);
// $F=frac($nu,$t,$c);
$F = charge_fraction($ii, $s, $meanQ);
if ($F < 0.001) {
continue;
}
echo "<div class=\"t2\">";
echo $ii;
echo "</div>";
echo "<div class=\"t3\">";
if ($F > 0.001) {
printf("%5.3f", $F);
} else {
echo "0.000";
}
echo "</div><div class=\"t4\">";
printf("%5.3f", ($thismass - $ii * $emass) / $ii);
echo "</div><div class=\"t5\">";
