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 | |
---|