[thread display] [new arrival display] [word search] [past log] [管理用]
| Subject | ： input file |
| Date | ： 2006/01/30(Mon) 20:16:24 |
| Contributor | ： Kentaro Harada < > |
!********************************************************************
!*
!* NSLS Xray-RING
!*
!********************************************************************
MOMENTUM=2.8GeV
;
MARK MKN =()
MKL =()
;
MONI PM =()
;
CAVI CAVI=(L=0.8 VOLT= 1.0MV HARM=30)
;
BEND B =(L = 2.7 ANGLE = .392699082 E1 = 0.5 E2=0.5)
;
QUAD
Q1 =(L = .45 K1 = -0.45*1.4650455999)
Q2 =(L = .8 K1 = 0.8*1.3717947451)
Q3 =(L = .45 K1 = -0.45*1.4109287101)
Q4 =(L = .225 K1 = 0.225*1.381258)
;
SEXT SF =(L =0 K2 = 0 )
SD =(L =0 K2 = 0 )
;
SEXT SD =(L =.2 K2 =-8.8021195157993 )
SF =(L =.2 K2 =5.5762274523103 )
;
DRIFT D1 =(L = 2.25 )
D2 =(L = 0.685 )
D3 =(L = 0.33 )
D4 =(L = 0.69)
D5 =(L = 0.8 )
D7 =(L = 0.7 )
D6 =(L = 0.15 )
D1R =(L = 1.85 )
;
LINE
NCELLH =(MKN D1 Q1 D2 Q2 D3 Q3 D4
B
D5 SD D7 SF D6 Q4
MKL)
NCELLR =(MKN D1R Q1 D2 Q2 D3 Q3 D4
B
D5 SD D7 SF D6 Q4
MKL)
NCELL =(NCELLH -NCELLH)
RING =(6*NCELL NCELLH -NCELLR CAVI NCELLR -NCELLH)
;
!
! dynamic aperture survey command 1
!
ON LOG;
ON RAD FLUC COD RFSW RADCOD;
FFS;
$FORM = "15.10";
PageWidth = 1999;
USE RING;
cell calc emit;
disp $$$;
w = EMITX;
sigx = Sqrt[Twiss["BX",1] * w];
sigy = Sqrt[Twiss["BY",1] * w];
Print[sigx * 1000, sigy * 1000];
SAVE ALL;
DAPWIDTH = 51;
FFS["cell calc emit"];
dp = 0.01 / SIGE;
cell calc;
zrange = Range[-1.5 * dp, 1.6 * dp, dp /2];
! Print[zrange];
range = {{0,100},{1,1},zrange};
resx = DynamicApertureSurvey[range,1000,Output->6];
Do[
Print[zrange[i],resx[2,2,i,2],(resx[2,2,i,2] - 1) * sigx * 1000 * 2];
,{i,1,Length[zrange]}];
range = {{1,1},{0,100},zrange};
resy = DynamicApertureSurvey[range,1000,Output->6];
Do[
Print[zrange[i],resy[2,2,i,2],(resy[2,2,i,2] - 1) * sigy * 1000 * 2];
,{i,1,Length[zrange]}];
stop;
stop;