//Minimal example
//author: Sebastian Thomas
//mail: sthomas@uni-mainz.de
//last change (ymd): 2020/01/21

OPTION, ECHO=TRUE; //Output
OPTION, ASCIIDUMP=TRUE;
OPTION, VERSION=22000;


REAL nstep=10000;
REAL timestep=5E-12;
REAL Edes=5E-3;
REAL gamma=(Edes+EMASS)/EMASS;
REAL beta=sqrt(1-(1/gamma^2));
REAl P0=gamma*beta*EMASS;



DIPOL1: SBEND, L = 0.6572973205894721, ANGLE = 0.6769836773, E1 = 0.6769836773, E2 = 0.0,
        GAP = 0.02, PSI = -1.570796327, ELEMEDGE = 0,
        DESIGNENERGY = 30.0;

CAVITY1: RFCAVITY, L = 0.1154, ELEMEDGE = 1.0,
        VOLT = 1.8248202,
        FREQ = 1300.0,
        LAG = 0.8195368035664573;

beamline_1: LINE = (DIPOL1, CAVITY1);


Dist1:DISTRIBUTION,TYPE=GAUSS, 
	SIGMAX=2.64773468e-4, sigmapx=8.214924502e-4, corrx=0, 
	SIGMAY=3.89799397e-4, sigmapy=4.836093778e-4, corry=0, 
	SIGMAZ=577.5e-6, sigmapZ=4e-3, corrZ=0,
	writetofile = true;

Fs1:FIELDSOLVER, FSTYPE=NONE, MX=1, MY=1, MT=1, PARFFTX=false, PARFFTY=false, PARFFTT=true, BCFFTY=open, BCFFTY=open, BCFFTT=open, BBOXINCR=1, GREENSF=STANDARD;

beam1:BEAM, PARTICLE=ELECTRON, BFREQ=1300.0, BCURRENT=0.001, NPART=2, GAMMA=gamma;

track, line=beamline_1, beam=beam1, DT=timestep, MAXSTEPS=nstep;

run, METHOD="PARALLEL-T", BEAM=beam1, fieldsolver=Fs1, DISTRIBUTION=DIST1;

ENDTRACK;
STOP;