[thread display] [new arrival display] [word search] [past log] [管理用]
| Subject | ： input file |
| Date | ： 2006/01/30(Mon) 20:17:23 |
| 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)
;
!
! trackparticles sample
!
ON LOG;
OFF ECHO;
ON RAD FLUC COD RFSW RADCOD;
FFS;
$FORM = "20.10";
PageWidth = 1999;
USE NCELL;
cell calc emit;
disp;
SAVE ALL;
rqk = Sqrt[LINE["K1","Q2"]/LINE["L", "Q2"]];
ql = LINE["L", "Q2"];
Print[rqk, ql];
qkmath = {{0,0},{0,0}};
qkmatv = {{0,0},{0,0}};
qkmath[1] = {Cos[rqk * ql],Sin[rqk * ql]/rqk};
qkmath[2] = {-rqk * Sin[rqk * ql], Cos[rqk * ql]};
qkmatv[1] = {Cosh[rqk * ql], Sinh[rqk * ql]/rqk};
qkmatv[2] = {rqk * Sinh[rqk * ql], Cosh[rqk * ql]};
Print[qkmath[1,1],qkmath[1,2],0,0];
Print[qkmath[2,1],qkmath[2,2],0,0];
Print[0,0,qkmatv[1,1],qkmatv[1,2]];
Print[0,0,qkmatv[2,1],qkmatv[2,2]];
Print[Det[qkmath], Det[qkmatv]];
before = {5,{{0.001,0,0,0},{0,0.001,0,0},{0,0,0.001,0},{0,0,0,0.001},{0,0,0,0},{0,0,0,0},{1,1,1,1}}};
after = TrackParticles[before,6];
after[2] = after[2] * 1000;
Print[after[2,1,1],after[2,1,2],after[2,1,3],after[2,1,4]];
Print[after[2,2,1],after[2,2,2],after[2,2,3],after[2,2,4]];
Print[after[2,3,1],after[2,3,2],after[2,3,3],after[2,3,4]];
Print[after[2,4,1],after[2,4,2],after[2,4,3],after[2,4,4]];
stop;
stop;