source: PSPA/madxPSPA/tests/test-aperture/test-aperture.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: 3.3 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: 04.02.13 12.08.26 +
7  +++++++++++++++++++++++++++++++++++++++++++
8title,'LHC Cell APERTURE exemple, BJ Aug2009';
9
10
11
12REAL CONST ANGLE.MB = TWOPI/8/(23*6+2*2*4);
13
14
15
16REAL CONST l.MB = 14.30;
17
18REAL CONST l.MQ = 3.10;
19
20
21
22REAL CONST drift1 = 3.67;
23
24REAL CONST drift2 = 1.02;
25
26REAL CONST drift3 = 1.74;
27
28
29
30len_half_cell = l.MQ+drift1+3*l.MB+2*drift2+drift3;
31
32value, len_half_cell;
33
34len_half_cell =              53.45 ;
35
36
37!!!! apertures
38
39a=0.022; b=0.018; h=0.022; v=0.022;
40
41!!!! tolerances
42
43th=0.001; tv=0.0005; tr=0.001;
44
45
46
47MB: SBEND, L := l.MB, angle = ANGLE.MB;
48
49MQ: QUADRUPOLE, L := l.MQ, apertype=rectellipse,aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
50
51
52
53kqf =  0.008703550815;
54
55kqd = -0.008705002933;
56
57
58
59LHCCELL: SEQUENCE, refer = ENTRY, L = 2*len_half_cell;
60
61 
62
63  MQ1: MQ, k1:= kqf, at= 0.;
64
65
66
67  MB1: MB,  at= l.MQ+drift1 ,
68
69            apertype=rectellipse, aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
70
71  MB2: MB,  at= l.MQ+drift1+l.mb+drift2,
72
73            apertype=rectellipse, aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
74
75  MB3: MB,  at= l.MQ+drift1+2*(l.mb+drift2),
76
77            apertype=rectellipse, aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
78
79
80
81  MQ2: MQ, k1 := kqd,  at=  len_half_cell;
82
83
84
85  MB4: MB,  at= len_half_cell + l.MQ+drift1 ,
86
87            apertype=rectellipse, aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
88
89  MB5: MB,  at= len_half_cell + l.MQ+drift1+l.mb+drift2,
90
91            apertype=rectellipse, aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
92
93  MB6: MB,  at= len_half_cell + l.MQ+drift1+2*(l.mb+drift2),
94
95            apertype=rectellipse, aperture = {a,b,h,v}, aper_tol = {th,tv,tr};
96
97
98
99endsequence;
100
101
102
103beam,particle=proton,energy=450.0;
104
105
106
107set, format="9.5f";
108
109use,period=LHCCELL;
110
111select, flag=twiss, clear;
112
113select, flag=twiss, column=name,s,l,betx,bety,alfx,alfy,dx,dpx,mux,muy,n1;
114
115twiss, sequence=LHCCELL;
116
117enter Twiss module
118++++++ info: Zero value of SIGT replaced by 1.
119++++++ info: Zero value of SIGE replaced by 1/1000.
120 
121iteration:   1 error:   0.000000E+00 deltap:   0.000000E+00
122orbit:   0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
123
124++++++ table: summ
125
126   length    orbit5      alfa   gammatr
127106.90000  -0.00000   0.00045  47.39218
128
129       q1       dq1   betxmax     dxmax
130  0.25000  -0.31791 177.45326   2.18723
131
132    dxrms    xcomax    xcorms        q2
133  1.69855   0.00000   0.00000   0.25000
134
135      dq2   betymax     dymax     dyrms
136 -0.31804 177.46816   0.00000   0.00000
137
138   ycomax    ycorms    deltap   synch_1
139  0.00000   0.00000   0.00000   0.00000
140
141  synch_2   synch_3   synch_4   synch_5
142  0.00000   0.00000   0.00000   0.00000
143
144
145!!! uses default beam and halo data (lhc at injection)
146
147select, flag=aperture, clear;
148
149select, flag=aperture, column=name,s,betx,bety,dx,x,y,n1,apertype,xtol,ytol,rtol;
150
151aperture,offsetelem="offset.tfs",file="ap.tfs";
152
153
154Processing apertures from lhccell$start:1 to lhccell$end:1...
155+++++++ deltap from TWISS            0
156
157 Reading offsets from tfs "offset.tfs"
158
159reference node: lhccell$start:1
160
161  Number of warnings: 0
162
163  ++++++++++++++++++++++++++++++++++++++++++++
164  + MAD-X 5.00.19 (64 bit) finished normally +
165  ++++++++++++++++++++++++++++++++++++++++++++
Note: See TracBrowser for help on using the repository browser.