[thread display] [new arrival display] [word search] [past log] [管理用]
| Subject | ： Re: その５ |
| Date | ： 2006/01/30(Mon) 20:14:06 |
| Contributor | ： Kentaro Harada < > |
!********************************************************************
!*
!* NSLS XRING
!*
!********************************************************************
MOMENTUM=2.5GeV
;
MARK MKN =()
MKL =()
;
MONI PM =()
;
CAVI CAVI=(VOLT= 1.12MV HARM = 30 )
;
BEND B =(L = 2.7 K1 = -0.0724921758 ANGLE = .39269908 E1 = 0.5 E2=0.5)
;
QUAD
Q1 =(L = .45 K1 = -0.675839592 )
Q2 =(L = .8 K1 = 1.069849888 )
Q3 =(L = .4132 K1 = -0.57926284872 )
Q4 =(L = .225 K1 = 0.2923738695 )
;
SEXT SF =(L =0.2 K2 = 0 )
SD =(L =0.2 K2 = 0 )
;
DRIFT D1 =(L = 2.25 )
D2 =(L = 0.685 )
D3 =(L = 0.3484 )
D4 =(L = 0.70825)
D5 =(L = 0.8 )
D7 =(L = 0.7 )
D6 =(L = 0.15 )
;
LINE
NCELLH =(MKN D1 Q1 D2 Q2 D3 Q3 D4
B
D5 SD D7 SF D6 Q4
MKL)
NCELL =(NCELLH -NCELLH)
RING =(4*NCELL CAVI 4*NCELL )
;
!
! chromaticity correction
!
ON LOG;
ON RAD FLUC COD RFSW RADCOD;
FFS;
$FORM = "15.10";
PageWidth = 1999;
USE RING;
cell calc;
{nx0, ny0}=Twiss[{"NX","NY"},"$$$"]/2/Pi;
fit nx nx0 2;
fit ny ny0 2;
DP=0.0001;
free S*;
go;
FITP 1;
type S*;
SAVE S*;
tunex = Table[{0,0},{41}];
tuney = Table[{0,0},{41}];
Do[
DP0 = (i - 21) * 0.001;
FFS["cell calc"];
tunex[i] = {DP0,Twiss["NX","$$$"]/2/Pi};
tuney[i] = {DP0,Twiss["NY","$$$"]/2/Pi};
! Print[i,DP0,tunex[i,2],tuney[i,2]];
,{i,1,41,1}];
Print[PolynomialFit[tunex,3]];
Print[PolynomialFit[tuney,3]];
stop;
stop;