forked from lhcopt/lhcmask
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsubmodule_01e_final.madx
67 lines (57 loc) · 2.69 KB
/
submodule_01e_final.madx
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
print, text="";
print, text="";
print, text=" --- Submodule 1e: final";
print, text=" ----------------------";
print, text="";
print, text="";
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++;
! STORE the optics
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++;
exec, crossing_disable;
if (par_verbose==1){
if (mylhcbeam<3){
use, sequence=lhcb1;
select, flag=twiss,clear;
select, flag=twiss, pattern="IP1", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
select, flag=twiss, pattern="IP2", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
select, flag=twiss, pattern="IP5", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
select, flag=twiss, pattern="IP8", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
twiss, sequence=lhcb1, file=twiss_crossing_disable.b1;
twiss, sequence=lhcb1, table=twisscrossingdisableb1; ! handle for further analysis
system, "cat twiss_crossing_disable.b1";
};
use, sequence=lhcb2;
select, flag=twiss,clear;
select, flag=twiss, pattern="IP1", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
select, flag=twiss, pattern="IP2", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
select, flag=twiss, pattern="IP5", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
select, flag=twiss, pattern="IP8", column=name,s,betx,bety,alfx,alfy,dx,dpx,mux,muy;
twiss, sequence=lhcb2, file=twiss_crossing_disable.b2;
twiss, sequence=lhcb2, table=twisscrossingdisableb2; ! handle for further analysis
system, "cat twiss_crossing_disable.b2";
};
! Switch on the crossing scheme
exec, crossing_restore;
if (par_verbose==1){
if (mylhcbeam<3){
use, sequence=lhcb1;
select,flag=twiss,clear;
select, flag=twiss, pattern="IP1", column=name,s,x,y,px,py;
select, flag=twiss, pattern="IP2", column=name,s,x,y,px,py;
select, flag=twiss, pattern="IP5", column=name,s,x,y,px,py;
select, flag=twiss, pattern="IP8", column=name,s,x,y,px,py;
twiss, sequence=lhcb1, file=twiss_crossing_enable.b1;
twiss, sequence=lhcb1, table=twisscrossingenableb1; ! handle for further analysis
system, "cat twiss_crossing_enable.b1";
};
use, sequence=lhcb2;
select, flag=twiss,clear;
select, flag=twiss, pattern="IP1", column=name,s,x,y,px,py;
select, flag=twiss, pattern="IP2", column=name,s,x,y,px,py;
select, flag=twiss, pattern="IP5", column=name,s,x,y,px,py;
select, flag=twiss, pattern="IP8", column=name,s,x,y,px,py;
twiss, sequence=lhcb2, file=twiss_crossing_enable.b2;
twiss, sequence=lhcb2, table=twisscrossingenableb2; ! handle for further analysis
system, "cat twiss_crossing_enable.b2";
};
if(mylhcbeam==1) {use, sequence=lhcb1;}; if(mylhcbeam>1) {use, sequence=lhcb2;};