option,echo TITLE,'Test run for MAD-8 package' ! Warning: ! This is not a real machine, do not try to make it in hardware ! Define constants for field errors set,gcutr,3.0 set, b1r, 1.e-2 set, b2r, 2.e-2 set, b3r, 3.e-2 set, b4r, 4.e-2 set, b5r, 5.e-2 set, b6r, 6.e-2 set, b7r, 7.e-2 set, b8r, 8.e-2 set, b9r, 9.e-2 set, b10r, 10.e-2 set, b11r, 11.e-2 set, a1r, 1.e-2 set, a2r, 2.e-2 set, a3r, 3.e-2 set, a4r, 4.e-2 set, a5r, 5.e-2 set, a6r, 6.e-2 set, a7r, 7.e-2 set, a8r, 8.e-2 set, a9r, 9.e-2 set, a10r, 10.e-2 set, a11r, 11.e-2 ! define beam for beam-beam elements Beam, particle = electron, energy = 50.0 ,& exn = 6.88E-6*4 ,& eyn = 3.75E-6*4 ,& npart=1.05e11, sige= 4.5e-4 ! define elements q1: quad l=0.34 k1=-1.58789 q2: quad l=0.50 k1=2.2453 q3: quad l=0.34 k1=-0.896051 q4: quad l=0.50 k1=1.80457 q5: quad l=0.17 k1=-1.11814 bd: sben l=1.44/2 angle=0.261799/2 qa: mult k1l=-0.4297*1.44 d1: drif l=3.076-0.3/24 d2: drif l=0.29 d3: drif l=0.21 d4: drif l=0.475 d5: drif l=2.094 d6: drif l=0.31 s0: mult k0l=1.e-4 k1l=4.e-6 k2l=3.e-8 k3l=4.e-10 lrad=10. s1: mult k1l=1.e-6 t1 s2: mult k2l=1.e-8 t2 s3: mult k3l=1.e-10 t3 rfc : rfcavity, type = aas, l = 0.3/24,volt=1.7/24, & shunt = 26.166, harmon = 432, tfill = 90 bb: beambeam sigx=1.e-3 sigy=2.e-3 xma=-3.e-3 yma=-4.e-3 charge=-1. ! define sequence ring: sequence, refer=centre rfc_1: rfc, at=.625e-2 q1_1: q1, at=3.246 s0_1: s0, at=3.706 s1_1: s1, at=3.706 q2_1: q2, at=4.246 q3_1: q3, at=4.876 bb_1: bb, at=5.046 bd_1: bd, at=5.881 qa_1: qa, at=6.241 bd_2: bd, at=6.601 q4_1: q4, at=7.686 s2_1: s2, at=8.226 s3_1: s3, at=10.32 q5_1: q5, at=10.715 q5_2: q5, at=10.885 s3_2: s3, at=11.28 s2_2: s2, at=13.374 q4_2: q4, at=13.914 bd_3: bd, at=14.999 qa_2: qa, at=15.359 bd_4: bd, at=15.719 bb_2: bb, at=16.554 q3_2: q3, at=16.724 q2_2: q2, at=17.354 s1_2: s1, at=17.894 s0_2: s0, at=17.894 q1_2: q1, at=18.354 rfc_2: rfc, at=21.59375 rfc_3: rfc, at=21.60625 q1_3: q1, at=24.846 s0_3: s0, at=25.306 s1_3: s1, at=25.306 q2_3: q2, at=25.846 q3_3: q3, at=26.476 bb_3: bb, at=26.646 bd_5: bd, at=27.481 qa_3: qa, at=27.841 bd_6: bd, at=28.201 q4_3: q4, at=29.286 s2_3: s2, at=29.826 s3_3: s3, at=31.92 q5_3: q5, at=32.315 q5_4: q5, at=32.485 s3_4: s3, at=32.88 s2_4: s2, at=34.974 q4_4: q4, at=35.514 bd_7: bd, at=36.599 qa_4: qa, at=36.959 bd_8: bd, at=37.319 bb_4: bb, at=38.154 q3_4: q3, at=38.324 q2_4: q2, at=38.954 s1_4: s1, at=39.494 s0_4: s0, at=39.494 q1_4: q1, at=39.954 rfc_4: rfc, at=43.19375 rfc_5: rfc, at=43.20625 q1_5: q1, at=46.446 s0_5: s0, at=46.906 s1_5: s1, at=46.906 q2_5: q2, at=47.446 q3_5: q3, at=48.076 bb_5: bb, at=48.246 bd_9: bd, at=49.081 qa_5: qa, at=49.441 bd_10: bd, at=49.801 q4_5: q4, at=50.886 s2_5: s2, at=51.426 s3_5: s3, at=53.52 q5_5: q5, at=53.915 q5_6: q5, at=54.085 s3_6: s3, at=54.48 s2_6: s2, at=56.574 q4_6: q4, at=57.114 bd_11: bd, at=58.199 qa_6: qa, at=58.559 bd_12: bd, at=58.919 bb_6: bb, at=59.754 q3_6: q3, at=59.924 q2_6: q2, at=60.554 s1_6: s1, at=61.094 s0_6: s0, at=61.094 q1_6: q1, at=61.554 rfc_6: rfc, at=64.79375 rfc_7: rfc, at=64.80625 q1_7: q1, at=68.046 s0_7: s0, at=68.506 s1_7: s1, at=68.506 q2_7: q2, at=69.046 q3_7: q3, at=69.676 bb_7: bb, at=69.846 bd_13: bd, at=70.681 qa_7: qa, at=71.041 bd_14: bd, at=71.401 q4_7: q4, at=72.486 s2_7: s2, at=73.026 s3_7: s3, at=75.12 q5_7: q5, at=75.515 q5_8: q5, at=75.685 s3_8: s3, at=76.08 s2_8: s2, at=78.174 q4_8: q4, at=78.714 bd_15: bd, at=79.799 qa_8: qa, at=80.159 bd_16: bd, at=80.519 bb_8: bb, at=81.354 q3_8: q3, at=81.524 q2_8: q2, at=82.154 s1_8: s1, at=82.694 s0_8: s0, at=82.694 q1_8: q1, at=83.154 rfc_8: rfc, at=86.39375 rfc_9: rfc, at=86.40625 q1_9: q1, at=89.646 s0_9: s0, at=90.106 s1_9: s1, at=90.106 q2_9: q2, at=90.646 q3_9: q3, at=91.276 bb_9: bb, at=91.446 bd_17: bd, at=92.281 qa_9: qa, at=92.641 bd_18: bd, at=93.001 q4_9: q4, at=94.086 s2_9: s2, at=94.626 s3_9: s3, at=96.72 q5_9: q5, at=97.115 q5_10: q5, at=97.285 s3_10: s3, at=97.68 s2_10: s2, at=99.774 q4_10: q4, at=1.00314e2 bd_19: bd, at=1.01399e2 qa_10: qa, at=1.01759e2 bd_20: bd, at=1.02119e2 bb_10: bb, at=1.02954e2 q3_10: q3, at=1.03124e2 q2_10: q2, at=1.03754e2 s1_10: s1, at=1.04294e2 s0_10: s0, at=1.04294e2 q1_10: q1, at=1.04754e2 rfc_10: rfc, at=1.0799375e2 rfc_11: rfc, at=1.0800625e2 q1_11: q1, at=1.11246e2 s0_11: s0, at=1.11706e2 s1_11: s1, at=1.11706e2 q2_11: q2, at=1.12246e2 q3_11: q3, at=1.12876e2 bb_11: bb, at=1.13046e2 bd_21: bd, at=1.13881e2 qa_11: qa, at=1.14241e2 bd_22: bd, at=1.14601e2 q4_11: q4, at=1.15686e2 s2_11: s2, at=1.16226e2 s3_11: s3, at=1.1832e2 q5_11: q5, at=1.18715e2 q5_12: q5, at=1.18885e2 s3_12: s3, at=1.1928e2 s2_12: s2, at=1.21374e2 q4_12: q4, at=1.21914e2 bd_23: bd, at=1.22999e2 qa_12: qa, at=1.23359e2 bd_24: bd, at=1.23719e2 bb_12: bb, at=1.24554e2 q3_12: q3, at=1.24724e2 q2_12: q2, at=1.25354e2 s1_12: s1, at=1.25894e2 s0_12: s0, at=1.25894e2 q1_12: q1, at=1.26354e2 rfc_12: rfc, at=1.2959375e2 rfc_13: rfc, at=1.2960625e2 q1_13: q1, at=1.32846e2 s0_13: s0, at=1.33306e2 s1_13: s1, at=1.33306e2 q2_13: q2, at=1.33846e2 q3_13: q3, at=1.34476e2 bb_13: bb, at=1.34646e2 bd_25: bd, at=1.35481e2 qa_13: qa, at=1.35841e2 bd_26: bd, at=1.36201e2 q4_13: q4, at=1.37286e2 s2_13: s2, at=1.37826e2 s3_13: s3, at=1.3992e2 q5_13: q5, at=1.40315e2 q5_14: q5, at=1.40485e2 s3_14: s3, at=1.4088e2 s2_14: s2, at=1.42974e2 q4_14: q4, at=1.43514e2 bd_27: bd, at=1.44599e2 qa_14: qa, at=1.44959e2 bd_28: bd, at=1.45319e2 bb_14: bb, at=1.46154e2 q3_14: q3, at=1.46324e2 q2_14: q2, at=1.46954e2 s1_14: s1, at=1.47494e2 s0_14: s0, at=1.47494e2 q1_14: q1, at=1.47954e2 rfc_14: rfc, at=1.5119375e2 rfc_15: rfc, at=1.5120625e2 q1_15: q1, at=1.54446e2 s0_15: s0, at=1.54906e2 s1_15: s1, at=1.54906e2 q2_15: q2, at=1.55446e2 q3_15: q3, at=1.56076e2 bb_15: bb, at=1.56246e2 bd_29: bd, at=1.57081e2 qa_15: qa, at=1.57441e2 bd_30: bd, at=1.57801e2 q4_15: q4, at=1.58886e2 s2_15: s2, at=1.59426e2 s3_15: s3, at=1.6152e2 q5_15: q5, at=1.61915e2 q5_16: q5, at=1.62085e2 s3_16: s3, at=1.6248e2 s2_16: s2, at=1.64574e2 q4_16: q4, at=1.65114e2 bd_31: bd, at=1.66199e2 qa_16: qa, at=1.66559e2 bd_32: bd, at=1.66919e2 bb_16: bb, at=1.67754e2 q3_16: q3, at=1.67924e2 q2_16: q2, at=1.68554e2 s1_16: s1, at=1.69094e2 s0_16: s0, at=1.69094e2 q1_16: q1, at=1.69554e2 rfc_16: rfc, at=1.7279375e2 rfc_17: rfc, at=1.7280625e2 q1_17: q1, at=1.76046e2 s0_17: s0, at=1.76506e2 s1_17: s1, at=1.76506e2 q2_17: q2, at=1.77046e2 q3_17: q3, at=1.77676e2 bb_17: bb, at=1.77846e2 bd_33: bd, at=1.78681e2 qa_17: qa, at=1.79041e2 bd_34: bd, at=1.79401e2 q4_17: q4, at=1.80486e2 s2_17: s2, at=1.81026e2 s3_17: s3, at=1.8312e2 q5_17: q5, at=1.83515e2 q5_18: q5, at=1.83685e2 s3_18: s3, at=1.8408e2 s2_18: s2, at=1.86174e2 q4_18: q4, at=1.86714e2 bd_35: bd, at=1.87799e2 qa_18: qa, at=1.88159e2 bd_36: bd, at=1.88519e2 bb_18: bb, at=1.89354e2 q3_18: q3, at=1.89524e2 q2_18: q2, at=1.90154e2 s1_18: s1, at=1.90694e2 s0_18: s0, at=1.90694e2 q1_18: q1, at=1.91154e2 rfc_18: rfc, at=1.9439375e2 rfc_19: rfc, at=1.9440625e2 q1_19: q1, at=1.97646e2 s0_19: s0, at=1.98106e2 s1_19: s1, at=1.98106e2 q2_19: q2, at=1.98646e2 q3_19: q3, at=1.99276e2 bb_19: bb, at=1.99446e2 bd_37: bd, at=2.00281e2 qa_19: qa, at=2.00641e2 bd_38: bd, at=2.01001e2 q4_19: q4, at=2.02086e2 s2_19: s2, at=2.02626e2 s3_19: s3, at=2.0472e2 q5_19: q5, at=2.05115e2 q5_20: q5, at=2.05285e2 s3_20: s3, at=2.0568e2 s2_20: s2, at=2.07774e2 q4_20: q4, at=2.08314e2 bd_39: bd, at=2.09399e2 qa_20: qa, at=2.09759e2 bd_40: bd, at=2.10119e2 bb_20: bb, at=2.10954e2 q3_20: q3, at=2.11124e2 q2_20: q2, at=2.11754e2 s1_20: s1, at=2.12294e2 s0_20: s0, at=2.12294e2 q1_20: q1, at=2.12754e2 rfc_20: rfc, at=2.1599375e2 rfc_21: rfc, at=2.1600625e2 q1_21: q1, at=2.19246e2 s0_21: s0, at=2.19706e2 s1_21: s1, at=2.19706e2 q2_21: q2, at=2.20246e2 q3_21: q3, at=2.20876e2 bb_21: bb, at=2.21046e2 bd_41: bd, at=2.21881e2 qa_21: qa, at=2.22241e2 bd_42: bd, at=2.22601e2 q4_21: q4, at=2.23686e2 s2_21: s2, at=2.24226e2 s3_21: s3, at=2.2632e2 q5_21: q5, at=2.26715e2 q5_22: q5, at=2.26885e2 s3_22: s3, at=2.2728e2 s2_22: s2, at=2.29374e2 q4_22: q4, at=2.29914e2 bd_43: bd, at=2.30999e2 qa_22: qa, at=2.31359e2 bd_44: bd, at=2.31719e2 bb_22: bb, at=2.32554e2 q3_22: q3, at=2.32724e2 q2_22: q2, at=2.33354e2 s1_22: s1, at=2.33894e2 s0_22: s0, at=2.33894e2 q1_22: q1, at=2.34354e2 rfc_22: rfc, at=2.3759375e2 rfc_23: rfc, at=2.3760625e2 q1_23: q1, at=2.40846e2 s0_23: s0, at=2.41306e2 s1_23: s1, at=2.41306e2 q2_23: q2, at=2.41846e2 q3_23: q3, at=2.42476e2 bb_23: bb, at=2.42646e2 bd_45: bd, at=2.43481e2 qa_23: qa, at=2.43841e2 bd_46: bd, at=2.44201e2 q4_23: q4, at=2.45286e2 s2_23: s2, at=2.45826e2 s3_23: s3, at=2.4792e2 q5_23: q5, at=2.48315e2 q5_24: q5, at=2.48485e2 s3_24: s3, at=2.4888e2 s2_24: s2, at=2.50974e2 q4_24: q4, at=2.51514e2 bd_47: bd, at=2.52599e2 qa_24: qa, at=2.52959e2 bd_48: bd, at=2.53319e2 bb_24: bb, at=2.54154e2 q3_24: q3, at=2.54324e2 q2_24: q2, at=2.54954e2 s1_24: s1, at=2.55494e2 s0_24: s0, at=2.55494e2 q1_24: q1, at=2.55954e2 rfc_24: rfc, at=2.5919375e2 end: marker, at=2.592e2 endsequence ! use part of the machine, make first plot use ring range=#s/q1_2 twiss save plot haxis=s vaxis=betx,bety colour=100 ! use full machine, add alignment and field errors use ring select,flag=error,class=quadrupole ealign dx=tgauss(2.5)*0.1e-4,dy=tgauss(2.5)*0.1e-4 select, error, clear; select, error, class= s0 efcomp, radius = 0.0010, order= 0, & dblsr= 1e-4* a1r , & dblnr= 1e-4* b1r , & dklsr(1)= 1e-4* a2r , & dklnr(1)= 1e-4* b2r , & dklsr(2)= 1e-4* a3r , & dklnr(2)= 1e-4* b3r , & dklsr(3)= 1e-4* a4r , & dklnr(3)= 1e-4* b4r , & dklsr(4)= 1e-4* a5r , & dklnr(4)= 1e-4* b5r , & dklsr(5)= 1e-4* a6r , & dklnr(5)= 1e-4* b6r , & dklsr(6)= 1e-4* a7r , & dklnr(6)= 1e-4* b7r , & dklsr(7)= 1e-4* a8r , & dklnr(7)= 1e-4* b8r , & dklsr(8)= 1e-4* a9r , & dklnr(8)= 1e-4* b9r , & dklsr(9)= 1e-4* a10r, & dklnr(9)= 1e-4* b10r, & dklsr(10)= 1e-4* a11r, & dklnr(10)= 1e-4* b11r; select, error, clear; select, error, class= s1 efcomp, radius = 0.010, order= 1, & dblsr= 1e-4* a1r , & dblnr= 1e-4* b1r , & dklsr(1)= 1e-4* a2r , & dklnr(1)= 1e-4* b2r , & dklsr(2)= 1e-4* a3r , & dklnr(2)= 1e-4* b3r , & dklsr(3)= 1e-4* a4r , & dklnr(3)= 1e-4* b4r , & dklsr(4)= 1e-4* a5r , & dklnr(4)= 1e-4* b5r , & dklsr(5)= 1e-4* a6r , & dklnr(5)= 1e-4* b6r , & dklsr(6)= 1e-4* a7r , & dklnr(6)= 1e-4* b7r , & dklsr(7)= 1e-4* a8r , & dklnr(7)= 1e-4* b8r , & dklsr(8)= 1e-4* a9r , & dklnr(8)= 1e-4* b9r ! make second plot twiss,couple,save plot haxis=s vaxis=betx,bety colour=100 stop