//FILE: FODOUMERrev.MADX
//MADX Example 1: FODO cell
//Author: L. Deniau (CERN), S. Bernal (U. of MD)
//Date: July 6th, 2020

BEAM, PARTICLE=ELECTRON, PC=10158E-06;

D : DRIFT, L=0.10836;
QF: QUADRUPOLE, L=0.05164/2, K1= 168.36;
QD: QUADRUPOLE, L=0.05164  , K1=-168.36;

FODO: LINE=(QF,D,QD,D,QF);
USE, PERIOD=FODO;

! Non periodic solution
TWISS, BETX=0.9, BETY=0.7, ALFX=0, ALFY=0;
MYPLT: PLOT, HAXIS=S, HMAX=0.35,
             VAXIS1={BETX,BETY}, VAXIS2={MUX,MUY},
             COLOUR=100, INTERPOLATE;

! Periodic solution
TWISS;
EXEC, MYPLT;

! Match phase advance to 90 degrees
MATCH, SEQUENCE=FODO;
VARY, NAME=QF->K1;
VARY, NAME=QD->K1;
CONSTRAINT, RANGE=#E, MUX=0.25;
CONSTRAINT, RANGE=#E, MUY=0.25;
LMDIF, TOLERANCE=1e-20;
ENDMATCH;

! Matched solution
TWISS;
EXEC, MYPLT;