!+ ! display command ! !- ON LOG; ON RAD FLUC COD RFSW RADCOD; FFS; $FORM = "15.10"; PageWidth = 1999; WriteGraphs[]:=( Print["WriteGraphs Called"]; FFS["draw bx & by & ex Q*",12]; FFS["draw nx & ny & epx Q*",13]; FFS["draw geo",14]; FFS["draw dx & dy Q*",15]; ); USE RECCELL; !! after optics roughly fixed !! fit AC1 for R56 ! cell calc emit; ! Emittance[]; ! CONVERGENCE=1E-10; ! fit EX 0; ! fit EPX 0; ! FitFunction:={(MomentumCompaction/.Emittance[])-0.000000}; ! free QA*; ! COUPLE QA11 QA21 1; ! COUPLE QA12 QA22 1; ! COUPLE QA13 QA23 1; ! go; ! VAR; ! type QA*; ! FIX Q*; ! Clear[FitFunction]; ! CONVERGENCE=1E-9; !!stop; !!stop; BXI 7; BYI 14; EXI 0; EPXI 0; AXI 0; AYI 0; ins calc; !! MKSTC -> MKAC1C bmax = 15; fit MKAC1C AX 0; fit MKAC1C AY 0; FFS["fit B BX 15"]; FFS["fit B BY 15"]; FitValue["B","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; FitValue["B","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; FFS["fit QA13 BX 15"]; FitValue["QA13","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; FFS["fit Q31 BX 15"]; FitValue["Q31","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; FFS["fit Q41 BY 15"]; FitValue["Q41","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; free Q22H, Q23, Q3*, Q4*; go; FIX Q*; SAVE; !! type Q*; !! disp; !! WriteGraphs[]; !!stop; !!stop; !! MKAC1C -> MKID bmax2 = 20; fit MKID AX 0; fit MKID AY 0; FFS["fit MKID BX 15"]; FFS["fit MKID BY 15"]; FitValue["MKID","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax2,bmax2,Null]]; FitValue["MKID","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax2,bmax2,Null]]; FFS["fit Q51 BY 15"]; FitValue["Q51","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; FFS["fit Q51 BX 15"]; FitValue["Q51","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; ! FFS["fit Q62 BY 15"]; ! FitValue["Q62","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; ! FFS["fit Q63 BX 15"]; ! FitValue["Q63","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax2,bmax2,Null]]; FFS["fit Q61 BX 15"]; FitValue["Q61","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax2,bmax2,Null]]; free Q5*, Q6*; go; FIX Q*; !SAVE; ! WriteGraphs[]; ! type Q*; ! disp; !stop; !stop; !! MKID -> MKAC2C bmax2 = 20; fit MKAC2C AX 0; fit MKAC2C AY 0; ! FFS["fit Q72 BX 15"]; ! FitValue["Q72","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; ! FFS["fit QA22 BY 15"]; ! FitValue["QA22","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax2,bmax2,Null]]; ! FFS["fit QA23 BX 15"]; ! FitValue["QA23","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; free Q7*; go; FIX Q*; WriteGraphs[]; ! type Q*; ! disp; !stop; !stop; fit AX 0; fit AY 0; fit BX 7; fit BY 14; ! FFS["fit Q84 BY 15"]; ! FitValue["Q84","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; ! FFS["fit Q82 BY 15"]; ! FitValue["Q82","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; FFS["fit Q11 BX 15"]; FitValue["Q11","BX",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; ! FFS["fit Q81 BY 15"]; ! FitValue["Q81","BY",dp_,vgoal_,vnow_]:=If[vnow<0,0,If[vnow>bmax,bmax,Null]]; free Q8*, Q1*, Q21; go; ! WriteGraphs[]; type Q*; !stop; !stop; WriteGraphs[]; ! disp a; disp; ! Print[LINE["NAME","B*"]]; ! Print[LINE["NAME","MK*"]]; ! Print[LINE["GEO","^^^"],LINE["GEO","$$$"]]; SAVE; stop; stop;