Research/DesignStudy/Meetings/2020/2020-03-10: Inj2.MADX

File Inj2.MADX, 1.9 KB (added by jpastern, 4 years ago)

Injection for Stage 2

Line 
1TITLE, 'LhARA Stage 2 Injection';
2
3k11:=33.6725;
4k21:=-23.186;
5k31:=-29.5275;
6k41:=5.43383;
7k51:=0.897643;
8k61:=-2.31316;
9k71:=-1.44269;
10k81:=-0.255517;
11k91:=0.0287971;
12k101:=1.84299;
13
14m1: marker;
15//Septum
16D1: SBend, L=0.9, angle=0.85;
17o1: Drift, L=0.8;
18Q1: Quadrupole, L=0.1, K1:=k11;
19o2: Drift, L=0.2;
20Q2: Quadrupole, L=0.1, K1:=k21;
21o3: Drift, L=0.4;
22D2: SBend, L=0.7, angle=0.64;
23o4: Drift, L=0.2;
24D3: SBend, L=0.7, angle=0.64;
25o5: Drift, L=0.2;
26D4: SBend, L=0.7, angle=0.64;
27o7: Drift, L=0.2;
28Q3: Quadrupole, L=0.1, K1:=k31;
29o81: Drift, L=0.2;
30Q4: Quadrupole, L=0.1, K1:=k41;
31o82: Drift, L=0.95;
32Q5: Quadrupole, L=0.1, K1:=k51;
33o9: Drift, L=0.2;
34D5: SBend, L=0.7, angle=-0.64;
35//passing through FFA
36o10: Drift, L=1.85;
37Q6: Quadrupole, L=0.1, K1:=k61;
38o10a: Drift, L=0.2;
39Q7: Quadrupole, L=0.1, K1:=k71;
40o10b: Drift, L=0.2;
41D6: SBend, L=0.7, angle=0.64;
42o11: Drift, L=0.2;
43D7: SBend, L=0.7, angle=0.64;
44o12: Drift, L=0.2;
45Q8: Quadrupole, L=0.1, K1:=k81;
46o13: Drift, L=0.2;
47Q9: Quadrupole, L=0.1, K1:=k91;
48o14: Drift, L=0.6;
49Q10: Quadrupole, L=0.1, K1:=k101;
50o15: Drift, L=1.0;
51//switching dipole
52D8: SBend, L=0.7, angle=0.4375;
53//o16: Drift, L=0.2;
54m2: marker;
55
56
57injection0:line=(m1, D1, o1, Q1, o2, Q2, o3, D2, o4, D3, o5, D4, o7, Q3, o81, Q4, o82, Q5, o9, D5, o10, Q6, o10a, Q7, o10b, D6, o11, D7, o12, Q8, o13, Q9, o14, Q10, o15, D8, m2);
58injection: line=(-injection0);
59
60Beam, particle = proton, sequence = injection, energy = 450,
61      NPART=1.05E11, sige=4.5e-4 ;
62
63
64//beam, particle=proton, sequence=injection, energy = 0.95327231, NPART=1.05E11, sige=4.5e-4, ex=4.08e-8, ey=4.08e-8;
65use, sequence = injection;
66
67
68select, flag = twiss, column=name,s,betx,alfx,mux,
69         bety,alfy, muy,dx,dpx,dy;
70twiss, betx=50,alfx=0.0,bety=50.0, alfy=0.0, Dx=0.0, DPx=0.0, save, centre, file=twissnew1i.out;
71
72plot,file=CCAPi, haxis=s, vaxis=betx,bety,dx, colour=100;
73plot,file=CCAPid, haxis=s, vaxis=dx, colour=100;
74
75
76
77