[4] | 1 | %% |
---|
| 2 | %[x, ATIndex, LostBeam] = getturns([.0 0, 0.0, 0, 0.01681, 0.0]', 1000, [1]); |
---|
| 3 | %[x, ATIndex, LostBeam] = getturns([.0 0, 0.0, 0, 0.01677, 0.0]', 1000, [1]); |
---|
| 4 | %[x, ATIndex, LostBeam] = getturns([.0, 0, 0.0, 0, 0.01661, 0.0]', 2000, [1]); |
---|
| 5 | [x, ATIndex, LostBeam] = getturns([1e-6, 0, 1e-6, 0, 0.02571, 0.0]', 1000, [1]); |
---|
| 6 | |
---|
| 7 | if LostBeam |
---|
| 8 | ATIndex |
---|
| 9 | else |
---|
| 10 | |
---|
| 11 | figure(11) |
---|
| 12 | plot(x(1,:,6),x(1,:,5)*100,'k'); grid on; hold on |
---|
| 13 | xlabel('ctau') |
---|
| 14 | ylabel('delta (%)') |
---|
| 15 | |
---|
| 16 | figure(2) |
---|
| 17 | subplot(2,2,1) |
---|
| 18 | plot(x(1,:,1)*1e3); grid on |
---|
| 19 | xlabel('turn number') |
---|
| 20 | ylabel('x (mm)') |
---|
| 21 | |
---|
| 22 | subplot(2,2,2) |
---|
| 23 | plot(x(1,:,3)*1e3); grid on |
---|
| 24 | xlabel('turn number') |
---|
| 25 | ylabel('z (mm)') |
---|
| 26 | |
---|
| 27 | subplot(2,2,3) |
---|
| 28 | plot(x(1,:,5)*100); grid on |
---|
| 29 | xlabel('turn number') |
---|
| 30 | ylabel('delta (%)') |
---|
| 31 | |
---|
| 32 | subplot(2,2,4) |
---|
| 33 | plot(x(1,:,1)*1e3, x(1,:,2)*1e3,'.'); grid on |
---|
| 34 | xlabel('x (mm)') |
---|
| 35 | ylabel('px (mrad)') |
---|
| 36 | |
---|
| 37 | end |
---|
| 38 | |
---|
| 39 | %% Compute alpha to all order |
---|
| 40 | physics_mcf; |
---|
| 41 | |
---|