option, beampipeRadius = 8.21 * mm,
        boxSize = 0.5 * m,
        tunnelRadius= 2 * m,
        beampipeThickness = 0.2 * cm,
        chordStepMinimum = 1e-10 * m,
        deltaIntersection = 1e-6 * m,
        deltaChord = 1e-8 * m,
        lengthSafety = 0 * m,
        thresholdCutCharged = 10 * GeV,
        thresholdCutPhotons = 1000 * MeV,
        randomSeed=-1,
        storeTrajectory=0,
        stopTracks=0;


beam,   particle="e-",
        energy=1496 * GeV,
        nparticles=1e+3,
        distrType="guineapig_bunch",
        distrFile="halo.dat";
!	distrFile="particles.in";

kod = 0.0e2;
lod = 2.0;
myscale = 5.0;
mybendfac = 0.083333333333;
boctf = 0.0;
boctd = 1.0e3;
coll_efact = 0.5;
bsign = -1.0;
coll_e0 = 5.0e2*coll_efact;
coll_brho = coll_e0*1.0e10/2.99792458e8;
emitxn = 0.68e-6;
emityn = 0.2e-7;
bleng = 0.9e-4;
esprd = 0.265e-2;
tbetx = 31.276142753063;
talfx = 0.799862640848e-2;
tdx = 0.0;
tdpx = 0.0;
tbety = 6.841720075737;
talfy = 0.035712458697;
lmm = 0.25;
lscale = 1.0;
etascale = 1.0;
lqq = 60.0*lscale;
l_linq = 0.9779;
a_linq = 0.635e-2;
qlinfact = 1.0/(a_linq*coll_brho)/myscale/myscale;
l_ebm = 2.0;
a_ebm = 0.01;
qebmfact = 1.0/(a_ebm*coll_brho)/myscale/myscale;
l_q_ecol = 1.0;
a_q_ecol = 0.0254;
qecfact = 1.0/(a_q_ecol*coll_brho)/myscale/myscale;
l_q_ecmi = 1.0;
a_q_ecmi = 0.015;
qmifact = 1.0/(a_q_ecmi*coll_brho)/myscale/myscale;
l_s_ecol = 1.0;
a_s_ecol = 0.0254;
lbtot_ec = lqq-2.0*l_s_ecol-15.0*lmm;
secfact = 2.0/(a_s_ecol*a_s_ecol*coll_brho);
l_s_dglg = 0.5;
a_s_dglg = 0.01;
sdogfact = 2.0/(a_s_dglg*a_s_dglg*coll_brho);
l_bend = lbtot_ec/12.0;
a_ecol = 0.0254;
bd300 = 7.811696283042*coll_efact;
bd301 = -8.573353230276*coll_efact;
bd302 = -bd301;
bd303 = bd301;
bd304 = -bd301;
bd305 = bd301;
bd306 = -bd301;
bd307 = bd301;
bd308 = -bd301;
bd309 = bd301;
bd310 = -bd301;
bd311 = bd301;
bf1ec = 6.910540869978*coll_efact;
bd1ec = -8.086730540558*coll_efact;
bf2ec = 6.892954704438*coll_efact;
bf3ec = 3.05604346908*coll_efact;
ld3bc = 23.26703;
ld4bc = 0.385153;
bf1bc = 53.8414;
bd1bc = -43.86097;
bf2bc = 14.16959;
bfdb1 = 3.86014*coll_efact;
bddb1 = -5.32672*coll_efact;
bddb2 = -2.35471*coll_efact;
bddb3 = -1.23565*coll_efact;
lmi = 19.8249;
bb1 = 0.44648899*etascale/(lscale*lscale)*bsign*coll_efact;
bb2 = -0.09911677*etascale/(lscale*lscale)*bsign*coll_efact;
bb3 = -bb1;
bb4 = -bb2;
bqec = 4.7354188/lscale*coll_efact;
bsyec1 = 0.018897;
bsxec1 = -0.935639e-2;
bsxec2 = -bsxec1;
bsyec2 = -bsyec1;
aoctecol = 0.015;
loctecol = 1.0;
ecolofac = 6.0*loctecol/(aoctecol*aoctecol*aoctecol*coll_brho);
boctec1 = -4.35185*coll_efact;
boctec2 = 4.30281*coll_efact;
boctec3 = -boctec1;
boctec4 = -boctec2;
bmi = 8.25416*coll_efact;
x_aper = 0.08e-3;
y_aper = 0.102e-3;
bfbcol = 3.67714*coll_efact;
bdbcol = -5.66162*coll_efact;
bcbcol = -5.05905*coll_efact;
bpole = 4.71383*coll_efact;
xf1 = -0.198257e-2*bsign;
xd1 = 0.221897e-2*bsign;
xf2 = -xf1;
xd2 = -xd1;
xmf1 = -0.345352e-2*bsign;
xmf2 = 0.800993e-3*bsign;
xmd = 0.642315e-3*bsign;
xsf1 = -xmf1;
xsf2 = -xmf2;
xsd = -xmd;
qebm0fact = 1.0/(a_ebm*coll_brho);
bsxdl1 = 0.0*bsign*coll_efact;
bsydl1 = 0.0*bsign*coll_efact;
bsxdl2 = 0.0*bsign*coll_efact;
bsydl2 = 0.0*bsign*coll_efact;
bsxdl3 = -bsxdl1;
bsydl3 = -bsydl1;
bsxdl4 = -bsxdl2;
bsydl4 = -bsydl2;
lfbc = 0.0;
lspace = 0.0;
lfbk = 0.0;
lws = 0.0;
lltcd4 = 0.7072719;
lltcd3 = 0.6119686;
lltcd2 = 15.48765;
lltcd1 = 3.551866;
lltcd0 = 30.69051;
kqltc4 = 0.02054873;
kqltc3 = -0.04256154;
kqltc2 = 0.01989328;
kqltc1 = -0.02502006;
lbtfd4 = 0.0525519;
lbtfd3 = 19.97658;
lbtfd2 = 19.9951;
lbtfd1 = 8.712693;
lbtfd0 = 0.0507269;
kbtfq4 = -0.666287e-2;
kbtfq3 = 0.5048062e-2;
kbtfq2 = 0.3044046e-2;
kbtfq1 = -0.01209993;
kmax = 1.5e2;
bxmax = 1.5e3;
bymax = 5.0e2;
lmax = 20.0;
lsca = 1.0;
tscal = 1.0;
nbfac = 5.0;
mfac = 3.4;
gama = 0.0;
efact = 1.0;
e0 = 1.5e3*efact;
cb = 33.356409519815;
brho = cb*e0/3.0;
ft2m = 0.3048;
lffb4 = 21.031;
lffb3 = 2.95e2*ft2m;
lffb3b = 60.0*ft2m;
lffb2 = 60.0*ft2m;
lffb1 = 60.0*ft2m;
lffsb1 = 60.0*ft2m;
lffsb1a = 60.0*ft2m;
lffsb1b = 60.0*ft2m;
bffb4 = -0.046*bsign*efact;
bffb3 = 0.14313*bsign*efact;
bffb3b = 0.14313*bsign*efact;
bffb1 = 0.043*bsign*efact;
bffb2 = 0.14313*bsign*efact;
bffsb1 = 0.043*bsign*efact;
bffsb1a = 0.043*bsign*efact;
bffsb1b = 0.043*bsign*efact;
affb4 = bffb4*lffb4/brho;
affb3 = bffb3*lffb3/brho;
affb3b = bffb3b*lffb3b/brho;
affb2 = bffb2*lffb2/brho;
affb1 = bffb1*lffb1/brho;
affsb1 = bffsb1*lffsb1/brho;
affsb1a = bffsb1a*lffsb1a/brho;
affsb1b = bffsb1b*lffsb1b/brho;
thetab4 = affb4/2.0/nbfac/mfac;
thetab3 = affb3/2.0/nbfac/mfac;
thetab3b = affb3b/2.0/nbfac/mfac;
thetab2 = affb2/2.0/nbfac/mfac;
thetab1 = affb1/2.0/nbfac/mfac;
thetasb1 = affsb1/2.0/nbfac/mfac;
thetasb1a = affsb1a/2.0/nbfac/mfac;
thetasb1b = affsb1b/2.0/nbfac/mfac;
bqmd5 = 3.953358864681*efact;
bqmd6 = -5.975913600923*efact;
bqmd7 = 6.335101647528*efact;
bqmd8 = -6.686048673284*efact;
bqmd9 = -8.816988086064*efact;
bqmd10 = 9.52435610138*efact;
bqmd11 = 4.689561065089;
bqmd12 = -5.753341;
bqmd13 = 4.140486830511;
bqmd14 = -3.615914513634;
ltrick = 0.0;
l12mm = 0.012;
l25cm = 0.25;
l50cm = 0.5;
l1m = 1.0;
l2m = 2.0;
aperq1 = l12mm/2.0;
qfact1 = 1.0/(aperq1*brho);
aperq2 = l12mm/2.0;
qfact2 = 1.0/(aperq2*brho);
aperq3 = l12mm/2.0;
qfact3 = 1.0/(aperq3*brho);
bqdf90 = 9.32501092109*efact;
bqdd90 = -bqdf90;
bqdf2 = 6.381144203992*efact;
bqdd2 = -4.782835045388*efact;
bsq1 = 0.0;
bsq2 = 0.0;
bsq3 = 0.0;
bsq4 = 0.0;
bqdd = -8.460739400782*efact;
bqdf = -bqdd;
lqf8 = 6.5*ft2m;
lqd7 = 6.5*ft2m;
lqd6c = 6.5*ft2m;
lqd6b = 6.5*ft2m;
lqf5b = 6.5*ft2m;
lqf5a = 6.5*ft2m;
lqd4b = 6.5*ft2m;
lqd4a = 6.5*ft2m;
lqf3b = 6.5*ft2m;
lqf3a = 6.5*ft2m;
lqd2 = 6.5*ft2m;
lqf1 = 13.11*ft2m;
lqd0 = 11.0*ft2m;
bqf8 = 8.080458001106;
bqd7 = 0.056565957538*efact;
bqd6c = -3.138988802384*efact;
bqd6b = -4.444636077355*efact;
bqf5b = 2.430643979095*efact;
bqf5a = 2.223200792307*efact;
bqd4b = -1.467820200265*efact;
bqd4a = -1.47391*efact;
bqf3b = 0.71978*efact;
bqf3a = 2.96367*efact;
bqd2 = -3.52328*efact;
bqf1 = 15.241*efact;
bqd0 = -25.7555*efact;
aqf8 = 0.8e-2;
aqd7 = 0.8e-2;
aqd6c = 0.8e-2;
aqd6b = 0.8e-2;
aqf5b = 0.8e-2;
aqf5a = 0.8e-2;
aqd4b = 0.8e-2;
aqd4a = 0.8e-2;
aqf3b = 0.016;
aqf3a = 0.016;
aqd2 = 0.016;
aqf1 = 0.034;
aqd0 = 0.02;
kqf8 = bqf8/aqf8/brho;
kqd7 = bqd7/aqd7/brho;
kqd6c = bqd6c/aqd6c/brho;
kqd6b = bqd6b/aqd6b/brho;
kqf5b = bqf5b/aqf5b/brho;
kqf5a = bqf5a/aqf5a/brho;
kqd4b = bqd4b/aqd4b/brho;
kqd4a = bqd4a/aqd4a/brho;
kqf3b = bqf3b/aqf3b/brho;
kqf3a = bqf3a/aqf3a/brho;
kqd2 = bqd2/aqd2/brho;
kqf1 = bqf1/aqf1/brho;
kqd0 = bqd0/aqd0/brho;
lsf6 = 3.28*ft2m;
lsf5 = 3.28*ft2m;
lsd4 = 3.28*ft2m;
lsf1 = 2.0*ft2m;
lsd0 = 2.0*ft2m;
bsf6 = -1.04*bsign*efact;
bsf5 = 0.3472*bsign*efact;
bsd4 = -1.4037*bsign*efact;
bsf1 = 15.87309*bsign*efact;
bsd0 = -19.76126*bsign*efact;
asf6 = 0.8e-2;
asf5 = 0.8e-2;
asd4 = 0.8e-2;
asf1 = 0.034;
asd0 = 0.02;
ksf6 = 2.0857;
ksf5 = -0.65055;
ksd4 = 2.6357;
ksf1 = -1.6311;
ksd0 = 5.9237;
ksf6d = 1.0/lsca/lsca/lsca*ksf6*mfac;
ksf5d = 1.0/lsca/lsca/lsca*ksf5*mfac;
ksd4d = 1.0/lsca/lsca/lsca*ksd4*mfac;
ksf1d = 1.0/lsca/lsca/lsca*ksf1*mfac;
ksd0d = 1.0/lsca/lsca/lsca*ksd0*mfac;
energy = e0;
emitx = emitxn/(energy/0.51099906e-3);
emity = emityn/(energy/0.51099906e-3);
sig66 = esprd*esprd;
qx = 4.161075063395;
qy = 6.86628460456;
qs = 0.0;
qx_ = -11.95662271795;
qy_ = 31.723218651627;
alfx0 = -1.95133;
alfy0 = 0.605865;
betx0 = 64.171;
bety0 = 18.244;
x0 = 0.0;
px0 = 0.0;
y0 = 0.0;
py0 = 0.0;
t0 = 0.0;
pt0 = 0.0;
kltcq4 = 0.04776785;
kltcq3 = 0.04579693;
kltcq2 = -0.1319859e-2;
kltcq1 = 0.04982816;
kltcq0 = 0.06739183;
betay0 = 18.2438;

qlinac: quadrupole, type=linq, l=l_linq/2.0*myscale, aperture=a_linq;
qebm: quadrupole, type=ebm, l=l_ebm/2.0*myscale, aperture=a_ebm;
qecol: quadrupole, type=q_ecol, l=l_q_ecol/2.0*myscale, aperture=a_q_ecol;
qmi: quadrupole, type=q_ecmi, l=l_q_ecmi/2.0*myscale, aperture=a_q_ecmi;
secol: sextupole, type=s_ecol, l=l_s_ecol/4.0*myscale, aperture=a_s_ecol;
s_dogleg: sextupole, l=l_s_dglg/4.0*myscale, aperture=a_s_dglg;
becol1: sbend, l=l_bend/2.0*myscale, e2=0.0, hgap=a_ecol, fint=0.5,
fintx=0.0;
becol2: sbend, l=l_bend/2.0*myscale, e1=0.0, hgap=a_ecol, fint=0.0,
fintx=0.5;

bndfact = becol1[l]/coll_brho*mybendfac;

b1_ebm: sbend, l=l_ebm/2.0*myscale, hgap=a_ebm;
b2_ebm: sbend, l=l_ebm/2.0*myscale, e1=0.0, hgap=a_ebm, fint=0.5,
fintx=0.0;

bebmfact = b1_ebm[l]/coll_brho;
thetf1 = bpole*xf1/b1_ebm[hgap]*bebmfact;
thetd1 = -bpole*xd1/b1_ebm[hgap]*bebmfact;
thetf2 = bpole*xf2/b1_ebm[hgap]*bebmfact;
thetd2 = -bpole*xd2/b1_ebm[hgap]*bebmfact;
thetmf1 = bpole*xmf1/b1_ebm[hgap]*bebmfact;
thetmf2 = bpole*xmf2/b1_ebm[hgap]*bebmfact;
thetmd = -bpole*xmd/b1_ebm[hgap]*bebmfact;
thetsf1 = bpole*xsf1/b1_ebm[hgap]*bebmfact;
thetsf2 = bpole*xsf2/b1_ebm[hgap]*bebmfact;
thetsd = -bpole*xsd/b1_ebm[hgap]*bebmfact;

sbdkick: drift, l=1.5*myscale;
dmm: drift, l=lmm*myscale;
bpmq: marker;
bpmmb: marker;
bpmfb: marker;
di1b1: drift, l=0.390440685;
dd301: drift, l=14.0;
dd301fb1: drift, l=dd301[l]-(lfbc/2.0+lspace+lfbk+lspace);
dd301fb2: drift, l=lfbc/2.0+lspace+lfbk/2.0;
dd301fb3: drift, l=lfbk/2.0+lspace;
dd302: drift, l=10.0-l_linq;
dd302fb1: drift, l=dd302[l]-(lfbc/2.0+lspace+lfbk+lspace);
dd302fb2: drift, l=lfbc/2.0+lspace+lfbk/2.0;
dd302fb3: drift, l=lfbk/2.0+lspace;
dd302ws1: drift, l=lws/2.0+lspace;
dd302ws2: drift, l=dd302[l]-(lws/2.0+lspace);
mmover: marker; ! type=laser_marker;
ws: marker; ! type=laser_marker;
fbcxy: marker;! type=fb_corr;
fbkxy: marker;! type=fb_kick;
diag3beg: marker;
diag3fit: marker;
diag3end: marker;
m21: marker;
qqd300: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd300*qlinfact, aperture=a_linq;
qqd301: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd301*qlinfact, aperture=a_linq;
qqd302: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd302*qlinfact, aperture=a_linq;
qqd303: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd303*qlinfact, aperture=a_linq;
qqd304: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd304*qlinfact, aperture=a_linq;
qqd305: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd305*qlinfact, aperture=a_linq;
qqd306: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd306*qlinfact, aperture=a_linq;
qqd307: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd307*qlinfact, aperture=a_linq;
qqd308: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd308*qlinfact, aperture=a_linq;
qqd309: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd309*qlinfact, aperture=a_linq;
qqd310: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd310*qlinfact, aperture=a_linq;
qqd311: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd311*qlinfact, aperture=a_linq;
d1ec: drift, l=0.2*myscale;
d2ec: drift, l=3.2528*myscale;
d3ec: drift, l=8.0*myscale;
d4ec: drift, l=8.0*myscale;
qf1ec: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bf1ec*qlinfact, aperture=a_linq;
qd1ec: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bd1ec*qlinfact, aperture=a_linq;
qf2ec: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bf2ec*qlinfact, aperture=a_linq;
qf3ec: qlinac, type=linq, l=l_linq/2.0*myscale, k1=bf3ec*qlinfact, aperture=a_linq;
dmm0: drift, l=lmm;
qebm0: quadrupole, type=ebm, l=l_ebm/2.0, aperture=a_ebm;
d1bc: drift, l=0.2;
d2bc: drift, l=0.25;
d3bc: drift, l=ld3bc;
d4bc: drift, l=ld4bc;
qf1bc: qebm0, type=ebm, l=l_ebm/2.0, k1=bf1bc*qebmfact, aperture=a_ebm;
qd1bc: qebm0, type=ebm, l=l_ebm/2.0, k1=bd1bc*qebmfact, aperture=a_ebm;
qf2bc: qebm0, type=ebm, l=l_ebm/2.0, k1=bf2bc*qebmfact, aperture=a_ebm;
dbd1: drift, l=15.82195*myscale;
dbd2: drift, l=17.5*myscale;
dbd4: drift, l=0.4*myscale;
qfdb1: qebm, type=ebm, l=l_ebm/2.0*myscale, k1=bfdb1*qebmfact, aperture=a_ebm;
qddb1: qebm, type=ebm, l=l_ebm/2.0*myscale, k1=bddb1*qebmfact, aperture=a_ebm;
qddb2: qebm, type=ebm, l=l_ebm/2.0*myscale, k1=bddb2*qebmfact, aperture=a_ebm;
qddb3: qebm, type=ebm, l=l_ebm/2.0*myscale, k1=bddb3*qebmfact, aperture=a_ebm;
engysp: rcol, type=ensp, material="Beryllium", xsize=3.51*mm, ysize=a_ecol, l=18*mm;
engyab: rcol, type=enab, material="Titanium", xsize=5.4*mm, ysize=a_ecol,  
l=0.7*m;
engyspm: marker;
engyabm: marker;
esm: marker;
dummy: marker;
trans: marker;

dutil: drift, l=2.0*myscale;
dmm2: drift, l=lmm/2.0*myscale;

dmm2esp: drift, l= dmm2[l]-engysp[l]/2;
dutileab: drift, l=dutil[l]-engyab[l]/2;

dsxl: drift, l=l_s_ecol*myscale;
ygirder: drift, l=l_s_ecol*myscale;
!octec1: multipole, l=1*nm, knl={0,0,0,boctec1*ecolofac}, ksl={0,0,0,0};
!octec2: multipole, l=1*nm, knl={0,0,0,boctec2*ecolofac}, ksl={0,0,0,0};
octec3: multipole, l=1*nm, knl={0,0,0,boctec3*ecolofac}, ksl={0,0,0,0};
!octec3 : marker;
octec4: multipole, l=1*nm, knl={0,0,0,boctec4*ecolofac}, ksl={0,0,0,0};
!octec4 : marker;
b1a: becol1, l=l_bend/2.0*myscale, angle=bb1*bndfact, e1=bb1*bndfact, e2=0.0, hgap=a_ecol, fint=0.5, fintx=0.0;
b1b: becol2, l=l_bend/2.0*myscale, angle=bb1*bndfact, e1=0.0, e2=bb1*bndfact, hgap=a_ecol, fint=0.0, fintx=0.5;
b2a: becol1, l=l_bend/2.0*myscale, angle=bb2*bndfact, e1=bb2*bndfact, e2=0.0, hgap=a_ecol, fint=0.5, fintx=0.0;
b2b: becol2, l=l_bend/2.0*myscale, angle=bb2*bndfact, e1=0.0, e2=bb2*bndfact, hgap=a_ecol, fint=0.0, fintx=0.5;
b3a: becol1, l=l_bend/2.0*myscale, angle=bb3*bndfact, e1=bb3*bndfact, e2=0.0, hgap=a_ecol, fint=0.5, fintx=0.0;
b3b: becol2, l=l_bend/2.0*myscale, angle=bb3*bndfact, e1=0.0, e2=bb3*bndfact, hgap=a_ecol, fint=0.0, fintx=0.5;
b4a: becol1, l=l_bend/2.0*myscale, angle=bb4*bndfact, e1=bb4*bndfact, e2=0.0, hgap=a_ecol, fint=0.5, fintx=0.0;
b4b: becol2, l=l_bend/2.0*myscale, angle=bb4*bndfact, e1=0.0, e2=bb4*bndfact, hgap=a_ecol, fint=0.0, fintx=0.5;
qfec: qecol, type=q_ecol, l=l_q_ecol/2.0*myscale, k1=bqec*qecfact, aperture=a_q_ecol;
qdec: qecol, type=q_ecol, l=l_q_ecol/2.0*myscale, k1=-bqec*qecfact, aperture=a_q_ecol;
sxec1: secol, type=s_ecol, l=l_s_ecol/4.0*myscale, k2=bsxec1*secfact,
aperture=a_s_ecol;
syec1: secol, type=s_ecol, l=l_s_ecol/4.0*myscale, k2=bsyec1*secfact,
aperture=a_s_ecol;
sxec2: secol, type=s_ecol, l=l_s_ecol/4.0*myscale, k2=bsxec2*secfact,
aperture=a_s_ecol;
syec2: secol, type=s_ecol, l=l_s_ecol/4.0*myscale, k2=bsyec2*secfact,
aperture=a_s_ecol;
dmi: drift, l=lmi*myscale;
dmi2: drift, l=lmi*myscale-4.0*sbdkick[l]-2.0*lmm*myscale;
dshort: drift, l=0.3*myscale;
qfmi: qmi, type=q_ecmi, l=l_q_ecmi/2.0*myscale, k1=bmi*qmifact, aperture=a_q_ecmi;
qdmi: qmi, type=q_ecmi, l=l_q_ecmi/2.0*myscale, k1=-bmi*qmifact, aperture=a_q_ecmi;
xsp: rcol, type=xspt, material="Beryllium", xsize=x_aper, ysize=a_ebm, l=0.180*m;
ysp: rcol, type=yspt, material="Beryllium", xsize=a_ebm, ysize=y_aper, l=0.180*m;
!xsp : marker;
!ysp : marker;
xab: ecol, type=xabt, material="Titanium", xsize=1.0e-3, ysize=1.0e-3,l=0.7*m;
yab: ecol, type=yabt, material="Titanium", xsize=1.0e-3, ysize=1.0e-3,l=0.7*m;
!pc: ecol, type=pct, material="Beryllium", xsize=0.6e-2, ysize=0.6e-2;
!xab : marker;
!yab : marker;
pc  : marker;
xsp1: xsp, type=xspt;
xsp2: xsp, type=xspt;
xsp3: xsp, type=xspt;
xsp4: xsp, type=xspt;
ysp1: ysp, type=yspt;
ysp2: ysp, type=yspt;
ysp3: ysp, type=yspt;
ysp4: ysp, type=yspt;
xsp1m : marker;
ysp1m : marker;
xsp2m : marker;
ysp2m : marker;
xsp3m : marker;
ysp3m : marker;
xsp4m : marker;
ysp4m : marker;

xab1: xab, type=xabt;
xab2: xab, type=xabt;
xab3: xab, type=xabt;
xab4: xab, type=xabt;
yab1: yab, type=yabt;
yab2: yab, type=yabt;
yab3: yab, type=yabt;
yab4: yab, type=yabt;
d1bcol: drift, l=35.16752;
d1bcolxab: drift, l=d1bcol[l]-xab[l];
d2bcol: drift, l=15.82195;
d3bcol: drift, l=0.43201;
d2bcolyab: drift, l=d2bcol[l]-yab[l];
d2bcolsp: drift, l=d2bcol[l]-ysp[l]-xsp[l];
qfbcol: qebm0, type=ebm, l=l_ebm/2.0, k1=bfbcol*qebm0fact, aperture=a_ebm;
qdbcol: qebm0, type=ebm, l=l_ebm/2.0, k1=bdbcol*qebm0fact, aperture=a_ebm;
qcbcol: qebm0, type=ebm, l=l_ebm/2.0, k1=bcbcol*qebm0fact, aperture=a_ebm;
bqfdl1a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetf1, k1=bpole*qebm0fact, e1=thetf1, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqfdl1b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetf1, k1=bpole*qebm0fact, e1=0.0, e2=thetf1, hgap=a_ebm, fint=0.5, fintx=0.0;
bqddl1a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetd1, k1=-bpole*qebm0fact, e1=thetd1, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqddl1b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetd1, k1=-bpole*qebm0fact, e1=0.0, e2=thetd1, hgap=a_ebm, fint=0.5, fintx=0.0;
bqfdl2a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetf2, k1=bpole*qebm0fact, e1=thetf2, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqfdl2b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetf2, k1=bpole*qebm0fact, e1=0.0, e2=thetf2, hgap=a_ebm, fint=0.5, fintx=0.0;
bqddl2a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetd2, k1=-bpole*qebm0fact, e1=thetd2, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqddl2b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetd2, k1=-bpole*qebm0fact, e1=0.0, e2=thetd2, hgap=a_ebm, fint=0.5, fintx=0.0;
bqmf1a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetmf1, k1=bpole*qebm0fact, e1=thetmf1, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqmf1b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetmf1, k1=bpole*qebm0fact, e1=0.0, e2=thetmf1, hgap=a_ebm, fint=0.5, fintx=0.0;
bqmf2a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetmf2, k1=bpole*qebm0fact, e1=thetmf2, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqmf2b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetmf2, k1=bpole*qebm0fact, e1=0.0, e2=thetmf2, hgap=a_ebm, fint=0.5, fintx=0.0;
bqmda: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetmd, k1=-bpole*qebm0fact, e1=thetmd, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqmdb: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetmd, k1=-bpole*qebm0fact, e1=0.0, e2=thetmd, hgap=a_ebm, fint=0.5, fintx=0.0;
bqsf1a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetsf1, k1=bpole*qebm0fact, e1=thetsf1, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqsf1b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetsf1, k1=bpole*qebm0fact, e1=0.0, e2=thetsf1, hgap=a_ebm, fint=0.5, fintx=0.0;
bqsf2a: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetsf2, k1=bpole*qebm0fact, e1=thetsf2, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqsf2b: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetsf2, k1=bpole*qebm0fact, e1=0.0, e2=thetsf2, hgap=a_ebm, fint=0.5, fintx=0.0;
bqsda: b1_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetsd, k1=-bpole*qebm0fact, e1=thetsd, e2=0.0, hgap=a_ebm, fint=0.5, fintx=0.0;
bqsdb: b2_ebm, type=ebm, l=l_ebm/2.0*myscale, angle=thetsd, k1=-bpole*qebm0fact, e1=0.0, e2=thetsd, hgap=a_ebm, fint=0.5, fintx=0.0;
qfdl: qebm, type=ebm, l=l_ebm/2.0*myscale, k1=bpole*qebmfact, aperture=a_ebm;
qddl: qebm, type=ebm, l=l_ebm/2.0*myscale, k1=-bpole*qebmfact, aperture=a_ebm;
sxdl1: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsxdl1*sdogfact, aperture=a_s_dglg;
sydl1: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsydl1*sdogfact, aperture=a_s_dglg;
sxdl2: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsxdl2*sdogfact, aperture=a_s_dglg;
sydl2: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsydl2*sdogfact, aperture=a_s_dglg;
sxdl3: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsxdl3*sdogfact, aperture=a_s_dglg;
sydl3: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsydl3*sdogfact, aperture=a_s_dglg;
sxdl4: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsxdl4*sdogfact, aperture=a_s_dglg;
sydl4: s_dogleg, l=l_s_dglg/4.0*myscale, k2=bsydl4*sdogfact, aperture=a_s_dglg;
begecol: marker;
endecol: marker;
begpcol: marker;
endpcol: marker;
ltcq0: quadrupole, l=0.4525417, k1=kltcq0;
ltcq1: quadrupole, l=5.0, k1=kqltc1;
ltcq2: quadrupole, l=5.0, k1=kqltc2;
ltcq3: quadrupole, l=5.0, k1=kqltc3;
ltcq4: quadrupole, l=5.0, k1=kqltc4;
ltcd0: drift, l=lltcd0;
ltcd1: drift, l=lltcd1;
ltcd2: drift, l=lltcd2;
ltcd3: drift, l=lltcd3;
ltcd4: drift, l=lltcd4;
btfq1: quadrupole, l=5.0, k1=kbtfq1;
btfq2: quadrupole, l=5.0, k1=kbtfq2;
btfq3: quadrupole, l=5.0, k1=kbtfq3;
btfq4: quadrupole, l=5.0, k1=kbtfq4;
btfd0: drift, l=lbtfd0;
btfd1: drift, l=lbtfd1;
btfd2: drift, l=lbtfd2;
btfd3: drift, l=lbtfd3;
btfd4: drift, l=lbtfd4;
qbds1: quadrupole, l=lsca*l50cm/2.0;
qbds2: quadrupole, l=lsca*l1m/2.0;
qbds3: quadrupole, l=lsca*l2m/2.0;
xcorfbck: drift, l=lsca*l25cm;
ycorfbck: drift, l=lsca*l25cm;
sffb4: sbend, l=lsca*lffb4/2.0/nbfac, angle=0.15927e-5;
sffb3: sbend, l=lsca*lffb3/2.0/nbfac, angle=-0.22699e-4;
sffb3b: sbend, l=lsca*lffb3b/2.0/nbfac, angle=-0.45857e-5;
sffb2: sbend, l=lsca*lffb2/2.0/nbfac, angle=thetab2, outR=5*cm;
sffb1: sbend, l=lsca*lffb1/2.0/nbfac, angle=-0.14e-5;
sffsb1: sbend, l=lsca*lffsb1/2.0/nbfac, angle=thetasb1;
sffsb1a: sbend, l=lsca*lffsb1a/2.0/nbfac, angle=thetasb1a;
sffsb1b: sbend, l=lsca*lffsb1b/2.0/nbfac, angle=thetasb1b;
mbffb4: matrix, kick6=0.188e-14*gama^3.0*thetab4^2.0/lffb4*2.0*nbfac;
mbffb3: matrix, kick6=0.188e-14*gama^3.0*thetab3^2.0/lffb3*2.0*nbfac;
mbffb3b: matrix, kick6=0.188e-14*gama^3.0*thetab3b^2.0/lffb3b*2.0*nbfac;
mbffb2: matrix, kick6=0.188e-14*gama^3.0*thetab2^2.0/lffb2*2.0*nbfac;
mbffb1: matrix, kick6=0.188e-14*gama^3.0*thetab1^2.0/lffb1*2.0*nbfac;
mbffsb1: matrix, kick6=0.188e-14*gama^3.0*thetasb1^2.0/lffsb1*2.0*nbfac;
mbffsb1a: matrix, kick6=0.188e-14*gama^3.0*thetasb1a^2.0/lffsb1a*2.0*nbfac;
mbffsb1b: matrix, kick6=0.188e-14*gama^3.0*thetasb1b^2.0/lffsb1b*2.0*nbfac;
qmd5: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd5*qfact3;
qmd6: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd6*qfact3;
qmd7: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd7*qfact3;
qmd8: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd8*qfact3;
qdf90: qbds2, l=lsca*l1m/2.0, k1=1.0/lsca/lsca*bqdf90*qfact2;
qdd90: qbds2, l=lsca*l1m/2.0, k1=1.0/lsca/lsca*bqdd90*qfact2;
qdf2: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqdf2*qfact3;
qdd2: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqdd2*qfact3;
sq1: qbds1, l=lsca*l50cm/2.0, k1=1.0/lsca/lsca*bsq1*qfact1, tilt=0.785398163397;
sq2: qbds1, l=lsca*l50cm/2.0, k1=1.0/lsca/lsca*bsq2*qfact1, tilt=0.785398163397;
sq3: qbds1, l=lsca*l50cm/2.0, k1=1.0/lsca/lsca*bsq3*qfact1, tilt=0.785398163397;
sq4: qbds1, l=lsca*l50cm/2.0, k1=1.0/lsca/lsca*bsq4*qfact1, tilt=0.785398163397;
qmd9: qbds2, l=lsca*l1m/2.0, k1=1.0/lsca/lsca*bqmd9*qfact2;
qmd10: qbds2, l=lsca*l1m/2.0, k1=1.0/lsca/lsca*bqmd10*qfact2;
qdd: qbds2, l=lsca*l1m/2.0, k1=1.0/lsca/lsca*bqdd*qfact2;
qdf: qbds2, l=lsca*l1m/2.0, k1=1.0/lsca/lsca*bqdf*qfact2;
qmd11: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd11*qfact3;
qmd12: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd12*qfact3;
qmd13: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd13*qfact3;
qmd14: qbds3, l=lsca*l2m/2.0, k1=1.0/lsca/lsca*bqmd14*qfact3;
qf8: quadrupole, l=lsca*lqf8/2.0, k1=1.0/lsca/lsca*kqf8;
qd7: quadrupole, l=lsca*lqd7/2.0, k1=1.0/lsca/lsca*kqd7;
qd6c: quadrupole, l=lsca*lqd6c/2.0, k1=-0.021346/lsca/lsca;
qd6b: quadrupole, l=lsca*lqd6b/2.0, k1=-0.032590718909/lsca/lsca;
qf5b: quadrupole, l=lsca*lqf5b/2.0, k1=1.0/lsca/lsca*kqf5b;
qf5a: quadrupole, l=lsca*lqf5a/2.0, k1=1.0/lsca/lsca*kqf5a;
qd4b: quadrupole, l=lsca*lqd4b/2.0, k1=1.0/lsca/lsca*kqd4b;
qd4a: quadrupole, l=lsca*lqd4a/2.0, k1=1.0/lsca/lsca*kqd4a;
qf3b: quadrupole, l=lsca*lqf3b/2.0, k1=1.0/lsca/lsca*kqf3b;
qf3a: quadrupole, l=lsca*lqf3a/2.0, k1=1.0/lsca/lsca*kqf3a;
qd2: quadrupole, l=lsca*lqd2/2.0, k1=1.0/lsca/lsca*kqd2;
qf1: quadrupole, l=lsca*lqf1/2.0, k1=1.0/lsca/lsca*kqf1;
qd0: quadrupole, l=lsca*lqd0/2.0+ltrick/2.0, k1=1.0/lsca/lsca*kqd0;
sf6: sextupole, l=lsca*lsf6/2.0, k2=ksf6d*tscal;
sf5: sextupole, l=lsca*lsf5/2.0, k2=ksf5d*tscal;
sd4: sextupole, l=lsca*lsd4/2.0, k2=ksd4d*tscal;
sf1: sextupole, l=lsca*lsf1/2.0, k2=ksf1d*tscal;
sd0: sextupole, l=lsca*lsd0/2.0, k2=ksd0d*tscal;
!sd0: drift, l=lsca*lsd0/2.0;

lmd1: drift, l=lsca*2.5;
lmd2: drift, l=lsca*1.803;
ld1sk: drift, l=lsca*11.705;
ld2a90: drift, l=lsca*14.508;
ld2sk: drift, l=lsca*2.313;
ld3sk: drift, l=lsca*6.758810956106;
ld4: drift, l=lsca*8.96504994685;
ld3a: drift, l=lsca*9.571810956106;
ldcorr: drift, l=lsca*0.25;
ld2b90: drift, l=lsca*10.695;
ldcorr2: drift, l=lsca*9.85;
ld2c90: drift, l=lsca*3.908;
ld4sk: drift, l=lsca*1.0;
lmd3: drift, l=lsca*6.420346326541;
ld2a: drift, l=lsca*8.39216;
ld1: drift, l=lsca*0.5;
ld2c: drift, l=lsca*6.89216;
ld2b: drift, l=lsca*7.89216;
lmd4: drift, l=lsca*6.968971864391;
lmd5: drift, l=lsca*4.244867662675;
lmd6: drift, l=lsca*7.2;
lmd7: drift, l=lsca*7.2;
lmd8: drift, l=lsca*6.9046;
oct1: multipole, l=1*nm, knl={0,0,0,kod1}, ksl={0,0,0,0};
oct2: multipole, l=1*nm, knl={0,0,0,kod2}, ksl={0,0,0,0};
octdrift: drift, l=lod;
d5odt: drift, l=18.0;
d5od: drift, l=lsca*1.4910001e2*ft2m-lod-18.0;
d4od: drift, l=lsca*35.49178*ft2m-lod;
d2od: drift, l=lsca*53.49556*ft2m-lod;
d1od: drift, l=lsca*10.0*ft2m-lod;
d8: drift, l=lsca*30.0*ft2m;
lx0: drift, l=lsca*1.0*ft2m;
d7: drift, l=lsca*31.806;
d6: drift, l=lsca*1.1572e2*ft2m;
d3: drift, l=lsca*94.90237*ft2m;
dec0: drift, l=lsca*2.0*ft2m;
d0: drift, l=lsca*14.1*ft2m-ltrick;
xcfbws9: xcorfbck, l=lsca*l25cm;
xcfbws10: xcorfbck, l=lsca*l25cm;
xcfbws11: xcorfbck, l=lsca*l25cm;
xcfbws12: xcorfbck, l=lsca*l25cm;
ycfbws9: ycorfbck, l=lsca*l25cm;
ycfbws10: ycorfbck, l=lsca*l25cm;
ycfbws11: ycorfbck, l=lsca*l25cm;
ycfbws12: ycorfbck, l=lsca*l25cm;
wsl9: marker;
wsl10: marker;
wsl11: marker;
wsl12: marker;
begeff1: marker;
ffin: marker;
ip: marker;
endeff1: marker;
coll1esp: marker;
coll2eab: marker;
coll3esp: marker;
coll4eab: marker;
coll5bdsp: marker;
coll6bfsp: marker;
coll7bfab: marker;
coll8bdab: marker;
coll9bdsp: marker;
coll10bfsp: marker;
coll11bfab: marker;
coll12bdab: marker;
coll13bdsp: marker;
coll14bfsp: marker;
coll15bfab: marker;
coll16bdab: marker;
coll17bdsp: marker;
coll18bfsp: marker;
coll19bfab: marker;
coll20bdab: marker;

qd0mark1: marker;
qd0mark2: marker;

!octf1: multipole, l=1*nm, knl={0,0,0,boctf}, ksl={0,0,0,0};
octd1: multipole, l=1*nm, knl={0,0,0,boctd}, ksl={0,0,0,0};
octf1 : marker;
!octd1 : marker;
octdoub: line=(oct1,octdrift,oct2);
d5: line=(d5od,octdrift,oct2,d5odt);
d4: line=(d4od,octdrift,oct1);
d2: line=(d2od,octdrift);
d1: line=(d1od,octdrift);
qqd300m: line=(qqd300,bpmq,bpmmb,mmover,qqd300);
qqd301m: line=(qqd301,bpmq,bpmmb,mmover,qqd301);
qqd302m: line=(qqd302,bpmq,bpmmb,mmover,qqd302);
qqd303m: line=(qqd303,bpmq,bpmmb,mmover,qqd303);
qqd304m: line=(qqd304,bpmq,bpmmb,mmover,qqd304);
qqd305m: line=(qqd305,bpmq,bpmmb,mmover,qqd305);
qqd306m: line=(qqd306,bpmq,bpmmb,mmover,qqd306);
qqd307m: line=(qqd307,bpmq,bpmmb,mmover,qqd307);
qqd308m: line=(qqd308,bpmq,bpmmb,mmover,qqd308);
qqd309m: line=(qqd309,bpmq,bpmmb,mmover,qqd309);
qqd310m: line=(qqd310,bpmq,bpmmb,mmover,qqd310);
qqd311m: line=(qqd311,bpmq,bpmmb,mmover,qqd311);
celd3: line=(qqd301,qqd301,dd302,qqd302,qqd302,dd302);
diag3: line=(diag3beg,di1b1,qqd300m,dd301fb1,fbcxy,dd301fb2,fbkxy,dd301fb3,
diag3fit,qqd301m,dd302,qqd302m,dd302,qqd303m,dd302,qqd304m,dd302fb1,fbcxy,
dd302fb2,fbkxy,dd302fb3,qqd305m,dd302ws1,ws,bpmq,dd302ws2,qqd306m,dd302,
qqd307m,dd302ws1,ws,bpmq,dd302ws2,qqd308m,dd302,qqd309m,dd302ws1,ws,bpmq,
dd302ws2,qqd310m,dd302,qqd311m,dd302ws1,ws,bpmq,diag3end);
qf1ecm: line=(qf1ec,bpmq,qf1ec);
qd1ecm: line=(qd1ec,bpmq,qd1ec);
qf2ecm: line=(qf2ec,bpmq,qf2ec);
qf3ecm: line=(qf3ec,bpmq,qf3ec);
ltecoll: line=(d1ec,qf3ecm,dmm,qf3ecm,d4ec,qf2ecm,dmm,qf2ecm,d3ec,qd1ecm,dmm,
qd1ecm,dmm,qd1ecm,d2ec,qf1ecm,dmm,qf1ecm,d1ec);
qf1bcm: line=(qf1bc,bpmq,qf1bc);
qd1bcm: line=(qd1bc,bpmq,qd1bc);
qf2bcm: line=(qf2bc,bpmq,qf2bc);
etbcoll: line=(d1bc,qf1bcm,dmm0,qf1bcm,d2bc,qd1bcm,dmm0,qd1bcm,dmm0,qd1bcm,
dmm0,qd1bcm,d3bc,qf2bcm,dmm0,qf2bcm,d4bc);
qfdb1m: line=(qfdb1,bpmq,qfdb1);
qddb1m: line=(qddb1,bpmq,qddb1);
qddb2m: line=(qddb2,bpmq,qddb2);
qddb3m: line=(qddb3,bpmq,qddb3);
btdogleg: line=(qddb1m,dbd1,qfdb1m,dbd2,qddb2m,dbd3,qddb3m,dbd4);
b1: line=(b1a,b1b);
b2: line=(b2a,b2b);
b3: line=(b3a,b3b);
b4: line=(b4a,b4b);

m0: marker;
m1: marker;
m2: marker;
m3: marker;
m4: marker;
m5: marker;
m6: marker;
m7: marker;
m8: marker;
m9: marker;
m10: marker;
m11: marker;
m12: marker;
m13: marker;

qfecm: line=(qfec,bpmq,qfec);
qdecm: line=(qdec,bpmq,qdec);
sxec1m: line=(4*sxec1);
syec1m: line=(4*syec1);
sxec2m: line=(4*sxec2);
syec2m: line=(4*syec2);
!ecolhcel(b01, sx, sy, spoiler, dmarker): line=(dutil,qfecm,dmm,sx,dmm,bline(
!b01),dmm2,spoiler,dmarker,dmm2,bline(b01),dmm,sy,dmm,qdecm,dutil);
mb2: line=(-b2);
mb1: line=(-b1);
mb4: line=(-b4);
mb3: line=(-b3);

!bline(b): line=(5*(b,dmm),b);
!ecol1: line=(ecolhcel(b1,dsxl,syec1m,dummy,dummy));
ecol1: line=(dutil,qfecm,dmm,dsxl,dmm,b1,dmm,b1,dmm,b1,dmm,b1,dmm,b1,dmm,
b1,dmm2,dummy,dummy,dmm2,b1,dmm,b1,dmm,b1,dmm,b1,dmm,b1,dmm,b1,dmm,
syec1m,dmm,qdecm,dutil);
!ecol2: line=(-ecolhcel(mb2,sxec1m,syec1m,engysp,coll1esp));
ecol2: line=-(dutil,qfecm,dmm,sxec1m,dmm,mb2,dmm,mb2,dmm,mb2,dmm,mb2,dmm,mb2,dmm,
mb2,dmm2,engysp,coll1esp,dmm2,mb2,dmm,mb2,dmm,mb2,dmm,mb2,dmm,mb2,dmm,mb2,dmm,
syec1m,dmm,qdecm,dutil);
!ecol3: line=(ecolhcel(b2,sxec1m,syec1m,dummy,dummy));
ecol3: line=(dutil,qfecm,dmm,sxec1m,dmm,b2,dmm,b2,dmm,b2,dmm,b2,dmm,b2,dmm,
b2,dmm2,dummy,dummy,dmm2,b2,dmm,b2,dmm,b2,dmm,b2,dmm,b2,dmm,b2,dmm,
syec1m,dmm,qdecm,dutil);
!ecol4: line=(-ecolhcel(mb1,dsxl,syec1m,dummy,dummy));
ecol4: line=-(dutil,qfecm,dmm,dsxl,dmm,mb1,dmm,mb1,dmm,mb1,dmm,mb1,dmm,mb1,dmm,
mb1,dmm2,dummy,dummy,dmm2,mb1,dmm,mb1,dmm,mb1,dmm,mb1,dmm,mb1,dmm,mb1,dmm,
syec1m,dmm,qdecm,dutil);
!ecol5: line=(ecolhcel(b3,ygirder,syec2m,dummy,dummy));
ecol5: line=(dutil,qfecm,dmm,ygirder,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,
b3,dmm2,dummy,dummy,dmm2,b3,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,
syec2m,dmm,qdecm,dutil);
!ecol6: line=(-ecolhcel(mb4,sxec2m,syec2m,engysp,coll3esp));
ecol6: line=-(dutileab,qfecm,dmm,sxec2m,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,
mb4,dmm2esp,engysp,engyspm,coll3esp,dmm2esp,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,
syec2m,dmm,qdecm,dutil);
!ecol7: line=(ecolhcel(b4,sxec2m,syec2m,dummy,dummy));
ecol7: line=(dutileab,qfecm,dmm,sxec2m,dmm,b4,dmm,b4,dmm,b4,dmm,b4,dmm,b4,dmm,
b4,dmm2,dummy,dummy,dmm2,b4,dmm,b4,dmm,b4,dmm,b4,dmm,b4,dmm,b4,dmm,
syec2m,dmm,qdecm,dutil);
!ecol8: line=(-ecolhcel(mb3,dsxl,syec2m,dummy,dummy));
ecol8: line=-(dutil,qfecm,dmm,dsxl,dmm,mb3,dmm,mb3,dmm,mb3,dmm,mb3,dmm,mb3,dmm,
mb3,dmm2,dummy,dummy,dmm2,mb3,dmm,mb3,dmm,mb3,dmm,mb3,dmm,mb3,dmm,mb3,dmm,
syec2m,dmm,qdecm,dutil);
qfmim: line=(qfmi,bpmq,qfmi);
qdmim: line=(qdmi,bpmq,qdmi);
minusi: line=(dshort,qfmim,dmi,qdmim,dshort,dshort,qdmim,dmi,qfmim,dshort,
dshort,qfmim,dmi,qdmim,dshort,dshort,qdmim,dmi2,2*sbdkick,dmm,sbdkick,dmm,
sbdkick,dmm,qfmim,dshort);
nlcecoll: line=(ecol5,octec3,ecol6,engyabm,engyab,coll4eab,ecol7,octec4,ecol8);
qfbcolm: line=(qfbcol,bpmq,qfbcol);
qdbcolm: line=(qdbcol,bpmq,qdbcol);
qcbcolm: line=(qcbcol,bpmq,qcbcol);
!bcol1(colf, cold, dmarker1, dmarker2): line=(d3bcol,pc,qcbcolm,d1bcol,colf,
!dmarker1,qfbcolm,d2bcol,cold,dmarker2);
!bcol2(colf, cold, dmarker1, dmarker2): line=(qdbcolm,cold,dmarker2,d2bcol,colf,
!dmarker1,qfbcolm,d1bcol,pc,qcbcolm,d3bcol);
!nlcbcoll: line=(bcol1(pc,pc,dummy,dummy),bcol2(xsp1,ysp1,coll6bfsp,coll5bdsp);
!,bcol1(xab1,yab1,coll7bfab,coll8bdab),bcol2(xsp2,ysp2,coll10bfsp,coll9bdsp),
!bcol1(xab2,yab2,coll11bfab,coll12bdab),bcol2(xsp3,ysp3,coll14bfsp,coll13bdsp);
!,bcol1(xab3,yab3,coll15bfab,coll16bdab),bcol2(xsp4,ysp4,coll18bfsp,
!coll17bdsp),bcol1(xab4,yab4,coll19bfab,coll20bdab));

nlcbcoll: line=(d3bcol,pc,qcbcolm,d1bcol,pc,dummy,qfbcolm,d2bcol,pc,dummy,
qdbcolm,ysp1m,ysp1,coll5bdsp,d2bcolsp,xsp1m,xsp1,coll6bfsp,
qfbcolm,d1bcol,pc,qcbcolm,d3bcol,d3bcol,pc,qcbcolm,
d1bcolxab,xab1m,xab1,coll7bfab,qfbcolm,d2bcolyab,yab1m,yab1,coll8bdab,
qdbcolm,ysp2m,ysp2,coll9bdsp,d2bcolsp,xsp2m,xsp2,coll10bfsp,
qfbcolm,d1bcol,pc,qcbcolm,d3bcol,d3bcol,pc,qcbcolm,
d1bcolxab,xab2m,xab2,coll11bfab,qfbcolm,d2bcolyab,yab2m,yab2,coll12bdab,
qdbcolm,ysp3m,ysp3,coll13bdsp,d2bcolsp,xsp3m,xsp3,coll14bfsp,
qfbcolm,d1bcol,pc,qcbcolm,d3bcol,d3bcol,pc,qcbcolm,
d1bcolxab,xab3m,xab3,coll15bfab,qfbcolm,d2bcolyab,yab3m,yab3,coll16bdab,
qdbcolm,ysp4m,ysp4,coll17bdsp,d2bcolsp,xsp4m,xsp4,coll18bfsp,
qfbcolm,d1bcol,pc,qcbcolm,d3bcol,d3bcol,pc,qcbcolm,
d1bcolxab,xab4m,xab4,coll19bfab,qfbcolm,d2bcolyab,yab4m,yab4,coll20bdab);

bqfdl1: line=(bqfdl1a,bpmq,bqfdl1b);
bqddl1: line=(bqddl1a,bpmq,bqddl1b);
bqfdl2: line=(bqfdl2a,bpmq,bqfdl2b);
bqddl2: line=(bqddl2a,bpmq,bqddl2b);
bqmf1: line=(bqmf1a,bpmq,bqmf1b);
bqmf2: line=(bqmf2a,bpmq,bqmf2b);
bqmd: line=(bqmda,bpmq,bqmdb);
bqsf1: line=(bqsf1a,bpmq,bqsf1b);
bqsf2: line=(bqsf2a,bpmq,bqsf2b);
bqsd: line=(bqsda,bpmq,bqsdb);
qfdlm: line=(qfdl,bpmq,qfdl);
qddlm: line=(qddl,bpmq,qddl);
sxdl1m: line=(4*sxdl1);
sydl1m: line=(4*sydl1);
sxdl2m: line=(4*sxdl2);
sydl2m: line=(4*sydl2);
sxdl3m: line=(4*sxdl3);
sydl3m: line=(4*sydl3);
sxdl4m: line=(4*sxdl4);
sydl4m: line=(4*sydl4);
spool: line=(2*dmm);
!dlbasic(sx, sy, qf1, qd, qf2): line=(sx,dmm,qf1,dmm,qf1,4*dmm,qd,dmm,qd,dmm,
!sy,dmm,qd,dmm,qd,4*dmm,qf2,dmm,qf2,dmm);
!dletain: line=(dlbasic(spool,spool,bqmf1,bqmd,bqmf2));
dletain: line=(spool,dmm,bqmf1,dmm,bqmf1,4*dmm,bqmd,dmm,bqmd,dmm,
spool,dmm,bqmd,dmm,bqmd,4*dmm,bqmf2,dmm,bqmf2,dmm);
!dlcell1: line=(dlbasic(sxdl1m,sydl1m,bqfdl1,bqddl1,bqfdl1));
dlcell1: line=(sxdl1m,dmm,bqfdl1,dmm,bqfdl1,4*dmm,bqddl1,dmm,bqddl1,dmm,
sydl1m,dmm,bqddl1,dmm,bqddl1,4*dmm,bqfdl1,dmm,bqfdl1,dmm);
!dlcell2: line=(dlbasic(sxdl2m,sydl2m,bqfdl1,bqddl1,bqfdl1));
dlcell2: line=(sxdl2m,dmm,bqfdl1,dmm,bqfdl1,4*dmm,bqddl1,dmm,bqddl1,dmm,
sydl2m,dmm,bqddl1,dmm,bqddl1,4*dmm,bqfdl1,dmm,bqfdl1,dmm);
!dl90deg: line=(dlbasic(spool,spool,qfdlm,qddlm,qfdlm));
dl90deg: line=(spool,dmm,qfdlm,dmm,qfdlm,4*dmm,qddlm,dmm,qddlm,dmm,
spool,dmm,qddlm,dmm,qddlm,4*dmm,qfdlm,dmm,qfdlm,dmm);
!dlcell3: line=(dlbasic(sxdl3m,sydl3m,bqfdl2,bqddl2,bqfdl2));
dlcell3: line=(sxdl3m,dmm,bqfdl2,dmm,bqfdl2,4*dmm,bqddl2,dmm,bqddl2,dmm,
sydl3m,dmm,bqddl2,dmm,bqddl2,4*dmm,bqfdl2,dmm,bqfdl2,dmm);
!dlcell4: line=(dlbasic(sxdl4m,sydl4m,bqfdl2,bqddl2,bqfdl2));
dlcell4: line=(sxdl4m,dmm,bqfdl2,dmm,bqfdl2,4*dmm,bqddl2,dmm,bqddl2,dmm,
sydl4m,dmm,bqddl2,dmm,bqddl2,4*dmm,bqfdl2,dmm,bqfdl2,dmm);
!dletaout: line=(dlbasic(spool,spool,bqsf2,bqsd,bqsf1));
dletaout: line=(spool,dmm,bqsf2,dmm,bqsf2,4*dmm,bqsd,dmm,bqsd,dmm,
spool,dmm,bqsd,dmm,bqsd,4*dmm,bqsf1,dmm,bqsf1,dmm);
dogleg: line=(dl90deg,dletain,dlcell1,dlcell2,2*dl90deg,dlcell3,dlcell4,
dletaout);
postcoll: line=(diag3,ltecoll,nlcecoll,etbcoll,nlcbcoll,btdogleg,dogleg);
ecoll: line=(begecol,postcoll,endecol);
pcol: line=(begpcol,postcoll,endpcol);

entr: line=(ltcq0,ltcd0,ltcq1,ltcd1,ltcq2,ltcd2,ltcq3,ltcd3,ltcq4,ltcd4);
btff: line=(btfd0,btfq1,btfd1,btfq2,btfd2,btfq3,btfd3,btfq4,btfd4);
mycolsys: line=(entr,trans,nlcecoll,trans,etbcoll,trans,nlcbcoll,trans,btff,
trans);

ksent  = 0.0;
!ksent  = -20.0;
sent : multipole, l=1*nm, knl={0,0,ksent }, ksl={0,0,0};
ksexit = -20.0;
sexit: multipole, l=1*nm, knl={0,0,ksexit}, ksl={0,0,0};
noentrcolsys: line=(sent ,trans,nlcecoll,trans,etbcoll, sexit, trans,nlcbcoll,trans,btff, trans);
ffb4: line=(sffb4,mbffb4,sffb4,mbffb4,sffb4,mbffb4,sffb4,mbffb4,sffb4,mbffb4);
ffb3: line=(sffb3,mbffb3,sffb3,mbffb3,sffb3,mbffb3,sffb3,mbffb3,sffb3,mbffb3);
ffb3b: line=(sffb3b,mbffb3b,sffb3b,mbffb3b,sffb3b,mbffb3b,sffb3b,mbffb3b,sffb3b,mbffb3b);
ffb2: line=(sffb2,mbffb2,sffb2,mbffb2,sffb2,mbffb2,sffb2,mbffb2,sffb2,mbffb2);
ffb1: line=(sffb1,mbffb1,sffb1,mbffb1,sffb1,mbffb1,sffb1,mbffb1,sffb1,mbffb1);
ffsb1: line=(sffsb1,mbffsb1,sffsb1,mbffsb1,sffsb1,mbffsb1,sffsb1,mbffsb1,sffsb1,mbffsb1);
ffsb1a: line=(sffsb1a,mbffsb1a,sffsb1a,mbffsb1a,sffsb1a,mbffsb1a,sffsb1a,mbffsb1a,sffsb1a,mbffsb1a);
ffsb1b: line=(sffsb1b,mbffsb1b,sffsb1b,mbffsb1b,sffsb1b,mbffsb1b,sffsb1b,mbffsb1b,sffsb1b,mbffsb1b);
qmd5m: line=(qmd5,bpmq,qmd5);
qmd6m: line=(qmd6,bpmq,qmd6);
qmd7m: line=(qmd7,bpmq,qmd7);
qmd8m: line=(qmd8,bpmq,qmd8);
qqdf90m: line=(qdf90,bpmq,qdf90);
qdd90m: line=(qdd90,bpmmb,qdd90);
qdf2m: line=(qdf2,bpmq,qdf2);
qdd2m: line=(qdd2,bpmq,qdd2);
qmd9m: line=(qmd9,bpmmb,qmd9);
qmd10m: line=(qmd10,bpmq,qmd10);
qddm: line=(qdd,bpmmb,qdd);
qdfm: line=(qdf,bpmq,qdf);
sq1m: line=(2*sq1);
sq2m: line=(2*sq2);
sq3m: line=(2*sq3);
sq4m: line=(2*sq4);
qmd11m: line=(qmd11,bpmq,qmd11);
qmd12m: line=(qmd12,bpmq,qmd12);
qf90m: line=(qf90,bpmq,qf90);
qd90m: line=(qd90,bpmq,qd90);
qmd13m: line=(qmd13,bpmq,qmd13);
qmd14m: line=(qmd14,bpmq,qmd14);
qmd15m: line=(qmd15,bpmq,qmd15);
qmd16m: line=(qmd16,bpmq,qmd16);
sk90c: line=(qdf90,ld2a90,qdd90,qdd90,ld2a90,qdf90);
bbdgn: line=(lmd1,qmd5m,lmd1,qmd6m,lmd1,qmd7m,lmd1,qmd8m,lmd2);
skcor: line=(sq1m,ld1sk,qdf90m,ld2a90,qdd90m,ld2sk,sq2m,ld3sk,qdf2m,ld4,
qdd2m,ld4,qdf2m,ld3a,qdd90m,ld2sk,sq3m,ldcorr,ycfbws9,ldcorr,xcfbws9,ld2b90,
qdf90m,ldcorr,ycfbws10,ldcorr2,xcfbws10,ld2c90,qmd9m,ld2sk,sq4m,ld4sk,wsl9,
bpmfb,ldcorr,ycfbws11,ldcorr,xcfbws11,lmd3,qmd10m,ld2a,qddm,ld1,wsl10,bpmfb,
ldcorr,xcfbws12,ldcorr,ycfbws12,ld2c,qdfm,ld2a,qddm,ld1,wsl11,bpmfb,ld2b,
qdfm,ld2a,qddm,ld1,wsl12,bpmfb);
ffm: line=(lmd4,qmd11,qmd11,lmd5,qmd12,qmd12,lmd6,qmd13,qmd13,lmd7,qmd14,
qmd14,lmd8,qf8);
ffpr: line=(ffin,qf8,d8,qd7,qd7,lx0,ffb4,ffb4,lx0,qd6c,qd6c,d7,qd6b,qd6b,lx0,
qd6b,qd6b,d6,sf6,sf6,lx0,qf5b,qf5b,lx0,qf5b,qf5b,d5,qf5a,qf5a,lx0,qf5a,qf5a,
lx0,sf5,sf5,d4,qd4b,qd4b,lx0,sd4,sd4,lx0,qd4a,qd4a,d3,ffb3,ffb3,lx0,qf3b,
qf3b,lx0,ffb3b,ffb3b,lx0,qf3a,qf3a,lx0,ffb2,ffb2,lx0,qd2,qd2,lx0,ffb1,ffb1,
lx0,ffsb1,ffsb1,lx0,ffsb1,ffsb1,lx0,ffsb1,ffsb1,lx0,ffsb1,ffsb1,lx0,ffsb1a,
ffsb1a,lx0,ffsb1b,ffsb1b,d2,sf1,sf1,m21,lx0,qf1,octf1,qf1,d1,dec0,lx0,sd0,
sd0,lx0,qd0mark1,qd0,octd1,qd0mark2,qd0,d0,ip);
testo: line=(ffin,qf8,d8,qd7,qd7,lx0,ffb4,ffb4,lx0,qd6c,qd6c,d7,qd6b,qd6b,
lx0,qd6b,qd6b,d6,sf6,sf6,lx0,qf5b,qf5b,lx0,qf5b,qf5b,d5,qf5a,qf5a,lx0,qf5a,
qf5a,lx0,sf5,sf5,d4,qd4b,qd4b,lx0,sd4,sd4,lx0,qd4a,qd4a,d3,ffb3,ffb3,lx0,
qf3b,qf3b,lx0,ffb3b,ffb3b,lx0,qf3a,qf3a,lx0,ffb2,ffb2,lx0,qd2,qd2,lx0,ffb1,
ffb1,lx0,ffsb1,ffsb1,lx0,ffsb1,ffsb1,lx0,ffsb1,ffsb1,lx0,ffsb1,ffsb1,lx0,
ffsb1a,ffsb1a,lx0,ffsb1b,ffsb1b,d2,sf1,sf1,m21);
eff1: line=(begeff1,bbdgn,skcor,ffm,ffpr,endeff1);
eff1s: line=(ffm,ffpr,endeff1);
beamlinei: line=(-beamline);
bds: line=(mycolsys,eff1s);

noentrbds: line=(noentrcolsys,eff1s);
shortff: line=(mycolsys,eff1s);
testo2: line=(mycolsys,ffm,testo);

ecol5temp: line=(m1,dutil,m2,qfecm,m3,dmm,m4,ygirder,m5,dmm,m6,b3,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,
b3,m7,dmm2,m8,dummy,dummy,dmm2,b3,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,b3,dmm,
syec2m,dmm,qdecm,dutil);

ecol6temp: line=-(m12,dutil,m11,qfecm,m10,dmm,m9,sxec2m,m8,dmm,m7,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,
mb4,dmm2,engysp,coll3esp,dmm2,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,mb4,dmm,
m6,syec2m,m5,dmm,m4,qdecm,m3,dutil,m2);

!temp: line = (entr,trans,ecol5,m1,octec3,ecol6temp);
dri: drift, l=50e-6*m;
testSextupole: sextupole, l=1*m, k2=87;
testSbend: sbend, l=1*m,angle=0.001;
temp : line = (dri,m1,testSbend,m2);
use,period=bds;
!sample,range=ip;
!sample,range=m0;
!sample,range=m1;
!sample,range=m2;
!sample,range=m3;
!sample,range=m4;
!sample,range=m5;
!sample,range=m6;
!sample,range=m7;
!sample,range=m8;
!sample,range=m9;
!sample,range=m10;
!sample,range=m11;
!sample,range=m12;
!sample,range=m13;
!dump,range=engyspm;
!dump,range=engyabm;
!dump,range=xsp1m;
!dump,range=xsp2m;
!dump,range=xsp3m;
!dump,range=xsp4m;
!dump,range=ysp1m;
!dump,range=ysp2m;
!dump,range=ysp3m;
!dump,range=ysp4m;
sample, range=bpmq;
sample, range=m21;
!sample, range=engyab;
sample, range=dummy;
sample, range=trans;
sample, range=pc;
sample, range=ffin;
sample, range=ip;
sample, range=coll3esp;
sample, range=coll4eab;
sample, range=coll5bdsp;
sample, range=coll6bfsp;
sample, range=coll7bfab;
sample, range=coll8bdab;
sample, range=coll9bdsp;
sample, range=coll10bfsp;
sample, range=coll11bfab;
sample, range=coll12bdab;
sample, range=coll13bdsp;
sample, range=coll14bfsp;
sample, range=coll15bfab;
sample, range=coll16bdab;
sample, range=coll17bdsp;
sample, range=coll18bfsp;
sample, range=coll19bfab;
sample, range=coll20bdab;
sample, range=qd0mark1;
sample, range=qd0mark2;

option,ngenerate=10000,backgroundScaleFactor=1e5,
!synchRadOn=0,srTrackPhotons=0,srLowX=0.5,srLowGamE=10*MeV,
physicsList="em_standard";

!twiss, srRescale=1;