OPTION, ECHO=FALSE;
OPTION, PSDUMPFREQ=50;
OPTION, TELL=TRUE;

// OPAL Input forLA2line
TFast = 5e-11;
TSlow = 5e-12;
Edes = 0.074649;
gamma = (Edes+PMASS)/PMASS;
beta = sqrt(1-(1/gamma^2));
gambet = gamma*beta;
P0 = gamma*beta*PMASS;
brho = (PMASS*1.0e9*gambet) / CLIGHT;
VALUE,{gamma,brho,Edes,beta,gambet};

rf = 50.6328e6;
NProt = 500000;

DISTRIB1: DISTRIBUTION,
SIGMAX = 0.00415344, SIGMAPX = 0.0395536*gambet,
SIGMAY = 0.00384516, SIGMAPY = 0.0395112*gambet,
SIGMAZ = 0.000180777, SIGMAPZ = 0.0243023*gambet,
CORRX = 0.654264, CORRY = 0.687577, CORRZ = 0.0889277,
R51 = 0, R52 = 0, R61 = 0, R62 = 0,
CUTOFFX = 5, CUTOFFPX = 5, CUTOFFY = 5,
CUTOFFPY = 5, CUTOFFLONG = 5, CUTOFFPZ = 5,
OFFSETX = 0, OFFSETY =0, OFFSETZ = 4.20585,
DISTRIBUTION = GAUSS, INPUTMOUNITS = NONE;

D16: DRIFT, L=0.00465, ELEMEDGE=4.20585;
KMA3i: ECOLLIMATOR, XSIZE=0.0035, YSIZE=0.0035, L=0.002, ELEMEDGE=4.2105;
D17: DRIFT, L=0.071, ELEMEDGE=4.2125;
KMA3o: ECOLLIMATOR, XSIZE=0.0051365, YSIZE=0.0051365, L=0.002, ELEMEDGE=4.2835;
D18: DRIFT, L=0.0045, ELEMEDGE=4.2855;
KMA4i: ECOLLIMATOR, XSIZE=0.0065, YSIZE=0.0065, L=0.002, ELEMEDGE=4.29;
D19: DRIFT, L=0.246, ELEMEDGE=4.292;
KMA4o: ECOLLIMATOR, XSIZE=0.0125, YSIZE=0.0125, L=0.002, ELEMEDGE=4.538;
D20: DRIFT, L=0.1212, ELEMEDGE=4.54;
SMA2: DRIFT, L=0.3876, ELEMEDGE=4.6612;
D21: DRIFT, L=0.1687, ELEMEDGE=5.0488;
MMAP9X: MONITOR, OUTFN="MMAP9X.h5", ELEMEDGE=5.2175;
MMAP10Y: MONITOR, OUTFN="MMAP10Y.h5", ELEMEDGE=5.2175;
D22: DRIFT, L=0.15475, ELEMEDGE=5.2175;
KMAV5: RCOLLIMATOR, XSIZE=0.018, YSIZE=0.04, L=0.0785, ELEMEDGE=5.37225;
D23: DRIFT, L=0.00525, ELEMEDGE=5.45075;
KMA5i: ECOLLIMATOR, XSIZE=0.015, YSIZE=0.015, L=0.002, ELEMEDGE=5.456;
D24: DRIFT, L=0.071, ELEMEDGE=5.458;
KMA5o: ECOLLIMATOR, XSIZE=0.01585, YSIZE=0.01585, L=0.002, ELEMEDGE=5.529;
D25: DRIFT, L=0.286, ELEMEDGE=5.531;
QMA4: QUADRUPOLE, K1=-5.47098, L=0.368, ELEMEDGE=5.817;
QMA4_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=5.817;
D26: DRIFT, L=0.202, ELEMEDGE=6.185;
QMA5: QUADRUPOLE, K1=4.60107, L=0.368, ELEMEDGE=6.387;
QMA5_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=6.387;
D27: DRIFT, L=0.3135, ELEMEDGE=6.755;
MMAP11X: MONITOR, OUTFN="MMAP11X.h5", ELEMEDGE=7.0685;
MMAP12Y: MONITOR, OUTFN="MMAP12Y.h5", ELEMEDGE=7.0685;
D28: DRIFT, L=0.1679, ELEMEDGE=7.0685;
KAMA1A: RCOLLIMATOR, XSIZE=0.045, YSIZE=0.027, L=0.004, ELEMEDGE=7.2364;
D29: DRIFT, L=0.2955, ELEMEDGE=7.2404;
AMS1A: DRIFT, L=0, ELEMEDGE=7.5359;
AMA1: SBEND, K0=0.583629, L=1.83999, ELEMEDGE=7.5359,
         //ANGLE=0.872665, HAPERT=45*2/1000,  K1=0,
         E1=0.425424, E2=0.425424,
  FMAPFN="1DPROFILE1-DEFAULT",
         DESIGNENERGY=7.4649e+07, GAP=32.5*2/1000;
AMA1_r0: ECOLLIMATOR, XSIZE=0.045, YSIZE=0.027, L=1.8997, ELEMEDGE=7.5359;
AMS1B: DRIFT, L=0, ELEMEDGE=9.4356;
D29A: DRIFT, L=0.292, ELEMEDGE=9.4356;
KAMA1B: RCOLLIMATOR, XSIZE=0.045, YSIZE=0.027, L=0.004, ELEMEDGE=9.7276;
KMA6: RCOLLIMATOR, XSIZE=0.019, YSIZE=0.027, L=0.045, ELEMEDGE=9.7316;
D30: DRIFT, L=0.1565, ELEMEDGE=9.7766;
MMAP13X: MONITOR, OUTFN="MMAP13X.h5", ELEMEDGE=9.9331;
MMAP14Y: MONITOR, OUTFN="MMAP14Y.h5", ELEMEDGE=9.9331;
D31: DRIFT, L=0.1175, ELEMEDGE=9.9331;
MMAHL01: ECOLLIMATOR, XSIZE=0.0409, YSIZE=0.0409, L=0.004, ELEMEDGE=10.0506;
D32: DRIFT, L=0.0364, ELEMEDGE=10.0546;
KMA7: RCOLLIMATOR, XSIZE=0.02, YSIZE=0.027, L=0.045, ELEMEDGE=10.091;
D32A: DRIFT, L=0.1154, ELEMEDGE=10.136;
QMA6: QUADRUPOLE, K1=-4.20864, L=0.368, ELEMEDGE=10.2514;
QMA6_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=10.2514;
D33: DRIFT, L=0.202, ELEMEDGE=10.6194;
QMA7: QUADRUPOLE, K1=3.31344, L=0.368, ELEMEDGE=10.8214;
QMA7_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=10.8214;
D34: DRIFT, L=0.4835, ELEMEDGE=11.1894;
MMAP15X: MONITOR, OUTFN="MMAP15X.h5", ELEMEDGE=11.6729;
MMAP16Y: MONITOR, OUTFN="MMAP16Y.h5", ELEMEDGE=11.6729;
D34A: DRIFT, L=0.085, ELEMEDGE=11.6729;
FMA1: RCOLLIMATOR, XSIZE=0.025, YSIZE=0.025, L=0.075, ELEMEDGE=11.7579;
D35: DRIFT, L=0.0775, ELEMEDGE=11.8329;
SMA3: DRIFT, L=0.238, ELEMEDGE=11.9104;
D36: DRIFT, L=0.085, ELEMEDGE=12.1484;
MMAHL02: ECOLLIMATOR, XSIZE=0.0409, YSIZE=0.0409, L=0.004, ELEMEDGE=12.2334;
D37: DRIFT, L=0.164, ELEMEDGE=12.2374;
QMA8: QUADRUPOLE, K1=3.31344, L=0.368, ELEMEDGE=12.4014;
QMA8_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=12.4014;
D38: DRIFT, L=0.202, ELEMEDGE=12.7694;
QMA9: QUADRUPOLE, K1=-4.20864, L=0.368, ELEMEDGE=12.9714;
QMA9_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=12.9714;
D39: DRIFT, L=0.3175, ELEMEDGE=13.3394;
MMAP17X: MONITOR, OUTFN="MMAP17X.h5", ELEMEDGE=13.6569;
MMAP18Y: MONITOR, OUTFN="MMAP18Y.h5", ELEMEDGE=13.6569;
D40: DRIFT, L=0.1625, ELEMEDGE=13.6569;
KAMA2A: RCOLLIMATOR, XSIZE=0.045, YSIZE=0.027, L=0.004, ELEMEDGE=13.8194;
D40A: DRIFT, L=0.332, ELEMEDGE=13.8234;
AMS2A: DRIFT, L=0, ELEMEDGE=14.1554;
AMA2: SBEND, K0=0.583878, L=1.83999, ELEMEDGE=14.1554,
         //ANGLE=0.872665, HAPERT=45*2/1000,  K1=0,
         E1=0.425424, E2=0.425424,
  FMAPFN="1DPROFILE1-DEFAULT",
         DESIGNENERGY=7.4649e+07, GAP=32.5*2/1000;
AMA2_r0: ECOLLIMATOR, XSIZE=0.045, YSIZE=0.027, L=1.8997, ELEMEDGE=14.1554;
AMS2B: DRIFT, L=0, ELEMEDGE=16.0551;
D40B: DRIFT, L=0.3, ELEMEDGE=16.0551;
KAMA2B: RCOLLIMATOR, XSIZE=0.045, YSIZE=0.027, L=0.004, ELEMEDGE=16.3551;
D41: DRIFT, L=0.0389, ELEMEDGE=16.3591;
MMAHL03: ECOLLIMATOR, XSIZE=0.0409, YSIZE=0.0409, L=0.004, ELEMEDGE=16.398;
D42: DRIFT, L=0.202, ELEMEDGE=16.402;
QMA10: QUADRUPOLE, K1=4.59089, L=0.368, ELEMEDGE=16.604;
QMA10_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=16.604;
D43: DRIFT, L=0.202, ELEMEDGE=16.972;
QMA11: QUADRUPOLE, K1=-5.27468, L=0.368, ELEMEDGE=17.174;
QMA11_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=17.174;
D44: DRIFT, L=0.3135, ELEMEDGE=17.542;
MMAP19X: MONITOR, OUTFN="MMAP19X.h5", ELEMEDGE=17.8555;
MMAP20Y: MONITOR, OUTFN="MMAP20Y.h5", ELEMEDGE=17.8555;
D44a: DRIFT, L=0.517, ELEMEDGE=17.8555;
SMA4: DRIFT, L=0.238, ELEMEDGE=18.3725;
D44b: DRIFT, L=0.1685, ELEMEDGE=18.6105;
KFR3: ECOLLIMATOR, XSIZE=0.045, YSIZE=0.045, L=0.004, ELEMEDGE=18.779;
MMAP21X: MONITOR, OUTFN="MMAP21X.h5", ELEMEDGE=18.783;
MMAP22Y: MONITOR, OUTFN="MMAP22Y.h5", ELEMEDGE=18.783;
D45: DRIFT, L=0.221, ELEMEDGE=18.783;
KMA8: ECOLLIMATOR, XSIZE=0.007, YSIZE=0.007, L=0.075, ELEMEDGE=19.004;
D46: DRIFT, L=0.147, ELEMEDGE=19.079;
D47: DRIFT, L=0.078, ELEMEDGE=19.226;
D49: DRIFT, L=0.9785, ELEMEDGE=19.304;
MMAP25X: MONITOR, OUTFN="MMAP25X.h5", ELEMEDGE=20.2825;
MMAP26Y: MONITOR, OUTFN="MMAP26Y.h5", ELEMEDGE=20.2825;
D50: DRIFT, L=0.1175, ELEMEDGE=20.2825;
MMAHL04: ECOLLIMATOR, XSIZE=0.0409, YSIZE=0.0409, L=0.004, ELEMEDGE=20.4;
D51: DRIFT, L=0.192, ELEMEDGE=20.404;
QMA12: QUADRUPOLE, K1=-4.70808, L=0.368, ELEMEDGE=20.596;
QMA12_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=20.596;
D52: DRIFT, L=0.202, ELEMEDGE=20.964;
QMA13: QUADRUPOLE, K1=4.56473, L=0.368, ELEMEDGE=21.166;
QMA13_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=21.166;
D53: DRIFT, L=0.251, ELEMEDGE=21.534;
D53A: DRIFT, L=0.4392, ELEMEDGE=21.785;
AMA3: DRIFT, L=1.29679, ELEMEDGE=22.2242;
D53B: DRIFT, L=0.3296, ELEMEDGE=23.521;
KAMA3B: RCOLLIMATOR, XSIZE=0.045, YSIZE=0.026, L=0.004, ELEMEDGE=23.8506;
D54: DRIFT, L=0.2325, ELEMEDGE=23.8546;
MMBP1X: MONITOR, OUTFN="MMBP1X.h5", ELEMEDGE=24.0871;
MMBP2Y: MONITOR, OUTFN="MMBP2Y.h5", ELEMEDGE=24.0871;
D56: DRIFT, L=0.16, ELEMEDGE=24.0871;
Myl1: DRIFT, L=0.000275, ELEMEDGE=24.2471;
D56A: DRIFT, L=0.065, ELEMEDGE=24.2474;
KMB0: DRIFT, L=0.075, ELEMEDGE=24.3124;
Myl2: DRIFT, L=0.2, ELEMEDGE=24.3874;
D56C: DRIFT, L=0.0325, ELEMEDGE=24.5874;
BMB1: ECOLLIMATOR, XSIZE=0.045, YSIZE=0.045, L=0.004, ELEMEDGE=24.6199;
D57: DRIFT, L=0.1922, ELEMEDGE=24.6239;
SMB1: DRIFT, L=0.3876, ELEMEDGE=24.8161;
D58: DRIFT, L=1.1639, ELEMEDGE=25.2037;
MMBP3X: MONITOR, OUTFN="MMBP3X.h5", ELEMEDGE=26.3676;
MMBP4Y: MONITOR, OUTFN="MMBP4Y.h5", ELEMEDGE=26.3676;
D58A: DRIFT, L=0.1175, ELEMEDGE=26.3676;
MMBHL01: ECOLLIMATOR, XSIZE=0.0409, YSIZE=0.0409, L=0.004, ELEMEDGE=26.4851;
D60: DRIFT, L=0.197, ELEMEDGE=26.4891;
QMB1: QUADRUPOLE, K1=-4.91772, L=0.368, ELEMEDGE=26.6861;
QMB1_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=26.6861;
D61: DRIFT, L=0.202, ELEMEDGE=27.0541;
QMB2: QUADRUPOLE, K1=4.38547, L=0.368, ELEMEDGE=27.2561;
QMB2_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=27.2561;
D62: DRIFT, L=1.6914, ELEMEDGE=27.6241;
SMB2: DRIFT, L=0.3876, ELEMEDGE=29.3155;
D63: DRIFT, L=0.1687, ELEMEDGE=29.7031;
MMBP5X: MONITOR, OUTFN="MMBP5X.h5", ELEMEDGE=29.8718;
MMBP6Y: MONITOR, OUTFN="MMBP6Y.h5", ELEMEDGE=29.8718;
D64: DRIFT, L=0.4093, ELEMEDGE=29.8718;
KMBV1A: ECOLLIMATOR, XSIZE=0.045, YSIZE=0.045, L=0.4, ELEMEDGE=30.2811;
D65: DRIFT, L=0.085, ELEMEDGE=30.6811;
D66: DRIFT, L=0.9, ELEMEDGE=30.7661;
D67: DRIFT, L=0.9, ELEMEDGE=31.6661;
MMBHL02: ECOLLIMATOR, XSIZE=0.0409, YSIZE=0.0409, L=0.004, ELEMEDGE=32.5661;
D69: DRIFT, L=0.196, ELEMEDGE=32.5701;
QMB3: QUADRUPOLE, K1=-3.67009, L=0.368, ELEMEDGE=32.7661;
QMB3_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=32.7661;
D70: DRIFT, L=0.202, ELEMEDGE=33.1341;
QMB4: QUADRUPOLE, K1=6.26124, L=0.368, ELEMEDGE=33.3361;
QMB4_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=33.3361;
D71: DRIFT, L=0.202, ELEMEDGE=33.7041;
QMB5: QUADRUPOLE, K1=-3.333, L=0.368, ELEMEDGE=33.9061;
QMB5_r0: ECOLLIMATOR, XSIZE=0.0465, YSIZE=0.0465, L=0.368, ELEMEDGE=33.9061;
D71A: DRIFT, L=0.3175, ELEMEDGE=34.2741;
MMBP7X: MONITOR, OUTFN="MMBP7X.h5", ELEMEDGE=34.5916;
MMBP8Y: MONITOR, OUTFN="MMBP8Y.h5", ELEMEDGE=34.5916;
D72: DRIFT, L=0.1687, ELEMEDGE=34.5916;
SMB3: DRIFT, L=0.3876, ELEMEDGE=34.7603;
D73: DRIFT, L=1.0442, ELEMEDGE=35.1479;
KMB1: ECOLLIMATOR, XSIZE=0.015, YSIZE=0.015, L=0.075, ELEMEDGE=36.1921;
D74A: DRIFT, L=0.044, ELEMEDGE=36.2671;
Myl3: DRIFT, L=0.000122, ELEMEDGE=36.3111;
D74B: DRIFT, L=0.02, ELEMEDGE=36.3112;
MMBP9X: MONITOR, OUTFN="MMBP9X.h5", ELEMEDGE=36.3312;
MMBP10Y: MONITOR, OUTFN="MMBP10Y.h5", ELEMEDGE=36.3312;
D74: DRIFT, L=0.043, ELEMEDGE=36.3312;
D75: DRIFT, L=0.015, ELEMEDGE=36.3742;
KMB2i: ECOLLIMATOR, XSIZE=0.004, YSIZE=0.004, L=0.004, ELEMEDGE=36.3892;
UPG2A: DRIFT, L=0.0335, ELEMEDGE=36.3932;
UPG2X: MONITOR, OUTFN="UPG2X.h5", ELEMEDGE=36.4267;
UPG2Y: MONITOR, OUTFN="UPG2Y.h5", ELEMEDGE=36.4267;
END: DRIFT, L=0.2, ELEMEDGE=36.4267;

BEAMLINE_B: LINE=(D16,KMA3i,D17,KMA3o,D18,KMA4i,D19,KMA4o,D20,SMA2,D21,MMAP9X,MMAP10Y,D22,KMAV5,D23,KMA5i,D24,KMA5o,D25,QMA4,QMA4_r0,D26,QMA5,QMA5_r0,D27,MMAP11X,MMAP12Y,D28,KAMA1A,D29,AMS1A,AMA1,AMA1_r0,AMS1B,D29A,KAMA1B,KMA6,D30,MMAP13X,MMAP14Y,D31,MMAHL01,D32,KMA7,D32A,QMA6,QMA6_r0,D33,QMA7,QMA7_r0,D34,MMAP15X,MMAP16Y,D34A,FMA1,D35,SMA3,D36,MMAHL02,D37,QMA8,QMA8_r0,D38,QMA9,QMA9_r0,D39,MMAP17X,MMAP18Y,D40,KAMA2A,D40A,AMS2A,AMA2,AMA2_r0,AMS2B,D40B,KAMA2B,D41,MMAHL03,D42,QMA10,QMA10_r0,D43,QMA11,QMA11_r0,D44,MMAP19X,MMAP20Y,D44a,SMA4,D44b,KFR3,MMAP21X,MMAP22Y,D45,KMA8,D46,D47,D49,MMAP25X,MMAP26Y,D50,MMAHL04,D51,QMA12,QMA12_r0,D52,QMA13,QMA13_r0,D53,D53A,AMA3,D53B,KAMA3B,D54,MMBP1X,MMBP2Y,D56,Myl1,D56A,KMB0,Myl2,D56C,BMB1,D57,SMB1,D58,MMBP3X,MMBP4Y,D58A,MMBHL01,D60,QMB1,QMB1_r0,D61,QMB2,QMB2_r0,D62,SMB2,D63,MMBP5X,MMBP6Y,D64,KMBV1A,D65,D66,D67,MMBHL02,D69,QMB3,QMB3_r0,D70,QMB4,QMB4_r0,D71,QMB5,QMB5_r0,D71A,MMBP7X,MMBP8Y,D72,SMB3,D73,KMB1,D74A,Myl3,D74B,MMBP9X,MMBP10Y,D74,D75,KMB2i,UPG2A,UPG2X,UPG2Y,END);

FS1:FIELDSOLVER, FSTYPE=NONE, MX=64, MY=64, MT=64, PARFFTX=true, PARFFTY=true, PARFFTT=true, BCFFTX=open, BCFFTY=open, BCFFTT=open, BBOXINCR=1, GREENSF=STANDARD;
BEAM_B: BEAM, PARTICLE=PROTON, PC=P0, NPART=NProt, BCURRENT=NProt*1.6e-19*rf, BFREQ=rf,  CHARGE=1;

STEP0= (5.631- 4.20585)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP0, DT=TSlow;//ZSTOP=5.631;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP1= (6.955- 5.631)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP1, DT=TFast;//ZSTOP=6.955;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP2= (7.3404- 6.955)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP2, DT=TSlow;//ZSTOP=7.3404;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP3= (9.8266- 7.3404)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP3, DT=TFast;//ZSTOP=9.8266;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP4= (10.2014- 9.8266)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP4, DT=TSlow;//ZSTOP=10.2014;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP5= (11.4394- 10.2014)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP5, DT=TFast;//ZSTOP=11.4394;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP6= (12.3014- 11.4394)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP6, DT=TSlow;//ZSTOP=12.3014;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP7= (13.5569- 12.3014)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP7, DT=TFast;//ZSTOP=13.5569;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP8= (13.9554- 13.5569)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP8, DT=TSlow;//ZSTOP=13.9554;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP9= (16.2551- 13.9554)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP9, DT=TFast;//ZSTOP=16.2551;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP10= (20.646- 16.2551)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP10, DT=TSlow;//ZSTOP=20.646;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP11= (23.9546- 20.646)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP11, DT=TFast;//ZSTOP=23.9546;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP12= (24.5699- 23.9546)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP12, DT=TSlow;//ZSTOP=24.5699;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP13= (26.3851- 24.5699)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP13, DT=TFast;//ZSTOP=26.3851;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP14= (29.2655- 26.3851)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP14, DT=TSlow;//ZSTOP=29.2655;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP15= (32.8661- 29.2655)/(beta*CLIGHT*TFast);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP15, DT=TFast;//ZSTOP=32.8661;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STEP16= (36.6267- 32.8661)/(beta*CLIGHT*TSlow);
SELECT, LINE=BEAMLINE_B;
TRACK, LINE=BEAMLINE_B, BEAM=BEAM_B, MAXSTEPS= STEP16, DT=TSlow;//ZSTOP=34.2741;
RUN, METHOD = "PARALLEL-T", BEAM=BEAM_B, FIELDSOLVER=FS1, DISTRIBUTION=DISTRIB1;
ENDTRACK;

STOP;