! FODO for MAD-X ! FODO --: with SIS cell length and SIS bending magnets !--------------------- DRIFT SPACES DEFINITION --------------------- dr: drift, L=1.191; !--------------------------------------------------------------------- !-------------------------- BENDING MAGNET --------------------------- bang=0.1; ALPHA := 15 * RADDEG; LL := 150 * RADDEG; PFR := 7.3* RADDEG; FI := 0.7; RB: SBEND,L=LL,ANGLE= ALPHA;!,E1= PFR,E2=PFR,HGAP= 0.035,FINT= FI; ksf = .0198492943e1; ! ------------------- horizontal quadrupole ---------------------------- qf: quadrupole, L= 4.0, K1= 0.05; qd: quadrupole, L= 4.0, K1= -0.05; hk: hkicker, hkick=0e-3; ! hk: hkicker, hkick=0e-3; mscbh: sextupole, l=1.1, k2:=0.0;!ksf; ! mb1: multipole; ! !--------------------------------------------------------------- ! --------------------------- LINE DEFINITION --------------------------- cella: line=(qf,dr,rb,dr,qd,mscbh,mb1,dr,rb,dr,hk); !------------------------------------------------------------------------ ! beam,particle=proton,energy=0.942563; beam,particle=proton,energy=5000.; ! beam,particle=proton,energy=0.942563,sequence=cella; ! beam; use, period=cella; show,beam%cella; show,beam; !---------------MATCHING---------------------------------------- ! match,sequence=cella; ! VARY,NAME=QF[K1],STEP=0.01; ! VARY,NAME=QD[K1],STEP=0.01; ! CONSTRAINT,sequence=cella,RANGE=#E,MUX=0.3525,MUY=0.29; ! MIGRAD,calls=2000,tolerance=1.e-15; ! endmatch; !--------------------------------------------------------- select,flag=makethin,clear; select,flag=makethin,pattern="q*",slice=100; makethin,sequence=cella; use, period=cella; SELECT,FLAG=ERROR,pattern="mb1*"; !value,bang; EFCOMP, DKS:={0,0,0,0,-1.}; select,flag=sectormap,clear; select,flag=sectormap,range=#e; select,flag=twiss,clear; select,flag=twiss,column=name,s,keyword,L,betx,bety,alfx,alfy,dx,dpx; ! twiss,rmatrix,file,sectormap; ! setplot,post=1,ascale=1.5, lscale=1.5,rscale=1.5, sscale=1.5, lwidth=3; ! PLOT, HAXIS=S,VAXIS1=BETX,VAXIS2=BETY,colour=100; TWISS; !sodd, distort1,multipole_order_range=-8,8,start_stop= 0.,27000., ! print_at_end; ptc_create_universe; ptc_create_layout,model=1,method=2,nst=1;!,exact; ! select_ptc_normal, dx=1, dpx=1; ! select_ptc_normal, dy=1, dpy=1; ! select_ptc_normal, q1=0, q2=0; ! select_ptc_normal, dq1=1, dq2=1; ! select_ptc_normal, dq1=2, dq2=2; ! select_ptc_normal, anhx=1,0,0; ! select_ptc_normal, anhy=1,0,0; ! select_ptc_normal, anhx=1,0,1; ! select_ptc_normal, anhy=1,0,1; ! select_ptc_normal, anhx=0,0,2; ! select_ptc_normal, anhy=0,0,2; ! select_ptc_normal, anhx=2,0,0; ! select_ptc_normal, anhy=2,0,0; ! select_ptc_normal, anhx=1,1,0; ! select_ptc_normal, anhy=1,1,0; ! select_ptc_normal, haml=1,0,0; ! select_ptc_normal, haml=3,0,0; ! select_ptc_normal, haml=1,2,0; select_ptc_normal, gnfu=-5,0,0; ptc_normal,closed_orbit,normal,maptable,icase=5,no=7; ptc_end; write, table=normal_results,file=ptc_normal_results.tfs; write, table=map_table, file=ptc_map_table.tfs; stop; 2.781078939E-02 1.242652155E-01 1.273392469E-01 4 0 0 0 3.975594422E-01 3.184195825E-01 5.093569874E-01 3 1 0 0 7.929397998E-02 -2.447512708E+00 2.448796846E+00 2 0 2 0 -3.822632583E+00 -3.061683214E+00 4.897593691E+00 2 0 1 1 -2.371110819E+00 6.119145973E-01 2.448796846E+00 2 0 0 2 -2.936297813E+00 -3.919767738E+00 4.897593691E+00 1 1 2 0 -3.677149401E-01 1.255352415E+00 1.308099371E+00 0 0 4 0 3.137025701E+00 4.187726490E+00 5.232397482E+00 0 0 3 1 -3.946908288E+00 3.344118576E+01 3.367329788E+01 5 0 0 0 8.912411114E+01 1.428431572E+02 1.683664894E+02 4 1 0 0 3.177189767E+02 1.115515615E+02 3.367329788E+02 3 2 0 0 3.903183478E+02 -1.006204404E+03 1.079257020E+03 3 0 2 0 -1.142600560E+03 -1.831296486E+03 2.158514040E+03 3 0 1 1 -1.075351803E+03 -9.172903158E+01 1.079257020E+03 3 0 0 2 -9.731118583E+02 -3.088076221E+03 3.237771059E+03 2 1 2 0 -6.109893431E+03 -2.145191830E+03 6.475542119E+03 2 1 1 1 -2.690006778E+03 1.801950323E+03 3.237771059E+03 2 1 0 2 -1.008591393E+03 1.405026800E+03 1.729553961E+03 1 0 4 0 2.079269273E+03 6.598359627E+03 6.918215845E+03 1 0 3 1 9.791356640E+03 3.437758531E+03 1.037732377E+04 1 0 2 2 5.747795990E+03 -3.850266449E+03 6.918215845E+03 1 0 1 3 9.112079845E+01 -1.727151963E+03 1.729553961E+03 1 0 0 4 7.274978916E-02 -2.925112408E-01 3.014222251E-01 4 0 1 0 -1.904937682E-01 -2.335968366E-01 3.014222251E-01 4 0 0 1 -5.043138164E-01 -1.095149898E+00 1.205688900E+00 3 1 1 0 -1.178854180E+00 -2.529595761E-01 1.205688900E+00 3 1 0 1 -1.617368289E+00 -8.092667656E-01 1.808533351E+00 2 2 1 0 -9.200874016E-01 1.699030806E+00 1.932166273E+00 2 0 3 0 2.424551176E+00 5.265068866E+00 5.796498820E+00 2 0 2 1 5.667487575E+00 1.216134513E+00 5.796498820E+00 2 0 1 2 1.457069026E+00 -1.268943010E+00 1.932166273E+00 2 0 0 3 6.879860300E-01 3.802596646E+00 3.864332546E+00 1 1 3 0 1.036759800E+01 5.187533699E+00 1.159299764E+01 1 1 2 1 4.216153600E-01 -4.535882823E-01 6.192752552E-01 0 0 5 0 -5.512630178E-01 -3.046909110E+00 3.096376276E+00 0 0 4 1 -5.538168039E+00 -2.771079021E+00 6.192752552E+00 0 0 3 2 0.07274978916 "HAMS" 4 0 1 0 -0.2925112408 "HAMA" 4 0 1 0 0.3014222251 "HAMC" 4 0 0 1 -0.1904937682 "HAMS" 4 0 0 1 -0.2335968366 "HAMA" 4 0 0 1 0.3014222251 "HAMC" 3 1 1 0 -0.5043138164 "HAMS" 3 1 1 0 -1.095149898 "HAMA" 3 1 1 0 1.2056889 "HAMC" 3 1 0 1 -1.17885418 "HAMS" 3 1 0 1 -0.2529595761 "HAMA" 3 1 0 1 1.2056889 "HAMC" 2 2 1 0 -1.617368289 "HAMS" 2 2 1 0 -0.8092667656 "HAMA" 2 2 1 0 1.808533351 "HAMC" 2 0 3 0 -0.9200874016 "HAMS" 2 0 3 0 1.699030806 "HAMA" 2 0 3 0 1.932166273 "HAMC" 2 0 2 1 2.424551176 "HAMS" 2 0 2 1 5.265068866 "HAMA" 2 0 2 1 5.79649882 "HAMC" 2 0 1 2 5.667487575 "HAMS" 2 0 1 2 1.216134513 "HAMA" 2 0 1 2 5.79649882 "HAMC" 2 0 0 3 1.457069026 "HAMS" 2 0 0 3 -1.26894301 "HAMA" 2 0 0 3 1.932166273 "HAMC" 1 1 3 0 0.68798603 "HAMS" 1 1 3 0 3.802596646 "HAMA" 1 1 3 0 3.864332546 "HAMC" 1 1 2 1 10.367598 "HAMS" 1 1 2 1 5.187533699 "HAMA" 1 1 2 1 11.59299764 "HAMC" 0 0 5 0 0.42161536 "HAMS" 0 0 5 0 -0.4535882823 "HAMA" 0 0 5 0 0.6192752552 "HAMC" 0 0 4 1 -0.5512630178 "HAMS" 0 0 4 1 -3.04690911 "HAMA" 0 0 4 1 3.096376276 "HAMC" 0 0 3 2 -5.538168039 "HAMS" 0 0 3 2 -2.771079021 "HAMA" 0 0 3 2 6.192752552 "HAML" -5 0 0 0 14