+++++++++++++++++++++++++++++++++++++++++++ + MAD-X 5.00.19 (64 bit, Darwin) + + DEBUG Version - use at own risk! + + Code Modification Date: 2012.12.11 + + Execution Time Stamp: 05.02.13 17.31.54 + +++++++++++++++++++++++++++++++++++++++++++ Option, warn,info,echo; Title, "LHC 6.503 with errors and two beam orbit correction using MADX"; // CALL IN SELECTION, SEQUENCE and OPTICS FILES; Option, -echo,-warn,-info; Beam, particle=proton, sequence=lhcb1, energy=7000.0, NPART=1.05E11, sige=4.5e-4; Beam, particle=proton, sequence=lhcb2, energy=7000.0, NPART=1.05E11, sige=4.5e-4, bv=-1; ! the 'bv' parameter indicates that beam2 rotates counter-clockwise // Removing classes of correctors select, flag=seqedit, clear; select, flag=seqedit, class=MU; select, flag=seqedit, class=MKI; select, flag=seqedit, class=MKI__192; select, flag=seqedit, class=MKI__193; select, flag=seqedit, class=MKA; select, flag=seqedit, class=MKD; select, flag=seqedit, class=MKQ; select, flag=seqedit, class=MKQA; select, flag=seqedit, class=MSDA; select, flag=seqedit, class=MSDB; select, flag=seqedit, class=MSDB2; select, flag=seqedit, class=MSDC; select, flag=seqedit, class=MSIA; select, flag=seqedit, class=MSIB; select, flag=seqedit, class=MCBV; select, flag=seqedit, class=MCBH; select, flag=seqedit, class=MCBCA; select, flag=seqedit, class=MCBCB; select, flag=seqedit, class=MCBCC; select, flag=seqedit, class=MCBCD; select, flag=seqedit, class=MCBCH; select, flag=seqedit, class=MCBCV; select, flag=seqedit, class=MCBYA; select, flag=seqedit, class=MCBYB; select, flag=seqedit, class=MCBYH; select, flag=seqedit, class=MCBYV; select, flag=seqedit, class=MCBWH; select, flag=seqedit, class=MCBWV; seqedit, sequence=lhcb1; remove, element=selected; seqedit, sequence=lhcb2; remove, element=selected; endedit; select, flag=seqedit, clear; select, flag=seqedit, class=MBXWH; select, flag=seqedit, class=MBXWS; select, flag=seqedit, class=MBXWT; select, flag=seqedit, class=MBAW; select, flag=seqedit, class=MBLW; select, flag=seqedit, class=MBWMD; seqedit, sequence=lhcb1; remove, element=selected; seqedit, sequence=lhcb1; remove, element=selected; seqedit, sequence=lhcb2; remove, element=selected; seqedit, sequence=lhcb2; remove, element=selected; endedit; // assign alignment errors in ring 1 and ring 2 at Q1 Use, period=lhcb2; Select, flag=ERROR, clear; Select, flag=ERROR, range=MQXA.1L1; ealign, dy:=0.100e-3, dx:=0.070e-3; Assigned alignment errors to 1 elements Select, flag=ERROR, clear; Select, flag=ERROR, range=MQXA.1R5; ealign, dy:=0.050e-3, dx:=0.035e-3; Assigned alignment errors to 1 elements Use, period=lhcb1; Select, flag=ERROR, clear; Select, flag=ERROR, range=MQXA.1L1; ealign, dy:=0.100e-3, dx:=0.070e-3; Assigned alignment errors to 1 elements Select, flag=ERROR, clear; Select, flag=ERROR, range=MQXA.1R5; ealign, dy:=0.050e-3, dx:=0.035e-3; Assigned alignment errors to 1 elements // enable only MCBX at Q1 USEKICK,CLASS=mcbx, status=OFF; USEKICK,CLASS=mcbxv, status=OFF; USEKICK,CLASS=mcbxh, status=OFF; USEKICK,CLASS=mcbxa, status=OFF; USEKICK,class=mbxws, status=OFF; USEKICK,class=mbxwt, status=OFF; USEKICK,class=mbaw, status=OFF; USEKICK,class=mbwmd, status=OFF; !USEKICK,RANGE=mcbx.3l1, status=ON; !USEKICK,RANGE=mcbx.3r1, status=ON; !USEKICK,RANGE=mcbx.3l5, status=ON; !USEKICK,RANGE=mcbx.3r5, status=ON; !USEKICK,RANGE=mcbx.2l1, status=ON; !USEKICK,RANGE=mcbx.2r1, status=ON; !USEKICK,RANGE=mcbx.2l5, status=ON; !USEKICK,RANGE=mcbx.2r5, status=ON; !USEKICK,RANGE=mcbx.1l1, status=ON; !USEKICK,RANGE=mcbx.1r1, status=ON; !USEKICK,RANGE=mcbx.1l5, status=ON; !USEKICK,RANGE=mcbx.1r5, status=ON; USEKICK,RANGE=mcbxv.3l1, status=ON; USEKICK,RANGE=mcbxv.3r1, status=ON; USEKICK,RANGE=mcbxv.2l1, status=ON; USEKICK,RANGE=mcbxv.2r1, status=ON; USEKICK,RANGE=mcbxv.1l1, status=ON; USEKICK,RANGE=mcbxv.1r1, status=ON; USEKICK,RANGE=mcbxv.3l5, status=ON; USEKICK,RANGE=mcbxv.2l5, status=ON; USEKICK,RANGE=mcbxv.1l5, status=ON; USEKICK,RANGE=mcbxv.1r5, status=ON; USEKICK,RANGE=mcbxh.3l1, status=ON; USEKICK,RANGE=mcbxh.3r1, status=ON; USEKICK,RANGE=mcbxh.2l1, status=ON; USEKICK,RANGE=mcbxh.2r1, status=ON; USEKICK,RANGE=mcbxh.1l1, status=ON; USEKICK,RANGE=mcbxh.1r1, status=ON; USEKICK,RANGE=mcbxh.3l5, status=ON; USEKICK,RANGE=mcbxh.2l5, status=ON; USEKICK,RANGE=mcbxh.1l5, status=ON; USEKICK,RANGE=mcbxh.1r5, status=ON; // calculate TWISS table for correction, beam 1 on_x1 = 0.0; on_x2 = 0.0; on_x5 = 0.0; on_x8 = 0.0; on_sep1 = 0.0; on_sep2 = 0.0; on_sep5 = 0.0; on_sep8 = 0.0; on_alice = 0.0; on_lhcb = 0.0; select, flag=twiss, clear; select, flag=twiss, column=name,s,x,px,y,py,betx,bety,alfx,alfy; TWISS,table=twb1,sequence=lhcb1; enter Twiss module iteration: 1 error: 1.298519E-04 deltap: 0.000000E+00 orbit: -2.809322E-05 8.613956E-05 7.403069E-05 5.714884E-06 0.000000E+00 0.000000E+00 iteration: 2 error: 3.572254E-07 deltap: 0.000000E+00 orbit: -2.814728E-05 8.594499E-05 7.406703E-05 5.662320E-06 0.000000E+00 0.000000E+00 ++++++ table: summ length orbit5 alfa gammatr 26658.8832 -0 0.0003213810617 55.78145773 q1 dq1 betxmax dxmax 64.31061633 2.264495632 4400.02675 2.824586493 dxrms xcomax xcorms q2 1.400840729 0.004315142148 0.0008355982316 59.31981837 dq2 betymax dymax dyrms 2.000703717 4416.620285 0.5220811909 0.07296361243 ycomax ycorms deltap synch_1 0.005046416637 0.0007531963962 0 0 synch_2 synch_3 synch_4 synch_5 0 0 0 0 !plot,table=twb1,title="Orbit ..L1..before",vmin=-0.005,vmax=0.005,noline=true,haxis=s,vaxis=x; coption,print=1; !CORRECT, FLAG=ring, ! MODE=micado,cond=0,NCORR= 1,ERROR=1E-6,PLANE=y, ! mlist='my1.tab',clist='cy1.tab', RESOUT=1,corzero=1; USEKICK,RANGE=mcbxv.1l5, status=ON; // calculate TWISS table for correction, beam 2 on_x1 = 0.0; on_x2 = 0.0; on_x5 = 0.0; on_x8 = 0.0; on_sep1 = 0.0; on_sep2 = 0.0; on_sep5 = 0.0; on_sep8 = 0.0; on_alice = 0.0; on_lhcb = 0.0; select, flag=twiss, clear; select, flag=twiss, column=name,s,x,px,y,py,betx,bety,alfx,alfy; TWISS,table=twb2,sequence=lhcb2; enter Twiss module iteration: 1 error: 9.575207E-05 deltap: 0.000000E+00 orbit: 3.027800E-05 -1.061548E-04 -5.427964E-05 2.545669E-05 0.000000E+00 0.000000E+00 iteration: 2 error: 1.482572E-07 deltap: 0.000000E+00 orbit: 3.018352E-05 -1.061182E-04 -5.422462E-05 2.539827E-05 0.000000E+00 0.000000E+00 ++++++ table: summ length orbit5 alfa gammatr 26658.8832 -0 0.0003209203541 55.82148281 q1 dq1 betxmax dxmax 64.30982579 2.082564444 4399.172405 2.979380862 dxrms xcomax xcorms q2 1.402695416 0.005881468084 0.0009328547335 59.32023759 dq2 betymax dymax dyrms 1.999899436 4406.690664 0.701180927 0.1101403134 ycomax ycorms deltap synch_1 0.004717597621 0.0008591074516 0 0 synch_2 synch_3 synch_4 synch_5 0 0 0 0 !plot,table=twb2,title="Orbit ..L2..before",vmin=-0.005,vmax=0.005,noline=true,haxis=s,vaxis=x; CORRECT, FLAG=ring, tworing,beam1tab=twb1,beam2tab=twb2, MODE=micado,cond=0,NCORR= 4,ERROR=1E-6,PLANE=x, mlist='mx12.tab',clist='cx12.tab', RESOUT=1,corzero=1; Want to use orbits from: twb1 and : twb2 for two beams orbit corrections ... Want to use orbit from: twb1 Want to use orbit from: twb2 mons and corrs (beam 1) : 544 0 mons and corrs (beam 2) : 545 0 mons and corrs (beam 1+2) : 2 24 1091 monitors and 24 correctors found in input 1091 monitors and 10 correctors enabled enter MICADO correction ... start MICADO correction with 4 correctors Back from micado 0 RMS before 9.588667e-01 RMS after 8.998220e-02 mcbxh.1r1:1 -1.115765e-03 mcbxh.2r1:1 0.000000e+00 mcbxh.3r1:1 0.000000e+00 mcbxh.3l5:1 0.000000e+00 mcbxh.2l5:1 0.000000e+00 mcbxh.1l5:1 5.766483e-04 mcbxh.1r5:1 -1.084434e-03 mcbxh.3l1:1 0.000000e+00 mcbxh.2l1:1 0.000000e+00 mcbxh.1l1:1 2.271659e-03 CORRECTION SUMMARY: rms before correction: 0.958867 mm rms after correction: 0.089982 mm ptp before correction: 10.060801 mm ptp after correction: 2.637795 mm Max strength: 2.271659e-03 should be less than 1.000000e+00 select, flag=twiss, clear; select, flag=twiss, column=name,s,x,px,y,py,betx,bety,alfx,alfy; TWISS,sequence=lhcb1; enter Twiss module iteration: 1 error: 1.295291E-04 deltap: 0.000000E+00 orbit: -2.731998E-05 -7.415126E-07 7.386753E-05 5.793822E-06 0.000000E+00 0.000000E+00 iteration: 2 error: 1.924924E-07 deltap: 0.000000E+00 orbit: -2.742854E-05 -6.552875E-07 7.390359E-05 5.661291E-06 0.000000E+00 0.000000E+00 ++++++ table: summ length orbit5 alfa gammatr 26658.8832 -0 0.0003213813328 55.7814342 q1 dq1 betxmax dxmax 64.30993506 2.001003829 4403.339372 2.858350735 dxrms xcomax xcorms q2 1.398919814 9.324377601e-05 9.689405509e-06 59.32005551 dq2 betymax dymax dyrms 1.999882205 4404.179411 0.52015192 0.07282880219 ycomax ycorms deltap synch_1 0.005045837468 0.0007529500667 0 0 synch_2 synch_3 synch_4 synch_5 0 0 0 0 !plot,title="Orbit ..L1..after",vmin=-0.005,vmax=0.005,noline=true,haxis=s,vaxis=x; select, flag=twiss, clear; select, flag=twiss, column=name,s,x,px,y,py,betx,bety,alfx,alfy; TWISS,sequence=lhcb2; enter Twiss module iteration: 1 error: 8.988769E-05 deltap: 0.000000E+00 orbit: 3.930653E-05 -8.982999E-07 -5.419899E-05 2.550666E-05 0.000000E+00 0.000000E+00 iteration: 2 error: 2.080341E-07 deltap: 0.000000E+00 orbit: 3.932110E-05 -8.382999E-07 -5.420823E-05 2.537128E-05 0.000000E+00 0.000000E+00 ++++++ table: summ length orbit5 alfa gammatr 26658.8832 -0 0.0003209217334 55.82136285 q1 dq1 betxmax dxmax 64.30994132 1.996885918 4403.835848 2.886331755 dxrms xcomax xcorms q2 1.398792227 0.0001242991045 1.235064602e-05 59.32005207 dq2 betymax dymax dyrms 2.002444267 4403.828262 0.7037915936 0.1103657192 ycomax ycorms deltap synch_1 0.004716170502 0.0008584970928 0 0 synch_2 synch_3 synch_4 synch_5 0 0 0 0 !plot,title="Orbit ..L2..after",vmin=-0.005,vmax=0.005,noline=true,haxis=s,vaxis=x; Number of warnings: 0 ++++++++++++++++++++++++++++++++++++++++++++ + MAD-X 5.00.19 (64 bit) finished normally + ++++++++++++++++++++++++++++++++++++++++++++