Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

test-ptc-twiss-old7.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: 16.6 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: 18.02.13 18.32.26 +
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=my_sect_table,! added 16:10
176
177	!column=name,pos,k1,r1,r36,t1,t216; ! added 16:10
178
179	!sectormap,sectortable=my_sect_table,sectorfile=my_sect_file; ! added 16:10
180
181
182
183	select,flag=twiss,column=name,s,betx,bety,dx,dy; ! dx, dy are dispersions
184
185	TWISS,centre,file=twiss_madx_ring; ! jluc : centre
186
187	enter Twiss module
188
189	iteration: 1 error: 6.606181E-04 deltap: 0.000000E+00
190	orbit: 3.261640E-04 2.597068E-06 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
191
192	iteration: 2 error: 1.896028E-07 deltap: 0.000000E+00
193	orbit: 3.260716E-04 2.596347E-06 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
194
195	++++++ table: summ
196
197	length orbit5 alfa gammatr
198	534.6 -0 0.0004430990969 47.50612067
199
200	q1 dq1 betxmax dxmax
201	1.254675989 1.051070031 181.4719986 2.200422953
202
203	dxrms xcomax xcorms q2
204	1.658847454 0.000332287113 0.000213664876 1.254674956
205
206	dq2 betymax dymax dyrms
207	1.046764172 181.4881854 0 0
208
209	ycomax ycorms deltap synch_1
210	0 0 0 0
211
212	synch_2 synch_3 synch_4 synch_5
213	0 0 0 0
214
215
216	SHOW,TWSSip; ! SHOW parameters (see "*.out" file)
217
218	command: beta0 module: control
219	parameter: betx double value: 177.7010914
220	parameter: alfx double value: -2.41257182
221	parameter: mux double value: 1.254675989
222	parameter: bety double value: 32.0344759
223	parameter: alfy double value: 0.4796139615
224	parameter: muy double value: 1.254674956
225	parameter: x double value: 0.0003260716242
226	parameter: px double value: 2.596346674e-06
227	parameter: y double value: 0
228	parameter: py double value: 0
229	parameter: t double value: -1.662096185e-05
230	parameter: pt double value: 0
231	parameter: dx double value: 2.174564191
232	parameter: dpx double value: 0.02962824424
233	parameter: dy double value: 0
234	parameter: dpy double value: 0
235	parameter: wx double value: 0
236	parameter: phix double value: 0
237	parameter: dmux double value: 0
238	parameter: wy double value: 0
239	parameter: phiy double value: 0
240	parameter: dmuy double value: 0
241	parameter: ddx double value: 0
242	parameter: ddpx double value: 0
243	parameter: ddy double value: 0
244	parameter: ddpy double value: 0
245	parameter: r11 double value: -0
246	parameter: r12 double value: 0
247	parameter: r21 double value: -0
248	parameter: r22 double value: 0
249	parameter: energy double value: 100
250
251
252	! added 20 nov 2008:
253
254	plot vaxis1=betx,vaxis2=dx,haxis=s;
255
256	Plot - default table plotted: twiss
257
258	GXPLOT-X11 1.50 initialized
259
260	plot number = 1
261
262
263	! Initialize PTC
264
265	ptc_create_universe;
266
267	Now PTC
268	ptc_create_layout,model=2,method=6,nst=10,exact;
269
270	MAD-X Beam Parameters
271	Energy : 0.100000E+03
272	Kinetic Energy : 0.990617E+02
273	Particle Rest Mass : 0.938272E+00
274	Momentum : 0.999956E+02
275	Setting MADx with
276	energy 100.00000000000000
277	method 6
278	Num. of steps 10
279	charge 1.0000000000000000
280	Length of machine: 534.59999999999991
281	The machine is a RING
282	------------------------------------ PTC Survey ------------------------------------
283	Before start: 0.0000000000000000 0.0000000000000000 0.0000000000000000
284	Before end: 0.0000000000000000 0.0000000000000000 0.0000000000000000
285	After start: 0.0000000000000000 0.0000000000000000 0.0000000000000000
286	After end: -40.580758467025781 0.0000000000000000 532.53346053079656
287	!ptc_create_layout,model=2,method=6,nst=1,exact; ! steps pb more apparent
288
289	!ptc_create_layout,model=2,method=6,nst=4,exact; ! looks like 2 is minimum to split the element at the middle...
290
291	ptc_align;
292
293
294
295	select,flag=ptc_twiss,column=name,s,beta11,beta21,
296
297	beta12,beta22,disp1,disp3,x,px,y,py;
298
299
300
301	ptc_twiss,closed_orbit,icase=5,file=twiss_ptc_ring,summary_file=twiss_ptc_summary,no=4;
302
303	********** State Summary **************
304	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
305	Rectangular Bend: input arc length (rho alpha)
306	Default integration method 6
307	Default integration steps 10
308	This is a proton
309	EXACT_MODEL = TRUE
310	TOTALPATH = 0
311	RADIATION = FALSE
312	STOCHASTIC = FALSE
313	ENVELOPE = FALSE
314	NOCAVITY = TRUE
315	TIME = TRUE
316	FRINGE = FALSE
317	PARA_IN = FALSE
318	ONLY_4D = TRUE
319	DELTA = TRUE
320	SPIN = FALSE
321	MODULATION = FALSE
322	RAMPING = FALSE
323	ACCELERATE = FALSE
324	Closed orbit: 3.2607158925373049E-004 2.5963462592936979E-006 0.0000000000000000 0.0000000000000000 0.0000000000000000
325	Deviation from symplecticity 0.5251354906476990E-11 %
326	PARENT LAYOUT NAME :No name assigned
327	NUMBER OF ORIGINAL LAYOUT ELEMENTS : 146
328	NUMBER OF THIN OBJECTS : 1180
329	TOTAL IDEAL LENGTH OF STRUCTURE : 534.60000000000196
330	TOTAL INTEGRATION LENGTH OF STRUCTURE (mad8 style survey) : 534.60000000000196
331	Deviation from symplecticity 0.5251354906476990E-11 %
332
333	++++++ table: ptc_twiss_summary
334
335	length alpha_c alpha_c_p alpha_c_p2
336	534.6 0.0004430990863 -0.0001937066866 -1000000
337
338	alpha_c_p3 eta_c gamma_tr q1
339	-1000000 0.0003550636492 47.50612124 0.2546759882
340
341	q2 dq1 dq2 qs
342	0.2546749559 1.05106564 1.046759439 0
343
344	beta_x_min beta_x_max beta_y_min beta_y_max
345	32.03433146 177.7089538 32.03308004 177.7231017
346
347	deltap orbit_x orbit_px orbit_y
348	0 0.0003260715893 2.596346259e-06 0
349
350	orbit_py orbit_pt orbit_-cT xcorms
351	0 0 0 0.0002136648593
352
353	ycorms pxcorms pycorms xcomax
354	0 3.396649456e-06 0 0.0003315825523
355
356	ycomax pxcomax pycomax
357	0 5.310270774e-06 0
358	ptc_twiss,closed_orbit,icase=5,slice_magnets,file=twiss_ptc_ring_slices,no=4;
359
360	********** State Summary **************
361	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
362	Rectangular Bend: input arc length (rho alpha)
363	Default integration method 6
364	Default integration steps 10
365	This is a proton
366	EXACT_MODEL = TRUE
367	TOTALPATH = 0
368	RADIATION = FALSE
369	STOCHASTIC = FALSE
370	ENVELOPE = FALSE
371	NOCAVITY = TRUE
372	TIME = TRUE
373	FRINGE = FALSE
374	PARA_IN = FALSE
375	ONLY_4D = TRUE
376	DELTA = TRUE
377	SPIN = FALSE
378	MODULATION = FALSE
379	RAMPING = FALSE
380	ACCELERATE = FALSE
381	Closed orbit: 3.2607158925373049E-004 2.5963462592936979E-006 0.0000000000000000 0.0000000000000000 0.0000000000000000
382	Deviation from symplecticity 0.5251354906476990E-11 %
383	PARENT LAYOUT NAME :No name assigned
384	NUMBER OF ORIGINAL LAYOUT ELEMENTS : 146
385	NUMBER OF THIN OBJECTS : 1180
386	TOTAL IDEAL LENGTH OF STRUCTURE : 534.60000000000196
387	TOTAL INTEGRATION LENGTH OF STRUCTURE (mad8 style survey) : 534.60000000000196
388	PARENT LAYOUT NAME :No name assigned
389	NUMBER OF ORIGINAL LAYOUT ELEMENTS : 146
390	NUMBER OF THIN OBJECTS : 1180
391	TOTAL IDEAL LENGTH OF STRUCTURE : 534.60000000000196
392	TOTAL INTEGRATION LENGTH OF STRUCTURE (mad8 style survey) : 534.60000000000196
393	Deviation from symplecticity 0.5251354906476990E-11 %
394
395	++++++ table: ptc_twiss_summary
396
397	length alpha_c alpha_c_p alpha_c_p2
398	534.6 0.0004430990863 -0.0001937066866 -1000000
399
400	alpha_c_p3 eta_c gamma_tr q1
401	-1000000 0.0003550636492 47.50612124 0.2546759882
402
403	q2 dq1 dq2 qs
404	0.2546749559 1.05106564 1.046759439 0
405
406	beta_x_min beta_x_max beta_y_min beta_y_max
407	31.29632882 181.4719986 31.2954739 181.4881854
408
409	deltap orbit_x orbit_px orbit_y
410	0 0.0003260715893 2.596346259e-06 0
411
412	orbit_py orbit_pt orbit_-cT xcorms
413	0 0 0 0.0002136648593
414
415	ycorms pxcorms pycorms xcomax
416	0 3.396649456e-06 0 0.0003315825523
417
418	ycomax pxcomax pycomax
419	0 5.310270774e-06 0
420
421
422	PTC_NORMAL,closed_orbit,maptable,icase=5,no=4;
423
424	********** State Summary **************
425	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
426	Rectangular Bend: input arc length (rho alpha)
427	Default integration method 6
428	Default integration steps 10
429	This is a proton
430	EXACT_MODEL = TRUE
431	TOTALPATH = 0
432	RADIATION = FALSE
433	STOCHASTIC = FALSE
434	ENVELOPE = FALSE
435	NOCAVITY = TRUE
436	TIME = TRUE
437	FRINGE = FALSE
438	PARA_IN = FALSE
439	ONLY_4D = TRUE
440	DELTA = TRUE
441	SPIN = FALSE
442	MODULATION = FALSE
443	RAMPING = FALSE
444	ACCELERATE = FALSE
445	Closed orbit: 3.2607158925373049E-004 2.5963462592936979E-006 0.0000000000000000 0.0000000000000000 0.0000000000000000
446
447
448	ptc_end;
449
450	Killing LayoutNo name assigned
451	NODE LAYOUT HAS BEEN KILLED
452	Layout killed
453
454
455	!write,table="map_table",file="map_table";
456
457	write,table="map_table",file="ring_matrix_at_end";
458
459
460
461	!=========== TRANSFER-LINE PARAMETERS ====================================
462
463
464
465	call, file = "../share/line_lattice/line.seq";
466
467	TITLE, s="e- Beam Delivery System 1 ->20 mr (ILC2005)";
468
469	OPTION,ECHO = false;!, VERIFY = true;
470
471	++++++ info: ygapsp1 redefined
472	++++++ info: ygapsp2 redefined
473	++++++ info: ygapsp3 redefined
474	++++++ info: ygapsp4 redefined
475	++++++ info: ygapsp5 redefined
476	Now PTC
477	MAD-X Beam Parameters
478	Energy : 0.250000E+03
479	Kinetic Energy : 0.249999E+03
480	Particle Rest Mass : 0.510999E-03
481	Momentum : 0.250000E+03
482	Setting MADx with
483	energy 250.00000000000000
484	method 6
485	Num. of steps 5
486	charge -1.0000000000000000
487	Length of machine: 534.59999999999991
488	The machine is a RING
489	------------------------------------ PTC Survey ------------------------------------
490	Before start: 0.0000000000000000 0.0000000000000000 0.0000000000000000
491	Before end: 0.0000000000000000 0.0000000000000000 0.0000000000000000
492	After start: 0.0000000000000000 0.0000000000000000 0.0000000000000000
493	After end: -40.580758467025781 0.0000000000000000 532.53346053079656
494	********** State Summary **************
495	MADTHICK=>KIND = 37 MATRIX-KICK-MATRIX
496	Rectangular Bend: input arc length (rho alpha)
497	Default integration method 6
498	Default integration steps 5
499	This is an electron (positron actually if charge=1)
500	EXACT_MODEL = TRUE
501	TOTALPATH = 0
502	RADIATION = FALSE
503	STOCHASTIC = FALSE
504	ENVELOPE = FALSE
505	NOCAVITY = TRUE
506	TIME = TRUE
507	FRINGE = FALSE
508	PARA_IN = FALSE
509	ONLY_4D = TRUE
510	DELTA = TRUE
511	SPIN = FALSE
512	MODULATION = FALSE
513	RAMPING = FALSE
514	ACCELERATE = FALSE
515	Closed orbit specified by the user!
516	Closed orbit: 8.8978067121092011E-005 7.2839183526605370E-007 5.2160503323142552E-005 -6.7480410626117705E-007 0.0000000000000000
517	Reduced SUMM Table (closed orbit not requested)
518
519	++++++ table: ptc_twiss_summary
520
521	length alpha_c alpha_c_p alpha_c_p2
522	534.6 -1000000 -1000000 -1000000
523
524	alpha_c_p3 eta_c gamma_tr q1
525	-1000000 -1000000 -1000000 1.254646403
526
527	q2 dq1 dq2 qs
528	1.254592531 -1000000 -1000000 -1000000
529
530	beta_x_min beta_x_max beta_y_min beta_y_max
531	-1000000 -1000000 -1000000 -1000000
532
533	deltap orbit_x orbit_px orbit_y
534	0 -1000000 -1000000 -1000000
535
536	orbit_py orbit_pt orbit_-cT xcorms
537	-1000000 -1000000 -1000000 -1000000
538
539	ycorms pxcorms pycorms xcomax
540	-1000000 -1000000 -1000000 -1000000
541
542	ycomax pxcomax pycomax
543	-1000000 -1000000 -1000000
544	Killing LayoutNo name assigned
545	Layout killed
546
547	Number of warnings: 0
548
549	++++++++++++++++++++++++++++++++++++++++++++
550	+ MAD-X 5.00.19 (64 bit) finished normally +
551	++++++++++++++++++++++++++++++++++++++++++++

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

Download in other formats:

Original Format