source: PSPA/madxPSPA/tests/test-ptc-normal/test-ptc-normal.ref @ 430

Last change on this file since 430 was 430, checked in by touze, 11 years ago

import madx-5.01.00
File size: 8.3 KB
Line 
1
2  +++++++++++++++++++++++++++++++++++++++++++
3  +    MAD-X 5.00.19  (64 bit, Linux)       +
4  +    DEBUG Version - use at own risk!     +
5  + Code Modification Date: 2012.12.11      +
6  + Execution Time Stamp: 16.01.13 13.51.45 +
7  +++++++++++++++++++++++++++++++++++++++++++
8   !                     FODO for MAD-X
9
10   ! FODO --: with SIS cell length and SIS bending magnets
11
12   !---------------------  DRIFT SPACES DEFINITION  ---------------------
13
14 
15
16   dr: drift, L=1.191;
17
18   
19
20   !---------------------------------------------------------------------
21
22
23
24   !--------------------------  BENDING MAGNET ---------------------------
25
26   bang=0.1;
27
28   ALPHA  := 15 * RADDEG;
29
30   LL     := 150 * RADDEG;
31
32   PFR    := 7.3* RADDEG;
33
34   FI     := 0.7;
35
36   RB: SBEND,L=LL,ANGLE= ALPHA;!,E1= PFR,E2=PFR,HGAP= 0.035,FINT= FI;
37
38ksf = .0198492943e1;
39
40   ! ------------------- horizontal quadrupole ----------------------------
41
42   qf:  quadrupole, L= 4.0,  K1=  0.05;
43
44   qd:  quadrupole, L= 4.0,  K1= -0.05;
45
46   hk:  hkicker, hkick=0e-3;
47
48!   hk:  hkicker, hkick=0e-3;
49
50mscbh: sextupole, l=1.1, k2:=0.0;!ksf;
51
52!
53
54mb1: multipole;
55
56!
57
58   !---------------------------------------------------------------
59
60
61
62   ! ---------------------------  LINE DEFINITION ---------------------------
63
64   cella: line=(qf,dr,rb,dr,qd,mscbh,mb1,dr,rb,dr,hk);
65
66   !------------------------------------------------------------------------
67
68
69
70!   beam,particle=proton,energy=0.942563;
71
72   beam,particle=proton,energy=5000.;
73
74!   beam,particle=proton,energy=0.942563,sequence=cella;
75
76!   beam;
77
78
79
80   use, period=cella;
81
82   show,beam%cella;
83
84   show,beam;
85
86command: beam  module: control
87parameter: particle   string: proton
88parameter: sequence   string: default_beam
89parameter: bunched   logical: true
90parameter: radiate   logical: false
91parameter: mass   double value:        0.938272013
92parameter: charge   double value:                  1
93parameter: energy   double value:               5000
94parameter: pc   double value:        4999.999912
95parameter: gamma   double value:        5328.945051
96parameter: ex   double value:                  1
97parameter: exn   double value:        21315.77983
98parameter: ey   double value:                  1
99parameter: eyn   double value:        21315.77983
100parameter: et   double value:                  1
101parameter: sigt   double value:                  0
102parameter: sige   double value:                  0
103parameter: kbunch   double value:                  1
104parameter: npart   double value:                  0
105parameter: bcurrent   double value:                  0
106parameter: freq0   double value:                  0
107parameter: circ   double value:                  0
108parameter: dtbyds   double value:                  0
109parameter: deltap   double value:                  0
110parameter: beta   double value:       0.9999999824
111parameter: alfa   double value:    3.521417482e-08
112parameter: u0   double value:                  0
113parameter: qs   double value:                  0
114parameter: arad   double value:    1.534698248e-18
115parameter: bv   double value:                  1
116parameter: pdamp   double array: 1.000000e+00 1.000000e+00 2.000000e+00
117parameter: n1min   double value:                 -1
118!---------------MATCHING----------------------------------------
119
120
121
122!   match,sequence=cella;
123
124!   VARY,NAME=QF[K1],STEP=0.01;
125
126!   VARY,NAME=QD[K1],STEP=0.01;
127
128!   CONSTRAINT,sequence=cella,RANGE=#E,MUX=0.3525,MUY=0.29;
129
130!   MIGRAD,calls=2000,tolerance=1.e-15;
131
132!   endmatch;
133
134
135
136!---------------------------------------------------------
137
138 
139
140    select,flag=makethin,clear;
141
142    select,flag=makethin,pattern="q*",slice=100;
143
144    makethin,sequence=cella;
145
146makethin: slicing sequence : cella
147
148
149    use, period=cella;
150
151 
152
153SELECT,FLAG=ERROR,pattern="mb1*";
154
155!value,bang;
156
157EFCOMP, DKS:={0,0,0,0,-1.};
158
159Assigned field errors to 1 elements
160   select,flag=sectormap,clear;
161
162   select,flag=sectormap,range=#e;
163
164   select,flag=twiss,clear;
165
166   select,flag=twiss,column=name,s,keyword,L,betx,bety,alfx,alfy,dx,dpx;
167
168
169
170
171
172
173
174!   twiss,rmatrix,file,sectormap;
175
176
177
178!   setplot,post=1,ascale=1.5, lscale=1.5,rscale=1.5, sscale=1.5, lwidth=3;
179
180!   PLOT, HAXIS=S,VAXIS1=BETX,VAXIS2=BETY,colour=100;
181
182
183
184TWISS;
185
186enter Twiss module
187++++++ info: Zero value of SIGT replaced by 1.
188++++++ info: Zero value of SIGE replaced by 1/1000.
189 
190iteration:   1 error:   0.000000E+00 deltap:   0.000000E+00
191orbit:   0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
192
193++++++ table: summ
194
195            length             orbit5               alfa            gammatr
196       19.09998776                 -0      0.07360632992        3.685890406
197
198                q1                dq1            betxmax              dxmax
199      0.3642635312      -0.5650132384        36.91316294        3.946011124
200
201             dxrms             xcomax             xcorms                 q2
202       2.728571506                  0                  0       0.2988958103
203
204               dq2            betymax              dymax              dyrms
205     -0.4328308235        34.07505989                  0                  0
206
207            ycomax             ycorms             deltap            synch_1
208                 0                  0                  0                  0
209
210           synch_2            synch_3            synch_4            synch_5
211                 0                  0                  0                  0
212!sodd, distort1,multipole_order_range=-8,8,start_stop= 0.,27000.,
213
214!      print_at_end;
215
216   ptc_create_universe;
217
218 Now PTC
219   ptc_create_layout,model=1,method=2,nst=1;!,exact;
220
221 MAD-X Beam Parameters
222                  Energy : 0.500000E+04
223          Kinetic Energy : 0.499906E+04
224      Particle Rest Mass : 0.938272E+00
225                Momentum : 0.500000E+04
226 Setting MADx with
227     energy           5000.0000000000000     
228     method                   2
229     Num. of steps            1
230     charge           1.0000000000000000     
231  Length of machine:    19.099987755982863     
232 The machine is a RING
233 ------------------------------------ PTC Survey ------------------------------------
234 Before start:    0.0000000000000000        0.0000000000000000        0.0000000000000000     
235 Before   end:    0.0000000000000000        0.0000000000000000        0.0000000000000000     
236 After  start:    0.0000000000000000        0.0000000000000000        0.0000000000000000     
237 After    end:    0.0000000000000000        0.0000000000000000        19.099987755983108     
238!       select_ptc_normal, dx=1, dpx=1;
239
240!       select_ptc_normal, dy=1, dpy=1;
241
242!       select_ptc_normal, q1=0, q2=0;
243
244!       select_ptc_normal, dq1=1, dq2=1;
245
246!       select_ptc_normal, dq1=2, dq2=2;
247
248!       select_ptc_normal, anhx=1,0,0;
249
250!       select_ptc_normal, anhy=1,0,0;
251
252!       select_ptc_normal, anhx=1,0,1;
253
254!       select_ptc_normal, anhy=1,0,1;
255
256!       select_ptc_normal, anhx=0,0,2;
257
258!       select_ptc_normal, anhy=0,0,2;
259
260!       select_ptc_normal, anhx=2,0,0;
261
262!       select_ptc_normal, anhy=2,0,0;
263
264!       select_ptc_normal, anhx=1,1,0;
265
266!       select_ptc_normal, anhy=1,1,0; 
267
268!       select_ptc_normal, haml=1,0,0;
269
270!       select_ptc_normal, haml=3,0,0;
271
272!       select_ptc_normal, haml=1,2,0;
273
274        select_ptc_normal, gnfu=-5,0,0;
275
276                       
277
278   ptc_normal,closed_orbit,normal,maptable,icase=5,no=7;
279
280 ************ State Summary ****************
281 MADTHICK=>KIND =   32  DRIFT-KICK-DRIFT       
282  Rectangular Bend: input arc length (rho alpha)
283  Default integration method     2
284  Default integration steps      1
285 This is a proton
286       EXACT_MODEL =  FALSE
287       TOTALPATH   =     0
288       RADIATION   =  FALSE
289       STOCHASTIC  =  FALSE
290       ENVELOPE    =  FALSE
291       NOCAVITY    =  TRUE
292       TIME        =  TRUE
293       FRINGE      =  FALSE
294       PARA_IN     =  FALSE
295       ONLY_4D     =  TRUE
296       DELTA       =  TRUE
297       SPIN        =  FALSE
298       MODULATION  =  FALSE
299       RAMPING     =  FALSE
300       ACCELERATE  =  FALSE
301 Closed orbit:    0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000     
302 Normal Form Type 1
303Deviation from symplecticity  0.6217248937900877E-12 %
304   ptc_end;
305
306 Killing LayoutNo name assigned                                                                                                       
307 Layout killed
308   write, table=normal_results,file=ptc_normal_results.tfs;
309
310   write, table=map_table, file=ptc_map_table.tfs;
311
312    stop;
313
314
315  Number of warnings: 0
316
317  ++++++++++++++++++++++++++++++++++++++++++++
318  + MAD-X 5.00.19 (64 bit) finished normally +
319  ++++++++++++++++++++++++++++++++++++++++++++
Note: See TracBrowser for help on using the repository browser.