source: PSPA/madxPSPA/tests/test-twiss-6/test-twiss-6.madx @ 430

Last change on this file since 430 was 430, checked in by touze, 11 years ago

import madx-5.01.00
File size: 3.9 KB
Line 
1! Threader example
2OPTION, -ECHO, -WARN, INFO;
3en:=2;
4beam,energy:=en,ex=1.e-6,ey=1.e-3,sigt=1e-3,sige=1e-3;
5// element definitions
6mb:rbend, l=14.2, angle:=bang;!,tilt=pi/2;
7mq:quadrupole, l:=3.1,apertype=ellipse,aperture={1,2};
8qft:mq, l:=0.31, k1:=kqf;!,tilt=-pi/4;
9qf.1:mq, l:=3.1, k1:=kqf;
10qf.2:mq, l:=3.1, k1:=kqf;
11qf.3:mq, l:=3.1, k1:=kqf;
12qf.4:mq, l:=3.1, k1:=kqf;
13qf.5:mq, l:=3.1, k1:=kqf;
14qd.1:mq, l:=3.1, k1:=kqd;
15qd.2:mq, l:=3.1, k1:=kqd;
16qd.3:mq, l:=3.1, k1:=kqd;
17qd.4:mq, l:=3.1, k1:=kqd;
18qd.5:mq, l:=3.1, k1:=kqd;
19bph:hmonitor, l:=l.bpm;
20bpv:vmonitor, l:=l.bpm;
21cbh:hkicker;
22cbv:vkicker;
23cbh.1:cbh, kick:=acbh1;
24cbh.2:cbh, kick:=acbh2;
25cbh.3:cbh, kick:=acbh3;
26cbh.4:cbh, kick:=acbh4;
27cbh.5:cbh, kick:=acbh5;
28cbv.1:cbv, kick:=acbv1;
29cbv.2:cbv, kick:=acbv2;
30cbv.3:cbv, kick:=acbv3;
31cbv.4:cbv, kick:=acbv4;
32cbv.5:cbv, kick:=acbv5;
33mscbh:sextupole, l:=1.1, k2:=ksf;
34!mscbh:multipole, knl:={0,0,0,ksf},tilt=-pi/8;
35mscbv:sextupole, l:=1.1, k2:=ksd;
36!mscbv:octupole, l:=1.1, k3:=ksd,tilt=-pi/8;
37
38// sequence declaration
39
40fivecell: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;
114end:marker, at=5.346000e+02;
115endsequence;
116
117OPTION, ECHO, WARN, INFO;
118
119// forces and other constants
120
121l.bpm:=.3;
122bang:=.509998807401e-2;
123kqf:=.872651312e-2;
124kqd:=-.872777242e-2;
125ksf:=.0198492943;
126ksd:=-.039621283;
127acbv1:=1.e-3*gauss(1.);
128acbh1:=1.e-3*gauss(1.);
129acbv2:=1.e-3;
130acbh2:=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;
137use, period=fivecell;
138select, flag=twiss, range=#s, column=name,x, px, y, py;
139select, flag=twiss, range=#e, column=name,x, px, y, py;
140!coguess,x=1.e-3;
141!twiss;!,keeporbit;
142acbv1:=1.e-3;
143acbh1:=1.e-3;
144acbv2:=1.e-3;
145acbh2:=1.e-3;
146coguess,x=0;
147twiss, 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;
153option,threader;
154twiss, file = twiss.2;
Note: See TracBrowser for help on using the repository browser.