{Scott's 2016-04-05 200MeV lattice + Holger's 2016-04-13 fieldmaps with Scott's adjustments}

#R=5;

#lqf  = +1.3647784895686557e-01;
#lqd  = +1.2446026779894361e-01;
#lfd  = +7.0000000000000007e-02;
#ldf  = +1.2000000000000000e-01;

#thqf=2*asin(0.5*lqf/R);
#thqd=2*asin(0.5*lqd/R);
#thfd=2*asin(0.5*lfd/R);
#thdf=2*asin(0.5*ldf/R);

{Drift Length Angle}
DF #R*thdf# #thdf*-1#
QF #R*thqf# #thqf*-1#
FD #R*thfd# #thfd*-1#
QD #R*thqd# #thqd*-1#
HDF #DF.Length/2# #DF.Angle/2#
HFD #FD.Length/2# #FD.Angle/2#
HDFS #HDF.Length#
QFS #QF.Length#
HFDS #HFD.Length#
QDS #QD.Length#

#scalF = +1.010282587908758e+00;
#scalD = +1.005830514157469e+00;
#dispF = +4.782740318800096e-04;
#dispD = -2.541293642894889e-04;

{TrilinearFieldMap File Strength Xrel MirrorY=1}
FMF Cbeta\Iron_2016-04-13\160413-FieldMap_QfOn3_polar.txt #scalF*-1# #dispF-R#
FMD Cbeta\Iron_2016-04-13\160413-FieldMap_QdOn3_polar.txt #scalD*-1# #dispD-R#

CellFM: HDF,QF,HFD,FMF,FMD,HFD,QD,HDF;

#xf   = -7.9568051765194405e-03 +2.523940959531643e-03 +1.169211897450673e-04 +dispF; /**/
#xd   = +2.0336097668811705e-02 -1.849083154885980e-03 -5.380253526274470e-04 +dispD; // Wrong from here on

{TrilinearFieldMap File Strength Xrel MirrorY=1}
FMFS #FMF.File# #scalF*-1# #dispF-R-xf#
FMDS #FMD.File# #scalD*-1# #dispD-R-xd#

CellStraightFM: HDFS,QFS,HFDS,FMFS,FMDS,HFDS,QDS,HDFS;

#EmassMeV=0.510998928;
{MatchScan Estart Egoal Estep Species=Electron AllowUnstable=1}
Match #200-EmassMeV#MeV #54.5-EmassMeV#MeV 48.5MeV
MatchFine #Match.Estart# #Match.Egoal# 1MeV
MatchStraight #Match.Estart# #Match.Egoal# #Match.Estep#

{Match-Aperture} 
MatchEnd

CellFM,Match,MatchFine,CellFM,MatchEnd,CellFM;
#//CellStraightFM,MatchStraight,CellStraightFM,MatchEnd,CellStraightFM;
#//CellFM,Particles,10*CellFM,MatchEnd,CellFM;
#//4*CellFM;