Option, TFS=FALSE;
Option, ECHO=FALSE;
Option, PSDUMPFREQ=1;

Title,string="Proton";

Edes=0.5;
gamma=(Edes+PMASS)/PMASS;
beta=sqrt(1-(1/gamma^2));
gambet=gamma*beta;
P0 = gamma*beta*PMASS;
brho = (PMASS*1.0e9*gambet) / CLIGHT;

AND1:SBEND, L=2.855796, ELEMEDGE=1.0, ALPHA=0.0, K0=2.0, FMAPFN='FAND1.T7', DESIGNENERGY=500.0, EXITANGLE=0.0;
Q3: QUADRUPOLE, L=0.01, K1=0.0, ELEMEDGE=4.0;


l1: Line = (AND1,Q3);

Dist1:DISTRIBUTION, DISTRIBUTION=gauss,
                    sigmax=  0.005, sigmapx=0.0, corrx=0.0,
                    sigmay=  0.005, sigmapy=0.0, corry=0.0,
                    t = 0.0, sigmat= 1.0e-13, pt=P0*1e9, sigmapt=0.0,corrt=0.0;


Fs1:FIELDSOLVER, FSTYPE=None, MX=32, MY=32, MT=2048, 
                 PARFFTX=false, PARFFTY=false, PARFFTT=true,
                 BCFFTX=open, BCFFTY=open, BCFFTT=open, 
                 BBOXINCR=1, GREENSF=INTEGRATED;

beam1: BEAM, PARTICLE=PROTON, pc=P0, NPART=1e5, CHARGE=1, SPACECHARGE=false;

Select, Line=l1;

track,line=l1, beam=beam1, MAXSTEPS=100000, DT=1.0e-11, ZSTOP=4.1;
 run, method = "PARALLEL-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1;
endtrack;
Stop;
Quit;