1 | ! Threader example |
---|
2 | OPTION, -ECHO, -WARN, INFO; |
---|
3 | en:=2; |
---|
4 | beam,energy:=en,ex=1.e-6,ey=1.e-3,sigt=1e-3,sige=1e-3; |
---|
5 | // element definitions |
---|
6 | mb:rbend, l=14.2, angle:=bang;!,tilt=pi/2; |
---|
7 | mq:quadrupole, l:=3.1,apertype=ellipse,aperture={1,2}; |
---|
8 | qft:mq, l:=0.31, k1:=kqf;!,tilt=-pi/4; |
---|
9 | qf.1:mq, l:=3.1, k1:=kqf; |
---|
10 | qf.2:mq, l:=3.1, k1:=kqf; |
---|
11 | qf.3:mq, l:=3.1, k1:=kqf; |
---|
12 | qf.4:mq, l:=3.1, k1:=kqf; |
---|
13 | qf.5:mq, l:=3.1, k1:=kqf; |
---|
14 | qd.1:mq, l:=3.1, k1:=kqd; |
---|
15 | qd.2:mq, l:=3.1, k1:=kqd; |
---|
16 | qd.3:mq, l:=3.1, k1:=kqd; |
---|
17 | qd.4:mq, l:=3.1, k1:=kqd; |
---|
18 | qd.5:mq, l:=3.1, k1:=kqd; |
---|
19 | bph:hmonitor, l:=l.bpm; |
---|
20 | bpv:vmonitor, l:=l.bpm; |
---|
21 | cbh:hkicker; |
---|
22 | cbv:vkicker; |
---|
23 | cbh.1:cbh, kick:=acbh1; |
---|
24 | cbh.2:cbh, kick:=acbh2; |
---|
25 | cbh.3:cbh, kick:=acbh3; |
---|
26 | cbh.4:cbh, kick:=acbh4; |
---|
27 | cbh.5:cbh, kick:=acbh5; |
---|
28 | cbv.1:cbv, kick:=acbv1; |
---|
29 | cbv.2:cbv, kick:=acbv2; |
---|
30 | cbv.3:cbv, kick:=acbv3; |
---|
31 | cbv.4:cbv, kick:=acbv4; |
---|
32 | cbv.5:cbv, kick:=acbv5; |
---|
33 | mscbh:sextupole, l:=1.1, k2:=ksf; |
---|
34 | !mscbh:multipole, knl:={0,0,0,ksf},tilt=-pi/8; |
---|
35 | mscbv:sextupole, l:=1.1, k2:=ksd; |
---|
36 | !mscbv:octupole, l:=1.1, k3:=ksd,tilt=-pi/8; |
---|
37 | |
---|
38 | // sequence declaration |
---|
39 | |
---|
40 | fivecell:sequence, refer=centre, l=534.6; |
---|
41 | qf.1:qf.1, at=1.550000e+00; |
---|
42 | qft:qft, at=3.815000e+00; |
---|
43 | ! mscbh:mscbh, at=3.815000e+00; |
---|
44 | cbh.1:cbh.1, at=4.365000e+00; |
---|
45 | cbh,at=4.365000e+00; |
---|
46 | mb:mb, at=1.262000e+01; |
---|
47 | mb:mb, at=2.828000e+01; |
---|
48 | mb:mb, at=4.394000e+01; |
---|
49 | ! bpv:bpv, at=5.246000e+01; |
---|
50 | qd.1:qd.1, at=5.501000e+01; |
---|
51 | mscbv:mscbv, at=5.727500e+01; |
---|
52 | cbv.1:cbv.1, at=5.782500e+01; |
---|
53 | cbv,at=5.782500e+01; |
---|
54 | mb:mb, at=6.608000e+01; |
---|
55 | mb:mb, at=8.174000e+01; |
---|
56 | mb:mb, at=9.740000e+01; |
---|
57 | bph:bph, at=1.059200e+02; |
---|
58 | qf.2:qf.2, at=1.084700e+02; |
---|
59 | mscbh:mscbh, at=1.107350e+02; |
---|
60 | cbh.2:cbh.2, at=1.112850e+02; |
---|
61 | mb:mb, at=1.195400e+02; |
---|
62 | mb:mb, at=1.352000e+02; |
---|
63 | mb:mb, at=1.508600e+02; |
---|
64 | bpv:bpv, at=1.593800e+02; |
---|
65 | qd.2:qd.2, at=1.619300e+02; |
---|
66 | mscbv:mscbv, at=1.641950e+02; |
---|
67 | cbv.2:cbv.2, at=1.647450e+02; |
---|
68 | mb:mb, at=1.730000e+02; |
---|
69 | mb:mb, at=1.886600e+02; |
---|
70 | mb:mb, at=2.043200e+02; |
---|
71 | bph:bph, at=2.128400e+02; |
---|
72 | qf.3:qf.3, at=2.153900e+02; |
---|
73 | mscbh:mscbh, at=2.176550e+02; |
---|
74 | cbh.3:cbh.3, at=2.182050e+02; |
---|
75 | mb:mb, at=2.264600e+02; |
---|
76 | mb:mb, at=2.421200e+02; |
---|
77 | mb:mb, at=2.577800e+02; |
---|
78 | bpv:bpv, at=2.663000e+02; |
---|
79 | qd.3:qd.3, at=2.688500e+02; |
---|
80 | mscbv:mscbv, at=2.711150e+02; |
---|
81 | cbv.3:cbv.3, at=2.716650e+02; |
---|
82 | mb:mb, at=2.799200e+02; |
---|
83 | mb:mb, at=2.955800e+02; |
---|
84 | mb:mb, at=3.112400e+02; |
---|
85 | bph:bph, at=3.197600e+02; |
---|
86 | qf.4:qf.4, at=3.223100e+02; |
---|
87 | mscbh:mscbh, at=3.245750e+02; |
---|
88 | cbh.4:cbh.4, at=3.251250e+02; |
---|
89 | mb:mb, at=3.333800e+02; |
---|
90 | mb:mb, at=3.490400e+02; |
---|
91 | mb:mb, at=3.647000e+02; |
---|
92 | bpv:bpv, at=3.732200e+02; |
---|
93 | qd.4:qd.4, at=3.757700e+02; |
---|
94 | mscbv:mscbv, at=3.780350e+02; |
---|
95 | cbv.4:cbv.4, at=3.785850e+02; |
---|
96 | mb:mb, at=3.868400e+02; |
---|
97 | mb:mb, at=4.025000e+02; |
---|
98 | mb:mb, at=4.181600e+02; |
---|
99 | bph:bph, at=4.266800e+02; |
---|
100 | qf.5:qf.5, at=4.292300e+02; |
---|
101 | mscbh:mscbh, at=4.314950e+02; |
---|
102 | cbh.5:cbh.5, at=4.320450e+02; |
---|
103 | mb:mb, at=4.403000e+02; |
---|
104 | mb:mb, at=4.559600e+02; |
---|
105 | mb:mb, at=4.716200e+02; |
---|
106 | bpv:bpv, at=4.801400e+02; |
---|
107 | qd.5:qd.5, at=4.826900e+02; |
---|
108 | mscbv:mscbv, at=4.849550e+02; |
---|
109 | cbv.5:cbv.5, at=4.855050e+02; |
---|
110 | mb:mb, at=4.937600e+02; |
---|
111 | mb:mb, at=5.094200e+02; |
---|
112 | mb:mb, at=5.250800e+02; |
---|
113 | bph:bph, at=5.336000e+02; |
---|
114 | end:marker, at=5.346000e+02; |
---|
115 | endsequence; |
---|
116 | |
---|
117 | OPTION, ECHO, WARN, INFO; |
---|
118 | |
---|
119 | // forces and other constants |
---|
120 | |
---|
121 | l.bpm:=.3; |
---|
122 | bang:=.509998807401e-2; |
---|
123 | kqf:=.872651312e-2; |
---|
124 | kqd:=-.872777242e-2; |
---|
125 | ksf:=.0198492943; |
---|
126 | ksd:=-.039621283; |
---|
127 | acbv1:=1.e-3*gauss(1.); |
---|
128 | acbh1:=1.e-3*gauss(1.); |
---|
129 | acbv2:=1.e-3; |
---|
130 | acbh2:=1.e-3; |
---|
131 | !select,flag=save,class=variable,pattern="bh.*"; |
---|
132 | !select,flag=save,clear; |
---|
133 | !save,beam,file=save; |
---|
134 | ! action |
---|
135 | |
---|
136 | !threader,vector=4e-3,4e-3; |
---|
137 | use, period=fivecell; |
---|
138 | select, flag=twiss, range=#s, column=name,x, px, y, py; |
---|
139 | select, flag=twiss, range=#e, column=name,x, px, y, py; |
---|
140 | !coguess,x=1.e-3; |
---|
141 | !twiss;!,keeporbit; |
---|
142 | acbv1:=1.e-3; |
---|
143 | acbh1:=1.e-3; |
---|
144 | acbv2:=1.e-3; |
---|
145 | acbh2:=1.e-3; |
---|
146 | coguess,x=0; |
---|
147 | twiss, file = twiss.1;!,keeporbit; |
---|
148 | |
---|
149 | !acbv1:=1.e-4; |
---|
150 | !acbh1:=1.e-4; |
---|
151 | !acbv2:=1.e-4; |
---|
152 | !acbh2:=1.e-4; |
---|
153 | option,threader; |
---|
154 | twiss, file = twiss.2; |
---|