//Test Input for monitor output
//author: Sebastian Thomas
//mail: sthomas@uni-mainz.de
//last change (ymd): 2019/12/18

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

//SYSTEM, "rm monitor1.loss";
//SYSTEM, "rm monitor2.loss";
//SYSTEM, "rm monitor4.loss";

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


quaddub1: quadrupole, elemedge = 0.27, l = 0.08, k1 = -0.20242810945;
quaddub2: quadrupole, elemedge = 0.4, l = 0.08, k1 = -0.0453881472;
quaddub3: quadrupole, elemedge = 0.75, l = 0.08, k1 = 0.29538547175;
quaddub4: quadrupole, elemedge = 0.88, l = 0.08, k1 = -0.25782907015; 

moni1:monitor, elemedge = 0.01, outfn = "monitor1";
moni2:monitor, elemedge = 1.00, outfn = "monitor2";


testline: LINE = (moni1,quaddub1,quaddub2,quaddub3,moni2);


Dist1:DISTRIBUTION, TYPE=GAUSS, 
	SIGMAX=2.64773468e-4, sigmapx=8.214924502e-4, corrx=0, //correlation x tp px (R12)
	SIGMAY=3.89799397e-4, sigmapy=4.836093778e-4, corry=0, //R34
	SIGMAZ=577.5e-6, sigmapZ=4e-3, corrZ=0; //R56

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

/////////////////////////////////////////////////////////////////
/// If you change the NPART from 1000 to 4000 the monitor{}.loss file will append 3000 particles and add a new header in line 1002

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

///////////////////////////////////////////////////////////////////

SELECT, Line=testline;

track,line=testline, beam=beam1, DT=timestep, MAXSTEPS=nstep;
	run, METHOD="PARALLEL-T", BEAM=beam1, fieldsolver=Fs1, DISTRIBUTION=DIST1;
ENDTRACK;
STOP;