Context Navigation

test-ptc-twiss-old3.ref

Last change on this file was 430, checked in by touze, 11 years ago
import madx-5.01.00
File size: 13.1 KB

Line
1
2	+++++++++++++++++++++++++++++++++++++++++++
3	+ MAD-X 5.00.19 (64 bit, Darwin) +
4	+ DEBUG Version - use at own risk! +
5	+ Code Modification Date: 2012.12.11 +
6	+ Execution Time Stamp: 15.02.13 11.54.18 +
7	+++++++++++++++++++++++++++++++++++++++++++
8	!TITLE,'Test input for PTC_TWISS';
9
10
11
12	!=========== RING PARAMETERS ===========================================
13
14	call, file = "../share/ring_lattice/ring.seq";
15
16	eg := 100;
17
18	bg := eg/pmass;
19
20	en := 3.75e-06;
21
22	epsx := en/bg;
23
24	epsy := en/bg;
25
26	beam, particle = proton, energy = eg ,
27
28	sigt= 0.077 ,
29
30	sige= 1.1e-4 ,
31
32	npart= 1.05e11 ,
33
34	exn=4en, eyn=4en ,
35
36	kbunch = 10,
37
38	et = 0.002, bv = 1,
39
40	ex=epsx, ey=epsy;
41
42	value,epsx,epsy;
43
44	epsx = 3.518520049e-08 ;
45	epsy = 3.518520049e-08 ;
46	option,-echo;
47
48
49
50	use,period=fivecell;
51
52	// alignment errors
53
54
55
56	select,flag=error,clear;
57
58	select,flag=error,pattern="q.*",range=qf.1;
59
60	ealign,dx=0.0001;
61
62	Assigned alignment errors to 1 elements
63	!select,flag=error,class=quadrupole,range=qf.4/qf.5;
64
65	!ealign,dy=0.0002,dtheta=0.0003;
66
67
68
69	/*
70
71	select,flag=error,clear;
72
73
74
75	// field errors
76
77	gcutr=3.0;
78
79	b1r=1.e-4;
80
81	b2r=2.e-4;
82
83	b3r=3.e-4;
84
85	b4r=4.e-4;
86
87	b5r=5.e-4;
88
89	a1r=1.e-4;
90
91	a2r=2.e-4;
92
93	a3r=3.e-4;
94
95	a4r=4.e-4;
96
97	a5r=5.e-4;
98
99
100
101	Select, flag=error, clear = true;
102
103	select, flag=error, pattern="q.\..";
104
105	efcomp, order:=1, radius:=0.010,
106
107	dknr={0,1e-1*b2r},
108
109	dksr={0,1e-1*a2r};
110
111	esave;
112
113	*/
114
115
116
117	acbh1=1e-6;
118
119	acbh2=1e-6;
120
121	acbh3=1e-6;
122
123	acbh4=1e-6;
124
125	acbh5=1e-6;
126
127
128
129	/*
130
131	acbv1=1e-6;
132
133	acbv2=1e-6;
134
135	acbv3=1e-6;
136
137	acbv4=1e-6;
138
139	acbv5=1e-6;
140
141	*/
142
143
144
145	!select,flag=sectormap,class=drift;
146
147
148
149	!select,flag=my_table,class=quadrupole,
150
151	!column=name,parent,keyword,s,betx,bety,dx,dy,k1l,
152
153	!re11,re12,re13,re14;
154
155
156
157	!twiss,centre,rmatrix,
158
159	!sectormap,sectorfile=my_sectorfile,
160
161	!table=my_table,file=twiss_fv9;
162
163
164
165
166
167	! At the end of the ring the data_block is produced
168
169
170
171	SAVEBETA, label=TWSSip, place=#E,sequence=fivecell;
172
173
174
175	select,flag=twiss,column=name,s,betx,bety,dx,dy;
176
177
178
179	TWISS,file=twiss_madx_ring;
180
181	enter Twiss module
182
183	iteration: 1 error: 6.606181E-04 deltap: 0.000000E+00
184	orbit: 3.261640E-04 2.597068E-06 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
185
186	iteration: 2 error: 1.896028E-07 deltap: 0.000000E+00
187	orbit: 3.260716E-04 2.596347E-06 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
188
189	++++++ table: summ
190
191	length orbit5 alfa gammatr
192	534.6 -0 0.0004430990969 47.50612067
193
194	q1 dq1 betxmax dxmax
195	1.254675989 1.051070031 177.7089538 2.177407271
196
197	dxrms xcomax xcorms q2
198	1.658847454 0.0003315825458 0.000213664876 1.254674956
199
200	dq2 betymax dymax dyrms
201	1.046764172 177.7231016 0 0
202
203	ycomax ycorms deltap synch_1
204	0 0 0 0
205
206	synch_2 synch_3 synch_4 synch_5
207	0 0 0 0
208
209
210	SHOW,TWSSip; ! SHOW parameters (see "*.out" file)
211
212	command: beta0 module: control
213	parameter: betx double value: 177.7010914
214	parameter: alfx double value: -2.41257182
215	parameter: mux double value: 1.254675989
216	parameter: bety double value: 32.0344759
217	parameter: alfy double value: 0.4796139615
218	parameter: muy double value: 1.254674956
219	parameter: x double value: 0.0003260716242
220	parameter: px double value: 2.596346674e-06
221	parameter: y double value: 0
222	parameter: py double value: 0
223	parameter: t double value: -1.662096185e-05
224	parameter: pt double value: 0
225	parameter: dx double value: 2.174564191
226	parameter: dpx double value: 0.02962824424
227	parameter: dy double value: 0
228	parameter: dpy double value: 0
229	parameter: wx double value: 0
230	parameter: phix double value: 0
231	parameter: dmux double value: 0
232	parameter: wy double value: 0
233	parameter: phiy double value: 0
234	parameter: dmuy double value: 0
235	parameter: ddx double value: 0
236	parameter: ddpx double value: 0
237	parameter: ddy double value: 0
238	parameter: ddpy double value: 0
239	parameter: r11 double value: -0
240	parameter: r12 double value: 0
241	parameter: r21 double value: -0
242	parameter: r22 double value: 0
243	parameter: energy double value: 100
244
245
246
247
248	! Initialize PTC
249
250	ptc_create_universe;
251
252	Now PTC
253	ptc_create_layout,model=2,method=6,nst=10,exact;
254
255	MAD-X Beam Parameters
256	Energy : 0.100000E+03
257	Kinetic Energy : 0.990617E+02
258	Particle Rest Mass : 0.938272E+00
259	Momentum : 0.999956E+02
260	Setting MADx with
261	energy 100.000000000000
262	method 6
263	Num. of steps 10
264	charge 1.00000000000000
265	Length of machine: 534.600000000000
266	The machine is a RING
267	------------------------------------ PTC Survey ------------------------------------
268	Before start: 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
269	Before end: 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
270	After start: 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
271	After end: -40.5807584670258 0.000000000000000E+000 532.533460530797
272	ptc_align;
273
274
275
276	select,flag=ptc_twiss,column=name,s,beta11,beta21,
277
278	beta12,beta22,disp1,disp3,x,px,y,py;
279
280
281
282	ptc_twiss,closed_orbit,icase=5,file=twiss_ptc_ring,no=2;
283
284	********** State Summary **************
285	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
286	Rectangular Bend: input arc length (rho alpha)
287	Default integration method 6
288	Default integration steps 10
289	This is a proton
290	EXACT_MODEL = TRUE
291	TOTALPATH = 0
292	RADIATION = FALSE
293	STOCHASTIC = FALSE
294	ENVELOPE = FALSE
295	NOCAVITY = TRUE
296	TIME = TRUE
297	FRINGE = FALSE
298	PARA_IN = FALSE
299	ONLY_4D = TRUE
300	DELTA = TRUE
301	SPIN = FALSE
302	MODULATION = FALSE
303	RAMPING = FALSE
304	ACCELERATE = FALSE
305	Closed orbit: 3.260715892537305E-004 2.596346259293698E-006 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
306	Deviation from symplecticity 0.5251354906476990E-11 %
307	PARENT LAYOUT NAME :No name assigned
308	NUMBER OF ORIGINAL LAYOUT ELEMENTS : 146
309	NUMBER OF THIN OBJECTS : 1180
310	TOTAL IDEAL LENGTH OF STRUCTURE : 534.600000000002
311	TOTAL INTEGRATION LENGTH OF STRUCTURE (mad8 style survey) : 534.600000000002
312	Deviation from symplecticity 0.5251354906476990E-11 %
313
314	++++++ table: ptc_twiss_summary
315
316	length alpha_c alpha_c_p alpha_c_p2
317	534.6 0.0004430990863 -0.0001937066866 -1000000
318
319	alpha_c_p3 eta_c gamma_tr q1
320	-1000000 0.0003550636492 47.50612124 0.2546759882
321
322	q2 dq1 dq2 qs
323	0.2546749559 1.05106564 1.046759439 0
324
325	beta_x_min beta_x_max beta_y_min beta_y_max
326	32.03433146 177.7089538 32.03308004 177.7231017
327
328	deltap orbit_x orbit_px orbit_y
329	0 0.0003260715893 2.596346259e-06 0
330
331	orbit_py orbit_pt orbit_-cT xcorms
332	0 0 0 0.0002136648593
333
334	ycorms pxcorms pycorms xcomax
335	0 3.396649456e-06 0 0.0003315825523
336
337	ycomax pxcomax pycomax
338	0 5.310270774e-06 0
339
340
341	PTC_NORMAL,closed_orbit,maptable,icase=5,no=2;
342
343	********** State Summary **************
344	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
345	Rectangular Bend: input arc length (rho alpha)
346	Default integration method 6
347	Default integration steps 10
348	This is a proton
349	EXACT_MODEL = TRUE
350	TOTALPATH = 0
351	RADIATION = FALSE
352	STOCHASTIC = FALSE
353	ENVELOPE = FALSE
354	NOCAVITY = TRUE
355	TIME = TRUE
356	FRINGE = FALSE
357	PARA_IN = FALSE
358	ONLY_4D = TRUE
359	DELTA = TRUE
360	SPIN = FALSE
361	MODULATION = FALSE
362	RAMPING = FALSE
363	ACCELERATE = FALSE
364	Closed orbit: 3.260715892537305E-004 2.596346259293698E-006 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
365
366
367	ptc_end;
368
369	Killing LayoutNo name assigned
370	NODE LAYOUT HAS BEEN KILLED
371	Layout killed
372
373
374	!write,table="map_table",file="map_table";
375
376	write,table="map_table",file="ring_matrix_at_end";
377
378
379
380	!=========== TRANSFER-LINE PARAMETERS ====================================
381
382
383
384	call, file = "../share/line_lattice/line.seq";
385
386	TITLE, s="e- Beam Delivery System 1 ->20 mr (ILC2005)";
387
388	OPTION,ECHO = false;!, VERIFY = true;
389
390	++++++ info: ygapsp1 redefined
391	++++++ info: ygapsp2 redefined
392	++++++ info: ygapsp3 redefined
393	++++++ info: ygapsp4 redefined
394	++++++ info: ygapsp5 redefined
395	Now PTC
396	MAD-X Beam Parameters
397	Energy : 0.250000E+03
398	Kinetic Energy : 0.249999E+03
399	Particle Rest Mass : 0.510999E-03
400	Momentum : 0.250000E+03
401	Setting MADx with
402	energy 250.000000000000
403	method 6
404	Num. of steps 5
405	charge -1.00000000000000
406	Length of machine: 534.600000000000
407	The machine is a RING
408	------------------------------------ PTC Survey ------------------------------------
409	Before start: 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
410	Before end: 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
411	After start: 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000
412	After end: -40.5807584670258 0.000000000000000E+000 532.533460530797
413	********** State Summary **************
414	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
415	Rectangular Bend: input arc length (rho alpha)
416	Default integration method 6
417	Default integration steps 5
418	This is an electron (positron actually if charge=1)
419	EXACT_MODEL = TRUE
420	TOTALPATH = 0
421	RADIATION = FALSE
422	STOCHASTIC = FALSE
423	ENVELOPE = FALSE
424	NOCAVITY = TRUE
425	TIME = TRUE
426	FRINGE = FALSE
427	PARA_IN = FALSE
428	ONLY_4D = TRUE
429	DELTA = TRUE
430	SPIN = FALSE
431	MODULATION = FALSE
432	RAMPING = FALSE
433	ACCELERATE = FALSE
434	Closed orbit specified by the user!
435	Deviation from symplecticity 0.5251354906476990E-11 %
436	Reduced SUMM Table (closed orbit not requested)
437
438	++++++ table: ptc_twiss_summary
439
440	length alpha_c alpha_c_p alpha_c_p2
441	534.6 -1000000 -1000000 -1000000
442
443	alpha_c_p3 eta_c gamma_tr q1
444	-1000000 -1000000 -1000000 1.254494202
445
446	q2 dq1 dq2 qs
447	1.254632429 -1000000 -1000000 -1000000
448
449	beta_x_min beta_x_max beta_y_min beta_y_max
450	-1000000 -1000000 -1000000 -1000000
451
452	deltap orbit_x orbit_px orbit_y
453	0 -1000000 -1000000 -1000000
454
455	orbit_py orbit_pt orbit_-cT xcorms
456	-1000000 -1000000 -1000000 -1000000
457
458	ycorms pxcorms pycorms xcomax
459	-1000000 -1000000 -1000000 -1000000
460
461	ycomax pxcomax pycomax
462	-1000000 -1000000 -1000000
463	Killing LayoutNo name assigned
464	Layout killed
465
466	Number of warnings: 0
467
468	++++++++++++++++++++++++++++++++++++++++++++
469	+ MAD-X 5.00.19 (64 bit) finished normally +
470	++++++++++++++++++++++++++++++++++++++++++++

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format