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

Title,string="Proton";

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

value,{gamma,brho,Edes,beta,gambet};

dt1 = 3.3355e-12;
dt2 = 1.0e-12;;
cdt = dt2;

// L:        physical element lenght (real)
// KS:       field scaling factor (real)
// ELEMEDGE: physical start of the element on the floor (real)

Q1: QUADRUPOLE, L=0.20, K1=14.08450704, ELEMEDGE=0.10;
Q3: QUADRUPOLE, L=0.01, K1=0.0, ELEMEDGE=1.5276;


l1: Line = (Q1,Q3);



/*
 We can specity ONE distribution

 Note: sigmapt is in eV 

*/
qb=0.0;
rf=1498.956e6;
v0=beta*CLIGHT;
T1s = 36.6e-12;
lz  = T1s*v0;
value,{v0,lz,qb,rf};

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= 0.00040, pt=P0*1e9, sigmapt=0.0,corrt=0.0;

/*
   We can specity ONE field solver object
*/

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=cdt, ZSTOP=1.527647;
 run, method = "PARALLEL-T", beam=beam1, fieldsolver=Fs1, distribution=Dist1;
endtrack;
Stop;
Quit;