1 | beam, particle='proton', e=0.9468, ex=4.08e-8, ey=4.08e-8;
|
---|
2 |
|
---|
3 | //all lengths in m
|
---|
4 | //drift at the target-capture interface
|
---|
5 | o0: drift, l=0.1;
|
---|
6 | //length representing the distance from the physical start of the Gabor lens to the start of the hard edge focusing area
|
---|
7 | o1: drift, l=0.15;
|
---|
8 |
|
---|
9 | //gl1: solenoid, l=0.857, b=1.2868;
|
---|
10 | !gl1: solenoid, l=0.857, ks=2.290287;
|
---|
11 | //gl2: solenoid, l=0.857, b=0.6671;
|
---|
12 | !gl2: solenoid, l=0.857, ks=1.187325;
|
---|
13 | //drift for the first RF cavity
|
---|
14 | !o2: drift, l=0.5;
|
---|
15 | //gl3: solenoid, l=0.857, b=0.8139;
|
---|
16 | !gl3: solenoid, l=0.857, ks=1.448605;
|
---|
17 | //drift for the energy selection;
|
---|
18 | !o3: drift, l=2.55;
|
---|
19 | //gl4: solenoid, l=0.857, b=1.3555;
|
---|
20 | !gl4: solenoid, l=0.857, ks=2.412561;
|
---|
21 | !o4: drift, l=2.9;
|
---|
22 | //gl5: solenoid, l=0.857, b=1.3351;
|
---|
23 | !gl5: solenoid, l=0.857, ks=2.376253;
|
---|
24 |
|
---|
25 |
|
---|
26 | o0: drift, l=0.10;
|
---|
27 | o1: drift, l=0.15;
|
---|
28 | gl1: solenoid, l=0.857, ks=2.290287; ! anode radius of 3.65cm for all lenses
|
---|
29 | o2: drift, l=0.15;
|
---|
30 | o2a: drift, l=0.15;
|
---|
31 | gl2: solenoid, l=0.857, ks=1.187325;
|
---|
32 | o3: drift, l=0.15;
|
---|
33 | //RF cavity 1
|
---|
34 | s1rf1: drift, l=0.5;
|
---|
35 | o5: drift, l=0.15;
|
---|
36 | gl3: solenoid, l=0.857, ks=1.30567;
|
---|
37 | o6: drift, l=0.15;
|
---|
38 | o7: drift, l=2.55;
|
---|
39 | s1rf2: drift, l=0.5;
|
---|
40 | o8: drift, l=0.15;
|
---|
41 | gl4: solenoid, l=0.857, ks=2.50424;
|
---|
42 | o9: drift, l=0.15;
|
---|
43 | o10: drift, l=0.1;
|
---|
44 | o11: drift, l=0.15;
|
---|
45 | gl5: solenoid, l=0.857, ks=0.411532;
|
---|
46 | o12: drift, l=0.15;
|
---|
47 | //place for octupole now mounted on the flange of gl5
|
---|
48 | oct: drift, l=0.1;
|
---|
49 | //RF cavity 2
|
---|
50 | s1rf2: drift, l=0.5;
|
---|
51 | //the extraction dipole is positioned after the RF2
|
---|
52 |
|
---|
53 | capture: line=(o0,o1,gl1,o2,o2a,gl2,o3);
|
---|
54 | energyselection: line=(s1rf1,o5,gl3,o6,o7);
|
---|
55 | matching: line=(o8,gl4,o9,o10,o11,gl5,o12, oct, s1rf2);
|
---|
56 | ccapv4: line=(capture, energyselection, matching);
|
---|
57 |
|
---|
58 | // Assumed Initial Parameters
|
---|
59 | //init_betx := 1.36e-4;
|
---|
60 | //init_bety := 1.36e-4;
|
---|
61 |
|
---|
62 | //beam, particle=proton, sequence=ccapv6, energy = 0.95327231, NPART=1.05E11, sige=4.5e-4, ex=4.08e-8, ey=4.08e-8;
|
---|
63 | Beam, particle = proton, sequence = injection, energy = 450,
|
---|
64 | NPART=1.05E11, sige=4.5e-4 ;
|
---|
65 | use, sequence = ccapv4;
|
---|
66 |
|
---|
67 | //select, flag = twiss;
|
---|
68 | //twiss, sequence=ccapv4, betx=init_betx,bety=init_bety, save, file=full.twiss;
|
---|
69 |
|
---|
70 | select, flag = twiss, column=name,s,betx,alfx,mux,
|
---|
71 | bety,alfy, muy,dx,dpx,dy;
|
---|
72 | twiss, betx=0.002468,alfx=0.022969,bety=0.002468, alfy=0.022969, Dx=0.0, DPx=0.0, save, centre, file=twissnewi0.out;
|
---|
73 |
|
---|
74 | plot,file=injmatch0, haxis=s, vaxis=betx,bety, colour=100;
|
---|
75 | //plot,file=s1optics, haxis=s, vaxis=alfx,alfy, colour=100;
|
---|
76 |
|
---|