[430] | 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; |
---|