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;