{CornellD = Like CornellA but based on slide 3 of Cornell-Doublet-Magnets-Oct16-2014.pptx
NB: Curvature is opposite to cell lattices.  To perform mirror-image, invert Angle and Dipole but *not Quad.  Also invert x and x' for particles.}

#eRMS=0e-6; // Set to zero for no misalignment errors
#eRMSquadRel=0e-2;
#eSysSext=0;
#eSysOct=0;

{Multipole Length Angle Dipole Quad Fringe=1.5cm AlignMode=Integral ErrorX=#eRMS# ErrorY=#eRMS# ErrorZ=#eRMS# ErrorRelQuad=#eRMSquadRel# Sextupole=#eSysSext# Octupole=#eSysOct#}	   
QF 8.94cm 0.024218889 0.21 -39.1
BD 12.05cm -0.087050742 -0.56 25.42

{Drift Length}
D1 3cm
D2 12.138cm
HD2 #D2.Length/2#

Cell: HD2,QF,D1,BD,HD2;

{MatchScan Estart Egoal Estep Species}
Match 286MeV 76MeV 70MeV Electron

{Match-Aperture} 
MatchEnd

#RingCells=50*2;
{Multipole Length Angle Dipole Quad Fringe AlignMode=Integral ErrorX=#eRMS# ErrorY=#eRMS# ErrorZ=#eRMS# ErrorRelQuad=#eRMSquadRel# Sextupole=#eSysSext# Octupole=#eSysOct#}
#Acells=10*2; fcubic=0.2;
#for (n=1;n<=Acells;n++)
#{
#	f=n/(Acells+1);
#	f=(1-fcubic)*f+fcubic*(3*f*f-2*f*f*f);
QF#n# #QF.Length# #f*QF.Angle# #f*QF.Dipole# #QF.Quad# #QF.Fringe#
BD#n# #BD.Length# #f*BD.Angle# #f*BD.Dipole# #BD.Quad# #BD.Fringe#
Cell#n#: HD2,QF#n#,D1,BD#n#,HD2;
#}

STA: Cell1#for (n=2;n<=Acells;n++) {#,Cell#n}#;
ATS: Cell#Acells# #for (n=Acells-1;n>=1;n--) {#,Cell#n}#;

Arc: #RingCells/2-Acells/2#*Cell;

QFS #QF.Length# 0 0 #QF.Quad# #QF.Fringe#
BDS #BD.Length# 0 0 #BD.Quad# #BD.Fringe#

StraightCell: HD2,QFS,D1,BDS,HD2;
#Scells=14*2;
Straight: #Scells#*StraightCell;

{Drift Length}
SplitterSpace 4m

{BMPImport File PixelSize Zhere Xhere}
CornellPlan 'cornell_area.bmp' 6.85cm 27.5m 13.55m

Ring1: Arc,ATS,Straight,STA,Arc;
Ring: CornellPlan,SplitterSpace,Ring1,SplitterSpace;

global_bore=20m
dt=5e-12

{Particle Species=Electron Energy X vx vz=1 Trail=-1}	   
P286 286MeV 0.00963161348 0.016255949	   
P216 216MeV -0.0028542648 -0.00644622481	   
P146 146MeV -0.00904905799 -0.0351376209	   
P76 76MeV -0.00640387947 -0.0742182004	   
PS286 286MeV	   
PS216 216MeV	   
PS146 146MeV	   
PS76 76MeV	 

ArcParticles:P76,P146,P216,P286;
StraightParticles:PS76,PS146,PS216,PS286;
 
{Beam Energy betax alphax betay alphay x RotZX TransEmitNorm=500e-6 Type=Twiss Species=Electron Population=1000 TransDist=Waterbag LongDist=Uniform}	   
B286 286MeV 0.601716 -0.823071 1.84473 2.3703 0.00963161348 0.016255949	   
B216 216MeV 0.479319 -0.890997 0.90815 1.56454 -0.0028542648 -0.00644622481	   
B146 146MeV 0.340456 -1.02423 0.484006 1.28572 -0.00904905799 -0.0351376209	   
B76 76MeV 0.319468 -2.97958 0.24484 1.6294 -0.00640387947 -0.0742182004	 

ArcBeams:B76,B146,B216,B286;

RingOnly: CornellPlan,SplitterSpace,Ring1,SplitterSpace;
RingParticles: CornellPlan,SplitterSpace,ArcParticles,Ring1,MatchEnd,SplitterSpace;
RingBeams: CornellPlan,SplitterSpace,ArcBeams,Ring1,MatchEnd,SplitterSpace;

#//Cell,Match,Cell,MatchEnd,Cell;
RingParticles;
#//RingOnly;
#//RingBeams;
#//StraightParticles,STA,Arc;