1 | !TITLE,'Test input for MAD-X'; |
---|
2 | |
---|
3 | eg := 100; |
---|
4 | bg := eg/pmass; |
---|
5 | en := 3.75e-06; |
---|
6 | epsx := en/bg; |
---|
7 | epsy := en/bg; |
---|
8 | beam, particle = proton, energy = eg , |
---|
9 | sigt= 0.077 , |
---|
10 | sige= 1.1e-4 , |
---|
11 | npart= 1.05e11 , |
---|
12 | exn=4*en, eyn=4*en , |
---|
13 | kbunch = 10, |
---|
14 | et = 0.002, bv = -1, |
---|
15 | ex=epsx, ey=epsy; |
---|
16 | value,epsx,epsy; |
---|
17 | option,-echo,-warn; |
---|
18 | call file=fv9.opt; // string without blanks does not need ".." |
---|
19 | call file="fv9.seq"; |
---|
20 | option,echo,warn; |
---|
21 | |
---|
22 | use,period=fivecell; |
---|
23 | SAVEBETA, LABEL=BETA1,sequence=fivecell, PLACE =partstart; |
---|
24 | TWISS; |
---|
25 | |
---|
26 | BETAXX1=BETA1->BETX; |
---|
27 | ALFAXX1=BETA1->ALFX; |
---|
28 | DXX1=BETA1->DX; |
---|
29 | DXXP1=BETA1->DPX; |
---|
30 | BETAYY1=BETA1->BETY; |
---|
31 | ALFAYY1=BETA1->ALFY; |
---|
32 | MUXX=BETA1->MUX; |
---|
33 | MUYY=BETA1->MUY; |
---|
34 | |
---|
35 | MUXX1:=MUXX/360.;! changed on request of Oliver Bruening |
---|
36 | MUYY1:=MUYY/360.;! changed on request of Oliver Bruening |
---|
37 | !MUXX1:=MUXX*TWOPI/360.; |
---|
38 | !MUYY1:=MUYY*TWOPI/360.; |
---|
39 | value, muxx,muyy,muxx1,muyy1; |
---|
40 | |
---|
41 | R111:=cos(MUXX1)+ALFAXX1*sin(MUXX1); |
---|
42 | R112:=BETAXX1*sin(MUXX1); |
---|
43 | R116:=(1.-R111)*DXX1-R112*DXXP1; |
---|
44 | R121:=-sin(MUXX1)/BETAXX1*(1+ALFAXX1^2); |
---|
45 | R122:=cos(MUXX1)-ALFAXX1*sin(MUXX1); |
---|
46 | R126:=(1.-R122)*DXXP1-R121*DXX1; |
---|
47 | R151:=R121*R116-R111*R126; |
---|
48 | R152:=R122*R116-R112*R126; |
---|
49 | R133:=cos(MUYY1)+ALFAYY1*sin(MUYY1); |
---|
50 | R134:=BETAYY1*sin(MUYY1); |
---|
51 | R143:=-sin(MUYY1)/BETAYY1*(1+ALFAYY1^2); |
---|
52 | R144:=cos(MUYY1)-ALFAYY1*sin(MUYY1); |
---|
53 | value,DXX1,DXXP1,R111,R112,R116,R121,R122,R126; |
---|
54 | |
---|
55 | MATSUP1 : MATRIX,RM11=R111,RM12=R112,RM21=R121,RM22=R122, |
---|
56 | RM16=R116,RM26=R126,RM51=R151,RM52=R152, |
---|
57 | RM33=R133,RM34=R134,RM43=R143,RM44=R144, |
---|
58 | RM55=1.,RM66=1.; |
---|
59 | |
---|
60 | SEQEDIT,SEQUENCE=fivecell; |
---|
61 | INSTALL, element=MATSUP1, at=0.0,from=partstart; |
---|
62 | ENDEDIT; |
---|
63 | |
---|
64 | use,period=fivecell; |
---|
65 | select,flag=twiss,class=mq,pattern="d.*", |
---|
66 | column=name,parent,keyword,s,betx,bety,dx,dy,k1l; |
---|
67 | twiss,file; |
---|
68 | stop; |
---|